I am trying calculate the energy of the spin model for Jxy=1 and Jxy=-1.

Why that energy does not depend on the sign of the Jxy (worm , dirloop_sse , fulldiag)?

Best wishes.