diff --git a/src/check.c b/src/check.c index a97cc681..a6e07c2b 100644 --- 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: