/dports/science/chemps2/CheMPS2-1.8.10/CheMPS2/ |
H A D | TensorL.cpp | 75 int dimRup = bk_up ->gCurrentDim( index, NRup, TwoSRup, IRup ); in create_right() local 110 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, Tdown, &dimLup, &add, storag… in create_right() 127 …dgemm_( &trans, ¬rans, &dimRup, &dimLdown, &dimLup, &alpha, Tup, &dimLup, Opart, &dimLup, &set,… in create_right() 129 …rans, ¬rans, &dimRup, &dimRdown, &dimLdown, &one, workmem, &dimRup, Tdown, &dimLdown, &one, sto… in create_right() 164 int dimRup = bk_up ->gCurrentDim( index + 1, NR, TwoSR, IR ); in create_left() local 168 assert( dimRup == dimRdown ); in create_left() 169 if ( dimRup > 0 ){ in create_left() 181 …dgemm_( ¬rans, &trans, &dimLup, &dimLdown, &dimRup, &alpha, Tup, &dimLup, Tdown, &dimLdown, &ad… in create_left() 185 if (( dimRup > 0 ) && ( dimRdown > 0 )){ in create_left() 198 …dgemm_( ¬rans, ¬rans, &dimLup, &dimRdown, &dimRup, &alpha, Tup, &dimLup, Opart, &dimRup, &se… in create_left()
|
H A D | Tensor3RDM.cpp | 93 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in a1() 120 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in a1() 169 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in b1() 196 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in b1() 245 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in c1() 272 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in c1() 321 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in d1() 385 …mm_( &trans, ¬rans, &dimRup, &dimRdown, &dimL, &alpha, Tup, &dimL, Tdown, &dimL, &beta, storage… in extra1() 424 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in extra2() 463 … &trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, workmem, &dimLup, &beta, sto… in extra3() [all …]
|
H A D | HeffDiagrams5.cpp | 48 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); local 97 int size = dimRup * dimRdown; 177 int size = dimRup * dimRdown; 256 int size = dimRup * dimRdown; 328 int size = dimRup * dimRdown; 380 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); local 427 int size = dimRup * dimRdown; 904 int size = dimRup * dimRdown; 1149 int size = dimRup * dimRdown; 1540 int size = dimRup * dimRdown; [all …]
|
H A D | ThreeDM.cpp | 1609 if ( dimRup > 0 ){ in diagram10() 1728 if ( dimRup > 0 ){ in diagram12() 1846 if ( dimRup > 0 ){ in diagram14() 2025 if ( dimRup > 0 ){ in diagram17() 2144 if ( dimRup > 0 ){ in diagram19() 2262 if ( dimRup > 0 ){ in diagram21() 3706 if (( dimLup > 0 ) && ( dimRup > 0 )){ in fill_tens_46_48() 3757 if (( dimLup > 0 ) && ( dimRup > 0 )){ in fill_53_54() 3803 if (( dimLup > 0 ) && ( dimRup > 0 )){ in fill_55_to_60() 3849 if (( dimLup > 0 ) && ( dimRup > 0 )){ in fill_61() [all …]
|
H A D | HeffDiagrams4.cpp | 691 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram4B1and4B2spin0() local 891 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram4B1and4B2spin1() local 1103 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram4B3and4B4spin0() local 3071 int size = dimRup * dimRdown; 3114 int size = dimRup * dimRdown; 3162 int size = dimRup * dimRdown; 3238 int size = dimRup * dimRdown; 3339 int size = dimRup * dimRdown; 3382 int size = dimRup * dimRdown; 3431 int size = dimRup * dimRdown; [all …]
|
H A D | Excitation.cpp | 129 assert( dimRup == dimRdown ); in neighbours() 133 int size = dimLup * dimRup; in neighbours() 340 …dgemm_( ¬rans, &trans, &dimLup, &dimRup, &dimRdown, &factor, block_down, &dimLup, block_right, … in second_left() 375 …dgemm_( ¬rans, &trans, &dimLup, &dimRup, &dimRdown, &one, block_down, &dimLup, block_right, &di… in third_left() 378 int size = dimLup * dimRup; in third_left() 409 assert( dimRup == dimRdown ); in first_right() 482 assert( dimRup == dimRdown ); in second_right() 551 assert( dimRup == dimRdown ); in third_right() 563 int size = dimLup * dimRup; in third_right() 707 …dgemm_( ¬rans, &trans, &dimLup, &dimRup, &dimRdown, &one, workmem1, &dimLup, block_right,… in third_middle() [all …]
|
H A D | TwoDM.cpp | 697 …dgemm_(¬rans,&trans,&dimL,&dimRup,&dimRdown,&alpha,Tdown,&dimL,Lblock,&dimRup,&beta,workmem,&di… in doD2() 701 int length = dimL * dimRup; in doD2() 742 int length = dimL * dimRup; in doD3() 866 int length = dimL * dimRup; in doD6() 940 if ((dimLup>0) && (dimRup>0)){ in doD8() 966 int length = dimLup * dimRup; in doD8() 1065 if ((dimLup>0) && (dimRup>0)){ in doD12() 1093 int length = dimLup * dimRup; in doD12() 1121 if ((dimLup>0) && (dimRup>0)){ in doD13() 1149 int length = dimLup * dimRup; in doD13() [all …]
|
H A D | TensorO.cpp | 100 int dimRup = bk_up->gCurrentDim( index, NR, TwoSR, IR ); in create_right() local 138 …&trans, ¬rans, &dimRup, &dimRdown, &dimLup, &alpha, Tup, &dimLup, Tdown, &dimLdown, &beta, stor… in create_right() 179 int dimRup = bk_up->gCurrentDim( index + 1, NR, TwoSR, IR ); in create_left() local 181 if (( dimRup > 0 ) && ( dimRdown > 0 ) && ( dimRup == dimRdown )){ in create_left() 189 …dgemm_( ¬rans, &trans, &dimLup, &dimLdown, &dimRup, &alpha, Tup, &dimLup, Tdown, &dimLdown, &be… in create_left()
|
H A D | HeffDiagrams3.cpp | 312 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram3C() local 347 …dgemm_(¬ra,&trans,&dimLup,&dimRup,&dimRdown,&alpha,temp,&dimLup,Lblock,&dimRup,&beta,memHeff+de… in addDiagram3C() 493 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram3Kand3F() local 531 int size = dimRup * dimRdown; in addDiagram3Kand3F() 590 int size = dimRup * dimRdown; in addDiagram3Kand3F() 601 …dgemm_(¬r,&tran,&dimL,&dimRup,&dimRdown,&factor,memS+denS->gKappa2index(memSkappa),&dimL,temp,&… in addDiagram3Kand3F() 626 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram3Land3G() local 664 int size = dimRup * dimRdown; in addDiagram3Land3G() 723 int size = dimRup * dimRdown; in addDiagram3Land3G() 760 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram3J() local [all …]
|
H A D | TensorX.cpp | 312 int dimRup = bk_up->gCurrentDim(index+1, NRup, TwoSRup, IRup); in addTermQLLeft() local 315 if ((dimRup>0) && (dimRdown>0)){ in addTermQLLeft() 332 int dimRupdown = dimRup * dimRdown; in addTermQLLeft() 350 …dgemm_(¬r,¬r,&dimL,&dimRdown,&dimRup,&factor,BlockTup,&dimL,ptr,&dimRup,&beta,workmemLR,&dim… in addTermQLLeft() 415 …int dimRup = bk_up->gCurrentDim(index+1, sector_nelec_up[ikappa], sector_spin_up[ikappa], sect… in addTermALeft() local 418 if ((dimRup>0) && (dimRdown>0)){ in addTermALeft() 428 …dgemm_(¬r,¬r,&dimL,&dimRdown,&dimRup,&factor,BlockTup,&dimL,BlockA,&dimRup,&beta,workmemLR,&… in addTermALeft() 619 int dimRup = bk_up->gCurrentDim(index+1,NR,TwoSRup, IR); in addTermDLeft() local 622 if ((dimRup>0) && (dimRdown>0)){ in addTermDLeft() 634 …dgemm_(¬r, &trans, &dimL, &dimRdown, &dimRup, &factor, BlockTup, &dimL, BlockD, &dimRdown, &bet… in addTermDLeft()
|
H A D | HeffDiagrams2.cpp | 1076 int dimRup = denBK->gCurrentDim(theindex+2,NR, TwoSR,IR); in addDiagram2e1and2e2() local 1087 …dgemm_(¬rans,¬rans,&dimL,&dimRup,&dimRdown,&alpha,memS+denS->gKappa2index(memSkappa),&dimL,B… in addDiagram2e1and2e2() 1104 int dimRup = denBK->gCurrentDim(theindex+2,NR, TwoSR,IR); in addDiagram2e1and2e2() local 1116 …dgemm_(¬rans,&trans,&dimL,&dimRup,&dimRdown,&alpha,memS+denS->gKappa2index(memSkappa),&dimL,Blo… in addDiagram2e1and2e2() 1139 int dimRup = denBK->gCurrentDim(theindex+2,NR, TwoSR,IR); in addDiagram2f1and2f2() local 1150 …dgemm_(¬rans,¬rans,&dimL,&dimRup,&dimRdown,&alpha,memS+denS->gKappa2index(memSkappa),&dimL,B… in addDiagram2f1and2f2() 1167 int dimRup = denBK->gCurrentDim(theindex+2,NR, TwoSR,IR); in addDiagram2f1and2f2() local 1179 …dgemm_(¬rans,&trans,&dimL,&dimRup,&dimRdown,&alpha,memS+denS->gKappa2index(memSkappa),&dimL,Blo… in addDiagram2f1and2f2() 1422 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR); in addDiagram2e3spin1() local 1446 …dgemm_(¬r,¬r,&dimL,&dimRup,&dimRdown,&alpha,memS+denS->gKappa2index(memSkappa),&dimL,Dblock,… in addDiagram2e3spin1() [all …]
|
H A D | DMRGtechnics.cpp | 574 int dimRup = Exc_BKs[state_number]->gCurrentDim(index+2,NR,TwoSR,IR); in calcVeffTilde() local 581 int dimBlock = dimLup * dimRup; in calcVeffTilde() 591 …dgemm_(¬rans,¬rans,&dimLdown,&dimRup,&dimLup,&alpha,Opart,&dimLdown,SupPart,&dimLup,&beta,wo… in calcVeffTilde() 610 …dgemm_(¬rans,&trans,&dimLdown,&dimRdown,&dimRup,&alpha,workmem,&dimLdown,Opart,&dimRdown,&beta,… in calcVeffTilde()
|