From: Alexander Schmidt <alex@treefish.org>
Date: Thu, 28 Nov 2013 08:47:10 +0000 (+0100)
Subject: Initializin fields wirh randmom angle.
X-Git-Url: http://git.treefish.org/~alex/phys/heatbath.git/commitdiff_plain/e6bdf191ec43615a0e5d05a5095336150a15417b?hp=96b86c0d9c02d0fcc9b650fb3e5069a8ec013223

Initializin fields wirh randmom angle.
---

diff --git a/sim-metro.hpp b/sim-metro.hpp
index b67f8a7..a91a4b5 100644
--- a/sim-metro.hpp
+++ b/sim-metro.hpp
@@ -58,7 +58,7 @@ void sim::_newParas() {
   *log << "SIM: Resetting fields." << endl << flush;
 
   for (int ix=0; ix<LSIZE2; ix++)
-    conf[ix].phi = 1;
+    conf[ix].phi = polar( 1., gsl_rng_uniform(rangsl) * 2*M_PI );
 }
 
 int sim::updatePhi(const int& x) {
diff --git a/sim-nonmetro.hpp b/sim-nonmetro.hpp
index 7acab4d..b72104f 100644
--- a/sim-nonmetro.hpp
+++ b/sim-nonmetro.hpp
@@ -86,7 +86,7 @@ void sim::_newParas() {
   *log << "SIM: Resetting fields." << endl << flush;
 
   for (int ix=0; ix<LSIZE2; ix++)
-    conf[ix].phi = 1;
+    conf[ix].phi = polar( 1., gsl_rng_uniform(rangsl) * 2*M_PI );
 }
 
 #endif