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.

1 files changed, 4 insertions, 0 deletions

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);
 	}