Vectorized '_set_cell_slope()' in MinMod.

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.