From 84189867ef5257df5e93b6a0bce333609623f001 Mon Sep 17 00:00:00 2001 From: Antonin Portelli Date: Tue, 2 Oct 2018 13:05:05 +0100 Subject: [PATCH] Hadrons: eigenpack converter with RB grids (to be generalised) --- Hadrons/Utilities/EigenPackCast.cc | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/Hadrons/Utilities/EigenPackCast.cc b/Hadrons/Utilities/EigenPackCast.cc index 3d74a3aa..e9be2c15 100644 --- a/Hadrons/Utilities/EigenPackCast.cc +++ b/Hadrons/Utilities/EigenPackCast.cc @@ -16,15 +16,17 @@ void convert(const std::string outFilename, const std::string inFilename, typedef typename FOut::vector_type VTypeOut; typedef typename FIn::vector_type VTypeIn; - std::shared_ptr gIn, gOut; + std::shared_ptr fgIn, fgOut, gIn, gOut; auto dim = GridDefaultLatt(); unsigned int nd = dim.size(); auto simdOut = GridDefaultSimd(nd, VTypeOut::Nsimd()); auto simdIn = GridDefaultSimd(nd, VTypeIn::Nsimd()); - gOut.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdOut, GridDefaultMpi())); - gIn.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdIn, GridDefaultMpi())); - + fgOut.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdOut, GridDefaultMpi())); + gOut.reset(SpaceTimeGrid::makeFourDimRedBlackGrid(fgOut.get())); + fgIn.reset(SpaceTimeGrid::makeFourDimGrid(dim, simdIn, GridDefaultMpi())); + gIn.reset(SpaceTimeGrid::makeFourDimRedBlackGrid(fgOut.get())); + FOut bufOut(gOut.get()); FIn bufIn(gIn.get()), testIn(gIn.get());