...

[phys/heatbath.git] / sim-1mr-.hpp

diff --git a/sim-1mr-.hpp b/sim-1mr-.hpp
index 00e15beed03479c6404a7a18b6ac96a70f18f66c..01701b959f5651b3cfe7d9d25a7ccbb55b484c03 100644 (file)
--- a/sim-1mr-.hpp
+++ b/sim-1mr-.hpp
@@ -53,10 +53,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;
 }