mirror of
https://github.com/paboyle/Grid.git
synced 2025-04-09 21:50:45 +01:00
Remove some old code in Wilson MG
This commit is contained in:
Daniel Richtmann
parent
39558cce52
commit
6cf635d61c
@ -498,225 +498,6 @@ public:
|
|||||||
|
|
||||||
};
|
};
|
||||||
|
|
||||||
#if 0
|
|
||||||
int main (int argc, char ** argv)
|
|
||||||
{
|
|
||||||
Grid_init(&argc,&argv);
|
|
||||||
|
|
||||||
const int Ls=params.Ls;
|
|
||||||
|
|
||||||
GridCartesian * UGrid = SpaceTimeGrid::makeFourDimGrid(GridDefaultLatt(), GridDefaultSimd(Nd,vComplex::Nsimd()),GridDefaultMpi());
|
|
||||||
GridRedBlackCartesian * UrbGrid = SpaceTimeGrid::makeFourDimRedBlackGrid(UGrid);
|
|
||||||
|
|
||||||
///////////////////////////////////////////////////
|
|
||||||
// Construct a coarsened grid; utility for this?
|
|
||||||
///////////////////////////////////////////////////
|
|
||||||
std::vector<int> block ({2,2,2,2});
|
|
||||||
const int nbasis= 32;
|
|
||||||
|
|
||||||
std::vector<int> clatt = GridDefaultLatt();
|
|
||||||
for(int d=0;d<clatt.size();d++){
|
|
||||||
clatt[d] = clatt[d]/block[d];
|
|
||||||
}
|
|
||||||
GridCartesian *Coarse4d = SpaceTimeGrid::makeFourDimGrid(clatt, GridDefaultSimd(Nd,vComplex::Nsimd()),GridDefaultMpi());;
|
|
||||||
|
|
||||||
std::vector<int> seeds4({1,2,3,4});
|
|
||||||
std::vector<int> seeds5({5,6,7,8});
|
|
||||||
std::vector<int> cseeds({5,6,7,8});
|
|
||||||
GridParallelRNG RNG5(FGrid); RNG5.SeedFixedIntegers(seeds5);
|
|
||||||
GridParallelRNG RNG4(UGrid); RNG4.SeedFixedIntegers(seeds4);
|
|
||||||
GridParallelRNG CRNG(Coarse5d);CRNG.SeedFixedIntegers(cseeds);
|
|
||||||
|
|
||||||
Gamma g5(Gamma::Algebra::Gamma5);
|
|
||||||
|
|
||||||
LatticeFermion src(FGrid); gaussian(RNG5,src);// src=src+g5*src;
|
|
||||||
LatticeFermion result(FGrid); result=zero;
|
|
||||||
LatticeFermion ref(FGrid); ref=zero;
|
|
||||||
LatticeFermion tmp(FGrid);
|
|
||||||
LatticeFermion err(FGrid);
|
|
||||||
LatticeGaugeField Umu(UGrid);
|
|
||||||
LatticeGaugeField UmuDD(UGrid);
|
|
||||||
LatticeColourMatrix U(UGrid);
|
|
||||||
LatticeColourMatrix zz(UGrid);
|
|
||||||
|
|
||||||
FieldMetaData header;
|
|
||||||
std::string file("./ckpoint_lat.4000");
|
|
||||||
NerscIO::readConfiguration(Umu,header,file);
|
|
||||||
|
|
||||||
|
|
||||||
if ( params.domaindecompose ) {
|
|
||||||
Lattice<iScalar<vInteger> > coor(UGrid);
|
|
||||||
zz=zero;
|
|
||||||
for(int mu=0;mu<Nd;mu++){
|
|
||||||
LatticeCoordinate(coor,mu);
|
|
||||||
U = PeekIndex<LorentzIndex>(Umu,mu);
|
|
||||||
U = where(mod(coor,params.domainsize)==(Integer)0,zz,U);
|
|
||||||
PokeIndex<LorentzIndex>(UmuDD,U,mu);
|
|
||||||
}
|
|
||||||
} else {
|
|
||||||
UmuDD = Umu;
|
|
||||||
}
|
|
||||||
// SU3::ColdConfiguration(RNG4,Umu);
|
|
||||||
// SU3::TepidConfiguration(RNG4,Umu);
|
|
||||||
// SU3::HotConfiguration(RNG4,Umu);
|
|
||||||
// Umu=zero;
|
|
||||||
|
|
||||||
RealD mass=params.mq;
|
|
||||||
RealD M5=1.8;
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Building g5R5 hermitian DWF operator" <<std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
DomainWallFermionR Ddwf(Umu,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5);
|
|
||||||
DomainWallFermionR DdwfDD(UmuDD,*FGrid,*FrbGrid,*UGrid,*UrbGrid,mass,M5);
|
|
||||||
|
|
||||||
typedef Aggregation<vSpinColourVector,vTComplex,nbasis> Subspace;
|
|
||||||
typedef CoarsenedMatrix<vSpinColourVector,vTComplex,nbasis> CoarseOperator;
|
|
||||||
typedef CoarseOperator::CoarseVector CoarseVector;
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Calling Aggregation class to build subspace" <<std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
MdagMLinearOperator<DomainWallFermionR,LatticeFermion> HermDefOp(Ddwf);
|
|
||||||
Subspace Aggregates(Coarse5d,FGrid);
|
|
||||||
// Aggregates.CreateSubspace(RNG5,HermDefOp,nbasis);
|
|
||||||
assert ( (nbasis & 0x1)==0);
|
|
||||||
int nb=nbasis/2;
|
|
||||||
std::cout<<GridLogMessage << " nbasis/2 = "<<nb<<std::endl;
|
|
||||||
// Aggregates.CreateSubspace(RNG5,HermDefOp,nb);
|
|
||||||
Aggregates.CreateSubspaceLanczos(RNG5,HermDefOp,nb);
|
|
||||||
for(int n=0;n<nb;n++){
|
|
||||||
G5R5(Aggregates.subspace[n+nb],Aggregates.subspace[n]);
|
|
||||||
std::cout<<GridLogMessage<<n<<" subspace "<<norm2(Aggregates.subspace[n+nb])<<" "<<norm2(Aggregates.subspace[n]) <<std::endl;
|
|
||||||
}
|
|
||||||
for(int n=0;n<nbasis;n++){
|
|
||||||
std::cout<<GridLogMessage << "vec["<<n<<"] = "<<norm2(Aggregates.subspace[n]) <<std::endl;
|
|
||||||
}
|
|
||||||
|
|
||||||
// for(int i=0;i<nbasis;i++){
|
|
||||||
// result = Aggregates.subspace[i];
|
|
||||||
// Aggregates.subspace[i]=result+g5*result;
|
|
||||||
// }
|
|
||||||
result=zero;
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Building coarse representation of Indef operator" <<std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
Gamma5R5HermitianLinearOperator<DomainWallFermionR,LatticeFermion> HermIndefOp(Ddwf);
|
|
||||||
Gamma5R5HermitianLinearOperator<DomainWallFermionR,LatticeFermion> HermIndefOpDD(DdwfDD);
|
|
||||||
CoarsenedMatrix<vSpinColourVector,vTComplex,nbasis> LDOp(*Coarse5d);
|
|
||||||
LDOp.CoarsenOperator(FGrid,HermIndefOp,Aggregates);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Testing some coarse space solvers " <<std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
CoarseVector c_src (Coarse5d);
|
|
||||||
CoarseVector c_res (Coarse5d);
|
|
||||||
gaussian(CRNG,c_src);
|
|
||||||
c_res=zero;
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Solving posdef-CG on coarse space "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
MdagMLinearOperator<CoarseOperator,CoarseVector> PosdefLdop(LDOp);
|
|
||||||
ConjugateGradient<CoarseVector> CG(1.0e-6,100000);
|
|
||||||
// CG(PosdefLdop,c_src,c_res);
|
|
||||||
|
|
||||||
// std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// std::cout<<GridLogMessage << "Solving indef-MCR on coarse space "<< std::endl;
|
|
||||||
// std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// HermitianLinearOperator<CoarseOperator,CoarseVector> HermIndefLdop(LDOp);
|
|
||||||
// ConjugateResidual<CoarseVector> MCR(1.0e-6,100000);
|
|
||||||
//MCR(HermIndefLdop,c_src,c_res);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Building deflation preconditioner "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
|
|
||||||
MultiGridPreconditioner <vSpinColourVector,vTComplex,nbasis,DomainWallFermionR> Precon (Aggregates, LDOp,
|
|
||||||
HermIndefOp,Ddwf,
|
|
||||||
HermIndefOp,Ddwf);
|
|
||||||
|
|
||||||
MultiGridPreconditioner <vSpinColourVector,vTComplex,nbasis,DomainWallFermionR> PreconDD(Aggregates, LDOp,
|
|
||||||
HermIndefOp,Ddwf,
|
|
||||||
HermIndefOpDD,DdwfDD);
|
|
||||||
TrivialPrecon<LatticeFermion> simple;
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Testing smoother efficacy"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// Precon.SmootherTest(src);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Testing DD smoother efficacy"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// PreconDD.SmootherTest(src);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Testing SAP smoother efficacy"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// PreconDD.SAP(src,result);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Unprec CG "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
|
|
||||||
// TrivialPrecon<LatticeFermion> simple;
|
|
||||||
// ConjugateGradient<LatticeFermion> fCG(1.0e-8,100000);
|
|
||||||
// fCG(HermDefOp,src,result);
|
|
||||||
// exit(0);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Testing GCR on indef matrix "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// PrecGeneralisedConjugateResidual<LatticeFermion> UPGCR(1.0e-8,100000,simple,8,128);
|
|
||||||
// UPGCR(HermIndefOp,src,result);
|
|
||||||
|
|
||||||
|
|
||||||
/// Get themax eval
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage <<" Applying power method to find spectral range "<<std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
Precon.PowerMethod(src);
|
|
||||||
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Building a two level DDPGCR "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
// PrecGeneralisedConjugateResidual<LatticeFermion> PGCRDD(1.0e-8,100000,PreconDD,8,128);
|
|
||||||
// result=zero;
|
|
||||||
// std::cout<<GridLogMessage<<"checking norm src "<<norm2(src)<<std::endl;
|
|
||||||
// PGCRDD(HermIndefOp,src,result);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Building a two level PGCR "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
PrecGeneralisedConjugateResidual<LatticeFermion> PGCR(1.0e-8,100000,Precon,8,8);
|
|
||||||
std::cout<<GridLogMessage<<"checking norm src "<<norm2(src)<<std::endl;
|
|
||||||
result=zero;
|
|
||||||
PGCR(HermIndefOp,src,result);
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Red Black Prec CG "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
SchurDiagMooeeOperator<DomainWallFermionR,LatticeFermion> HermOpEO(Ddwf);
|
|
||||||
ConjugateGradient<LatticeFermion> pCG(1.0e-8,10000);
|
|
||||||
|
|
||||||
LatticeFermion src_o(FrbGrid);
|
|
||||||
LatticeFermion result_o(FrbGrid);
|
|
||||||
pickCheckerboard(Odd,src_o,src);
|
|
||||||
result_o=zero;
|
|
||||||
|
|
||||||
pCG(HermOpEO,src_o,result_o);
|
|
||||||
|
|
||||||
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "Done "<< std::endl;
|
|
||||||
std::cout<<GridLogMessage << "**************************************************"<< std::endl;
|
|
||||||
Grid_finalize();
|
|
||||||
}
|
|
||||||
|
|
||||||
#else
|
|
||||||
int main (int argc, char ** argv)
|
int main (int argc, char ** argv)
|
||||||
{
|
{
|
||||||
Grid_init(&argc,&argv);
|
Grid_init(&argc,&argv);
|
||||||
@ -883,4 +664,3 @@ int main (int argc, char ** argv)
|
|||||||
|
|
||||||
Grid_finalize();
|
Grid_finalize();
|
||||||
}
|
}
|
||||||
#endif
|
|
||||||
|
|||||||
Loading…
x
Reference in New Issue
Block a user