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

Slight clean up

Browse Source
This commit is contained in:
Peter Boyle 2023-10-06 10:47:32 -04:00
parent e7020017c5
commit 2525ad4623
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
Grid/algorithms/iterative/AdefGeneric.h
View File

@ -161,12 +161,12 @@ class TwoLevelCG : public LinearFunction<Field>
axpy(tmp,-1.0,src,mmp); axpy(tmp,-1.0,src,mmp);
RealD mmpnorm = sqrt(norm2(mmp)); RealD mmpnorm = sqrt(norm2(mmp));
RealD psinorm = sqrt(norm2(x)); RealD xnorm = sqrt(norm2(x));
RealD srcnorm = sqrt(norm2(src)); RealD srcnorm = sqrt(norm2(src));
RealD tmpnorm = sqrt(norm2(tmp)); RealD tmpnorm = sqrt(norm2(tmp));
RealD true_residual = tmpnorm/srcnorm; RealD true_residual = tmpnorm/srcnorm;
std::cout<<GridLogMessage<<"HDCG: true residual is "<<true_residual std::cout<<GridLogMessage<<"HDCG: true residual is "<<true_residual
<<" solution "<<psinorm<<" source "<<srcnorm<<std::endl; <<" solution "<<xnorm<<" source "<<srcnorm<<std::endl;
return; return;
} }
@ -326,7 +326,7 @@ public:
// Can just inherit existing M3 // Can just inherit existing M3
// Simple vstart - do nothing // Simple vstart - do nothing
virtual void Vstart(Field & x,const Field & src){}; virtual void Vstart(Field & x,const Field & src){ x=src; };
// Override PcgM1 // Override PcgM1
virtual void PcgM1(Field & in, Field & out) virtual void PcgM1(Field & in, Field & out)