Fix a sign error in the definition of A and B in kondo_preprocessv1.2.1

The operators A and B introduced by kondo_preprocess had the wrong sign.
This bug does not affect the beta function for the hierarchical Kondo model.

2 files changed, 5 insertions, 1 deletions

diff --git a/src/definitions.cpp b/src/definitions.cpp
index d32537d..982f7a7 100644
--- a/src/definitions.cpp
+++ b/src/definitions.cpp
@@ -17,7 +17,7 @@ limitations under the License.
 #ifndef DEFINITIONS_GCC
 #define DEFINITIONS_GCC
 
-#define VERSION "1.2"
+#define VERSION "1.2.1"
 
 // number of entries in a configuration file
 #define ARG_COUNT 10
diff --git a/src/kondo.c b/src/kondo.c
index c86b4b3..39d74bb 100644
--- a/src/kondo.c
+++ b/src/kondo.c
@@ -1039,6 +1039,8 @@ int kondo_resolve_ABh(char* str, Polynomial* output, Fields_Table fields){
 	polynomial_conjugate(poly_conjugate);
 	polynomial_multiply_scalar(poly_conjugate, pauli_mat.matrix[a][b]);
 	polynomial_prod_chain(psi[a],&poly_conjugate,fields);
+	// correct sign: psi[a]^+ should be on the left of psi[b]^-
+	polynomial_multiply_Qscalar(poly_conjugate,quot(-1,1));
 	// add to poly
 	polynomial_concat_noinit(poly_conjugate, output);
       }
@@ -1235,6 +1237,8 @@ int kondo_polynomial_vector(int offset, int index, Polynomial (*polys)[3], Field
 	  polynomial_conjugate(poly_conjugate);
 	  polynomial_multiply_scalar(poly_conjugate, pauli_mat.matrix[a][b]);
 	  polynomial_prod_chain(psi[a],&poly_conjugate,fields);
+	  // correct sign: psi[a]^+ should be on the left of psi[b]^-
+	  polynomial_multiply_Qscalar(poly_conjugate,quot(-1,1));
 	  // add to polys[j]
 	  polynomial_concat_noinit(poly_conjugate, (*polys)+i);
 	}