diff --git a/src/USER-MISC/compute_pressure_grem.cpp b/src/USER-MISC/compute_pressure_grem.cpp index bcacb43cd2f069bee4c09a1594b2a3e6bb2fe76d..c00e5ebdfdafe5d31db4b61fe1e0f4afe3d52b13 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);