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

uses evec4d now

Browse Source
This commit is contained in:
Felix Erben 2019-01-18 15:58:10 +00:00
parent 2568504821
commit f0f1ba0307
1 changed files with 51 additions and 49 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -8,59 +8,60 @@
BEGIN_HADRONS_NAMESPACE BEGIN_HADRONS_NAMESPACE
/****************************************************************************** /******************************************************************************
* * Class to generate Distillation $\varrho$ and $\varphi$ vectors * * Class to generate Distillation $\varrho$ and $\varphi$ vectors *
* ******************************************************************************/ ******************************************************************************/
/* /*
* template <typename FImpl> template <typename FImpl>
* class DistilVectorsFromPerambulator class DistilVectorsFromPerambulator
* { {
* public: public:
* FERM_TYPE_ALIASES(FImpl,); FERM_TYPE_ALIASES(FImpl,);
* SOLVER_TYPE_ALIASES(FImpl,); SOLVER_TYPE_ALIASES(FImpl,);
* public: public:
* DistilVectorsFromPerambulator(FMat &action); DistilVectorsFromPerambulator(FMat &action);
* virtual ~DistilVectorsFromPerambulator(void) = default; virtual ~DistilVectorsFromPerambulator(void) = default;
* void makeRho(FermionField &rhoOut, void makeRho(FermionField &rhoOut,
* const FermionField &evec, const Complex &noises); const FermionField &evec, const Complex &noises);
* void makePhi(FermionField &phiOut, void makePhi(FermionField &phiOut,
* const FermionField &evec, const SpinVector &Perambulators); const FermionField &evec, const SpinVector &Perambulators);
* private: private:
* GridBase *fGrid_, *frbGrid_, *gGrid_; GridBase *fGrid_, *frbGrid_, *gGrid_;
* bool is5d_; bool is5d_;
* FermionField src_o_, sol_e_, sol_o_, tmp_, tmp5_; FermionField src_o_, sol_e_, sol_o_, tmp_, tmp5_;
* }; };
*
* */ */
/****************************************************************************** /******************************************************************************
* * * A2AVectorsSchurDiagTwo template implementation * * A2AVectorsSchurDiagTwo template implementation *
* * ******************************************************************************/ ******************************************************************************/
/* /*
*
* template <typename FImpl> template <typename FImpl>
* A2AVectorsSchurDiagTwo<FImpl>::A2AVectorsSchurDiagTwo(FMat &action) A2AVectorsSchurDiagTwo<FImpl>::A2AVectorsSchurDiagTwo(FMat &action)
* : action_(action) : action_(action)
* , fGrid_(action_.FermionGrid()) , fGrid_(action_.FermionGrid())
* , frbGrid_(action_.FermionRedBlackGrid()) , frbGrid_(action_.FermionRedBlackGrid())
* , gGrid_(action_.GaugeGrid()) , gGrid_(action_.GaugeGrid())
* , src_o_(frbGrid_) , src_o_(frbGrid_)
* , sol_e_(frbGrid_) , sol_e_(frbGrid_)
* , sol_o_(frbGrid_) , sol_o_(frbGrid_)
* , tmp_(frbGrid_) , tmp_(frbGrid_)
* , tmp5_(fGrid_) , tmp5_(fGrid_)
* , op_(action_) , op_(action_)
* {} {}
*
* template <typename FImpl> template <typename FImpl>
* void A2AVectorsSchurDiagTwo<FImpl>::makeRho(FermionField &rhoOut, const FermionField &evec, const Complex &noises) void A2AVectorsSchurDiagTwo<FImpl>::makeRho(FermionField &rhoOut, const FermionField &evec, const Complex &noises)
* { {
*
* LatticeSpinColourVector tmp2(grid4d); LatticeSpinColourVector tmp2(grid4d);
* LatticeSpinColourVector tmp3d(grid3d); LatticeSpinColourVector tmp3d(grid3d);
* LatticeColourVector tmp3d_nospin(grid3d); LatticeColourVector tmp3d_nospin(grid3d);
* LatticeColourVector tmp_nospin(grid4d); LatticeColourVector tmp_nospin(grid4d);
* LatticeColourVector evec3d(grid3d);
for (int inoise = 0; inoise < nnoise; inoise++) { for (int inoise = 0; inoise < nnoise; inoise++) {
for (int dk = 0; dk < LI; dk++) { for (int dk = 0; dk < LI; dk++) {
@ -74,7 +75,8 @@ BEGIN_HADRONS_NAMESPACE
if( it >= Ntfirst && it < Ntfirst + Ntlocal ) { if( it >= Ntfirst && it < Ntfirst + Ntlocal ) {
for (int ik = dk; ik < nvec; ik += LI){ for (int ik = dk; ik < nvec; ik += LI){
for (int is = ds; is < Ns; is += Ns){ //at the moment, full spin dilution is enforced for (int is = ds; is < Ns; is += Ns){ //at the moment, full spin dilution is enforced
tmp3d_nospin = eig[it].evec[ik] * noises[inoise][it][ik]()(is)(); //noises do not have to be a spin vector ExtractSliceLocal(evec3d,eig4d.evec[ik],0,it,3);
tmp3d_nospin = evec3d * noises[inoise][it][ik]()(is)(); //noises do not have to be a spin vector
tmp3d=zero; tmp3d=zero;
pokeSpin(tmp3d,tmp3d_nospin,is); pokeSpin(tmp3d,tmp3d_nospin,is);
tmp2=zero; tmp2=zero;