From 47bba94dcd9e7894472d59292738a9aae35cc49b Mon Sep 17 00:00:00 2001
From: sniklas142 <niklas.schnierstein@googlemail.com>
Date: Fri, 23 Aug 2019 13:31:19 +0200
Subject: [PATCH] added func to add ou modes to modes of field
---
cpp/full_code/ornstein_uhlenbeck_process.cpp | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/cpp/full_code/ornstein_uhlenbeck_process.cpp b/cpp/full_code/ornstein_uhlenbeck_process.cpp
index 42cdfdcb..101093e9 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.cpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.cpp
@@ -174,9 +174,6 @@ void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_replace(
src->cval(cindex,1,1) = this->ou_field->cval(cindex,1,1);
src->cval(cindex,2,1) = this->ou_field->cval(cindex,2,1);
}
- // else {
- // assert(src->cval(cindex,0,0) == this->ou_field->cval(cindex,0,0));
- // }
}
@@ -197,7 +194,13 @@ void ornstein_uhlenbeck_process<rnumber,be>::add_to_field_sharp(
if (k2 <= this->ou_kmax_squ && k2 >= this->ou_kmin_squ)
{
+ src->cval(cindex,0,0) += this->ou_field->cval(cindex,0,0);
+ src->cval(cindex,1,0) += this->ou_field->cval(cindex,1,0);
+ src->cval(cindex,2,0) += this->ou_field->cval(cindex,2,0);
+ src->cval(cindex,0,1) += this->ou_field->cval(cindex,0,1);
+ src->cval(cindex,1,1) += this->ou_field->cval(cindex,1,1);
+ src->cval(cindex,2,1) += this->ou_field->cval(cindex,2,1);
}
}
--
GitLab