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

Hadrons: precision fixes

Browse Source
This commit is contained in:
Antonin Portelli 2018-08-28 16:58:23 +01:00
parent f3a0158213
commit fb3596f968
2 changed files with 7 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

12
Hadrons/EigenPack.hpp
View File

@ -56,7 +56,7 @@ public:
VecRecord(void): index(0), eval(0.) {} VecRecord(void): index(0), eval(0.) {}
}; };
public: public:
std::vector<Real> eval; std::vector<RealD> eval;
std::vector<F> evec; std::vector<F> evec;
PackRecord record; PackRecord record;
public: public:
@ -119,7 +119,7 @@ protected:
} }
template <typename T> template <typename T>
void basicRead(std::vector<T> &evec, std::vector<double> &eval, void basicRead(std::vector<T> &evec, std::vector<RealD> &eval,
const std::string filename, const unsigned int size) const std::string filename, const unsigned int size)
{ {
ScidacReader binReader; ScidacReader binReader;
@ -144,7 +144,7 @@ protected:
} }
template <typename T> template <typename T>
void basicReadSingle(T &evec, double &eval, const std::string filename, void basicReadSingle(T &evec, RealD &eval, const std::string filename,
const unsigned int index) const unsigned int index)
{ {
ScidacReader binReader; ScidacReader binReader;
@ -166,7 +166,7 @@ protected:
template <typename T> template <typename T>
void basicWrite(const std::string filename, std::vector<T> &evec, void basicWrite(const std::string filename, std::vector<T> &evec,
const std::vector<double> &eval, const unsigned int size) const std::vector<RealD> &eval, const unsigned int size)
{ {
ScidacWriter binWriter(evec[0]._grid->IsBoss()); ScidacWriter binWriter(evec[0]._grid->IsBoss());
XmlWriter xmlWriter("", "eigenPackPar"); XmlWriter xmlWriter("", "eigenPackPar");
@ -190,7 +190,7 @@ protected:
template <typename T> template <typename T>
void basicWriteSingle(const std::string filename, T &evec, void basicWriteSingle(const std::string filename, T &evec,
const double eval, const unsigned int index) const RealD eval, const unsigned int index)
{ {
ScidacWriter binWriter(evec._grid->IsBoss()); ScidacWriter binWriter(evec._grid->IsBoss());
XmlWriter xmlWriter("", "eigenPackPar"); XmlWriter xmlWriter("", "eigenPackPar");
@ -215,7 +215,7 @@ class CoarseEigenPack: public EigenPack<FineF>
public: public:
typedef CoarseF CoarseField; typedef CoarseF CoarseField;
public: public:
std::vector<Real> evalCoarse; std::vector<RealD> evalCoarse;
std::vector<CoarseF> evecCoarse; std::vector<CoarseF> evecCoarse;
public: public:
CoarseEigenPack(void) = default; CoarseEigenPack(void) = default;

2
Hadrons/Modules/MContraction/A2AMesonField.hpp
View File

@ -94,7 +94,7 @@ private:
bool hasPhase_{false}; bool hasPhase_{false};
std::string momphName_; std::string momphName_;
std::vector<Gamma::Algebra> gamma_; std::vector<Gamma::Algebra> gamma_;
std::vector<std::vector<double>> mom_; std::vector<std::vector<Real>> mom_;
std::vector<std::pair<unsigned int, unsigned int>> nodeFile_; std::vector<std::pair<unsigned int, unsigned int>> nodeFile_;
}; };