diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_Advect.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_Advect.cc | 68 |
1 files changed, 34 insertions, 34 deletions
diff --git a/ML_BSSN/src/ML_BSSN_Advect.cc b/ML_BSSN/src/ML_BSSN_Advect.cc index a519fdf..e03d073 100644 --- a/ML_BSSN/src/ML_BSSN_Advect.cc +++ b/ML_BSSN/src/ML_BSSN_Advect.cc @@ -1036,14 +1036,14 @@ static void ML_BSSN_Advect_Body(const cGH* restrict const cctkGH, const int dir, double xx = x[index]; int origin = fabs(xx) < 1e-8; - PDupwindNthAnti2At11 = origin ? -2 * PDupwindNthAnti1At12 : -2 * At12[index] / xx; - PDupwindNthSymm2At11 = origin ? -2 * PDupwindNthSymm1At12 : -2 * At12[index] / xx; + PDupwindNthAnti2At11 = 0; + PDupwindNthSymm2At11 = 0; PDupwindNthAnti2At12 = origin ? PDupwindNthAnti1At11 - PDupwindNthAnti1At22 : (At11[index] - At22[index]) / xx; PDupwindNthSymm2At12 = origin ? PDupwindNthSymm1At11 - PDupwindNthSymm1At22 : (At11[index] - At22[index]) / xx; - PDupwindNthAnti2At13 = origin ? -PDupwindNthAnti1At23 : - At23[index] / xx; - PDupwindNthSymm2At13 = origin ? -PDupwindNthSymm1At23 : - At23[index] / xx; - PDupwindNthAnti2At22 = origin ? 2 * PDupwindNthAnti1At12 : 2 * At12[index] / xx; - PDupwindNthSymm2At22 = origin ? 2 * PDupwindNthSymm1At12 : 2 * At12[index] / xx; + PDupwindNthAnti2At13 = 0; + PDupwindNthSymm2At13 = 0; + PDupwindNthAnti2At22 = 0; + PDupwindNthSymm2At22 = 0; PDupwindNthAnti2At23 = origin ? PDupwindNthAnti1At13 : At13[index] / xx; PDupwindNthSymm2At23 = origin ? PDupwindNthSymm1At13 : At13[index] / xx; PDupwindNthAnti2At33 = 0; @@ -2059,34 +2059,34 @@ extern "C" void ML_BSSN_Advect(CCTK_ARGUMENTS) "ML_BSSN::ML_trace_curvrhs"}; GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Advect", 14, groups); - switch (fdOrder) - { - case 2: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 2, 2, 2); - break; - } - - case 4: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 3, 3, 3); - break; - } - - case 6: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 4, 4, 4); - break; - } - - case 8: - { - GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 5, 5, 5); - break; - } - default: - CCTK_BUILTIN_UNREACHABLE(); - } + //switch (fdOrder) + //{ + // case 2: + // { + // GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 2, 2, 2); + // break; + // } + // + // case 4: + // { + // GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 3, 3, 3); + // break; + // } + // + // case 6: + // { + // GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 4, 4, 4); + // break; + // } + // + // case 8: + // { + // GenericFD_EnsureStencilFits(cctkGH, "ML_BSSN_Advect", 5, 5, 5); + // break; + // } + // default: + // CCTK_BUILTIN_UNREACHABLE(); + //} GenericFD_LoopOverInterior(cctkGH, ML_BSSN_Advect_Body); |