Still one bug

This commit is contained in:
Miikka Vaisala
2019-08-07 19:10:39 +08:00
parent 7fdbd76aa2
commit 0bb568642f

View File

@@ -335,11 +335,13 @@ compute_gradient(const int i, const int j, const int k, const ModelScalar* arr)
return (ModelVector){derx(i, j, k, arr), dery(i, j, k, arr), derz(i, j, k, arr)}; return (ModelVector){derx(i, j, k, arr), dery(i, j, k, arr), derz(i, j, k, arr)};
} }
#if LUPWD
static inline ModelVector static inline ModelVector
compute_upwind(const int i, const int j, const int k, const ModelScalar* arr) compute_upwind(const int i, const int j, const int k, const ModelScalar* arr)
{ {
return (ModelVector){der6x_upwd(i, j, k, arr), der6y_upwd(i, j, k, arr), der6z_upwd(i, j, k, arr)}; return (ModelVector){der6x_upwd(i, j, k, arr), der6y_upwd(i, j, k, arr), der6z_upwd(i, j, k, arr)};
} }
#endif
static inline ModelMatrix static inline ModelMatrix
compute_hessian(const int i, const int j, const int k, const ModelScalar* arr) compute_hessian(const int i, const int j, const int k, const ModelScalar* arr)