working

cpp/full_code/ornstein_uhlenbeck_process.hpp

@@ -32,7 +32,7 @@ class ornstein_uhlenbeck_process{
         std::vector<double> aom; //absolute k value of modes
 
         std::random_device rd{};
-        std::mt19937 gen{rd()};
+        std::mt19937 gen{rd()}; //32-bit, maybe fix for double prec TODO
         std::normal_distribution<double> d{0,1};
 
         ornstein_uhlenbeck_process(
@@ -58,7 +58,7 @@ class ornstein_uhlenbeck_process{
         inline double gamma(double kabs)
         {
             //return kabs*sqrt(this->energy(kabs));
-            return 0.3333;
+            return kolmogorov_constant*sqrt(kolmogorov_constant/this->ou_energy_amplitude);
         }