From: Alexander Schmidt 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?hp=eb8194d0b056a8a616c23532966faafbc3f24be4 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( 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( 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( cos(theta), sin(theta) ) + * polar(1.0, theta) + V / M; }