From e57ae1ce3f7e4be96e972cd1ab9f83f667e23b25 Mon Sep 17 00:00:00 2001
From: Axel Kohlmeyer <akohlmey@gmail.com>
Date: Wed, 16 Nov 2016 12:45:13 -0500
Subject: [PATCH] compute scaled kinetic energy tensor without destroying the
original data
---
src/USER-MISC/compute_pressure_grem.cpp | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/USER-MISC/compute_pressure_grem.cpp b/src/USER-MISC/compute_pressure_grem.cpp
index bcacb43cd2..c00e5ebdfd 100644
--- a/src/USER-MISC/compute_pressure_grem.cpp
+++ b/src/USER-MISC/compute_pressure_grem.cpp
@@ -120,14 +120,14 @@ void ComputePressureGrem::compute_vector()
// invoke temperature if it hasn't been already
- double *ke_tensor;
+ double ke_tensor[6];
if (keflag) {
if (temperature->invoked_vector != update->ntimestep)
temperature->compute_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) {
inv_volume = 1.0 / (domain->xprd * domain->yprd * domain->zprd);
--
GitLab