From: Alexander Schmidt <alex@treefish.org>
Date: Mon, 25 Nov 2013 12:29:01 +0000 (+0100)
Subject: Simplified rand theta update syntax.
X-Git-Url: http://git.treefish.org/~alex/phys/heatbath.git/commitdiff_plain/4b4dcf9a4b5a913c1574a9cc1aa479f1aa35a003

Simplified rand theta update syntax.
---

diff --git a/sim-1mr+.hpp b/sim-1mr+.hpp
index 8902558..4c70b47 100644
--- 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 beadd7d..d03ef22 100644
--- 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 0cbce5e..c3b6888 100644
--- 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;
 }