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

Adding report at the end of the DWF HMC tests

Browse Source
This commit is contained in:
Guido Cossu 2016-12-12 04:21:34 +00:00
parent 2bd4233919
commit 7bc2065113
2 changed files with 30 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
tests/hmc/Test_hmc_EODWFRatioLsVectorised_Binary.cc
View File

@ -63,7 +63,8 @@ public:
{ {
typedef DomainWallVec5dImplR ImplPolicy; typedef DomainWallVec5dImplR ImplPolicy;
typedef ScaledShamirFermion<ImplPolicy> FermionAction; //typedef ScaledShamirFermion<ImplPolicy> FermionAction;
typedef DomainWallFermion<ImplPolicy> FermionAction;
typedef typename FermionAction::FermionField FermionField; typedef typename FermionAction::FermionField FermionField;
const int Ls = 8; const int Ls = 8;
@ -83,15 +84,21 @@ public:
LatticeGaugeField U(UGrid); LatticeGaugeField U(UGrid);
// Gauge action // Gauge action
double beta = 4.0; double beta = 5.6;
WilsonGaugeActionR Waction(beta); WilsonGaugeActionR Waction(beta);
Real mass = 0.04; Real mass = 0.04;
Real pv = 1.0; Real pv = 1.0;
RealD M5 = 1.5; RealD M5 = 1.5;
RealD scale = 2.0; RealD scale = 2.0;
/*
// Scaled Shamir
FermionAction DenOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,mass,M5,scale); FermionAction DenOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,mass,M5,scale);
FermionAction NumOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,pv,M5,scale); FermionAction NumOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,pv,M5,scale);
*/
// Std Domain wall fermion
FermionAction DenOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,mass,M5);
FermionAction NumOp(U,*FGrid,*FrbGrid,*sUGrid,*sUrbGrid,pv,M5);
double StoppingCondition = 1.0e-8; double StoppingCondition = 1.0e-8;
double MaxCGIterations = 10000; double MaxCGIterations = 10000;
@ -99,7 +106,7 @@ public:
TwoFlavourEvenOddRatioPseudoFermionAction<ImplPolicy> Nf2(NumOp, DenOp,CG,CG); TwoFlavourEvenOddRatioPseudoFermionAction<ImplPolicy> Nf2(NumOp, DenOp,CG,CG);
// Set smearing (true/false), default: false // Set smearing (true/false), default: false
Nf2.is_smeared = true; Nf2.is_smeared = false;
// Collect actions // Collect actions
// here an example of 2 level integration // here an example of 2 level integration
@ -140,8 +147,17 @@ public:
UGrid, Nsmear, Stout); UGrid, Nsmear, Stout);
/////////////////// ///////////////////
Run(argc, argv, Checkpoint, SmearingPolicy); NumOp.ZeroCounters();
//Run(argc, argv, Checkpoint); // no smearing DenOp.ZeroCounters();
//Run(argc, argv, Checkpoint, SmearingPolicy);
Run(argc, argv, Checkpoint); // no smearing
std::cout << GridLogMessage << "Numerator report, Pauli-Villars term : " << std::endl;
NumOp.Report();
std::cout << GridLogMessage << "Denominator report, Dw(m) term (includes CG) : " << std::endl;
DenOp.Report();
}; };
}; };
} }

9
tests/hmc/Test_hmc_EODWFRatio_Binary.cc
View File

@ -135,8 +135,17 @@ public:
UGrid, Nsmear, Stout); UGrid, Nsmear, Stout);
/////////////////// ///////////////////
NumOp.ZeroCounters();
DenOp.ZeroCounters();
Run(argc, argv, Checkpoint, SmearingPolicy); Run(argc, argv, Checkpoint, SmearingPolicy);
//Run(argc, argv, Checkpoint); // no smearing //Run(argc, argv, Checkpoint); // no smearing
std::cout << GridLogMessage << "Numerator report, Pauli-Villars term : " << std::endl;
NumOp.Report();
std::cout << GridLogMessage << "Denominator report, Dw(m) term (includes CG) : " << std::endl;
DenOp.Report();
}; };
}; };
} }