Initializin fields wirh randmom angle.

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++)
   *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) {
 }
 
 int sim::updatePhi(const int& x) {

diff --git a/sim-nonmetro.hpp b/sim-nonmetro.hpp
index 7acab4df05de57ce0384990f30f1223038449a38..b72104f8342803be998034a0bc9a781b1e28e275 100644 (file)
--- 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++)
   *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
 }
 
 #endif