portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-13 01:05:36 +00:00
Code Releases Activity

Merge branch 'feature/distil' of github.com:mmphys/Grid into feature/distil

Browse Source
This commit is contained in:
Michael Marshall 2019-01-30 21:16:30 +00:00
commit dae7b30b92
2 changed files with 14 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
Hadrons/Modules/MDistil/DistilVectors.hpp
View File

@ -96,11 +96,16 @@ void TDistilVectors<FImpl>::setup(void)
{ {
//auto &noise = envGet(std::vector<std::vector<std::vector<SpinVector>>>, par().noise); //auto &noise = envGet(std::vector<std::vector<std::vector<SpinVector>>>, par().noise);
auto &noise = envGet(std::vector<Complex>, par().noise); auto &noise = envGet(std::vector<Complex>, par().noise);
int nnoise=par().nnoise;
int LI=par().LI;
int Ns=par().Ns;
int Nt_inv=par().Nt_inv;
envCreate(std::vector<FermionField>, getName() + "_rho", 1, envCreate(std::vector<FermionField>, getName() + "_rho", 1,
noise.size(), envGetGrid(FermionField)); nnoise*LI*Ns*Nt_inv, envGetGrid(FermionField));
envCreate(std::vector<FermionField>, getName() + "_phi", 1, envCreate(std::vector<FermionField>, getName() + "_phi", 1,
noise.size(), envGetGrid(FermionField)); nnoise*LI*Ns*Nt_inv, envGetGrid(FermionField));
GridCartesian * grid4d = env().getGrid(); GridCartesian * grid4d = env().getGrid();
@ -218,6 +223,8 @@ void TDistilVectors<FImpl>::execute(void)
} }
} }
std::cout << "size rho" << rho.size() << std::endl;
std::cout << "size phi" << phi.size() << std::endl;
} }

7
tests/hadrons/Test_hadrons_distil.cc
View File

@ -134,12 +134,13 @@ void test_MesonField(Application &application)
// DistilVectors parameters // DistilVectors parameters
MContraction::A2AMesonField::Par A2AMesonFieldPar; MContraction::A2AMesonField::Par A2AMesonFieldPar;
A2AMesonFieldPar.left="DistilVecs_phi"; A2AMesonFieldPar.left="DistilVecs_phi";
A2AMesonFieldPar.right="DistilVecs_rho"; //A2AMesonFieldPar.right="DistilVecs_rho";
A2AMesonFieldPar.right="DistilVecs_phi";
A2AMesonFieldPar.output="DistilFields"; A2AMesonFieldPar.output="DistilFields";
A2AMesonFieldPar.gammas="all"; A2AMesonFieldPar.gammas="all";
A2AMesonFieldPar.mom={"0 0 0"}; A2AMesonFieldPar.mom={"0 0 0"};
A2AMesonFieldPar.cacheBlock=4; A2AMesonFieldPar.cacheBlock=2;
A2AMesonFieldPar.block=16; A2AMesonFieldPar.block=4;
application.createModule<MContraction::A2AMesonField>("DistilMesonField",A2AMesonFieldPar); application.createModule<MContraction::A2AMesonField>("DistilMesonField",A2AMesonFieldPar);
} }