yij implemented for all single excitations

cthree-40 · cthree-40 · commit 27b025910ce9 · 2020-08-29T17:03:19.000-04:00
diff --git a/mpi_source/citruncate.c b/mpi_source/citruncate.c
@@ -809,9 +809,9 @@ void generate_actvx(int nrep, struct occstr str, int str_docc, int str_actv,
                 xlist[*numx].fo[0] = scr[j];
                 xlist[*numx].fo[1] = 0;
                 
-                xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
-                                                       scr[j], elec);
-                
+                //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
+                //                                       scr[j], elec);
+                xlist[*numx].permx = str.yij[i][scr[j]];
                 (*numx)++;
             }
         }
@@ -919,9 +919,10 @@ void generate_actv2virtx(int nrep, struct occstr str, int str_docc, int str_actv
                     xlist[*numx].fo[0] = scr[j];
                     xlist[*numx].fo[1] = 0;
                     
-                    xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
-                                                       scr[j], elec);;
-                                    
+//                    xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
+//                                                       scr[j], elec);;
+                    xlist[*numx].permx = str.yij[i][scr[j]];
+
                     (*numx)++;
                 }
             }
@@ -999,8 +1000,9 @@ void generate_actv2virtx(int nrep, struct occstr str, int str_docc, int str_actv
                     xlist[*numx].fo[0] = scr[i];
                     xlist[*numx].fo[1] = 0;
                     
-                    xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],
-                                                       scr[i], elec);;
+//                    xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],
+//                                                       scr[i], elec);;
+                    xlist[*numx].permx = str.yij[str_docc + str_actv + j][scr[i]];
                     (*numx)++;
                 }
             }
@@ -1698,8 +1700,10 @@ void generate_docc2actvx(int nrep, struct occstr str, int str_docc, int str_actv
                     xlist[*numx].fo[0] = scr[j];
                     xlist[*numx].fo[1] = 0;
                     
-                    xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
-                                                           scr[j], elec);;
+                    //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
+                    //                                       scr[j], elec);
+                    xlist[*numx].permx = str.yij[i][scr[j]];
+
                     (*numx)++;
                 }
             }
@@ -1766,6 +1770,7 @@ void generate_docc2actvx(int nrep, struct occstr str, int str_docc, int str_actv
 
                     xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
                                                            scr[j], elec);;
+                    xlist[*numx].permx = str.yij[i][scr[j]];
                     (*numx)++;
                 }
             }
@@ -3186,8 +3191,9 @@ void generate_docc2virtx(int nrep, struct occstr str, int str_docc, int str_actv
                     xlist[*numx].fo[0] = scr[j];
                     xlist[*numx].fo[1] = 0;
 
-                    xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
-                                                           scr[j], elec);;
+//                    xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
+//                                                           scr[j], elec);;
+                    xlist[*numx].permx = str.yij[i][scr[j]];
                     (*numx)++;
                 }
             }
@@ -3266,8 +3272,9 @@ void generate_docc2virtx(int nrep, struct occstr str, int str_docc, int str_actv
                     xlist[*numx].fo[0] = scr[i];
                     xlist[*numx].fo[1] = 0;
 
-                    xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],
-                                                           scr[i], elec);;
+//                    xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],
+//                                                           scr[i], elec);;
+                    xlist[*numx].permx = str.yij[str_docc + str_actv + j][scr[i]];
                     (*numx)++;
                 }
             }
@@ -3354,8 +3361,10 @@ void generate_doccx(int nrep, struct occstr str, int str_docc, int str_actv,
                 xlist[*numx].fo[0] = scr[j];
                 xlist[*numx].fo[1] = 0;
 
-                xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
-                                                       scr[j], elec);;
+                //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],
+                //                                       scr[j], elec);
+                
+                xlist[*numx].permx = str.yij[i][scr[j]];
                 (*numx)++;
             }
         }
@@ -3452,7 +3461,8 @@ void generate_virtx(int nrep, struct occstr str, int str_docc, int str_actv,
                 xlist[*numx].fo[1] = 0;
 
                 xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[i],
-                                                       scr[j], elec);;
+                                                       scr[j], elec);
+                xlist[*numx].permx = str.yij[str_docc + str_actv + i][scr[j]];
                 (*numx)++;
             }
         }

Original file line number	Diff line number	Diff line change
`@@ -809,9 +809,9 @@ void generate_actvx(int nrep, struct occstr str, int str_docc, int str_actv,`
`809`	`809`	`xlist[*numx].fo[0] = scr[j];`
`810`	`810`	`xlist[*numx].fo[1] = 0;`
`811`	`811`
`812`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`813`		`- scr[j], elec);`
`814`		`-`
	`812`	`+ //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
	`813`	`+ // scr[j], elec);`
	`814`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
