1
0
mirror of https://github.com/paboyle/Grid.git synced 2024-09-20 01:05:38 +01:00

Tests on frontier with real speed up . 3.5x on 16^3 at mq=0.01

This commit is contained in:
Peter Boyle 2023-10-20 19:25:39 -04:00
parent b6ad1bafc7
commit 9c9c42d0df
2 changed files with 8 additions and 8 deletions

View File

@ -144,10 +144,10 @@ int main (int argc, char ** argv)
{ {
Grid_init(&argc,&argv); Grid_init(&argc,&argv);
const int Ls=16; const int Ls=24;
const int nbasis = 40; const int nbasis = 40;
const int cb = 0 ; const int cb = 0 ;
RealD mass=0.01; RealD mass=0.00078;
RealD M5=1.8; RealD M5=1.8;
RealD b=1.5; RealD b=1.5;
RealD c=0.5; RealD c=0.5;
@ -219,10 +219,10 @@ int main (int argc, char ** argv)
//////////////////////////////////////////////////////////// ////////////////////////////////////////////////////////////
LittleDiracOperator LittleDiracOp(geom,FrbGrid,Coarse5d); LittleDiracOperator LittleDiracOp(geom,FrbGrid,Coarse5d);
bool load=true; bool load=false;
if ( load ) { if ( load ) {
LoadBasis(Aggregates,"Subspace.scidac"); LoadBasis(Aggregates,"/lustre/orion/phy157/proj-shared/phy157_dwf/paboyle/Subspace.scidac");
LoadOperator(LittleDiracOp,"LittleDiracOp.scidac"); LoadOperator(LittleDiracOp,"/lustre/orion/phy157/proj-shared/phy157_dwf/paboyle/LittleDiracOp.scidac");
} else { } else {
Aggregates.CreateSubspaceChebyshev(RNG5,HermOpEO,nbasis, Aggregates.CreateSubspaceChebyshev(RNG5,HermOpEO,nbasis,
95.0,0.1, 95.0,0.1,
@ -235,8 +235,8 @@ int main (int argc, char ** argv)
100, 100,
0.0); 0.0);
LittleDiracOp.CoarsenOperator(FineHermOp,Aggregates); LittleDiracOp.CoarsenOperator(FineHermOp,Aggregates);
SaveBasis(Aggregates,"Subspace.scidac"); SaveBasis(Aggregates,"/lustre/orion/phy157/proj-shared/phy157_dwf/paboyle/Subspace.scidac");
SaveOperator(LittleDiracOp,"LittleDiracOp.scidac"); SaveOperator(LittleDiracOp,"/lustre/orion/phy157/proj-shared/phy157_dwf/paboyle/LittleDiracOp.scidac");
} }
// Try projecting to one hop only // Try projecting to one hop only

View File

@ -257,7 +257,7 @@ int main (int argc, char ** argv)
// Build a coarse lanczos // Build a coarse lanczos
////////////////////////////////////////// //////////////////////////////////////////
// Chebyshev<CoarseVector> IRLCheby(0.01,44.0,201); // 1 iter // Chebyshev<CoarseVector> IRLCheby(0.01,44.0,201); // 1 iter
Chebyshev<CoarseVector> IRLCheby(0.006,44.0,301); // 1 iter Chebyshev<CoarseVector> IRLCheby(0.005,44.0,401); // 1 iter
FunctionHermOp<CoarseVector> IRLOpCheby(IRLCheby,CoarseOp); FunctionHermOp<CoarseVector> IRLOpCheby(IRLCheby,CoarseOp);
PlainHermOp<CoarseVector> IRLOp (CoarseOp); PlainHermOp<CoarseVector> IRLOp (CoarseOp);
int Nk=160; int Nk=160;