Skip to content
Snippets Groups Projects
Commit 805bfda2 authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

compute spectra for joint acc vel stats as well

parent d9a4bea9
Branches
Tags
No related merge requests found
Pipeline #51428 passed
......@@ -154,17 +154,6 @@ int joint_acc_vel_stats<rnumber>::work_on_current_iteration(void)
max_acc_estimate[3] = max_acceleration_estimate;
max_vel_estimate[3] = max_velocity_estimate;
acc->compute_rspace_stats(
stat_group,
"acceleration",
this->iteration / this->niter_out,
max_acc_estimate);
vel->compute_rspace_stats(
stat_group,
"velocity",
this->iteration / this->niter_out,
max_vel_estimate);
/// compute joint PDF
joint_rspace_PDF(
acc, vel,
......@@ -174,6 +163,20 @@ int joint_acc_vel_stats<rnumber>::work_on_current_iteration(void)
max_acc_estimate,
max_vel_estimate);
/// compute real space stats and spectra
acc->compute_stats(
kk,
stat_group,
"acceleration",
this->iteration / this->niter_out,
max_acceleration_estimate / sqrt(3));
vel->compute_stats(
kk,
stat_group,
"velocity",
this->iteration / this->niter_out,
max_velocity_estimate / sqrt(3));
delete vel;
return EXIT_SUCCESS;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment