From c99209369fbdb86ac26cd05b5656d6ab5dee531f Mon Sep 17 00:00:00 2001 From: Cristian C Lalescu <Cristian.Lalescu@ds.mpg.de> Date: Thu, 12 Nov 2015 17:32:01 +0100 Subject: [PATCH] fix call to rspace_stats --- bfps/NavierStokes.py | 13 +++++++++++-- 1 file changed, 11 insertions(+), 2 deletions(-) diff --git a/bfps/NavierStokes.py b/bfps/NavierStokes.py index 1279e7e4..9f08dde7 100644 --- a/bfps/NavierStokes.py +++ b/bfps/NavierStokes.py @@ -102,6 +102,7 @@ class NavierStokes(bfps.fluid_base.fluid_particle_base): double *vorticity_moments = fftw_alloc_real(10*4); ptrdiff_t *hist_velocity = new ptrdiff_t[histogram_bins*4]; ptrdiff_t *hist_vorticity = new ptrdiff_t[histogram_bins*4]; + double max_estimates[4]; {0} *ksample_values; fs->compute_velocity(fs->cvorticity); if (fs->cd->myrank == 0) @@ -120,16 +121,24 @@ class NavierStokes(bfps.fluid_base.fluid_particle_base): fs->cospectrum(fs->cvelocity, fs->cvelocity, spec_velocity); fs->cospectrum(fs->cvorticity, fs->cvorticity, spec_vorticity); fs->ift_velocity(); + max_estimates[0] = max_velocity_estimate/sqrt(3); + max_estimates[1] = max_estimates[0]; + max_estimates[2] = max_estimates[0]; + max_estimates[3] = max_velocity_estimate; fs->compute_rspace_stats(fs->rvelocity, velocity_moments, hist_velocity, - max_velocity_estimate, + max_estimates, histogram_bins); fs->ift_vorticity(); + max_estimates[0] = max_velocity_estimate/sqrt(3); + max_estimates[1] = max_estimates[0]; + max_estimates[2] = max_estimates[0]; + max_estimates[3] = max_velocity_estimate; fs->compute_rspace_stats(fs->rvorticity, vorticity_moments, hist_vorticity, - max_vorticity_estimate, + max_estimates, histogram_bins); if (fs->cd->myrank == 0) {{ -- GitLab