@@ -146,17 +146,17 @@ void Spatial_Derivative::Spatial_Derivative_WCNS_X()
146
146
{
147
147
for (int iVar = 0 ; iVar < num_of_prim_vars; iVar++)
148
148
{
149
- rhs[1 ][j][iVar] = (-11.0 /12 * fluxVector[0 ][j][iVar] + 17.0 /24 * fluxVector[1 ][j][iVar]
149
+ rhs[1 ][j][iVar] =- (-11.0 /12 * fluxVector[0 ][j][iVar] + 17.0 /24 * fluxVector[1 ][j][iVar]
150
150
+ 3.0 /8 * fluxVector[2 ][j][iVar] - 5.0 /24 * fluxVector[3 ][j][iVar]
151
151
+ 1.0 /24 * fluxVector[4 ][j][iVar]) / ds;
152
152
153
- rhs[2 ][j][iVar] = (- 1.0 /24 * fluxVector[0 ][j][iVar] - 9.0 /8 * fluxVector[1 ][j][iVar]
153
+ rhs[2 ][j][iVar] =- (- 1.0 /24 * fluxVector[0 ][j][iVar] - 9.0 /8 * fluxVector[1 ][j][iVar]
154
154
+ 9.0 /8 * fluxVector[2 ][j][iVar] - 1.0 /24 * fluxVector[3 ][j][iVar]) / ds;
155
155
156
- rhs[ied ][j][iVar] = (- 1.0 /24 * fluxVector[ied + 1 ][j][iVar] - 9.0 /8 * fluxVector[ied ][j][iVar]
156
+ rhs[ied ][j][iVar] =- (- 1.0 /24 * fluxVector[ied + 1 ][j][iVar] - 9.0 /8 * fluxVector[ied ][j][iVar]
157
157
+ 9.0 /8 * fluxVector[ied - 1 ][j][iVar] - 1.0 /24 * fluxVector[ied - 2 ][j][iVar]) / ds;
158
158
159
- rhs[ied + 1 ][j][iVar] = (-11.0 /12 * fluxVector[ied + 1 ][j][iVar] + 17.0 /24 * fluxVector[ied ][j][iVar]
159
+ rhs[ied + 1 ][j][iVar] =- (-11.0 /12 * fluxVector[ied + 1 ][j][iVar] + 17.0 /24 * fluxVector[ied ][j][iVar]
160
160
+ 3.0 /8 * fluxVector[ied - 1 ][j][iVar] - 5.0 /24 * fluxVector[ied - 2 ][j][iVar]
161
161
+ 1.0 /24 * fluxVector[ied - 3 ][j][iVar]) / ds;
162
162
}
@@ -192,17 +192,17 @@ void Spatial_Derivative::Spatial_Derivative_WCNS_Y()
192
192
{
193
193
for (int iVar = 0 ; iVar < num_of_prim_vars; iVar++)
194
194
{
195
- rhs[i][1 ][iVar] = (-11.0 /12 * fluxVector[i][0 ][iVar] + 17.0 /24 * fluxVector[i][1 ][iVar]
195
+ rhs[i][1 ][iVar] =- (-11.0 /12 * fluxVector[i][0 ][iVar] + 17.0 /24 * fluxVector[i][1 ][iVar]
196
196
+ 3.0 /8 * fluxVector[i][2 ][iVar] - 5.0 /24 * fluxVector[i][3 ][iVar]
197
197
+ 1.0 /24 * fluxVector[i][4 ][iVar]) / ds;
198
198
199
- rhs[i][2 ][iVar] = (- 1.0 /24 * fluxVector[i][0 ][iVar] - 9.0 /8 * fluxVector[i][1 ][iVar]
199
+ rhs[i][2 ][iVar] =- (- 1.0 /24 * fluxVector[i][0 ][iVar] - 9.0 /8 * fluxVector[i][1 ][iVar]
200
200
+ 9.0 /8 * fluxVector[i][2 ][iVar] - 1.0 /24 * fluxVector[i][3 ][iVar]) / ds;
201
201
202
- rhs[i][jed ][iVar] = (- 1.0 /24 * fluxVector[i][jed + 1 ][iVar] - 9.0 /8 * fluxVector[i][jed ][iVar]
202
+ rhs[i][jed ][iVar] =- (- 1.0 /24 * fluxVector[i][jed + 1 ][iVar] - 9.0 /8 * fluxVector[i][jed ][iVar]
203
203
+ 9.0 /8 * fluxVector[i][jed - 1 ][iVar] - 1.0 /24 * fluxVector[i][jed - 2 ][iVar]) / ds;
204
204
205
- rhs[i][jed + 1 ][iVar] = (-11.0 /12 * fluxVector[i][jed + 1 ][iVar] + 17.0 /24 * fluxVector[i][jed ][iVar]
205
+ rhs[i][jed + 1 ][iVar] =- (-11.0 /12 * fluxVector[i][jed + 1 ][iVar] + 17.0 /24 * fluxVector[i][jed ][iVar]
206
206
+ 3.0 /8 * fluxVector[i][jed - 1 ][iVar] - 5.0 /24 * fluxVector[i][jed - 2 ][iVar]
207
207
+ 1.0 /24 * fluxVector[i][jed - 3 ][iVar]) / ds;
208
208
}
0 commit comments