From b98974a8b3b120c6530e005e50fcb3cc5ec6ad96 Mon Sep 17 00:00:00 2001
From: Niklas Schnierstein <niklass@antofalla.lmp.ds.mpg.de>
Date: Tue, 3 Dec 2019 11:05:39 +0100
Subject: [PATCH] calc gamma via eddy turnover time
---
cpp/full_code/ornstein_uhlenbeck_process.cpp | 2 +-
cpp/full_code/ornstein_uhlenbeck_process.hpp | 2 +-
2 files changed, 2 insertions(+), 2 deletions(-)
diff --git a/cpp/full_code/ornstein_uhlenbeck_process.cpp b/cpp/full_code/ornstein_uhlenbeck_process.cpp
index a3c0285b..3a095047 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.cpp
+++ b/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);
diff --git a/cpp/full_code/ornstein_uhlenbeck_process.hpp b/cpp/full_code/ornstein_uhlenbeck_process.hpp
index 6c6bfb33..4987fe2e 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.hpp
+++ b/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.);
}
--
GitLab