calc gamma via eddy turnover time

cpp/full_code/ornstein_uhlenbeck_process.cpp

@@ -181,7 +181,7 @@ void ornstein_uhlenbeck_process<rnumber,be>::initialize_B()
 template <class rnumber, field_backend be>
 void ornstein_uhlenbeck_process<rnumber, be>::let_converge(void)
 {
-  double time_step = 1./this->gamma(1)/1000.
+  double time_step = 1./this->gamma(1)/1000.;
   for(int i=0; i<5000; i++)
   {
     this->step(time_step);

cpp/full_code/ornstein_uhlenbeck_process.hpp

@@ -59,7 +59,7 @@ class ornstein_uhlenbeck_process{
             //return pow(this->ou_energy_amplitude/this->kolmogorov_constant,-0.5)
             //        *pow(kabs,-2.0/3.0);
             
-            return this->ou_gamma*pow(kabs,-2/3);
+            return this->ou_gamma*pow(this->epsilon,1./3.)*pow(kabs,2./3.);
                   
         }