portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-04-09 21:50:45 +01:00
Code Releases Activity

-Checkerboard fixes for Lanczos

Browse Source
This commit is contained in:
Christopher Kelly 2016-04-06 13:50:56 -04:00
parent 60d965f79e
commit af9c8d1372
1 changed files with 2 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
lib/algorithms/iterative/ImplicitlyRestartedLanczos.h
View File

@ -564,6 +564,7 @@ until convergence
for(int j=k1-1; j<k2+1; ++j){ for(int j=k1-1; j<k2+1; ++j){
for(int k=0; k<Nm; ++k){ for(int k=0; k<Nm; ++k){
B[j].checkerboard = evec[k].checkerboard;
B[j] += Qt[k+Nm*j] * evec[k]; B[j] += Qt[k+Nm*j] * evec[k];
} }
} }
@ -592,6 +593,7 @@ until convergence
for(int j = 0; j<Nk; ++j){ for(int j = 0; j<Nk; ++j){
for(int k = 0; k<Nk; ++k){ for(int k = 0; k<Nk; ++k){
B[j].checkerboard = evec[k].checkerboard;
B[j] += Qt[k+j*Nm] * evec[k]; B[j] += Qt[k+j*Nm] * evec[k];
} }
// std::cout << "norm(B["<<j<<"])="<<norm2(B[j])<<std::endl; // std::cout << "norm(B["<<j<<"])="<<norm2(B[j])<<std::endl;