diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-12-27 17:55:56 +0100 |
---|---|---|
committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-12-27 17:55:56 +0100 |
commit | f5a7c761dafe64caf81cf007f0e8bff5475401b9 (patch) | |
tree | 71d42353c1695318d8276da110c7d9860d73d16b /ML_ADM/src/ML_ADM_RHS.c | |
parent | a66116f18b5ae7e373b9e00cdd2acff8a8e6afbf (diff) |
Use only 2 time levels for evolved variables
Diffstat (limited to 'ML_ADM/src/ML_ADM_RHS.c')
-rw-r--r-- | ML_ADM/src/ML_ADM_RHS.c | 36 |
1 files changed, 18 insertions, 18 deletions
diff --git a/ML_ADM/src/ML_ADM_RHS.c b/ML_ADM/src/ML_ADM_RHS.c index 397e892..57b55c6 100644 --- a/ML_ADM/src/ML_ADM_RHS.c +++ b/ML_ADM/src/ML_ADM_RHS.c @@ -34,15 +34,15 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[ /* Declare predefined quantities */ - CCTK_REAL p1o12dx = INITVALUE; - CCTK_REAL p1o12dy = INITVALUE; - CCTK_REAL p1o12dz = INITVALUE; - CCTK_REAL p1o144dxdy = INITVALUE; - CCTK_REAL p1o144dxdz = INITVALUE; - CCTK_REAL p1o144dydz = INITVALUE; - CCTK_REAL pm1o12dx2 = INITVALUE; - CCTK_REAL pm1o12dy2 = INITVALUE; - CCTK_REAL pm1o12dz2 = INITVALUE; + CCTK_REAL qp1o12dx = INITVALUE; + CCTK_REAL qp1o12dy = INITVALUE; + CCTK_REAL qp1o12dz = INITVALUE; + CCTK_REAL qp1o144dxdy = INITVALUE; + CCTK_REAL qp1o144dxdz = INITVALUE; + CCTK_REAL qp1o144dydz = INITVALUE; + CCTK_REAL qpm1o12dx2 = INITVALUE; + CCTK_REAL qpm1o12dy2 = INITVALUE; + CCTK_REAL qpm1o12dz2 = INITVALUE; if (verbose > 1) { @@ -73,15 +73,15 @@ void ML_ADM_RHS_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, CCTK_REAL normal[ hdzi = 0.5 * dzi; /* Initialize predefined quantities */ - p1o12dx = INV(dx)/12.; - p1o12dy = INV(dy)/12.; - p1o12dz = INV(dz)/12.; - p1o144dxdy = (INV(dx)*INV(dy))/144.; - p1o144dxdz = (INV(dx)*INV(dz))/144.; - p1o144dydz = (INV(dy)*INV(dz))/144.; - pm1o12dx2 = -pow(dx,-2)/12.; - pm1o12dy2 = -pow(dy,-2)/12.; - pm1o12dz2 = -pow(dz,-2)/12.; + qp1o12dx = INV(dx)/12.; + qp1o12dy = INV(dy)/12.; + qp1o12dz = INV(dz)/12.; + qp1o144dxdy = (INV(dx)*INV(dy))/144.; + qp1o144dxdz = (INV(dx)*INV(dz))/144.; + qp1o144dydz = (INV(dy)*INV(dz))/144.; + qpm1o12dx2 = -pow(dx,-2)/12.; + qpm1o12dy2 = -pow(dy,-2)/12.; + qpm1o12dz2 = -pow(dz,-2)/12.; /* Loop over the grid points */ LC_LOOP3 (ML_ADM_RHS, |