fixed typo in check.c

issp-center-dev · Apr 14, 2024 · 442e147 · 442e147
1 parent 262091b
commit 442e147
Showing 1 changed file with 5 additions and 4 deletions.
diff --git a/src/check.c b/src/check.c
@@ -172,16 +172,17 @@ int check(struct BindStruct *X){
     //comb_2 = Binomial(NCond,Nup-u_loc)
     //comb_3 = Binomial(NCond,Ndown+u_loc-NLocSpn)
     NLocSpn  = X->Def.NLocSpn;
-    NCond    = X->Def.Nsite-NLocSpn;
+    NCond    = X->Def.Ne-NLocSpn;
+    int NsCond   = X->Def.Nsite-NLocSpn;
     comb_1   = pow(2,NLocSpn);//Tpow is not defined
     comb_sum = 0;
     for(int tmp_Nup=0;tmp_Nup<=NCond;tmp_Nup++){
-      comb_2     = Binomial(NCond,tmp_Nup,comb,NCond);
-      comb_3     = Binomial(NCond,NCond-tmp_Nup,comb,NCond);
+      comb_2     = Binomial(NsCond,tmp_Nup,comb,NsCond);
+      comb_3     = Binomial(NsCond,NCond-tmp_Nup,comb,NsCond);
       comb_sum  += comb_1*comb_2*comb_3;
       printf("tmp %ld %ld %ld\n",comb_sum,comb_2,comb_3);
     }
-    printf("Ne=%d NLocSpn=%d NCond=%d comb_1=%ld comb_2=%ld comb_3=%ld comb_sum=%ld\n",X->Def.Ne,NLocSpn,NCond,comb_1,comb_2,comb_3,comb_sum);
+    //printf("Ne=%d NLocSpn=%d NsCond %d NCond=%d comb_1=%ld comb_2=%ld comb_3=%ld comb_sum=%ld\n",X->Def.Ne,NLocSpn,NsCond,NCond,comb_1,comb_2,comb_3,comb_sum);
     break;
 
   case KondoGC: