modify diago_cg_test.cpp

root · root · commit 5199f78e4636 · 2022-03-25T15:54:50.000+08:00
diff --git a/source/src_pw/test/diago_cg_test.cpp b/source/src_pw/test/diago_cg_test.cpp
@@ -75,8 +75,9 @@ class DiagoCGPrepare
         {
             for (int j = 0; j < npw; j++)
             {
-                // psiguess(i,j) = ev(j,i)*(1+u(p)/10.);
-                psiguess(i, j) = ev(j, i) * u(p) / 10.;
+		double rand = static_cast<double>(u(p))/10.;
+                // psiguess(i,j) = ev(j,i)*(1+rand);
+                psiguess(i, j) = ev(j, i) * rand;
             }
         }
         // run cg

Original file line number	Diff line number	Diff line change
`@@ -75,8 +75,9 @@ class DiagoCGPrepare`
`75`	`75`	`{`
`76`	`76`	`for (int j = 0; j < npw; j++)`
`77`	`77`	`{`
`78`		`- // psiguess(i,j) = ev(j,i)*(1+u(p)/10.);`
`79`		`- psiguess(i, j) = ev(j, i) * u(p) / 10.;`
	`78`	`+ double rand = static_cast<double>(u(p))/10.;`
	`79`	`+ // psiguess(i,j) = ev(j,i)*(1+rand);`
	`80`	`+ psiguess(i, j) = ev(j, i) * rand;`
`80`	`81`	`}`
`81`	`82`	`}`
`82`	`83`	`// run cg`