Initializin fields wirh randmom angle.

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

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