...

diff --git a/sim-1mr+.hpp b/sim-1mr+.hpp
index 4c70b4718c5f6740adcdbafedac14d56b067fcde..2b141498f03d8e29e610f0939f6536da3ba1ea57 100644 (file)
--- a/sim-1mr+.hpp
+++ b/sim-1mr+.hpp
@@ -44,7 +44,6 @@ sim::sim(o815 *_O815) : o815::sim( _O815,
 void sim::updatePhi (const int& x) 
 {
   complex<double> V=0;
 void sim::updatePhi (const int& x) 
 {
   complex<double> V=0;

-  const double theta = gsl_rng_uniform(rangsl) * 2*M_PI;

 
   for (int nu=0; nu<4; nu++)
     V += conf[ (*nb)[x*4+nu] ].phi;
 
   for (int nu=0; nu<4; nu++)
     V += conf[ (*nb)[x*4+nu] ].phi;


@@ -59,7 +58,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)
                        + 2 * real( conf[x].phi * conj(V) ) );
   
   conf[x].phi = sqrt(std::log( 1./(1-r) )) / sqrt(M)

-    * polar(1.0,theta)
+    * polar( 1.0, arg(conf[x].phi) )

     + V / M;
 }
 
     + V / M;
 }