portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2025-11-04 05:54:32 +00:00
Code Releases Activity

meson fields are now the same

Browse Source
This commit is contained in:
ferben
2019-01-30 18:03:34 +00:00
parent c3273eff20
commit 7b66197534
2 changed files with 14 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -97,10 +97,15 @@ void TDistilVectors<FImpl>::setup(void)
//auto &noise = envGet(std::vector<std::vector<std::vector<SpinVector>>>, 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,
noise.size(), envGetGrid(FermionField));
nnoise*LI*Ns*Nt_inv, envGetGrid(FermionField));
envCreate(std::vector<FermionField>, getName() + "_phi", 1,
noise.size(), envGetGrid(FermionField));
nnoise*LI*Ns*Nt_inv, envGetGrid(FermionField));
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
MContraction::A2AMesonField::Par A2AMesonFieldPar;
A2AMesonFieldPar.left="DistilVecs_phi";
A2AMesonFieldPar.right="DistilVecs_rho";
//A2AMesonFieldPar.right="DistilVecs_rho";
A2AMesonFieldPar.right="DistilVecs_phi";
A2AMesonFieldPar.output="DistilFields";
A2AMesonFieldPar.gammas="all";
A2AMesonFieldPar.mom={"0 0 0"};
A2AMesonFieldPar.cacheBlock=4;
A2AMesonFieldPar.block=16;
A2AMesonFieldPar.cacheBlock=2;
A2AMesonFieldPar.block=4;
application.createModule<MContraction::A2AMesonField>("DistilMesonField",A2AMesonFieldPar);
}