diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-05-12 13:35:15 -0500 |
---|---|---|
committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2008-05-12 13:35:15 -0500 |
commit | efe228d96b57a11e9cd1d922b1f50703d1301853 (patch) | |
tree | adc1fa31300223e9512d82cdf52c70da30f4459a /ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | |
parent | c5fb2c853ffb3c9a5461c1807a835b73ecbcbd83 (diff) |
Update auto-generated code
Diffstat (limited to 'ML_BSSN_MP/src/ML_BSSN_MP_enforce.c')
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | 14 |
1 files changed, 13 insertions, 1 deletions
diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c index da4b3c4..fae01b8 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c @@ -23,7 +23,7 @@ #define CUB(x) ((x) * (x) * (x)) #define QAD(x) ((x) * (x) * (x) * (x)) -void ML_BSSN_MP_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[3], CCTK_REAL tangentA[3], CCTK_REAL tangentB[3], CCTK_INT min[3], CCTK_INT max[3], CCTK_INT n_subblock_gfs, CCTK_REAL *subblock_gfs[]) +void ML_BSSN_MP_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT const face, CCTK_REAL const normal[3], CCTK_REAL const tangentA[3], CCTK_REAL const tangentB[3], CCTK_INT const min[3], CCTK_INT const max[3], CCTK_INT const n_subblock_gfs, CCTK_REAL * const subblock_gfs[]) { DECLARE_CCTK_ARGUMENTS DECLARE_CCTK_PARAMETERS @@ -43,9 +43,15 @@ void ML_BSSN_MP_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL CCTK_REAL p1o144dxdy = INITVALUE; CCTK_REAL p1o144dxdz = INITVALUE; CCTK_REAL p1o144dydz = INITVALUE; + CCTK_REAL p1o2dx = INITVALUE; + CCTK_REAL p1o2dy = INITVALUE; + CCTK_REAL p1o2dz = INITVALUE; CCTK_REAL pm1o12dx2 = INITVALUE; CCTK_REAL pm1o12dy2 = INITVALUE; CCTK_REAL pm1o12dz2 = INITVALUE; + CCTK_REAL pm1o2dx = INITVALUE; + CCTK_REAL pm1o2dy = INITVALUE; + CCTK_REAL pm1o2dz = INITVALUE; if (verbose > 1) { @@ -82,9 +88,15 @@ void ML_BSSN_MP_enforce_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL p1o144dxdy = (INV(dx)*INV(dy))/144.; p1o144dxdz = (INV(dx)*INV(dz))/144.; p1o144dydz = (INV(dy)*INV(dz))/144.; + p1o2dx = khalf*INV(dx); + p1o2dy = khalf*INV(dy); + p1o2dz = khalf*INV(dz); pm1o12dx2 = -pow(dx,-2)/12.; pm1o12dy2 = -pow(dy,-2)/12.; pm1o12dz2 = -pow(dz,-2)/12.; + pm1o2dx = -(khalf*INV(dx)); + pm1o2dy = -(khalf*INV(dy)); + pm1o2dz = -(khalf*INV(dz)); /* Loop over the grid points */ #pragma omp parallel |