Skip to content
Snippets Groups Projects
Commit e57ae1ce authored by Axel Kohlmeyer's avatar Axel Kohlmeyer
Browse files

compute scaled kinetic energy tensor without destroying the original data

parent 1c68e42e
No related branches found
No related tags found
No related merge requests found
...@@ -120,14 +120,14 @@ void ComputePressureGrem::compute_vector() ...@@ -120,14 +120,14 @@ void ComputePressureGrem::compute_vector()
// invoke temperature if it hasn't been already // invoke temperature if it hasn't been already
double *ke_tensor; double ke_tensor[6];
if (keflag) { if (keflag) {
if (temperature->invoked_vector != update->ntimestep) if (temperature->invoked_vector != update->ntimestep)
temperature->compute_vector(); temperature->compute_vector();
ke_tensor = temperature->vector; ke_tensor = temperature->vector;
for (int i = 0; i < 6; ++i)
ke_tensor[i] = temperature->vector[i] / (*scale_grem);
} }
for (int i = 0; i < 6; i++)
ke_tensor[i] /= *scale_grem;
if (dimension == 3) { if (dimension == 3) {
inv_volume = 1.0 / (domain->xprd * domain->yprd * domain->zprd); inv_volume = 1.0 / (domain->xprd * domain->yprd * domain->zprd);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment