From f06df48595d34cf61f49f88d66e61637ac46113b Mon Sep 17 00:00:00 2001
From: Jose Agustin Arguedas Leiva <arguedas@telica.nld.ds.mpg.de>
Date: Wed, 12 Sep 2018 16:29:32 +0200
Subject: [PATCH] debugging: orientation does not change
---
bfps/cpp/full_code/NSVEcomplex_particles.cpp | 6 +++++-
bfps/cpp/particles/particles_system.hpp | 2 +-
2 files changed, 6 insertions(+), 2 deletions(-)
diff --git a/bfps/cpp/full_code/NSVEcomplex_particles.cpp b/bfps/cpp/full_code/NSVEcomplex_particles.cpp
index 4f0ccb50..364827cf 100644
--- a/bfps/cpp/full_code/NSVEcomplex_particles.cpp
+++ b/bfps/cpp/full_code/NSVEcomplex_particles.cpp
@@ -45,6 +45,7 @@ int NSVEcomplex_particles<rnumber>::initialize(void)
this->enable_inner = false;
current_particles_inner_computer.setEnable(enable_inner);
+
this->ps = particles_system_builder_with_p2p(
this->fs->cvelocity, // (field object)
this->fs->kk, // (kspace object, contains dkx, dky, dkz)
@@ -237,7 +238,10 @@ int NSVEcomplex_particles<rnumber>::read_parameters(void)
this->tracers0_smoothness = hdf5_tools::read_value<int>(parameter_file, "parameters/tracers0_smoothness");
this->enable_p2p = hdf5_tools::read_value<int>(parameter_file, "parameters/tracers0_enable_p2p");
this->enable_inner = hdf5_tools::read_value<int>(parameter_file, "parameters/tracers0_enable_inner");
- this->enable_vorticity_omega = hdf5_tools::read_value<int>(parameter_file, "parameters/tracers0_enable_vorticity_omega");
+ int tval = hdf5_tools::read_value<int>(parameter_file, "parameters/tracers0_enable_vorticity_omega");
+ this->enable_vorticity_omega = tval;
+ DEBUG_MSG("tracers0_enable_vorticity_omega = %d, this->enable_vorticity_omega = %d\n",
+ tval, this->enable_vorticity_omega);
this->cutoff = hdf5_tools::read_value<double>(parameter_file, "parameters/tracers0_cutoff");
this->inner_v0 = hdf5_tools::read_value<double>(parameter_file, "parameters/tracers0_inner_v0");
H5Fclose(parameter_file);
diff --git a/bfps/cpp/particles/particles_system.hpp b/bfps/cpp/particles/particles_system.hpp
index db2f33fa..e451f56a 100644
--- a/bfps/cpp/particles/particles_system.hpp
+++ b/bfps/cpp/particles/particles_system.hpp
@@ -278,7 +278,7 @@ public:
void completeLoopWithVorticity(const real_number dt,
const real_number particle_extra_rhs[]) final {
- TIMEZONE("particles_system::completeLoop");
+ TIMEZONE("particles_system::completeLoopWithVorticity");
compute();
compute_p2p();
compute_particles_inner(particle_extra_rhs);
--
GitLab