From 3b591888b5dad7cef02170743a92e2bf9c5831db Mon Sep 17 00:00:00 2001 From: Ian Jauslin Date: Tue, 14 Jul 2015 13:37:40 +0000 Subject: Fix a sign error in the definition of A and B in kondo_preprocess 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. --- src/definitions.cpp | 2 +- src/kondo.c | 4 ++++ 2 files changed, 5 insertions(+), 1 deletion(-) (limited to 'src') 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); } -- cgit v1.2.3-54-g00ecf