diff --git a/scripts/tcd/Limiter.py b/scripts/tcd/Limiter.py index ab2be32c620c7a433c2ba55fe1309a6d9239b18b..099719850f45bcf8aa306e7f7d4f9d229c1bf0e6 100644 --- a/scripts/tcd/Limiter.py +++ b/scripts/tcd/Limiter.py @@ -191,24 +191,11 @@ class MinMod(Limiter): Slope of the given cell. """ - slope = [] - for current_cell in range(len(projection[0])): - new_entry = sum( - projection[degree][current_cell] * (degree+0.5)**0.5 - for degree in range(1, len(projection))) - slope.append(new_entry) + root_vector = np.array([np.sqrt(degree+0.5) + for degree in range(len(projection))]) + slope = root_vector[1:] @ projection[1:] return slope[cell] - # # print(np.array(slope).shape) - # # print(slope) - # root_vector = np.array([np.sqrt(degree+0.5) - # for degree in range(len(projection))]) - # test = root_vector[1:] @ projection[1:] - # # print(test.shape) - # # print(np.isclose(test, slope, rtol=1e-16)) - # - # return test[cell] - class ModifiedMinMod(MinMod): """Class for modified minmod limiting function.