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

Begin fixes for single precision

Browse Source
This commit is contained in:
Michael Marshall 2019-03-25 20:40:05 +00:00
parent 625a97a466
commit 8502660023
3 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
Hadrons/Modules/MDistil/BContraction.hpp
View File

@ -215,7 +215,7 @@ void TBContraction<FImpl>::execute(void)
tmp111 = 0.5*(double)parity*(tmp111 + tmp222); // P_\pm * ... tmp111 = 0.5*(double)parity*(tmp111 + tmp222); // P_\pm * ...
diquark2 = factor23[0]*innerProduct(tmp22s,tmp333); diquark2 = factor23[0]*innerProduct(tmp22s,tmp333);
for (int is=0 ; is < 4 ; is++){ for (int is=0 ; is < 4 ; is++){
BField3(imom,ig,t,is,i1,i2,i3)+=(double)epsilon_sgn[ie]*tmp111()(is)()*diquark2; BField3(imom,ig,t,is,i1,i2,i3)+=static_cast<Real>(epsilon_sgn[ie])*tmp111()(is)()*diquark2;
} }
} }
} }

2
Hadrons/Modules/MDistil/Baryon2pt.hpp
View File

@ -169,7 +169,7 @@ void TBaryon2pt<FImpl>::execute(void)
Eigen::Tensor<Complex,3> B3L = B4L.chip(is,0); Eigen::Tensor<Complex,3> B3L = B4L.chip(is,0);
Eigen::Tensor<Complex,3> B3R = B4R.chip(is,0); Eigen::Tensor<Complex,3> B3R = B4R.chip(is,0);
Eigen::Tensor<Complex,0> C2 = B3L.contract(B3R,product_dims); Eigen::Tensor<Complex,0> C2 = B3L.contract(B3R,product_dims);
corr(imom,t) += (double)epsilon_sgn[pairs[ipair]]*C2(0); corr(imom,t) += static_cast<Real>(epsilon_sgn[pairs[ipair]])*C2(0);
} }
} }
} }

2
Hadrons/Modules/MDistil/Distil.hpp
View File

@ -146,7 +146,7 @@ inline void SliceShare( GridBase * gridLowDim, GridBase * gridHighDim, void * Bu
*************************************************************************************/ *************************************************************************************/
template<typename Field, typename GaugeField=LatticeGaugeFieldD> template<typename Field, typename GaugeField=LatticeGaugeField>
class LinOpPeardonNabla : public LinearOperatorBase<Field>, public LinearFunction<Field> { class LinOpPeardonNabla : public LinearOperatorBase<Field>, public LinearFunction<Field> {
typedef typename GaugeField::vector_type vCoeff_t; typedef typename GaugeField::vector_type vCoeff_t;
protected: // I don't really mind if _gf is messed with ... so make this public? protected: // I don't really mind if _gf is messed with ... so make this public?