X-Git-Url: http://git.treefish.org/~alex/phys/heatbath.git/blobdiff_plain/dca762861c66d8782be31811822812d347ee4577..274a66cdb22fd36ea8cbbc373783d848f1ba253e:/sim-1mr-.hpp

diff --git a/sim-1mr-.hpp b/sim-1mr-.hpp
index 00e15be..a7cde92 100644
--- a/sim-1mr-.hpp
+++ b/sim-1mr-.hpp
@@ -43,7 +43,6 @@ sim::sim(o815 *_O815) : o815::sim( _O815,
 
 void sim::updatePhi (const int& x) 
 {
-  const double theta = gsl_rng_uniform(rangsl) * 2*M_PI;
   complex<double> V=0;
 
   for (int nu=0; nu<4; nu++)
@@ -53,10 +52,8 @@ void sim::updatePhi (const int& x)
 			- 1./M * norm(V)
 			+ 2 * real( conf[x].phi * conj(V) ) );
 
-  double oldarg = arg(conf[x].phi);
-
   conf[x].phi = sqrt(std::log( 1./(1-r) )) / sqrt(M)
-    * polar(1.0, oldarg)
+    * polar( 1.0, arg(conf[x].phi) )
     + V / M;
 }