portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-10-31 12:04:33 +00:00
Code Releases Activity

Make solver parameters match in all MR|GMRES solver tests

Browse Source
This commit is contained in:
Daniel Richtmann
2018-02-08 17:33:10 +01:00
parent 9f79a87102
commit 13ae371ef8
7 changed files with 10 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tests/solver/Test_staggered_fgmres_prec.cc
View File

@@ -58,14 +58,14 @@ int main (int argc, char ** argv)
volume=volume*latt_size[mu]; volume=volume*latt_size[mu];
} }
RealD mass=0.1; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
TrivialPrecon<FermionField> simple; TrivialPrecon<FermionField> simple;
FlexibleGeneralisedMinimalResidual<FermionField> FGMRES(1.0e-6, 50, simple, 25); FlexibleGeneralisedMinimalResidual<FermionField> FGMRES(1.0e-8, 10000, simple, 25);
FGMRES(HermOp,src,result); FGMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();

4
tests/solver/Test_staggered_gmres_unprec.cc
View File

@@ -58,11 +58,11 @@ int main (int argc, char ** argv)
volume=volume*latt_size[mu]; volume=volume*latt_size[mu];
} }
RealD mass=0.1; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
GeneralisedMinimalResidual<FermionField> GMRES(1.0e-6, 50, 25); GeneralisedMinimalResidual<FermionField> GMRES(1.0e-8, 10000, 25);
GMRES(HermOp,src,result); GMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();

4
tests/solver/Test_staggered_mr_unprec.cc
View File

@@ -58,11 +58,11 @@ int main (int argc, char ** argv)
volume=volume*latt_size[mu]; volume=volume*latt_size[mu];
} }
RealD mass=0.1; RealD mass=0.5;
ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass); ImprovedStaggeredFermionR Ds(Umu,Umu,Grid,RBGrid,mass);
MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds); MdagMLinearOperator<ImprovedStaggeredFermionR,FermionField> HermOp(Ds);
MinimalResidual<FermionField> MR(1.0e-6,20000,0.8); MinimalResidual<FermionField> MR(1.0e-8,10000,0.8);
MR(HermOp,src,result); MR(HermOp,src,result);
Grid_finalize(); Grid_finalize();

2
tests/solver/Test_wilson_cagmres_unprec.cc
View File

@@ -58,7 +58,7 @@ int main (int argc, char ** argv)
WilsonFermionR Dw(Umu,Grid,RBGrid,mass); WilsonFermionR Dw(Umu,Grid,RBGrid,mass);
MdagMLinearOperator<WilsonFermionR,LatticeFermion> HermOp(Dw); MdagMLinearOperator<WilsonFermionR,LatticeFermion> HermOp(Dw);
CommunicationAvoidingGeneralisedMinimalResidual<LatticeFermion> CAGMRES(1.0e-8, 50, 25); CommunicationAvoidingGeneralisedMinimalResidual<LatticeFermion> CAGMRES(1.0e-8, 10000, 25);
CAGMRES(HermOp,src,result); CAGMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();

2
tests/solver/Test_wilson_fcagmres_prec.cc
View File

@@ -61,7 +61,7 @@ int main (int argc, char ** argv)
TrivialPrecon<LatticeFermion> simple; TrivialPrecon<LatticeFermion> simple;
FlexibleCommunicationAvoidingGeneralisedMinimalResidual<LatticeFermion> FCAGMRES(1.0e-8, 50, simple, 25); FlexibleCommunicationAvoidingGeneralisedMinimalResidual<LatticeFermion> FCAGMRES(1.0e-8, 10000, simple, 25);
FCAGMRES(HermOp,src,result); FCAGMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();

2
tests/solver/Test_wilson_fgmres_prec.cc
View File

@@ -61,7 +61,7 @@ int main (int argc, char ** argv)
TrivialPrecon<LatticeFermion> simple; TrivialPrecon<LatticeFermion> simple;
FlexibleGeneralisedMinimalResidual<LatticeFermion> FGMRES(1.0e-8, 50, simple, 25); FlexibleGeneralisedMinimalResidual<LatticeFermion> FGMRES(1.0e-8, 10000, simple, 25);
FGMRES(HermOp,src,result); FGMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();

2
tests/solver/Test_wilson_gmres_unprec.cc
View File

@@ -58,7 +58,7 @@ int main (int argc, char ** argv)
WilsonFermionR Dw(Umu,Grid,RBGrid,mass); WilsonFermionR Dw(Umu,Grid,RBGrid,mass);
MdagMLinearOperator<WilsonFermionR,LatticeFermion> HermOp(Dw); MdagMLinearOperator<WilsonFermionR,LatticeFermion> HermOp(Dw);
GeneralisedMinimalResidual<LatticeFermion> GMRES(1.0e-8, 50, 25); GeneralisedMinimalResidual<LatticeFermion> GMRES(1.0e-8, 10000, 25);
GMRES(HermOp,src,result); GMRES(HermOp,src,result);
Grid_finalize(); Grid_finalize();