portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2026-03-26 14:06:10 +00:00
Code Releases Activity

Reducing comments

Browse Source
This commit is contained in:
Chulwoo Jung
2025-12-09 14:23:22 +00:00
parent 88611659a3
commit df4c2a082b
2 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
Grid/algorithms/iterative/KrylovSchur.h
View File

@@ -349,7 +349,7 @@ class KrylovSchur {
startVec = u; // original code
start = Nk;
std::cout << GridLogMessage << "b after Arnoldi " << b << std::endl;
std::cout << GridLogDebug << "b after Arnoldi " << b << std::endl;
// checkKSDecomposition();
@@ -381,7 +381,7 @@ if(1){
for (int m=0;m<Nm;m++) temp(m,m) -= *shift;
Eigen::MatrixXcd RayleighS = temp.inverse(); // (B-tI)^-1
Eigen::MatrixXcd temp2 = RayleighS*temp;
std::cout << GridLogMessage << "Shift inverse check: shift= "<<*shift<<" "<< temp2 <<std::endl;
std::cout << GridLogDebug << "Shift inverse check: shift= "<<*shift<<" "<< temp2 <<std::endl;
temp2=RayleighS.adjoint(); //(B-tI)^-1*
Eigen::VectorXcd g = temp2*b; //g = (B-tI)^-1* * b

2
examples/Example_krylov_schur.cc
View File

@@ -540,7 +540,7 @@ int main (int argc, char ** argv)
// Hacked, really EvalImagSmall
KrylovSchur KrySchur (Dwilson, UGrid, resid,EvalImNormSmall);
// BlockKrylovSchur KrySchur (HermOp2, UGrid, Nu, resid,EvalNormSmall);
RealD shift=0.;
RealD shift=1.;
KrySchur(src[0], maxIter, Nm, Nk, Nstop,&shift);
std::cout << GridLogMessage << "evec.size= " << KrySchur.evecs.size()<< std::endl;