portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-13 01:05:36 +00:00
Code Releases Activity

Single precision compile fix. Soon deprecate single precision

Browse Source
This commit is contained in:
Peter Boyle 2019-07-18 22:02:31 +01:00
parent 9fa705c5a0
commit 0695f8cec2
2 changed files with 6 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
HMC/Mobius2p1fEOFA.cc
View File

@ -30,7 +30,9 @@ directory
/* END LEGAL */ /* END LEGAL */
#include <Grid/Grid.h> #include <Grid/Grid.h>
#ifdef GRID_DEFAULT_PRECISION_DOUBLE
#define MIXED_PRECISION #define MIXED_PRECISION
#endif
namespace Grid{ namespace Grid{
namespace QCD{ namespace QCD{
@ -348,8 +350,9 @@ int main(int argc, char **argv) {
#else #else
ExactOneFlavourRatioPseudoFermionAction<FermionImplPolicy> ExactOneFlavourRatioPseudoFermionAction<FermionImplPolicy>
EOFA(Strange_Op_L, Strange_Op_R, EOFA(Strange_Op_L, Strange_Op_R,
ActionCG,
ActionCG, ActionCG, ActionCG, ActionCG,
DerivativeCG, DerivativeCG, DerivativeCG, DerivativeCG,
OFRp, true); OFRp, true);
#endif #endif
Level1.push_back(&EOFA); Level1.push_back(&EOFA);

5
Hadrons/Modules/MNPR/Amputate.hpp
View File

@ -109,8 +109,6 @@ template <typename FImpl1, typename FImpl2>
std::vector<std::string> TAmputate<FImpl1, FImpl2>::getOutput(void) std::vector<std::string> TAmputate<FImpl1, FImpl2>::getOutput(void)
{ {
std::vector<std::string> output = {getName()}; std::vector<std::string> output = {getName()};
return output; return output;
} }
@ -123,7 +121,8 @@ SpinColourMatrix TAmputate<Fimpl1, Fimpl2>::invertspincolmat(SpinColourMatrix &s
for(int jc=0; jc<Nc; jc++){ for(int jc=0; jc<Nc; jc++){
for(int is=0; is<Ns; is++){ for(int is=0; is<Ns; is++){
for(int js=0; js<Ns; js++){ for(int js=0; js<Ns; js++){
scmat_2d(Ns*ic+is,Ns*jc+js) = scmat()(is,js)(ic,jc); auto z =scmat()(is,js)(ic,jc);;
scmat_2d(Ns*ic+is,Ns*jc+js) = std::complex<double>(real(z),imag(z));
}} }}
}} }}
Eigen::MatrixXcd scmat_2d_inv = scmat_2d.inverse(); Eigen::MatrixXcd scmat_2d_inv = scmat_2d.inverse();