portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-06-19 16:27:05 +01:00
Code Releases Activity

block Lanczos construction is added.

Browse Source
This commit is contained in:
Yong-Chull Jang
2017-12-03 23:55:22 -05:00
parent 2c35c89b92
commit 5139eaf491
3 changed files with 711 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

17
tests/lanczos/Test_dwf_block_lanczos.cc
View File

@ -2,7 +2,7 @@
Grid physics library, www.github.com/paboyle/Grid
Source file: ./tests/Test_dwf_lanczos.cc
Source file: ./tests/Test_dwf_block_lanczos.cc
Copyright (C) 2015
@ -78,6 +78,7 @@ int main (int argc, char ** argv)
// SchurDiagMooeeOperator<DomainWallFermionR,LatticeFermion> HermOp(Ddwf);
const int Nstop = 30;
const int Nu = 4;
const int Nk = 60;
const int Np = 60;
const int Nm = Nk+Np;
@ -92,19 +93,21 @@ int main (int argc, char ** argv)
// ChebyshevLanczos<LatticeFermion> Cheb(9.,1.,0.,20);
// Cheb.csv(std::cout);
// exit(-24);
ImplicitlyRestartedLanczos<FermionField> IRL(HermOp,Cheb,Nstop,Nk,Nm,resid,MaxIt);
ImplicitlyRestartedBlockLanczos<FermionField> IRBL(HermOp,Cheb,Nstop,Nu,Nk,Nm,resid,MaxIt);
std::vector<RealD> eval(Nm);
FermionField src(FrbGrid);
gaussian(RNG5rb,src);
std::vector<FermionField> src(Nu,FrbGrid);
for ( int i=0; i<Nu; ++i ) gaussian(RNG5rb,src[i]);
std::vector<FermionField> evec(Nm,FrbGrid);
for(int i=0;i<1;i++){
std::cout << GridLogMessage <<i<<" / "<< Nm<< " grid pointer "<<evec[i]._grid<<std::endl;
for(int i=0;i<1;++i){
clog << i <<" / "<< Nm <<" grid pointer "<< evec[i]._grid << std::endl;
};
int Nconv;
IRL.calc(eval,evec,src,Nconv);
IRBL.calc(eval,evec,src,Nconv);
Grid_finalize();