Simplified rand theta update syntax.

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

diff --git a/sim-1mr+.hpp b/sim-1mr+.hpp
index 8902558d655f75977ec826423df565e092bdb236..4c70b4718c5f6740adcdbafedac14d56b067fcde 100644 (file)
--- a/sim-1mr+.hpp
+++ b/sim-1mr+.hpp
@@ -59,7 +59,7 @@ void sim::updatePhi (const int& x)
                        + 2 * real( conf[x].phi * conj(V) ) );
   
   conf[x].phi = sqrt(std::log( 1./(1-r) )) / sqrt(M)
-    * complex<double>( cos(theta), sin(theta) )
+    * polar(1.0,theta)
     + V / M;
 }