diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_matter_constraints.c')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_matter_constraints.c | 36 |
1 files changed, 18 insertions, 18 deletions
diff --git a/ML_BSSN/src/ML_BSSN_matter_constraints.c b/ML_BSSN/src/ML_BSSN_matter_constraints.c index 9874d39..b8b41a3 100644 --- a/ML_BSSN/src/ML_BSSN_matter_constraints.c +++ b/ML_BSSN/src/ML_BSSN_matter_constraints.c @@ -34,15 +34,15 @@ void ML_BSSN_matter_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, C /* 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_BSSN_matter_constraints_Body(cGH *cctkGH, CCTK_INT dir, CCTK_INT face, C 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_BSSN_matter_constraints, |