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

Messagign

Browse Source
This commit is contained in:
paboyle 2017-10-25 23:53:19 +01:00
parent 3d63b4894e
commit e4d461cb03
1 changed files with 7 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
tests/lanczos/Test_dwf_compressed_lanczos.cc
View File

@ -26,9 +26,9 @@
// The following are now decoupled from the Lanczos and deal with grids. // The following are now decoupled from the Lanczos and deal with grids.
// Safe to replace functionality // Safe to replace functionality
///////////////////////////////////////////////////////////////////////////// /////////////////////////////////////////////////////////////////////////////
#include <Grid/algorithms/iterative/BlockImplicitlyRestartedLanczos/BlockedGrid.h> #include "BlockedGrid.h"
#include <Grid/algorithms/iterative/BlockImplicitlyRestartedLanczos/FieldBasisVector.h> #include "FieldBasisVector.h"
#include <Grid/algorithms/iterative/BlockImplicitlyRestartedLanczos/BlockProjector.h> #include "BlockProjector.h"
#include "FieldVectorIO.h" #include "FieldVectorIO.h"
#include "Params.h" #include "Params.h"
@ -431,6 +431,7 @@ void CoarseGridLanczos(BlockProjector<Field>& pr,RealD alpha2,RealD beta,int Npo
auto result = src_orig; auto result = src_orig;
// undeflated solve // undeflated solve
std::cout << GridLogMessage << " Undeflated solve "<<std::endl;
result = zero; result = zero;
CG(HermOp, src_orig, result); CG(HermOp, src_orig, result);
// if (UCoarseGrid->IsBoss()) // if (UCoarseGrid->IsBoss())
@ -438,6 +439,7 @@ void CoarseGridLanczos(BlockProjector<Field>& pr,RealD alpha2,RealD beta,int Npo
// CG.ResHistory.clear(); // CG.ResHistory.clear();
// deflated solve with all eigenvectors // deflated solve with all eigenvectors
std::cout << GridLogMessage << " Deflated solve with all evectors"<<std::endl;
result = zero; result = zero;
pr.deflate(coef,eval2,Nstop2,src_orig,result); pr.deflate(coef,eval2,Nstop2,src_orig,result);
CG(HermOp, src_orig, result); CG(HermOp, src_orig, result);
@ -446,6 +448,7 @@ void CoarseGridLanczos(BlockProjector<Field>& pr,RealD alpha2,RealD beta,int Npo
// CG.ResHistory.clear(); // CG.ResHistory.clear();
// deflated solve with non-blocked eigenvectors // deflated solve with non-blocked eigenvectors
std::cout << GridLogMessage << " Deflated solve with non-blocked evectors"<<std::endl;
result = zero; result = zero;
pr.deflate(coef,eval1,Nstop1,src_orig,result); pr.deflate(coef,eval1,Nstop1,src_orig,result);
CG(HermOp, src_orig, result); CG(HermOp, src_orig, result);
@ -454,6 +457,7 @@ void CoarseGridLanczos(BlockProjector<Field>& pr,RealD alpha2,RealD beta,int Npo
// CG.ResHistory.clear(); // CG.ResHistory.clear();
// deflated solve with all eigenvectors and original eigenvalues from proj // deflated solve with all eigenvectors and original eigenvalues from proj
std::cout << GridLogMessage << " Deflated solve with all eigenvectors and original eigenvalues from proj"<<std::endl;
result = zero; result = zero;
pr.deflate(coef,eval3,Nstop2,src_orig,result); pr.deflate(coef,eval3,Nstop2,src_orig,result);
CG(HermOp, src_orig, result); CG(HermOp, src_orig, result);