From 42e99d9303ba3727862198b250b51729a916b409 Mon Sep 17 00:00:00 2001 From: "Daniel R. Reynolds" Date: Wed, 25 Sep 2024 09:58:31 -0500 Subject: [PATCH] more braces --- .../arkode/CXX_manyvector/ark_sod_lsrk.hpp | 30 +++++++++++++++---- 1 file changed, 25 insertions(+), 5 deletions(-) diff --git a/examples/arkode/CXX_manyvector/ark_sod_lsrk.hpp b/examples/arkode/CXX_manyvector/ark_sod_lsrk.hpp index 7139ad0ac8..014b24bbd7 100644 --- a/examples/arkode/CXX_manyvector/ark_sod_lsrk.hpp +++ b/examples/arkode/CXX_manyvector/ark_sod_lsrk.hpp @@ -148,11 +148,11 @@ class EulerData const sunrealtype* my, const sunrealtype* mz, const sunrealtype* et, const long int& i) { - for (int l = 0; l < STSIZE; l++) this->w1d[l][0] = rho[i - 3 + l]; - for (int l = 0; l < STSIZE; l++) this->w1d[l][1] = mx[i - 3 + l]; - for (int l = 0; l < STSIZE; l++) this->w1d[l][2] = my[i - 3 + l]; - for (int l = 0; l < STSIZE; l++) this->w1d[l][3] = mz[i - 3 + l]; - for (int l = 0; l < STSIZE; l++) this->w1d[l][4] = et[i - 3 + l]; + for (int l = 0; l < STSIZE; l++) { this->w1d[l][0] = rho[i - 3 + l]; } + for (int l = 0; l < STSIZE; l++) { this->w1d[l][1] = mx[i - 3 + l]; } + for (int l = 0; l < STSIZE; l++) { this->w1d[l][2] = my[i - 3 + l]; } + for (int l = 0; l < STSIZE; l++) { this->w1d[l][3] = mz[i - 3 + l]; } + for (int l = 0; l < STSIZE; l++) { this->w1d[l][4] = et[i - 3 + l]; } } // Utility routine to pack 1-dimensional data for locations near the @@ -164,25 +164,45 @@ class EulerData const sunrealtype* et, const long int& i) { for (int l = 0; l < 3; l++) + { this->w1d[l][0] = (i < (3 - l)) ? rho[2 - (i + l)] : rho[i - 3 + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l][1] = (i < (3 - l)) ? mx[2 - (i + l)] : mx[i - 3 + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l][2] = (i < (3 - l)) ? my[2 - (i + l)] : my[i - 3 + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l][3] = (i < (3 - l)) ? mz[2 - (i + l)] : mz[i - 3 + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l][4] = (i < (3 - l)) ? et[2 - (i + l)] : et[i - 3 + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l + 3][0] = (i > (nx - l - 1)) ? rho[i + l - 3] : rho[i + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l + 3][1] = (i > (nx - l - 1)) ? mx[i + l - 3] : mx[i + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l + 3][2] = (i > (nx - l - 1)) ? my[i + l - 3] : my[i + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l + 3][3] = (i > (nx - l - 1)) ? mz[i + l - 3] : mz[i + l]; + } for (int l = 0; l < 3; l++) + { this->w1d[l + 3][4] = (i > (nx - l - 1)) ? et[i + l - 3] : et[i + l]; + } } // Equation of state -- compute and return pressure,