portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-10 07:55:35 +00:00
Code Releases Activity

Namespace

Browse Source
This commit is contained in:
paboyle 2018-01-14 22:01:17 +00:00
parent fdcbe0a0d1
commit 9aa34dc803
1 changed files with 11 additions and 11 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

22
lib/qcd/utils/GaugeFix.h
View File

@ -1,4 +1,4 @@
/************************************************************************************* /*************************************************************************************
grid` physics library, www.github.com/paboyle/Grid grid` physics library, www.github.com/paboyle/Grid
@ -22,18 +22,18 @@ Author: Peter Boyle <paboyle@ph.ed.ac.uk>
51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
See the full license in the file "LICENSE" in the top level distribution directory See the full license in the file "LICENSE" in the top level distribution directory
*************************************************************************************/ *************************************************************************************/
/* END LEGAL */ /* END LEGAL */
//#include <Grid/Grid.h> //#include <Grid/Grid.h>
#ifndef GRID_QCD_GAUGE_FIX_H #ifndef GRID_QCD_GAUGE_FIX_H
#define GRID_QCD_GAUGE_FIX_H #define GRID_QCD_GAUGE_FIX_H
namespace Grid {
namespace QCD { NAMESPACE_BEGIN(Grid);
template <class Gimpl> template <class Gimpl>
class FourierAcceleratedGaugeFixer : public Gimpl { class FourierAcceleratedGaugeFixer : public Gimpl {
public: public:
INHERIT_GIMPL_TYPES(Gimpl); INHERIT_GIMPL_TYPES(Gimpl);
typedef typename Gimpl::GaugeLinkField GaugeMat; typedef typename Gimpl::GaugeLinkField GaugeMat;
@ -60,7 +60,7 @@ class FourierAcceleratedGaugeFixer : public Gimpl {
Real trG; Real trG;
std::vector<GaugeMat> U(Nd,grid); std::vector<GaugeMat> U(Nd,grid);
GaugeMat dmuAmu(grid); GaugeMat dmuAmu(grid);
for(int i=0;i<maxiter;i++){ for(int i=0;i<maxiter;i++){
for(int mu=0;mu<Nd;mu++) U[mu]= PeekIndex<LorentzIndex>(Umu,mu); for(int mu=0;mu<Nd;mu++) U[mu]= PeekIndex<LorentzIndex>(Umu,mu);
@ -154,8 +154,8 @@ class FourierAcceleratedGaugeFixer : public Gimpl {
Fp = psqMax*one/psq; Fp = psqMax*one/psq;
/* /*
static int once; static int once;
if ( once == 0 ) { if ( once == 0 ) {
std::cout << " Fp " << Fp <<std::endl; std::cout << " Fp " << Fp <<std::endl;
once ++; once ++;
}*/ }*/
@ -188,6 +188,6 @@ class FourierAcceleratedGaugeFixer : public Gimpl {
} }
}; };
} NAMESPACE_END(Grid);
}
#endif #endif