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

Moving output order

Browse Source
This commit is contained in:
Guido Cossu 2016-12-14 02:02:42 +00:00
parent 2fb92dbc6e
commit af0ccdd8e9
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
tests/solver/Test_dwf_cg_prec_LsVec.cc
View File

@ -57,6 +57,9 @@ int main(int argc, char** argv) {
GridCartesian* FGrid = SpaceTimeGrid::makeFiveDimDWFGrid(Ls, UGrid); GridCartesian* FGrid = SpaceTimeGrid::makeFiveDimDWFGrid(Ls, UGrid);
GridRedBlackCartesian* FrbGrid = SpaceTimeGrid::makeFiveDimDWFRedBlackGrid(Ls, UGrid); GridRedBlackCartesian* FrbGrid = SpaceTimeGrid::makeFiveDimDWFRedBlackGrid(Ls, UGrid);
std::cout << GridLogMessage << "Lattice dimensions: " << GridDefaultLatt()
<< " Ls: " << Ls << std::endl;
std::vector<int> seeds4({1, 2, 3, 4}); std::vector<int> seeds4({1, 2, 3, 4});
std::vector<int> seeds5({5, 6, 7, 8}); std::vector<int> seeds5({5, 6, 7, 8});
GridParallelRNG RNG5(FGrid); GridParallelRNG RNG5(FGrid);
@ -64,17 +67,17 @@ int main(int argc, char** argv) {
GridParallelRNG RNG4(UGrid); GridParallelRNG RNG4(UGrid);
RNG4.SeedFixedIntegers(seeds4); RNG4.SeedFixedIntegers(seeds4);
std::cout << GridLogMessage << "Generating source field" << std::endl; std::cout << GridLogMessage << "Generating random ferrmion field" << std::endl;
LatticeFermion src(FGrid); LatticeFermion src(FGrid);
random(RNG5, src); random(RNG5, src);
LatticeFermion result(FGrid); LatticeFermion result(FGrid);
result = zero; result = zero;
LatticeGaugeField Umu(UGrid); LatticeGaugeField Umu(UGrid);
std::cout << GridLogMessage << "Generating random gauge field" << std::endl;
SU3::HotConfiguration(RNG4, Umu); SU3::HotConfiguration(RNG4, Umu);
std::cout << GridLogMessage << "Lattice dimensions: " << GridDefaultLatt()
<< " Ls: " << Ls << std::endl;
std::vector<LatticeColourMatrix> U(4, UGrid); std::vector<LatticeColourMatrix> U(4, UGrid);
for (int mu = 0; mu < Nd; mu++) { for (int mu = 0; mu < Nd; mu++) {