X-Git-Url: http://git.treefish.org/~alex/phys/u1casc.git/blobdiff_plain/7664dbeaba693ea47d5fd405c20f5db109d4da13..a0c10794ee1276a61a221a398f96df89658149dd:/u1casc-flux/sim.hpp diff --git a/u1casc-flux/sim.hpp b/u1casc-flux/sim.hpp index 0da2b8a..496a28d 100644 --- a/u1casc-flux/sim.hpp +++ b/u1casc-flux/sim.hpp @@ -22,19 +22,20 @@ public: double kappa[2], lambda[2], Mu[2], beta; int wArg(int x, int i); double WF(int nom, int denom, int iflav); - + int iPlaq(int i, int j); + double I(int p); + void _newParas(); + private: + void _resetConfig(); void _makeSweep(); - void _newParas(); gsl_rng* rangsl; int lineSweep(); int lineUpdate(int x, int vardir); int cubeUpdate(int x, int orient); void getCube(int cube[3], int orient); int plaqUpdate(int x, int mu, int nu, int i); - int iPlaq(int i, int j); int lpUpdate(int x, int mu, int i); - double I(int p); static double wkern(double t, void *params); double wGsl(int n, int iflav); struct bcache { @@ -92,7 +93,7 @@ sim::sim(o815 *_O815) : o815::sim( _O815, } void sim::_makeSweep() -{ +{ lineSweep(); for( int x=0; xparaQ)["beta"]; Mu[0] = (*O815->paraQ)["muone"]; Mu[1] = (*O815->paraQ)["mutwo"]; +} - memset(confMem, 0, sizeof(int)*lsize4*(2*4+2*4+6)); +void sim::_resetConfig() { + memset(confMem, 0, confmemSize); } int sim::lineSweep()