diff options
author | Peter Diener <diener@diener-3.lsu.edu> | 2008-07-01 14:21:21 -0500 |
---|---|---|
committer | Peter Diener <diener@diener-3.lsu.edu> | 2008-07-01 14:21:21 -0500 |
commit | 650e98e017aeecbff39cf9bdb0691aeb5f3e13ab (patch) | |
tree | 3deb750f996876321818833a487c31d3de1b8c81 /ML_BSSN_MP | |
parent | e4cf3ce9b6b4bf007ca1134db90b1beb94bd672a (diff) | |
parent | 9c1c9c787b183dba5ef41fbfb9028dbbf8b986e5 (diff) |
Merge commit 'origin/master'
Diffstat (limited to 'ML_BSSN_MP')
-rw-r--r-- | ML_BSSN_MP/param.ccl | 8 | ||||
-rw-r--r-- | ML_BSSN_MP/schedule.ccl | 34 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_ADMBaseBoundary.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_RHS.c | 35 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_RHSBoundary.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_boundary.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_constraints.c | 23 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_constraints_boundary.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_convertToADMBase.c | 2 | ||||
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_enforce.c | 2 |
13 files changed, 32 insertions, 86 deletions
diff --git a/ML_BSSN_MP/param.ccl b/ML_BSSN_MP/param.ccl index 4e5eb94..d00cc0e 100644 --- a/ML_BSSN_MP/param.ccl +++ b/ML_BSSN_MP/param.ccl @@ -116,6 +116,14 @@ KEYWORD my_boundary_condition "my_boundary_condition" } "none" restricted: +KEYWORD calculate_ADMBase_variables_at "calculate_ADMBase_variables_at" +{ + "MoL_PostStep" :: "MoL_PostStep" + "CCTK_EVOL" :: "CCTK_EVOL" + "CCTK_ANALYSIS" :: "CCTK_ANALYSIS" +} "MoL_PostStep" + +restricted: CCTK_INT ML_BSSN_MP_MaxNumEvolvedVars "Number of evolved variables used by this thorn" ACCUMULATOR-BASE=MethodofLines::MoL_Num_Evolved_Vars { 25:25 :: "Number of evolved variables used by this thorn" diff --git a/ML_BSSN_MP/schedule.ccl b/ML_BSSN_MP/schedule.ccl index ed8a8e0..92e1efa 100644 --- a/ML_BSSN_MP/schedule.ccl +++ b/ML_BSSN_MP/schedule.ccl @@ -220,14 +220,10 @@ if (CCTK_EQUALS(my_rhs_boundary_condition, "radiative")) } "ML_BSSN_MP_RHSBoundary" } - -if (CCTK_EQUALS(evolution_method, "ML_BSSN_MP")) +schedule ML_BSSN_MP_enforce IN MoL_PostStep BEFORE ML_BSSN_MP_BoundConds { - schedule ML_BSSN_MP_enforce IN MoL_PostStep BEFORE ML_BSSN_MP_BoundConds - { - LANG: C - } "ML_BSSN_MP_enforce" -} + LANG: C +} "ML_BSSN_MP_enforce" if (CCTK_EQUALS(my_boundary_condition, "Minkowski")) @@ -238,25 +234,21 @@ if (CCTK_EQUALS(my_boundary_condition, "Minkowski")) } "ML_BSSN_MP_boundary" } - -if (CCTK_EQUALS(evolution_method, "ML_BSSN_MP")) +schedule ML_BSSN_MP_convertToADMBase IN ML_BSSN_convertToADMBaseGroup { - schedule ML_BSSN_MP_convertToADMBase IN MoL_PostStep AFTER (ML_BSSN_MP_ApplyBCs ML_BSSN_MP_enforce) - { - LANG: C - SYNC: ADMBase::curv - SYNC: ADMBase::dtlapse - SYNC: ADMBase::dtshift - SYNC: ADMBase::lapse - SYNC: ADMBase::metric - SYNC: ADMBase::shift - } "ML_BSSN_MP_convertToADMBase" -} + LANG: C + SYNC: ADMBase::curv + SYNC: ADMBase::dtlapse + SYNC: ADMBase::dtshift + SYNC: ADMBase::lapse + SYNC: ADMBase::metric + SYNC: ADMBase::shift +} "ML_BSSN_MP_convertToADMBase" if (CCTK_EQUALS(my_boundary_condition, "Minkowski")) { - schedule ML_BSSN_MP_ADMBaseBoundary IN MoL_PostStep AFTER ML_BSSN_MP_convertToADMBase + schedule ML_BSSN_MP_ADMBaseBoundary IN ML_BSSN_convertToADMBaseGroup AFTER ML_BSSN_MP_convertToADMBase { LANG: C } "ML_BSSN_MP_ADMBaseBoundary" diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_ADMBaseBoundary.c b/ML_BSSN_MP/src/ML_BSSN_MP_ADMBaseBoundary.c index 5df7be8..adfb523 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_ADMBaseBoundary.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_ADMBaseBoundary.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_ADMBaseBoundary_Body(cGH const * const cctkGH, CCTK_INT const di #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_ADMBaseBoundary, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c b/ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c index 89254fd..eca7993 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_Minkowski.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_Minkowski_Body(cGH const * const cctkGH, CCTK_INT const dir, CCT #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_Minkowski, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_RHS.c b/ML_BSSN_MP/src/ML_BSSN_MP_RHS.c index 0f31fcd..f1f449a 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_RHS.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_RHS.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_RHS, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; @@ -169,10 +169,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33alpha = INITVALUE; CCTK_REAL PDstandardNth12alpha = INITVALUE; CCTK_REAL PDstandardNth13alpha = INITVALUE; - CCTK_REAL PDstandardNth21alpha = INITVALUE; CCTK_REAL PDstandardNth23alpha = INITVALUE; - CCTK_REAL PDstandardNth31alpha = INITVALUE; - CCTK_REAL PDstandardNth32alpha = INITVALUE; CCTK_REAL PDstandardNth1At11 = INITVALUE; CCTK_REAL PDstandardNth2At11 = INITVALUE; CCTK_REAL PDstandardNth3At11 = INITVALUE; @@ -208,10 +205,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33beta1 = INITVALUE; CCTK_REAL PDstandardNth12beta1 = INITVALUE; CCTK_REAL PDstandardNth13beta1 = INITVALUE; - CCTK_REAL PDstandardNth21beta1 = INITVALUE; CCTK_REAL PDstandardNth23beta1 = INITVALUE; - CCTK_REAL PDstandardNth31beta1 = INITVALUE; - CCTK_REAL PDstandardNth32beta1 = INITVALUE; CCTK_REAL PDstandardNth1beta2 = INITVALUE; CCTK_REAL PDstandardNth2beta2 = INITVALUE; CCTK_REAL PDstandardNth3beta2 = INITVALUE; @@ -220,10 +214,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33beta2 = INITVALUE; CCTK_REAL PDstandardNth12beta2 = INITVALUE; CCTK_REAL PDstandardNth13beta2 = INITVALUE; - CCTK_REAL PDstandardNth21beta2 = INITVALUE; CCTK_REAL PDstandardNth23beta2 = INITVALUE; - CCTK_REAL PDstandardNth31beta2 = INITVALUE; - CCTK_REAL PDstandardNth32beta2 = INITVALUE; CCTK_REAL PDstandardNth1beta3 = INITVALUE; CCTK_REAL PDstandardNth2beta3 = INITVALUE; CCTK_REAL PDstandardNth3beta3 = INITVALUE; @@ -232,10 +223,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33beta3 = INITVALUE; CCTK_REAL PDstandardNth12beta3 = INITVALUE; CCTK_REAL PDstandardNth13beta3 = INITVALUE; - CCTK_REAL PDstandardNth21beta3 = INITVALUE; CCTK_REAL PDstandardNth23beta3 = INITVALUE; - CCTK_REAL PDstandardNth31beta3 = INITVALUE; - CCTK_REAL PDstandardNth32beta3 = INITVALUE; CCTK_REAL PDstandardNth1gt11 = INITVALUE; CCTK_REAL PDstandardNth2gt11 = INITVALUE; CCTK_REAL PDstandardNth3gt11 = INITVALUE; @@ -244,10 +232,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt11 = INITVALUE; CCTK_REAL PDstandardNth12gt11 = INITVALUE; CCTK_REAL PDstandardNth13gt11 = INITVALUE; - CCTK_REAL PDstandardNth21gt11 = INITVALUE; CCTK_REAL PDstandardNth23gt11 = INITVALUE; - CCTK_REAL PDstandardNth31gt11 = INITVALUE; - CCTK_REAL PDstandardNth32gt11 = INITVALUE; CCTK_REAL PDstandardNth1gt12 = INITVALUE; CCTK_REAL PDstandardNth2gt12 = INITVALUE; CCTK_REAL PDstandardNth3gt12 = INITVALUE; @@ -256,10 +241,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt12 = INITVALUE; CCTK_REAL PDstandardNth12gt12 = INITVALUE; CCTK_REAL PDstandardNth13gt12 = INITVALUE; - CCTK_REAL PDstandardNth21gt12 = INITVALUE; CCTK_REAL PDstandardNth23gt12 = INITVALUE; - CCTK_REAL PDstandardNth31gt12 = INITVALUE; - CCTK_REAL PDstandardNth32gt12 = INITVALUE; CCTK_REAL PDstandardNth1gt13 = INITVALUE; CCTK_REAL PDstandardNth2gt13 = INITVALUE; CCTK_REAL PDstandardNth3gt13 = INITVALUE; @@ -268,10 +250,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt13 = INITVALUE; CCTK_REAL PDstandardNth12gt13 = INITVALUE; CCTK_REAL PDstandardNth13gt13 = INITVALUE; - CCTK_REAL PDstandardNth21gt13 = INITVALUE; CCTK_REAL PDstandardNth23gt13 = INITVALUE; - CCTK_REAL PDstandardNth31gt13 = INITVALUE; - CCTK_REAL PDstandardNth32gt13 = INITVALUE; CCTK_REAL PDstandardNth1gt22 = INITVALUE; CCTK_REAL PDstandardNth2gt22 = INITVALUE; CCTK_REAL PDstandardNth3gt22 = INITVALUE; @@ -280,10 +259,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt22 = INITVALUE; CCTK_REAL PDstandardNth12gt22 = INITVALUE; CCTK_REAL PDstandardNth13gt22 = INITVALUE; - CCTK_REAL PDstandardNth21gt22 = INITVALUE; CCTK_REAL PDstandardNth23gt22 = INITVALUE; - CCTK_REAL PDstandardNth31gt22 = INITVALUE; - CCTK_REAL PDstandardNth32gt22 = INITVALUE; CCTK_REAL PDstandardNth1gt23 = INITVALUE; CCTK_REAL PDstandardNth2gt23 = INITVALUE; CCTK_REAL PDstandardNth3gt23 = INITVALUE; @@ -292,10 +268,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt23 = INITVALUE; CCTK_REAL PDstandardNth12gt23 = INITVALUE; CCTK_REAL PDstandardNth13gt23 = INITVALUE; - CCTK_REAL PDstandardNth21gt23 = INITVALUE; CCTK_REAL PDstandardNth23gt23 = INITVALUE; - CCTK_REAL PDstandardNth31gt23 = INITVALUE; - CCTK_REAL PDstandardNth32gt23 = INITVALUE; CCTK_REAL PDstandardNth1gt33 = INITVALUE; CCTK_REAL PDstandardNth2gt33 = INITVALUE; CCTK_REAL PDstandardNth3gt33 = INITVALUE; @@ -304,10 +277,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33gt33 = INITVALUE; CCTK_REAL PDstandardNth12gt33 = INITVALUE; CCTK_REAL PDstandardNth13gt33 = INITVALUE; - CCTK_REAL PDstandardNth21gt33 = INITVALUE; CCTK_REAL PDstandardNth23gt33 = INITVALUE; - CCTK_REAL PDstandardNth31gt33 = INITVALUE; - CCTK_REAL PDstandardNth32gt33 = INITVALUE; CCTK_REAL PDstandardNth1phi = INITVALUE; CCTK_REAL PDstandardNth2phi = INITVALUE; CCTK_REAL PDstandardNth3phi = INITVALUE; @@ -316,10 +286,7 @@ void ML_BSSN_MP_RHS_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_INT CCTK_REAL PDstandardNth33phi = INITVALUE; CCTK_REAL PDstandardNth12phi = INITVALUE; CCTK_REAL PDstandardNth13phi = INITVALUE; - CCTK_REAL PDstandardNth21phi = INITVALUE; CCTK_REAL PDstandardNth23phi = INITVALUE; - CCTK_REAL PDstandardNth31phi = INITVALUE; - CCTK_REAL PDstandardNth32phi = INITVALUE; CCTK_REAL PDstandardNth1trK = INITVALUE; CCTK_REAL PDstandardNth2trK = INITVALUE; CCTK_REAL PDstandardNth3trK = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_RHSBoundary.c b/ML_BSSN_MP/src/ML_BSSN_MP_RHSBoundary.c index 0307fef..5a18cf0 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_RHSBoundary.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_RHSBoundary.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_RHSBoundary_Body(cGH const * const cctkGH, CCTK_INT const dir, C #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_RHSBoundary, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_boundary.c b/ML_BSSN_MP/src/ML_BSSN_MP_boundary.c index 3c591e6..13a0363 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_boundary.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_boundary.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_boundary_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_boundary, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_constraints.c b/ML_BSSN_MP/src/ML_BSSN_MP_constraints.c index 020e675..833d501 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_constraints.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_constraints.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_constraints, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; @@ -170,10 +170,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt11 = INITVALUE; CCTK_REAL PDstandardNth12gt11 = INITVALUE; CCTK_REAL PDstandardNth13gt11 = INITVALUE; - CCTK_REAL PDstandardNth21gt11 = INITVALUE; CCTK_REAL PDstandardNth23gt11 = INITVALUE; - CCTK_REAL PDstandardNth31gt11 = INITVALUE; - CCTK_REAL PDstandardNth32gt11 = INITVALUE; CCTK_REAL PDstandardNth1gt12 = INITVALUE; CCTK_REAL PDstandardNth2gt12 = INITVALUE; CCTK_REAL PDstandardNth3gt12 = INITVALUE; @@ -182,10 +179,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt12 = INITVALUE; CCTK_REAL PDstandardNth12gt12 = INITVALUE; CCTK_REAL PDstandardNth13gt12 = INITVALUE; - CCTK_REAL PDstandardNth21gt12 = INITVALUE; CCTK_REAL PDstandardNth23gt12 = INITVALUE; - CCTK_REAL PDstandardNth31gt12 = INITVALUE; - CCTK_REAL PDstandardNth32gt12 = INITVALUE; CCTK_REAL PDstandardNth1gt13 = INITVALUE; CCTK_REAL PDstandardNth2gt13 = INITVALUE; CCTK_REAL PDstandardNth3gt13 = INITVALUE; @@ -194,10 +188,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt13 = INITVALUE; CCTK_REAL PDstandardNth12gt13 = INITVALUE; CCTK_REAL PDstandardNth13gt13 = INITVALUE; - CCTK_REAL PDstandardNth21gt13 = INITVALUE; CCTK_REAL PDstandardNth23gt13 = INITVALUE; - CCTK_REAL PDstandardNth31gt13 = INITVALUE; - CCTK_REAL PDstandardNth32gt13 = INITVALUE; CCTK_REAL PDstandardNth1gt22 = INITVALUE; CCTK_REAL PDstandardNth2gt22 = INITVALUE; CCTK_REAL PDstandardNth3gt22 = INITVALUE; @@ -206,10 +197,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt22 = INITVALUE; CCTK_REAL PDstandardNth12gt22 = INITVALUE; CCTK_REAL PDstandardNth13gt22 = INITVALUE; - CCTK_REAL PDstandardNth21gt22 = INITVALUE; CCTK_REAL PDstandardNth23gt22 = INITVALUE; - CCTK_REAL PDstandardNth31gt22 = INITVALUE; - CCTK_REAL PDstandardNth32gt22 = INITVALUE; CCTK_REAL PDstandardNth1gt23 = INITVALUE; CCTK_REAL PDstandardNth2gt23 = INITVALUE; CCTK_REAL PDstandardNth3gt23 = INITVALUE; @@ -218,10 +206,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt23 = INITVALUE; CCTK_REAL PDstandardNth12gt23 = INITVALUE; CCTK_REAL PDstandardNth13gt23 = INITVALUE; - CCTK_REAL PDstandardNth21gt23 = INITVALUE; CCTK_REAL PDstandardNth23gt23 = INITVALUE; - CCTK_REAL PDstandardNth31gt23 = INITVALUE; - CCTK_REAL PDstandardNth32gt23 = INITVALUE; CCTK_REAL PDstandardNth1gt33 = INITVALUE; CCTK_REAL PDstandardNth2gt33 = INITVALUE; CCTK_REAL PDstandardNth3gt33 = INITVALUE; @@ -230,10 +215,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33gt33 = INITVALUE; CCTK_REAL PDstandardNth12gt33 = INITVALUE; CCTK_REAL PDstandardNth13gt33 = INITVALUE; - CCTK_REAL PDstandardNth21gt33 = INITVALUE; CCTK_REAL PDstandardNth23gt33 = INITVALUE; - CCTK_REAL PDstandardNth31gt33 = INITVALUE; - CCTK_REAL PDstandardNth32gt33 = INITVALUE; CCTK_REAL PDstandardNth1phi = INITVALUE; CCTK_REAL PDstandardNth2phi = INITVALUE; CCTK_REAL PDstandardNth3phi = INITVALUE; @@ -242,10 +224,7 @@ void ML_BSSN_MP_constraints_Body(cGH const * const cctkGH, CCTK_INT const dir, C CCTK_REAL PDstandardNth33phi = INITVALUE; CCTK_REAL PDstandardNth12phi = INITVALUE; CCTK_REAL PDstandardNth13phi = INITVALUE; - CCTK_REAL PDstandardNth21phi = INITVALUE; CCTK_REAL PDstandardNth23phi = INITVALUE; - CCTK_REAL PDstandardNth31phi = INITVALUE; - CCTK_REAL PDstandardNth32phi = INITVALUE; CCTK_REAL PDstandardNth1trK = INITVALUE; CCTK_REAL PDstandardNth2trK = INITVALUE; CCTK_REAL PDstandardNth3trK = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_constraints_boundary.c b/ML_BSSN_MP/src/ML_BSSN_MP_constraints_boundary.c index d074866..5bef2a7 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_constraints_boundary.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_constraints_boundary.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_constraints_boundary_Body(cGH const * const cctkGH, CCTK_INT con #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_constraints_boundary, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c index 51ceb9d..588f7f1 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBase.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_convertFromADMBase_Body(cGH const * const cctkGH, CCTK_INT const #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_convertFromADMBase, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c index d09e406..0ec22d6 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_convertFromADMBaseGamma.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_convertFromADMBaseGamma_Body(cGH const * const cctkGH, CCTK_INT #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_convertFromADMBaseGamma, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_convertToADMBase.c b/ML_BSSN_MP/src/ML_BSSN_MP_convertToADMBase.c index 447367e..7f4b0dc 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_convertToADMBase.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_convertToADMBase.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_convertToADMBase_Body(cGH const * const cctkGH, CCTK_INT const d #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_convertToADMBase, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c index d6628bb..49d98fd 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c +++ b/ML_BSSN_MP/src/ML_BSSN_MP_enforce.c @@ -102,7 +102,7 @@ void ML_BSSN_MP_enforce_Body(cGH const * const cctkGH, CCTK_INT const dir, CCTK_ #pragma omp parallel LC_LOOP3 (ML_BSSN_MP_enforce, i,j,k, min[0],min[1],min[2], max[0],max[1],max[2], - cctk_lsh[0],cctk_lsh[1],cctk_lsh[2]) + cctk_lssh[CCTK_LSSH_IDX(0,0)],cctk_lssh[CCTK_LSSH_IDX(0,1)],cctk_lssh[CCTK_LSSH_IDX(0,2)]) { int index = INITVALUE; int subblock_index = INITVALUE; |