diff options
Diffstat (limited to 'ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc')
-rw-r--r-- | ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc | 46 |
1 files changed, 1 insertions, 45 deletions
diff --git a/ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc b/ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc index 17fd8c1..1fc6d96 100644 --- a/ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc +++ b/ML_BSSN_UPW/src/ML_BSSN_UPW_Advect.cc @@ -53,9 +53,6 @@ extern "C" void ML_BSSN_UPW_Advect_SelectBCs(CCTK_ARGUMENTS) ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_UPW::ML_shiftrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_UPW::ML_shiftrhs."); - ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_UPW::ML_Thetarhs","flat"); - if (ierr < 0) - CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_UPW::ML_Thetarhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_UPW::ML_trace_curvrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_UPW::ML_trace_curvrhs."); @@ -280,8 +277,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC gt33rhsL = vec_load(gt33rhs[index]); CCTK_REAL_VEC phiL = vec_load(phi[index]); CCTK_REAL_VEC phirhsL = vec_load(phirhs[index]); - CCTK_REAL_VEC ThetaL = vec_load(Theta[index]); - CCTK_REAL_VEC ThetarhsL = vec_load(Thetarhs[index]); CCTK_REAL_VEC trKL = vec_load(trK[index]); CCTK_REAL_VEC trKrhsL = vec_load(trKrhs[index]); CCTK_REAL_VEC Xt1L = vec_load(Xt1[index]); @@ -373,9 +368,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC PDupwindNth1phi; CCTK_REAL_VEC PDupwindNth2phi; CCTK_REAL_VEC PDupwindNth3phi; - CCTK_REAL_VEC PDupwindNth1Theta; - CCTK_REAL_VEC PDupwindNth2Theta; - CCTK_REAL_VEC PDupwindNth3Theta; CCTK_REAL_VEC PDupwindNth1trK; CCTK_REAL_VEC PDupwindNth2trK; CCTK_REAL_VEC PDupwindNth3trK; @@ -455,9 +447,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const PDupwindNth1phi = PDupwindNthfdOrder21(&phi[index]); PDupwindNth2phi = PDupwindNthfdOrder22(&phi[index]); PDupwindNth3phi = PDupwindNthfdOrder23(&phi[index]); - PDupwindNth1Theta = PDupwindNthfdOrder21(&Theta[index]); - PDupwindNth2Theta = PDupwindNthfdOrder22(&Theta[index]); - PDupwindNth3Theta = PDupwindNthfdOrder23(&Theta[index]); PDupwindNth1trK = PDupwindNthfdOrder21(&trK[index]); PDupwindNth2trK = PDupwindNthfdOrder22(&trK[index]); PDupwindNth3trK = PDupwindNthfdOrder23(&trK[index]); @@ -536,9 +525,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const PDupwindNth1phi = PDupwindNthfdOrder41(&phi[index]); PDupwindNth2phi = PDupwindNthfdOrder42(&phi[index]); PDupwindNth3phi = PDupwindNthfdOrder43(&phi[index]); - PDupwindNth1Theta = PDupwindNthfdOrder41(&Theta[index]); - PDupwindNth2Theta = PDupwindNthfdOrder42(&Theta[index]); - PDupwindNth3Theta = PDupwindNthfdOrder43(&Theta[index]); PDupwindNth1trK = PDupwindNthfdOrder41(&trK[index]); PDupwindNth2trK = PDupwindNthfdOrder42(&trK[index]); PDupwindNth3trK = PDupwindNthfdOrder43(&trK[index]); @@ -617,9 +603,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const PDupwindNth1phi = PDupwindNthfdOrder61(&phi[index]); PDupwindNth2phi = PDupwindNthfdOrder62(&phi[index]); PDupwindNth3phi = PDupwindNthfdOrder63(&phi[index]); - PDupwindNth1Theta = PDupwindNthfdOrder61(&Theta[index]); - PDupwindNth2Theta = PDupwindNthfdOrder62(&Theta[index]); - PDupwindNth3Theta = PDupwindNthfdOrder63(&Theta[index]); PDupwindNth1trK = PDupwindNthfdOrder61(&trK[index]); PDupwindNth2trK = PDupwindNthfdOrder62(&trK[index]); PDupwindNth3trK = PDupwindNthfdOrder63(&trK[index]); @@ -698,9 +681,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const PDupwindNth1phi = PDupwindNthfdOrder81(&phi[index]); PDupwindNth2phi = PDupwindNthfdOrder82(&phi[index]); PDupwindNth3phi = PDupwindNthfdOrder83(&phi[index]); - PDupwindNth1Theta = PDupwindNthfdOrder81(&Theta[index]); - PDupwindNth2Theta = PDupwindNthfdOrder82(&Theta[index]); - PDupwindNth3Theta = PDupwindNthfdOrder83(&Theta[index]); PDupwindNth1trK = PDupwindNthfdOrder81(&trK[index]); PDupwindNth2trK = PDupwindNthfdOrder82(&trK[index]); PDupwindNth3trK = PDupwindNthfdOrder83(&trK[index]); @@ -744,7 +724,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC JacPDupwindNth1gt23; CCTK_REAL_VEC JacPDupwindNth1gt33; CCTK_REAL_VEC JacPDupwindNth1phi; - CCTK_REAL_VEC JacPDupwindNth1Theta; CCTK_REAL_VEC JacPDupwindNth1trK; CCTK_REAL_VEC JacPDupwindNth1Xt1; CCTK_REAL_VEC JacPDupwindNth1Xt2; @@ -770,7 +749,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC JacPDupwindNth2gt23; CCTK_REAL_VEC JacPDupwindNth2gt33; CCTK_REAL_VEC JacPDupwindNth2phi; - CCTK_REAL_VEC JacPDupwindNth2Theta; CCTK_REAL_VEC JacPDupwindNth2trK; CCTK_REAL_VEC JacPDupwindNth2Xt1; CCTK_REAL_VEC JacPDupwindNth2Xt2; @@ -796,7 +774,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC JacPDupwindNth3gt23; CCTK_REAL_VEC JacPDupwindNth3gt33; CCTK_REAL_VEC JacPDupwindNth3phi; - CCTK_REAL_VEC JacPDupwindNth3Theta; CCTK_REAL_VEC JacPDupwindNth3trK; CCTK_REAL_VEC JacPDupwindNth3Xt1; CCTK_REAL_VEC JacPDupwindNth3Xt2; @@ -867,9 +844,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth1phi = kmadd(J11L,PDupwindNth1phi,kmadd(J21L,PDupwindNth2phi,kmul(J31L,PDupwindNth3phi))); - JacPDupwindNth1Theta = - kmadd(J11L,PDupwindNth1Theta,kmadd(J21L,PDupwindNth2Theta,kmul(J31L,PDupwindNth3Theta))); - JacPDupwindNth1trK = kmadd(J11L,PDupwindNth1trK,kmadd(J21L,PDupwindNth2trK,kmul(J31L,PDupwindNth3trK))); @@ -945,9 +919,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth2phi = kmadd(J12L,PDupwindNth1phi,kmadd(J22L,PDupwindNth2phi,kmul(J32L,PDupwindNth3phi))); - JacPDupwindNth2Theta = - kmadd(J12L,PDupwindNth1Theta,kmadd(J22L,PDupwindNth2Theta,kmul(J32L,PDupwindNth3Theta))); - JacPDupwindNth2trK = kmadd(J12L,PDupwindNth1trK,kmadd(J22L,PDupwindNth2trK,kmul(J32L,PDupwindNth3trK))); @@ -1023,9 +994,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth3phi = kmadd(J13L,PDupwindNth1phi,kmadd(J23L,PDupwindNth2phi,kmul(J33L,PDupwindNth3phi))); - JacPDupwindNth3Theta = - kmadd(J13L,PDupwindNth1Theta,kmadd(J23L,PDupwindNth2Theta,kmul(J33L,PDupwindNth3Theta))); - JacPDupwindNth3trK = kmadd(J13L,PDupwindNth1trK,kmadd(J23L,PDupwindNth2trK,kmul(J33L,PDupwindNth3trK))); @@ -1082,8 +1050,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth1phi = PDupwindNth1phi; - JacPDupwindNth1Theta = PDupwindNth1Theta; - JacPDupwindNth1trK = PDupwindNth1trK; JacPDupwindNth1Xt1 = PDupwindNth1Xt1; @@ -1134,8 +1100,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth2phi = PDupwindNth2phi; - JacPDupwindNth2Theta = PDupwindNth2Theta; - JacPDupwindNth2trK = PDupwindNth2trK; JacPDupwindNth2Xt1 = PDupwindNth2Xt1; @@ -1186,8 +1150,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const JacPDupwindNth3phi = PDupwindNth3phi; - JacPDupwindNth3Theta = PDupwindNth3Theta; - JacPDupwindNth3trK = PDupwindNth3trK; JacPDupwindNth3Xt1 = PDupwindNth3Xt1; @@ -1227,9 +1189,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const Xt3rhsL = kadd(Xt3rhsL,kmadd(beta1L,JacPDupwindNth1Xt3,kmadd(beta2L,JacPDupwindNth2Xt3,kmul(beta3L,JacPDupwindNth3Xt3)))); - ThetarhsL = kadd(ThetarhsL,IfThen(conformalMethod == - 2,kmadd(beta1L,JacPDupwindNth1Theta,kmadd(beta2L,JacPDupwindNth2Theta,kmul(beta3L,JacPDupwindNth3Theta))),ToReal(0))); - trKrhsL = kadd(trKrhsL,kmadd(beta1L,JacPDupwindNth1trK,kmadd(beta2L,JacPDupwindNth2trK,kmul(beta3L,JacPDupwindNth3trK)))); @@ -1298,7 +1257,6 @@ static void ML_BSSN_UPW_Advect_Body(cGH const * restrict const cctkGH, int const vec_store_nta_partial(gt23rhs[index],gt23rhsL); vec_store_nta_partial(gt33rhs[index],gt33rhsL); vec_store_nta_partial(phirhs[index],phirhsL); - vec_store_nta_partial(Thetarhs[index],ThetarhsL); vec_store_nta_partial(trKrhs[index],trKrhsL); vec_store_nta_partial(Xt1rhs[index],Xt1rhsL); vec_store_nta_partial(Xt2rhs[index],Xt2rhsL); @@ -1340,11 +1298,9 @@ extern "C" void ML_BSSN_UPW_Advect(CCTK_ARGUMENTS) "ML_BSSN_UPW::ML_metricrhs", "ML_BSSN_UPW::ML_shift", "ML_BSSN_UPW::ML_shiftrhs", - "ML_BSSN_UPW::ML_Theta", - "ML_BSSN_UPW::ML_Thetarhs", "ML_BSSN_UPW::ML_trace_curv", "ML_BSSN_UPW::ML_trace_curvrhs"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_UPW_Advect", 20, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_UPW_Advect", 18, groups); switch(fdOrder) { |