From ba71d3acd4304dc96c6b34df374bf26423b75000 Mon Sep 17 00:00:00 2001
From: Niklas Schnierstein <niklass@antofalla.lmp.ds.mpg.de>
Date: Wed, 13 Nov 2019 10:33:11 +0100
Subject: [PATCH] symmetrize ou field

---
 cpp/full_code/ornstein_uhlenbeck_process.cpp | 2 +-
 cpp/full_code/ou_vorticity_equation.cpp      | 6 +++---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/cpp/full_code/ornstein_uhlenbeck_process.cpp b/cpp/full_code/ornstein_uhlenbeck_process.cpp
index f3955ecc..a73e48e4 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.cpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.cpp
@@ -105,9 +105,9 @@ void ornstein_uhlenbeck_process<rnumber,be>::step(
         }
     });
 
+    this->ou_field->symmetrize();
     this->calc_ou_vorticity();
 
-    // this->ou_field->symmetrize(); ???
 }
 
 // template <class rnumber,field_backend be>
diff --git a/cpp/full_code/ou_vorticity_equation.cpp b/cpp/full_code/ou_vorticity_equation.cpp
index 2b341ca7..f3bfe8c4 100644
--- a/cpp/full_code/ou_vorticity_equation.cpp
+++ b/cpp/full_code/ou_vorticity_equation.cpp
@@ -38,7 +38,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[1], this->v[0]);
+    //this->impose_forcing(this->v[1], this->v[0]);
 
     this->omega_nonlin(1);
     this->kk->CLOOP_K2(
@@ -60,7 +60,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[2], this->v[0]);
+    //this->impose_forcing(this->v[2], this->v[0]);
 
     this->omega_nonlin(2);
     // store old vorticity
@@ -83,7 +83,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[0], this->v[1]);
+    //this->impose_forcing(this->v[0], this->v[1]);
 
     this->kk->template force_divfree<rnumber>(this->cvorticity->get_cdata());
     this->cvorticity->symmetrize();
-- 
GitLab