`815`	`815`	`(*numx)++;`
`816`	`816`	`}`
`817`	`817`	`}`
`@@ -919,9 +919,10 @@ void generate_actv2virtx(int nrep, struct occstr str, int str_docc, int str_actv`
`919`	`919`	`xlist[*numx].fo[0] = scr[j];`
`920`	`920`	`xlist[*numx].fo[1] = 0;`
`921`	`921`
`922`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`923`		`- scr[j], elec);;`
`924`		`-`
	`922`	`+// xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
	`923`	`+// scr[j], elec);;`
	`924`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
	`925`	`+`
`925`	`926`	`(*numx)++;`
`926`	`927`	`}`
`927`	`928`	`}`
`@@ -999,8 +1000,9 @@ void generate_actv2virtx(int nrep, struct occstr str, int str_docc, int str_actv`
`999`	`1000`	`xlist[*numx].fo[0] = scr[i];`
`1000`	`1001`	`xlist[*numx].fo[1] = 0;`
`1001`	`1002`
`1002`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],`
`1003`		`- scr[i], elec);;`
	`1003`	`+// xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],`
	`1004`	`+// scr[i], elec);;`
	`1005`	`+ xlist[*numx].permx = str.yij[str_docc + str_actv + j][scr[i]];`
`1004`	`1006`	`(*numx)++;`
`1005`	`1007`	`}`
`1006`	`1008`	`}`
`@@ -1698,8 +1700,10 @@ void generate_docc2actvx(int nrep, struct occstr str, int str_docc, int str_actv`
`1698`	`1700`	`xlist[*numx].fo[0] = scr[j];`
`1699`	`1701`	`xlist[*numx].fo[1] = 0;`
`1700`	`1702`
`1701`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`1702`		`- scr[j], elec);;`
	`1703`	`+ //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
	`1704`	`+ // scr[j], elec);`
	`1705`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
	`1706`	`+`
`1703`	`1707`	`(*numx)++;`
`1704`	`1708`	`}`
`1705`	`1709`	`}`
`@@ -1766,6 +1770,7 @@ void generate_docc2actvx(int nrep, struct occstr str, int str_docc, int str_actv`
`1766`	`1770`
`1767`	`1771`	`xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`1768`	`1772`	`scr[j], elec);;`
	`1773`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
`1769`	`1774`	`(*numx)++;`
`1770`	`1775`	`}`
`1771`	`1776`	`}`
`@@ -3186,8 +3191,9 @@ void generate_docc2virtx(int nrep, struct occstr str, int str_docc, int str_actv`
`3186`	`3191`	`xlist[*numx].fo[0] = scr[j];`
`3187`	`3192`	`xlist[*numx].fo[1] = 0;`
`3188`	`3193`
`3189`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`3190`		`- scr[j], elec);;`
	`3194`	`+// xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
	`3195`	`+// scr[j], elec);;`
	`3196`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
`3191`	`3197`	`(*numx)++;`
`3192`	`3198`	`}`
`3193`	`3199`	`}`
`@@ -3266,8 +3272,9 @@ void generate_docc2virtx(int nrep, struct occstr str, int str_docc, int str_actv`
`3266`	`3272`	`xlist[*numx].fo[0] = scr[i];`
`3267`	`3273`	`xlist[*numx].fo[1] = 0;`
`3268`	`3274`
`3269`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],`
`3270`		`- scr[i], elec);;`
	`3275`	`+// xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[j],`
	`3276`	`+// scr[i], elec);;`
	`3277`	`+ xlist[*numx].permx = str.yij[str_docc + str_actv + j][scr[i]];`
`3271`	`3278`	`(*numx)++;`
`3272`	`3279`	`}`
`3273`	`3280`	`}`
`@@ -3354,8 +3361,10 @@ void generate_doccx(int nrep, struct occstr str, int str_docc, int str_actv,`
`3354`	`3361`	`xlist[*numx].fo[0] = scr[j];`
`3355`	`3362`	`xlist[*numx].fo[1] = 0;`
`3356`	`3363`
`3357`		`- xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
`3358`		`- scr[j], elec);;`
	`3364`	`+ //xlist[*numx].permx = pindex_single_rep(str.istr, str.istr[i],`
	`3365`	`+ // scr[j], elec);`
	`3366`	`+`
	`3367`	`+ xlist[*numx].permx = str.yij[i][scr[j]];`
`3359`	`3368`	`(*numx)++;`
`3360`	`3369`	`}`
`3361`	`3370`	`}`
`@@ -3452,7 +3461,8 @@ void generate_virtx(int nrep, struct occstr str, int str_docc, int str_actv,`
`3452`	`3461`	`xlist[*numx].fo[1] = 0;`
`3453`	`3462`
`3454`	`3463`	`xlist[*numx].permx = pindex_single_rep(str.istr, str.virtx[i],`
`3455`		`- scr[j], elec);;`
	`3464`	`+ scr[j], elec);`
	`3465`	`+ xlist[*numx].permx = str.yij[str_docc + str_actv + i][scr[j]];`
`3456`	`3466`	`(*numx)++;`
`3457`	`3467`	`}`
`3458`	`3468`	`}`