Lines Matching refs:theindex
451 const int theindex = denT->gIndex(); in FillSite() local
452 const int DIM = max(denBK->gMaxDimAtBound(theindex), denBK->gMaxDimAtBound(theindex+1)); in FillSite()
460 set_2rdm_A_DMRG(theindex,theindex,theindex,theindex, 2*d1); in FillSite()
461 set_2rdm_B_DMRG(theindex,theindex,theindex,theindex,-2*d1); in FillSite()
471 for (int j_index=theindex+1; j_index<L; j_index++){ in FillSite()
472 if (denBK->gIrrep(j_index) == denBK->gIrrep(theindex)){ in FillSite()
478 const double d2 = doD2(denT, Ltens[theindex][j_index-theindex-1], workmem); in FillSite()
479 set_2rdm_A_DMRG(theindex,j_index,theindex,theindex, 2*d2); in FillSite()
480 set_2rdm_B_DMRG(theindex,j_index,theindex,theindex,-2*d2); in FillSite()
485 const int dimTriangle = L - theindex - 1; in FillSite()
500 … const double d3 = doD3(denT, S0tens[theindex][k_index-j_index][j_index-theindex-1], workmem); in FillSite()
501 set_2rdm_A_DMRG(theindex,theindex,j_index,k_index, 2*d3); in FillSite()
502 set_2rdm_B_DMRG(theindex,theindex,j_index,k_index,-2*d3); in FillSite()
510 … const double d4 = doD4(denT, F0tens[theindex][k_index-j_index][j_index-theindex-1], workmem); in FillSite()
511 … const double d5 = doD5(denT, F0tens[theindex][k_index-j_index][j_index-theindex-1], workmem); in FillSite()
512 … const double d6 = doD6(denT, F1tens[theindex][k_index-j_index][j_index-theindex-1], workmem); in FillSite()
513 set_2rdm_A_DMRG(theindex,j_index,k_index,theindex, -2*d4 - 2*d5 - 3*d6); in FillSite()
514 set_2rdm_B_DMRG(theindex,j_index,k_index,theindex, -2*d4 - 2*d5 + d6); in FillSite()
515 set_2rdm_A_DMRG(theindex,j_index,theindex,k_index, 4*d4 + 4*d5); in FillSite()
516 set_2rdm_B_DMRG(theindex,j_index,theindex,k_index, 2*d6); in FillSite()
522 for (int g_index=0; g_index<theindex; g_index++){ in FillSite()
523 if (denBK->gIrrep(g_index) == denBK->gIrrep(theindex)){ in FillSite()
529 const double d7 = doD7(denT, Ltens[theindex-1][theindex-g_index-1], workmem); in FillSite()
530 set_2rdm_A_DMRG(g_index,theindex,theindex,theindex, 2*d7); in FillSite()
531 set_2rdm_B_DMRG(g_index,theindex,theindex,theindex,-2*d7); in FillSite()
536 const int globalsize8to12 = theindex * ( L - 1 - theindex ); in FillSite()
539 const int g_index = gj_index % theindex; in FillSite()
540 const int j_index = ( gj_index / theindex ) + theindex + 1; in FillSite()
548 …const double d8 = doD8(denT, Ltens[theindex-1][theindex-g_index-1], Ltens[theindex][j_index-theind… in FillSite()
550 …doD9andD10andD11(denT, Ltens[theindex-1][theindex-g_index-1], Ltens[theindex][j_index-theindex-1],… in FillSite()
551 set_2rdm_A_DMRG(g_index,theindex,j_index,theindex, -4*d8-d9); in FillSite()
552 set_2rdm_A_DMRG(g_index,theindex,theindex,j_index, 2*d8 + d11); in FillSite()
553 set_2rdm_B_DMRG(g_index,theindex,j_index,theindex, d9 - 2*d10); in FillSite()
554 set_2rdm_B_DMRG(g_index,theindex,theindex,j_index, 2*d8 + 2*d10 - d11); in FillSite()
557 …const double d12 = doD12(denT, Ltens[theindex-1][theindex-g_index-1], Ltens[theindex][j_index-thei… in FillSite()
558 set_2rdm_A_DMRG(g_index,j_index,theindex,theindex, 2*d12); in FillSite()
559 set_2rdm_B_DMRG(g_index,j_index,theindex,theindex,-2*d12); in FillSite()
564 const int globalsize = theindex * upperboundTriangle; in FillSite()
567 const int g_index = gjk_index % theindex; in FillSite()
568 const int global = gjk_index / theindex; in FillSite()
574 const int cnt2 = j_index - theindex - 1; in FillSite()
576 …if (Irreps::directProd(I_g, denBK->gIrrep(theindex)) == Irreps::directProd(denBK->gIrrep(j_index),… in FillSite()
582 …const double d13 = doD13(denT, Ltens[theindex-1][theindex-g_index-1], S0tens[theindex][cnt1][cnt2]… in FillSite()
583 …const double d14 = doD14(denT, Ltens[theindex-1][theindex-g_index-1], S0tens[theindex][cnt1][cnt2]… in FillSite()
587 …d15 = doD15(denT, Ltens[theindex-1][theindex-g_index-1], S1tens[theindex][cnt1][cnt2], workmem, wo… in FillSite()
588 …d16 = doD16(denT, Ltens[theindex-1][theindex-g_index-1], S1tens[theindex][cnt1][cnt2], workmem, wo… in FillSite()
590 set_2rdm_A_DMRG(g_index,theindex,j_index,k_index, 2*d13 + 2*d14 + 3*d15 + 3*d16); in FillSite()
591 set_2rdm_A_DMRG(g_index,theindex,k_index,j_index, 2*d13 + 2*d14 - 3*d15 - 3*d16); in FillSite()
592 set_2rdm_B_DMRG(g_index,theindex,j_index,k_index,-2*d13 - 2*d14 + d15 + d16); in FillSite()
593 set_2rdm_B_DMRG(g_index,theindex,k_index,j_index,-2*d13 - 2*d14 - d15 - d16); in FillSite()
601 …const double d17 = doD17orD21(denT, Ltens[theindex-1][theindex-g_index-1], F0tens[theindex][cnt1][… in FillSite()
602 …const double d18 = doD18orD22(denT, Ltens[theindex-1][theindex-g_index-1], F0tens[theindex][cnt1][… in FillSite()
603 …const double d19 = doD19orD23(denT, Ltens[theindex-1][theindex-g_index-1], F1tens[theindex][cnt1][… in FillSite()
604 …const double d20 = doD20orD24(denT, Ltens[theindex-1][theindex-g_index-1], F1tens[theindex][cnt1][… in FillSite()
605 set_2rdm_A_DMRG(g_index,j_index,k_index,theindex, -2*d17 - 2*d18 - 3*d19 - 3*d20); in FillSite()
606 set_2rdm_A_DMRG(g_index,j_index,theindex,k_index, 4*d17 + 4*d18 ); in FillSite()
607 set_2rdm_B_DMRG(g_index,j_index,k_index,theindex, -2*d17 - 2*d18 + d19 + d20); in FillSite()
608 set_2rdm_B_DMRG(g_index,j_index,theindex,k_index, 2*d19 + 2*d20); in FillSite()
611 …const double d21 = doD17orD21(denT, Ltens[theindex-1][theindex-g_index-1], F0tens[theindex][cnt1][… in FillSite()
612 …const double d22 = doD18orD22(denT, Ltens[theindex-1][theindex-g_index-1], F0tens[theindex][cnt1][… in FillSite()
613 …const double d23 = doD19orD23(denT, Ltens[theindex-1][theindex-g_index-1], F1tens[theindex][cnt1][… in FillSite()
614 …const double d24 = doD20orD24(denT, Ltens[theindex-1][theindex-g_index-1], F1tens[theindex][cnt1][… in FillSite()
615 set_2rdm_A_DMRG(g_index,k_index,j_index,theindex, -2*d21 - 2*d22 - 3*d23 - 3*d24); in FillSite()
616 set_2rdm_A_DMRG(g_index,k_index,theindex,j_index, 4*d21 + 4*d22 ); in FillSite()
617 set_2rdm_B_DMRG(g_index,k_index,j_index,theindex, -2*d21 - 2*d22 + d23 + d24); in FillSite()
618 set_2rdm_B_DMRG(g_index,k_index,theindex,j_index, 2*d23 + 2*d24); in FillSite()
644 int theindex = denT->gIndex(); in doD1() local
648 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD1()
649 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD1()
651 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL,IL); in doD1()
652 int dimR = denBK->gCurrentDim(theindex+1,NL+2,TwoSL,IL); in doD1()
672 int theindex = denT->gIndex(); in doD2() local
676 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD2()
677 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD2()
681 int IRup = Irreps::directProd(IL, denBK->gIrrep(theindex)); in doD2()
683 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL,IL); in doD2()
684 int dimRdown = denBK->gCurrentDim(theindex+1,NL+2,TwoSL,IL); in doD2()
685 int dimRup = denBK->gCurrentDim(theindex+1,NL+1,TwoSR,IRup); in doD2()
717 int theindex = denT->gIndex(); in doD3() local
721 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD3()
722 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD3()
725 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL,IL); in doD3()
726 int dimRdown = denBK->gCurrentDim(theindex+1,NL, TwoSL,IL); in doD3()
727 int dimRup = denBK->gCurrentDim(theindex+1,NL+2,TwoSL,IL); in doD3()
757 int theindex = denT->gIndex(); in doD4() local
761 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD4()
762 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD4()
765 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL,IL); in doD4()
766 int dimR = denBK->gCurrentDim(theindex+1,NL+2,TwoSL,IL); in doD4()
795 int theindex = denT->gIndex(); in doD5() local
799 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD5()
800 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD5()
804 int IR = Irreps::directProd(IL,denBK->gIrrep(theindex)); in doD5()
805 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL,IL); in doD5()
806 int dimR = denBK->gCurrentDim(theindex+1,NL+1,TwoSR,IR); in doD5()
836 int theindex = denT->gIndex(); in doD6() local
840 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD6()
841 … for (int TwoSL = denBK->gTwoSmin(theindex,NL); TwoSL<= denBK->gTwoSmax(theindex,NL); TwoSL+=2){ in doD6()
846 int IR = Irreps::directProd(IL,denBK->gIrrep(theindex)); in doD6()
847 int dimL = denBK->gCurrentDim(theindex, NL, TwoSL, IL); in doD6()
848 int dimRup = denBK->gCurrentDim(theindex+1,NL+1,TwoSRup, IR); in doD6()
849 int dimRdown = denBK->gCurrentDim(theindex+1,NL+1,TwoSRdown,IR); in doD6()
883 int theindex = denT->gIndex(); in doD7() local
887 for (int NLup = denBK->gNmin(theindex); NLup<=denBK->gNmax(theindex); NLup++){ in doD7()
888 …for (int TwoSLup = denBK->gTwoSmin(theindex,NLup); TwoSLup<= denBK->gTwoSmax(theindex,NLup); TwoSL… in doD7()
891 int dimLup = denBK->gCurrentDim(theindex, NLup, TwoSLup, ILup); in doD7()
895 int IR = Irreps::directProd(ILup,denBK->gIrrep(theindex)); in doD7()
896 int dimLdown = denBK->gCurrentDim(theindex, NLup-1, TwoSLdown, IR); in doD7()
897 int dimR = denBK->gCurrentDim(theindex+1, NLup+1, TwoSLdown, IR); in doD7()
929 int theindex = denT->gIndex(); in doD8() local
933 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD8()
934 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD8()
937 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD8()
938 int dimRup = denBK->gCurrentDim(theindex+1, NL+2, TwoSLup, ILup); in doD8()
946 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD8()
947 int dimRdown = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, Idown); in doD8()
987 int theindex = denT->gIndex(); in doD9andD10andD11() local
989 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD9andD10andD11()
990 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD9andD10andD11()
993 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD9andD10andD11()
996 int IRup = Irreps::directProd(ILup, denBK->gIrrep(theindex)); in doD9andD10andD11()
998 int IRdown = Irreps::directProd(ILdown, denBK->gIrrep(theindex)); in doD9andD10andD11()
1005 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown); in doD9andD10andD11()
1006 int dimRup = denBK->gCurrentDim(theindex+1, NL+1, TwoSRup, IRup); in doD9andD10andD11()
1007 int dimRdown = denBK->gCurrentDim(theindex+1, NL, TwoSRdown, IRdown); in doD9andD10andD11()
1057 int theindex = denT->gIndex(); in doD12() local
1059 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD12()
1060 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD12()
1063 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD12()
1064 int dimRup = denBK->gCurrentDim(theindex+1, NL, TwoSLup, ILup); in doD12()
1071 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD12()
1072 int dimRdown = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, Idown); in doD12()
1112 int theindex = denT->gIndex(); in doD13() local
1114 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD13()
1115 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD13()
1118 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD13()
1119 int dimRup = denBK->gCurrentDim(theindex+1, NL+2, TwoSLup, ILup); in doD13()
1124 int IRdown = Irreps::directProd(ILdown, denBK->gIrrep(theindex)); in doD13()
1128 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown); in doD13()
1129 int dimRdown = denBK->gCurrentDim(theindex+1, NL, TwoSLup, IRdown); in doD13()
1168 int theindex = denT->gIndex(); in doD14() local
1170 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD14()
1171 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD14()
1174 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD14()
1179 int IRup = Irreps::directProd(ILup, denBK->gIrrep(theindex)); in doD14()
1183 int dimRup = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, IRup); in doD14()
1184 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD14()
1185 int dimRdown = denBK->gCurrentDim(theindex+1, NL-1, TwoSLdown, Idown); in doD14()
1224 int theindex = denT->gIndex(); in doD15() local
1226 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD15()
1227 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD15()
1230 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD15()
1231 int dimRup = denBK->gCurrentDim(theindex+1, NL+2, TwoSLup, ILup); in doD15()
1236 int IRdown = Irreps::directProd(ILdown, denBK->gIrrep(theindex)); in doD15()
1241 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown); in doD15()
1242 int dimRdown = denBK->gCurrentDim(theindex+1, NL, TwoSRdown, IRdown); in doD15()
1283 int theindex = denT->gIndex(); in doD16() local
1285 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD16()
1286 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD16()
1289 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD16()
1294 int IRup = Irreps::directProd(ILup, denBK->gIrrep(theindex)); in doD16()
1299 int dimRup = denBK->gCurrentDim(theindex+1, NL+1, TwoSRup, IRup); in doD16()
1300 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD16()
1301 int dimRdown = denBK->gCurrentDim(theindex+1, NL-1, TwoSLdown, Idown); in doD16()
1342 int theindex = denT->gIndex(); in doD17orD21() local
1344 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD17orD21()
1345 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD17orD21()
1348 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD17orD21()
1353 int IRdown = Irreps::directProd(ILdown, denBK->gIrrep(theindex)); in doD17orD21()
1357 int dimRup = denBK->gCurrentDim(theindex+1, NL, TwoSLup, ILup); in doD17orD21()
1358 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown); in doD17orD21()
1359 int dimRdown = denBK->gCurrentDim(theindex+1, NL, TwoSLup, IRdown); in doD17orD21()
1399 int theindex = denT->gIndex(); in doD18orD22() local
1401 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD18orD22()
1402 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD18orD22()
1405 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD18orD22()
1410 int IRup = Irreps::directProd(ILup, denBK->gIrrep(theindex)); in doD18orD22()
1414 int dimRup = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, IRup); in doD18orD22()
1415 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD18orD22()
1416 int dimRdown = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, Idown); in doD18orD22()
1458 int theindex = denT->gIndex(); in doD19orD23() local
1460 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD19orD23()
1461 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD19orD23()
1464 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD19orD23()
1469 int IRdown = Irreps::directProd(ILdown, denBK->gIrrep(theindex)); in doD19orD23()
1474 int dimRup = denBK->gCurrentDim(theindex+1, NL, TwoSLup, ILup); in doD19orD23()
1475 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown); in doD19orD23()
1476 int dimRdown = denBK->gCurrentDim(theindex+1, NL, TwoSRdown, IRdown); in doD19orD23()
1526 int theindex = denT->gIndex(); in doD20orD24() local
1528 for (int NL = denBK->gNmin(theindex); NL<=denBK->gNmax(theindex); NL++){ in doD20orD24()
1529 …for (int TwoSLup = denBK->gTwoSmin(theindex,NL); TwoSLup<= denBK->gTwoSmax(theindex,NL); TwoSLup+=… in doD20orD24()
1532 int dimLup = denBK->gCurrentDim(theindex, NL, TwoSLup, ILup); in doD20orD24()
1537 int IRup = Irreps::directProd(ILup, denBK->gIrrep(theindex)); in doD20orD24()
1542 int dimRup = denBK->gCurrentDim(theindex+1, NL+1, TwoSRup, IRup); in doD20orD24()
1543 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, Idown); in doD20orD24()
1544 int dimRdown = denBK->gCurrentDim(theindex+1, NL+1, TwoSLdown, Idown); in doD20orD24()