From c54d2b4aadfb3f45598362fa8b087307911d9224 Mon Sep 17 00:00:00 2001
From: Tobias Baetge <tobae@zardkuh.lmp.ds.mpg.de>
Date: Thu, 19 Dec 2019 17:30:45 +0100
Subject: [PATCH] bugfix: NSVE_Stokes_particles uses correct Loop in step()
---
TurTLE/test/test_turtle_NSVE_Stokes_particles.py | 3 ++-
cpp/full_code/NSVE_Stokes_particles.cpp | 2 +-
2 files changed, 3 insertions(+), 2 deletions(-)
diff --git a/TurTLE/test/test_turtle_NSVE_Stokes_particles.py b/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
index 825ad1c8..9428b3f8 100644
--- a/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
+++ b/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
@@ -63,7 +63,8 @@ def main():
sys.argv[1:])
f = h5py.File('quiescent_nsve_stokes_particles_particles.h5', 'r')
print(f['tracers0/momentum/0'][:])
- print(f['tracers0/momentum/100'][:])
+ print(f['tracers0/momentum/0'][:]-f['tracers0/momentum/100'][:])
+ print(f['tracers0/position/0'][:]-f['tracers0/position/100'][:])
return None
if __name__ == '__main__':
diff --git a/cpp/full_code/NSVE_Stokes_particles.cpp b/cpp/full_code/NSVE_Stokes_particles.cpp
index eb9b4084..4c90dbc2 100644
--- a/cpp/full_code/NSVE_Stokes_particles.cpp
+++ b/cpp/full_code/NSVE_Stokes_particles.cpp
@@ -88,7 +88,7 @@ int NSVE_Stokes_particles<rnumber>::step(void)
TIMEZONE("NSVE_Stokes_particles::step");
this->fs->compute_velocity(this->fs->cvorticity);
this->fs->cvelocity->ift();
- this->ps->completeLoop(this->dt);
+ this->ps->complete2ndOrderLoop(this->dt);
this->NSVE<rnumber>::step();
return EXIT_SUCCESS;
}
--
GitLab