Lines Matching refs:theindex
46 int theindex = denS->gIndex(); local
47 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
48 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
54 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
76 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
80 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
86 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
92 … int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
93 … int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
100 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
102 … double prefac = factor * ( Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
103 … + ((TwoJdown==0)?1:-1) * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex) );
104 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
116 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
156 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
160 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
166 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
172 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
173 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
180 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
182 … double prefac = factor1 * Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
183 … + factor2 * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex);
184 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
196 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
235 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
239 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
245 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
251 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
252 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
259 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
261 … double prefac = factor1 * Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
262 … + factor2 * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex);
263 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
275 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
307 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
311 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
317 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
323 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
324 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
331 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
333 … double prefac = factor * ( Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
334 … + ((TwoJ==0)?1:-1) * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex) );
335 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
347 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
378 int theindex = denS->gIndex(); local
379 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
380 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
385 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
406 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
410 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
416 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
422 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
423 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
430 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
432 … double prefac = factor * ( Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
433 … + ((TwoJ==0)?1:-1) * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex) );
434 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
446 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
483 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
487 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
493 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
499 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
500 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
507 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
509 … double prefac = factor1 * Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
510 … + factor2 * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex);
511 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
523 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
561 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
565 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
571 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
577 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
578 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
585 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
587 … double prefac = factor1 * Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
588 … + factor2 * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex);
589 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
601 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
634 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
638 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
644 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
650 … int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
651 … int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
658 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
660 … double prefac = factor * ( Prob->gMxElement(l_alpha, l_beta, theindex, theindex+1)
661 … + ((TwoJdown==0)?1:-1) * Prob->gMxElement(l_alpha, l_beta, theindex+1, theindex) );
662 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
674 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
707 int theindex = denS->gIndex(); local
708 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
709 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
715 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
737 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
741 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
747 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
753 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
754 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
761 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
763 … double prefac = factor1 * Prob->gMxElement(l_alpha,theindex,theindex+1,l_beta)
764 … + factor2 * Prob->gMxElement(l_alpha,theindex,l_beta,theindex+1);
765 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
777 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
810 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
814 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
820 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
826 … int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
827 … int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
834 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
836 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
837 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
838 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
850 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
883 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
887 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
893 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
899 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
900 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
907 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
909 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
910 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
911 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
923 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
953 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
957 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
963 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
969 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
970 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
977 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
979 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
980 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
981 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
993 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
1025 int theindex = denS->gIndex(); local
1026 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
1027 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
1033 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
1055 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1059 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1065 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1071 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1072 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1079 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1081 … double prefac = factor1 * Prob->gMxElement(l_alpha,theindex,theindex+1,l_beta)
1082 … + factor2 * Prob->gMxElement(l_alpha,theindex,l_beta,theindex+1);
1083 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1095 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1128 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1132 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1138 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1144 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1145 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1152 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1154 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
1155 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
1156 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1168 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1200 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1204 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1210 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1216 … int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1217 … int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1224 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1226 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
1227 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
1228 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1240 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1272 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1276 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1282 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1288 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1289 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1296 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1298 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex, theindex+1, l_beta)
1299 … + factor2 * Prob->gMxElement(l_alpha, theindex, l_beta, theindex+1);
1300 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1312 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1344 int theindex = denS->gIndex(); local
1345 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
1346 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
1352 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
1373 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1377 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1383 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1389 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
1390 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
1397 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1399 … double prefac = factor1 * Prob->gMxElement(l_alpha,theindex+1,theindex,l_beta)
1400 … + factor2 * Prob->gMxElement(l_alpha,theindex+1,l_beta,theindex);
1401 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
1413 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
1446 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1450 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1456 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1462 … int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
1463 … int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
1470 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1472 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1473 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1474 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
1486 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
1519 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1523 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1529 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1535 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
1536 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
1543 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1545 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1546 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1547 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
1559 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
1589 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1593 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1599 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1605 int dimLdown = denBK->gCurrentDim(theindex, NL-1, TwoSLdown, ILdown);
1606 int dimRdown = denBK->gCurrentDim(theindex+2, NR-1, TwoSRdown, IRdown);
1613 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1615 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1616 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1617 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR-1,TwoSRdown,IRdown,NR,TwoSR,IR);
1629 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL-1,TwoSLdown,ILdown,NL,TwoSL,IL);
1661 int theindex = denS->gIndex(); local
1662 int dimLup = denBK->gCurrentDim(theindex ,NL,TwoSL,IL);
1663 int dimRup = denBK->gCurrentDim(theindex+2,NR,TwoSR,IR);
1669 int IprodMID = Irreps::directProd(denBK->gIrrep(theindex),denBK->gIrrep(theindex+1));
1690 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1694 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1700 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1706 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1707 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1714 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1716 … double prefac = factor1 * Prob->gMxElement(l_alpha,theindex+1,theindex,l_beta)
1717 … + factor2 * Prob->gMxElement(l_alpha,theindex+1,l_beta,theindex);
1718 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1730 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1761 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1765 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1771 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1777 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1778 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1785 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1787 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1788 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1789 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1801 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1832 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1836 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1842 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1848 … int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1849 … int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1856 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1858 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1859 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1860 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1872 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);
1904 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1908 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1914 for (int l_alpha=0; l_alpha<theindex; l_alpha++){
1920 int dimLdown = denBK->gCurrentDim(theindex, NL+1, TwoSLdown, ILdown);
1921 int dimRdown = denBK->gCurrentDim(theindex+2, NR+1, TwoSRdown, IRdown);
1928 for (int l_beta=theindex+2; l_beta<Prob->gL(); l_beta++){
1930 … double prefac = factor1 * Prob->gMxElement(l_alpha, theindex+1, theindex, l_beta)
1931 … + factor2 * Prob->gMxElement(l_alpha, theindex+1, l_beta, theindex);
1932 … double * LblockR = Lright[l_beta-theindex-2]->gStorage(NR,TwoSR,IR,NR+1,TwoSRdown,IRdown);
1944 … double * LblockL = Lleft[theindex-1-l_alpha]->gStorage(NL,TwoSL,IL,NL+1,TwoSLdown,ILdown);