cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Fixed a float/long double conversion that some complainers... um compilers complain about

Browse Source
This commit is contained in:
jpekkila
2019-08-08 18:12:24 +03:00
parent af24f17d1f
commit 168cdc9109
1 changed files with 37 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files

51
src/standalone/model/model_rk3.cc
View File

@@ -282,14 +282,11 @@ der6x_upwd(const int i, const int j, const int k, const ModelScalar* arr)
{ {
ModelScalar inv_ds = get(AC_inv_dsx); ModelScalar inv_ds = get(AC_inv_dsx);
return ModelScalar(1.0/60.0)*inv_ds* ( return ModelScalar(1.0 / 60.0) * inv_ds *
-ModelScalar( 20.0)* arr[IDX(i, j, k)] (-ModelScalar(20.0) * arr[IDX(i, j, k)] +
+ModelScalar( 15.0)*(arr[IDX(i+1, j, k)] ModelScalar(15.0) * (arr[IDX(i + 1, j, k)] + arr[IDX(i - 1, j, k)]) -
+ arr[IDX(i-1, j, k)]) ModelScalar(6.0) * (arr[IDX(i + 2, j, k)] + arr[IDX(i - 2, j, k)]) +
-ModelScalar( 6.0)*(arr[IDX(i+2, j, k)] arr[IDX(i + 3, j, k)] + arr[IDX(i - 3, j, k)]);
+ arr[IDX(i-2, j, k)])
+ arr[IDX(i+3, j, k)]
+ arr[IDX(i-3, j, k)]);
} }
static inline ModelScalar static inline ModelScalar
@@ -297,14 +294,11 @@ der6y_upwd(const int i, const int j, const int k, const ModelScalar* arr)
{ {
ModelScalar inv_ds = get(AC_inv_dsy); ModelScalar inv_ds = get(AC_inv_dsy);
return ModelScalar(1.0/60.0)*inv_ds* ( return ModelScalar(1.0 / 60.0) * inv_ds *
-ModelScalar( 20.0)* arr[IDX(i, j, k)] (-ModelScalar(20.0) * arr[IDX(i, j, k)] +
+ModelScalar( 15.0)*(arr[IDX(i, j+1, k)] ModelScalar(15.0) * (arr[IDX(i, j + 1, k)] + arr[IDX(i, j - 1, k)]) -
+ arr[IDX(i, j-1, k)]) ModelScalar(6.0) * (arr[IDX(i, j + 2, k)] + arr[IDX(i, j - 2, k)]) +
-ModelScalar( 6.0)*(arr[IDX(i, j+2, k)] arr[IDX(i, j + 3, k)] + arr[IDX(i, j - 3, k)]);
+ arr[IDX(i, j-2, k)])
+ arr[IDX(i, j+3, k)]
+ arr[IDX(i, j-3, k)]);
} }
static inline ModelScalar static inline ModelScalar
@@ -312,14 +306,11 @@ der6z_upwd(const int i, const int j, const int k, const ModelScalar* arr)
{ {
ModelScalar inv_ds = get(AC_inv_dsz); ModelScalar inv_ds = get(AC_inv_dsz);
return ModelScalar(1.0/60.0)*inv_ds* ( return ModelScalar(1.0 / 60.0) * inv_ds *
-ModelScalar( 20.0)* arr[IDX(i, j, k )] (-ModelScalar(20.0) * arr[IDX(i, j, k)] +
+ModelScalar( 15.0)*(arr[IDX(i, j, k+1)] ModelScalar(15.0) * (arr[IDX(i, j, k + 1)] + arr[IDX(i, j, k - 1)]) -
+ arr[IDX(i, j, k-1)]) ModelScalar(6.0) * (arr[IDX(i, j, k + 2)] + arr[IDX(i, j, k - 2)]) +
-ModelScalar( 6.0)*(arr[IDX(i, j, k+2)] arr[IDX(i, j, k + 3)] + arr[IDX(i, j, k - 3)]);
+ arr[IDX(i, j, k-2)])
+ arr[IDX(i, j, k+3)]
+ arr[IDX(i, j, k-3)]);
} }
#endif #endif
@@ -339,7 +330,8 @@ compute_gradient(const int i, const int j, const int k, const ModelScalar* arr)
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 #endif
@@ -416,8 +408,6 @@ gradients(const ModelVectorData& data)
return (ModelMatrix){gradient(data.x), gradient(data.y), gradient(data.z)}; return (ModelMatrix){gradient(data.x), gradient(data.y), gradient(data.z)};
} }
/* /*
* ============================================================================= * =============================================================================
* Level 0.3 (Built-in functions available during the Stencil Processing Stage) * Level 0.3 (Built-in functions available during the Stencil Processing Stage)
@@ -571,9 +561,9 @@ contract(const ModelMatrix& mat)
ModelScalar ModelScalar
upwd_der6(const ModelVectorData& uu, const ModelScalarData& lnrho) upwd_der6(const ModelVectorData& uu, const ModelScalarData& lnrho)
{ {
ModelScalar uux = fabs(value(uu).x); ModelScalar uux = fabsl(value(uu).x);
ModelScalar uuy = fabs(value(uu).y); ModelScalar uuy = fabsl(value(uu).y);
ModelScalar uuz = fabs(value(uu).z); ModelScalar uuz = fabsl(value(uu).z);
return uux * lnrho.upwind.x + uuy * lnrho.upwind.y + uuz * lnrho.upwind.z; return uux * lnrho.upwind.x + uuy * lnrho.upwind.y + uuz * lnrho.upwind.z;
} }
#endif #endif
@@ -774,7 +764,6 @@ helical_forcing(ModelScalar magnitude, ModelVector k_force, ModelVector xx, Mode
ModelScalar real_comp_phase = cos_k_dot_x * cos_phi - sin_k_dot_x * sin_phi; ModelScalar real_comp_phase = cos_k_dot_x * cos_phi - sin_k_dot_x * sin_phi;
ModelScalar imag_comp_phase = cos_k_dot_x * sin_phi + sin_k_dot_x * cos_phi; ModelScalar imag_comp_phase = cos_k_dot_x * sin_phi + sin_k_dot_x * cos_phi;
ModelVector force = (ModelVector){ff_re.x * real_comp_phase - ff_im.x * imag_comp_phase, ModelVector force = (ModelVector){ff_re.x * real_comp_phase - ff_im.x * imag_comp_phase,
ff_re.y * real_comp_phase - ff_im.y * imag_comp_phase, ff_re.y * real_comp_phase - ff_im.y * imag_comp_phase,
ff_re.z * real_comp_phase - ff_im.z * imag_comp_phase}; ff_re.z * real_comp_phase - ff_im.z * imag_comp_phase};