Fixed wrong 1mr+ update.

[phys/heatbath.git] / sim-r.hpp

diff --git a/sim-r.hpp b/sim-r.hpp
index 0cbce5ec4dba439fac5b031785e15c26a4430699..c3b68889ffff0684c094de5957ab7f46c2963203 100644 (file)
--- a/sim-r.hpp
+++ b/sim-r.hpp
@@ -51,7 +51,7 @@ void sim::updatePhi (const int& x)
     V += conf[ (*nb)[x*4+nu] ].phi;
 
   conf[x].phi = sqrt(std::log( 1./(1-r) )) / sqrt(M)
-    * complex<double>( cos(theta), sin(theta) )
+    * polar(1.0, theta)
     + V / M;
 }