diff --git a/bfps/cpp/full_code/NSVEcomplex_particles.cpp b/bfps/cpp/full_code/NSVEcomplex_particles.cpp index 4f0ccb506865b74a0d055a276ddff37addc4d1d3..364827cf21241f38ae37d8794868c788da3d7bad 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 db2f33fa3716a8b22c96817a8c2724434c17834a..e451f56a56b00a050d9d940085ad50ccf564ea30 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);