Improved Gauge Action
More...
#include <fermiqcd_gauge_actions.h>
Inheritance diagram for ImprovedGaugeAction:

Static Public Member Functions | |
| static gauge_stats | heatbath (gauge_field &U, coefficients &coeff, int n_iter=1, string model="MILC") |
Improved Gauge Action
Example using the MILC improved action:
int ns=2, steps=10;
gauge_field U(lattice,nc);
coefficients gauge;
U.load("myfield.0000");
gauge["beta"]=6.0;
gauge["zeta"]=1.0; // MUST BE ONE
gauge["u_t"]=1.0;
gauge["u_s"]=1.0;
ImprovedGaugeAction::heatbath(U,gauge,steps,"MILC");
U.save("myfield.0001");
int ns=2, steps=10;
gauge_field U(lattice,nc);
coefficients gauge;
U.load("myfield.0000");
gauge["beta"]=6.0;
gauge["zeta"]=1.0; // CAN BE != ONE
gauge["u_t"]=1.0;
gauge["u_s"]=1.0;
ImprovedGaugeAction::heatbath(U,gauge,steps,"Morningstar");
U.save("myfield.0001");
1.4.1