diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_Dissipation.cc')
-rw-r--r-- | ML_CCZ4/src/ML_CCZ4_Dissipation.cc | 18 |
1 files changed, 0 insertions, 18 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_Dissipation.cc b/ML_CCZ4/src/ML_CCZ4_Dissipation.cc index 673b790..da28695 100644 --- a/ML_CCZ4/src/ML_CCZ4_Dissipation.cc +++ b/ML_CCZ4/src/ML_CCZ4_Dissipation.cc @@ -258,10 +258,8 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int /* Assign local copies of grid functions */ - CCTK_REAL_VEC AL CCTK_ATTRIBUTE_UNUSED = vec_load(A[index]); CCTK_REAL_VEC alphaL CCTK_ATTRIBUTE_UNUSED = vec_load(alpha[index]); CCTK_REAL_VEC alpharhsL CCTK_ATTRIBUTE_UNUSED = vec_load(alpharhs[index]); - CCTK_REAL_VEC ArhsL CCTK_ATTRIBUTE_UNUSED = vec_load(Arhs[index]); CCTK_REAL_VEC At11L CCTK_ATTRIBUTE_UNUSED = vec_load(At11[index]); CCTK_REAL_VEC At11rhsL CCTK_ATTRIBUTE_UNUSED = vec_load(At11rhs[index]); CCTK_REAL_VEC At12L CCTK_ATTRIBUTE_UNUSED = vec_load(At12[index]); @@ -413,9 +411,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int { case 2: { - PDdissipationNth1A = PDdissipationNthfdOrder21(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder22(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder23(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder21(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder22(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder23(&alpha[index]); @@ -496,9 +491,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int case 4: { - PDdissipationNth1A = PDdissipationNthfdOrder41(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder42(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder43(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder41(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder42(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder43(&alpha[index]); @@ -579,9 +571,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int case 6: { - PDdissipationNth1A = PDdissipationNthfdOrder61(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder62(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder63(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder61(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder62(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder63(&alpha[index]); @@ -662,9 +651,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int case 8: { - PDdissipationNth1A = PDdissipationNthfdOrder81(&A[index]); - PDdissipationNth2A = PDdissipationNthfdOrder82(&A[index]); - PDdissipationNth3A = PDdissipationNthfdOrder83(&A[index]); PDdissipationNth1alpha = PDdissipationNthfdOrder81(&alpha[index]); PDdissipationNth2alpha = PDdissipationNthfdOrder82(&alpha[index]); PDdissipationNth3alpha = PDdissipationNthfdOrder83(&alpha[index]); @@ -1284,9 +1270,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int alpharhsL = kadd(alpharhsL,kmadd(epsdiss1,JacPDdissipationNth1alpha,kmadd(epsdiss2,JacPDdissipationNth2alpha,kmul(epsdiss3,JacPDdissipationNth3alpha)))); - ArhsL = - kadd(ArhsL,kmadd(epsdiss1,JacPDdissipationNth1A,kmadd(epsdiss2,JacPDdissipationNth2A,kmul(epsdiss3,JacPDdissipationNth3A)))); - beta1rhsL = kadd(beta1rhsL,kmadd(epsdiss1,JacPDdissipationNth1beta1,kmadd(epsdiss2,JacPDdissipationNth2beta1,kmul(epsdiss3,JacPDdissipationNth3beta1)))); @@ -1308,7 +1291,6 @@ static void ML_CCZ4_Dissipation_Body(const cGH* restrict const cctkGH, const int /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,vecimin,vecimax); vec_store_nta_partial(alpharhs[index],alpharhsL); - vec_store_nta_partial(Arhs[index],ArhsL); vec_store_nta_partial(At11rhs[index],At11rhsL); vec_store_nta_partial(At12rhs[index],At12rhsL); vec_store_nta_partial(At13rhs[index],At13rhsL); |