mirror of
https://github.com/paboyle/Grid.git
synced 2024-11-14 01:35:36 +00:00
Mom filter refresh sRNG
This commit is contained in:
parent
a198d59381
commit
0f08364e4f
@ -61,7 +61,9 @@ int main (int argc, char ** argv)
|
|||||||
std::vector<int> seeds({1,2,3,4});
|
std::vector<int> seeds({1,2,3,4});
|
||||||
|
|
||||||
GridParallelRNG pRNG(&Grid);
|
GridParallelRNG pRNG(&Grid);
|
||||||
|
GridSerialRNG sRNG;
|
||||||
pRNG.SeedFixedIntegers(seeds);
|
pRNG.SeedFixedIntegers(seeds);
|
||||||
|
sRNG.SeedFixedIntegers(seeds);
|
||||||
|
|
||||||
typedef PeriodicGimplR Gimpl;
|
typedef PeriodicGimplR Gimpl;
|
||||||
typedef WilsonGaugeAction<Gimpl> GaugeAction;
|
typedef WilsonGaugeAction<Gimpl> GaugeAction;
|
||||||
@ -115,7 +117,7 @@ int main (int argc, char ** argv)
|
|||||||
|
|
||||||
integrator.setMomentumFilter(filter);
|
integrator.setMomentumFilter(filter);
|
||||||
|
|
||||||
integrator.refresh(U, pRNG); //doesn't actually change the gauge field
|
integrator.refresh(U, sRNG, pRNG); //doesn't actually change the gauge field
|
||||||
|
|
||||||
//Check the momentum is zero on the boundary
|
//Check the momentum is zero on the boundary
|
||||||
const auto &P = integrator.getMomentum();
|
const auto &P = integrator.getMomentum();
|
||||||
|
Loading…
Reference in New Issue
Block a user