Initializin fields wirh randmom angle.

[phys/heatbath.git] / sim-metro.hpp

diff --git a/sim-metro.hpp b/sim-metro.hpp
index b67f8a783b5d2c4e1b7f33107ddfcadcbb245d8e..a91a4b562614c6e5d72c1beebf21192220b9252d 100644 (file)
--- 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) {