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

WilsonMG: Adapt staggered GMRES/MR tests to "new" constructor

Browse Source
This commit is contained in:
Daniel Richtmann 2018-06-18 16:18:56 +02:00
parent ddcb53bce2
commit cc5d025ea4
No known key found for this signature in database
GPG Key ID: B33C490AF0772057
5 changed files with 20 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
tests/solver/Test_staggered_cagmres_unprec.cc
View File

@ -59,7 +59,10 @@ int main (int argc, char ** argv)
} }
RealD mass=0.5; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); RealD c1=9.0/8.0;
RealD c2=-1.0/24.0;
RealD u0=1.0;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass,c1,c2,u0);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
CommunicationAvoidingGeneralisedMinimalResidual<FermionField> CAGMRES(1.0e-8, 10000, 25); CommunicationAvoidingGeneralisedMinimalResidual<FermionField> CAGMRES(1.0e-8, 10000, 25);

5
tests/solver/Test_staggered_fcagmres_prec.cc
View File

@ -59,7 +59,10 @@ int main (int argc, char ** argv)
} }
RealD mass=0.5; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); RealD c1=9.0/8.0;
RealD c2=-1.0/24.0;
RealD u0=1.0;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass,c1,c2,u0);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);

5
tests/solver/Test_staggered_fgmres_prec.cc
View File

@ -59,7 +59,10 @@ int main (int argc, char ** argv)
} }
RealD mass=0.5; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); RealD c1=9.0/8.0;
RealD c2=-1.0/24.0;
RealD u0=1.0;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass,c1,c2,u0);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);

5
tests/solver/Test_staggered_gmres_unprec.cc
View File

@ -59,7 +59,10 @@ int main (int argc, char ** argv)
} }
RealD mass=0.5; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); RealD c1=9.0/8.0;
RealD c2=-1.0/24.0;
RealD u0=1.0;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass,c1,c2,u0);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
GeneralisedMinimalResidual<FermionField> GMRES(1.0e-8, 10000, 25); GeneralisedMinimalResidual<FermionField> GMRES(1.0e-8, 10000, 25);

5
tests/solver/Test_staggered_mr_unprec.cc
View File

@ -59,7 +59,10 @@ int main (int argc, char ** argv)
} }
RealD mass=0.5; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); RealD c1=9.0/8.0;
RealD c2=-1.0/24.0;
RealD u0=1.0;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass,c1,c2,u0);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
MinimalResidual<FermionField> MR(1.0e-8,10000,0.8); MinimalResidual<FermionField> MR(1.0e-8,10000,0.8);