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

changed default behaviour as discussed with antonin

Browse Source
This commit is contained in:
ferben 2019-11-15 11:00:25 +00:00
parent 25d2521d77
commit b1e8b5b5ce
1 changed files with 9 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Hadrons/Modules/MDistil/PerambFromSolve.hpp
View File

@ -41,7 +41,13 @@ BEGIN_MODULE_NAMESPACE(MDistil)
This module computes a perambulator from an already completed solve. This module computes a perambulator from an already completed solve.
Optionally, the number of eigenvectors used in the perambulator and the Optionally, the number of eigenvectors used in the perambulator and the
parameter LI can be chosen to be lower than the ones in the solve, allowing parameter LI can be chosen to be lower than the ones in the solve, allowing
for a study of the signal with different values of nvec. for a study of the signal with different values of nvec.
LI_reduced : value of LI actually used in the computation
nvec_reduced: value of nvec actually used in the computation
LI : value of LI used to compute the 'solve'
nvec : value of nvec used to compute the 'solve'
******************************************************************************/ ******************************************************************************/
class PerambFromSolvePar: Serializable class PerambFromSolvePar: Serializable
@ -107,8 +113,8 @@ void TPerambFromSolve<FImpl>::setup(void)
const bool full_tdil{ dp.TI == Nt }; const bool full_tdil{ dp.TI == Nt };
const int Nt_inv{ full_tdil ? 1 : dp.TI }; const int Nt_inv{ full_tdil ? 1 : dp.TI };
grid3d.reset( MakeLowerDimGrid( env().getGrid() ) ); grid3d.reset( MakeLowerDimGrid( env().getGrid() ) );
const int nvec_reduced{par().nvec_reduced.empty() ? dp.nvec : std::stoi(par().nvec_reduced)}; const int nvec_reduced{par().nvec_reduced};
const int LI_reduced{ par().LI_reduced.empty() ? dp.LI : std::stoi(par().LI_reduced)}; const int LI_reduced{ par().LI_reduced};
envCreate(PerambTensor, getName(), 1, Nt,nvec_reduced,LI_reduced,dp.nnoise,Nt_inv,dp.SI); envCreate(PerambTensor, getName(), 1, Nt,nvec_reduced,LI_reduced,dp.nnoise,Nt_inv,dp.SI);
envCreate(NoiseTensor, getName() + "_noise", 1, dp.nnoise, Nt, dp.nvec, Ns ); envCreate(NoiseTensor, getName() + "_noise", 1, dp.nnoise, Nt, dp.nvec, Ns );
envTmp(LatticeColourVector, "result3d_nospin",1,LatticeColourVector(grid3d.get())); envTmp(LatticeColourVector, "result3d_nospin",1,LatticeColourVector(grid3d.get()));