Simplified rand theta update syntax.

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;
 }
 

diff --git a/sim-1mr.hpp b/sim-1mr.hpp
index beadd7d55bdcc626f735110277899a05c6d15761..d03ef224b33e735201abf9771b9da489b7d97cb1 100644 (file)
--- a/sim-1mr.hpp
+++ b/sim-1mr.hpp
@@ -54,7 +54,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;
 }
 

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;
 }