diff options
| author | Erik Schnetter <schnetter@cct.lsu.edu> | 2009-04-28 13:38:57 -0500 |
|---|---|---|
| committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2009-04-28 13:38:57 -0500 |
| commit | 10ee905e681cedb97d8a2e4d4f7b4445d74ee7bc (patch) | |
| tree | 21cf8b52fb83a88acf0b72088244a8e881092560 /ML_BSSN/src/ML_BSSN_enforce.c | |
| parent | 51a068080ddd36d89943b754c46fd7b0e265c1ff (diff) | |
Regenerate code
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_enforce.c')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_enforce.c | 35 |
1 files changed, 13 insertions, 22 deletions
diff --git a/ML_BSSN/src/ML_BSSN_enforce.c b/ML_BSSN/src/ML_BSSN_enforce.c index 2b02dd2..79cd691 100644 --- a/ML_BSSN/src/ML_BSSN_enforce.c +++ b/ML_BSSN/src/ML_BSSN_enforce.c @@ -43,15 +43,12 @@ void ML_BSSN_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT 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 p1odx = INITVALUE; + CCTK_REAL p1ody = INITVALUE; + CCTK_REAL p1odz = 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) { @@ -88,15 +85,12 @@ void ML_BSSN_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT 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); + p1odx = INV(dx); + p1ody = INV(dy); + p1odz = 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 @@ -146,22 +140,19 @@ void ML_BSSN_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT /* Calculate temporaries and grid functions */ detgt = 1; - CCTK_REAL const T1000001 = INV(detgt); + gtu11 = INV(detgt)*(gt22L*gt33L - SQR(gt23L)); - gtu11 = T1000001*(gt22L*gt33L - SQR(gt23L)); + gtu21 = (gt13L*gt23L - gt12L*gt33L)*INV(detgt); - gtu21 = (gt13L*gt23L - gt12L*gt33L)*T1000001; + gtu31 = (-(gt13L*gt22L) + gt12L*gt23L)*INV(detgt); - gtu31 = (-(gt13L*gt22L) + gt12L*gt23L)*T1000001; + gtu22 = INV(detgt)*(gt11L*gt33L - SQR(gt13L)); - gtu22 = T1000001*(gt11L*gt33L - SQR(gt13L)); + gtu32 = (gt12L*gt13L - gt11L*gt23L)*INV(detgt); - gtu32 = (gt12L*gt13L - gt11L*gt23L)*T1000001; + gtu33 = INV(detgt)*(gt11L*gt22L - SQR(gt12L)); - gtu33 = T1000001*(gt11L*gt22L - SQR(gt12L)); - - trAt = At11L*gtu11 + At22L*gtu22 + - 2*(At12L*gtu21 + At13L*gtu31 + At23L*gtu32) + At33L*gtu33; + trAt = At11L*gtu11 + At22L*gtu22 + 2*(At12L*gtu21 + At13L*gtu31 + At23L*gtu32) + At33L*gtu33; At11L = At11L - gt11L*kthird*trAt; |