portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-07-18 08:03:27 +01:00
Code Releases Activity

Adding Lanczos parameter file in examples directory

Browse Source
This commit is contained in:
Chulwoo Jung
2026-04-06 13:35:43 -04:00
parent d4b685cf2d
commit 82a8f1342e
2 changed files with 22 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
examples/Example_krylov_schur.cc
+2 -1
View File
@@ -351,6 +351,7 @@ int main (int argc, char ** argv)
KrylovSchur KrySchur (Dwilson, UGrid, resid,EvalImNormSmall); KrylovSchur KrySchur (Dwilson, UGrid, resid,EvalImNormSmall);
KrySchur(src[0], maxIter, Nm, Nk, Nstop); KrySchur(src[0], maxIter, Nm, Nk, Nstop);
// KrySchur(src[0], maxIter, Nm, Nk, Nstop,&shift); // KrySchur(src[0], maxIter, Nm, Nk, Nstop,&shift);
std::cout << GridLogMessage << "KrylovSchur evec.size= " << KrySchur.evecs.size()<< std::endl;
#else #else
int Nblock=4; int Nblock=4;
Nblock=LanParams.Nblock; Nblock=LanParams.Nblock;
@@ -362,8 +363,8 @@ int main (int argc, char ** argv)
KrySchur(src, maxIter, Nm, Nk, Nstop,Nblock,true,if_verify); KrySchur(src, maxIter, Nm, Nk, Nstop,Nblock,true,if_verify);
// HarmonicBlockKrylovSchur KrySchur (Dwilson, UGrid, resid,shift,EvalImNormSmall); // HarmonicBlockKrylovSchur KrySchur (Dwilson, UGrid, resid,shift,EvalImNormSmall);
// KrySchur(src, maxIter, Nm, Nk, Nstop,Nblock,true); // KrySchur(src, maxIter, Nm, Nk, Nstop,Nblock,true);
std::cout << GridLogMessage << "BlockKrylovSchur evec.size= " << KrySchur.evecs.size()<< std::endl;
#endif #endif
std::cout << GridLogMessage << "evec.size= " << KrySchur.evecs.size()<< std::endl;
src[0]=KrySchur.evecs[0]; src[0]=KrySchur.evecs[0];
std::cout << GridLogMessage << "KrySchur.evecs= "<< KrySchur.evecs.size() <<std::endl; std::cout << GridLogMessage << "KrySchur.evecs= "<< KrySchur.evecs.size() <<std::endl;
examples/LanParams.xml
+20
View File
@@ -0,0 +1,20 @@
<?xml version="1.0"?>
<grid>
<LanczosParameters>
<mass>0</mass>
<mstep>-0.025</mstep>
<M5>1.8</M5>
<Ls>48</Ls>
<Nstop>80</Nstop>
<Nk>100</Nk>
<Np>100</Np>
<ReadEvec>0</ReadEvec>
<maxIter>1000</maxIter>
<Nblock>4</Nblock>
<verify>0</verify>
<resid>1e-10</resid>
<ChebyLow>1</ChebyLow>
<ChebyHigh>100</ChebyHigh>
<ChebyOrder>51</ChebyOrder>
</LanczosParameters>
</grid>