Merge pull request #62 from k-ido/debug_lanczos

Debug complex Lanczos mode for systems with exchange and/or interall interactions
issp-center-dev · Oct 1, 2024 · 0e79814 · 0e79814
2 parents 14f95a6 + ea0ca5c
commit 0e79814
Show file tree

Hide file tree

Showing 2 changed files with 8 additions and 7 deletions.
diff --git a/src/mVMC/locgrn.c b/src/mVMC/locgrn.c
@@ -288,6 +288,10 @@ double complex GreenFuncN(const int n, int *rsi, int *rsj, const double complex
 
   MakeProjCnt(projCntNew,eleNum);
   x = ProjRatio(projCntNew,eleProjCnt);
+  if (FlagRBM) {
+    MakeRBMCnt(rbmCntNew, eleNum);
+    x *= RBMRatio(rbmCntNew,rbmCnt);
+  }
 
   /* calculateNewPfM */
   for(qpidx=0;qpidx<NQPFull;qpidx++) {
@@ -304,7 +308,7 @@ double complex GreenFuncN(const int n, int *rsi, int *rsj, const double complex
     eleNum[rsi[k]] = 0;
   }
 
-  return x*z/ip;
+  return conj(x*z/ip);
 }
 
 //

diff --git a/src/mVMC/lslocgrn.c b/src/mVMC/lslocgrn.c
@@ -626,12 +626,9 @@ double complex calHCACA1(const int ri, const int rj, const int rk, const int rl,
 
   ReleaseWorkSpaceInt();
   ReleaseWorkSpaceComplex();
-  if (FlagRBM) {
-    return e*conj(z*ipNew/ip);
-  }
-  else {
-    return e*z*ipNew/ip;
-  }
+
+  return e*conj(z*ipNew/ip);
+
 }
 
 /* calculate <psi| H C_is A_js C_kt A_lt|x>/<psi|x> for <psi|CACA|x>/<psi|x>=0 */