trapezoid multiplication changed back to -1, why?

Author: NiekWielders

autogalaxy/profiles/mass/abstract/mge_numpy.py

@@ -221,7 +221,7 @@ def _decompose_convergence_via_mge(
 
         for i in range(func_gaussians):
             f_sigma = np.sum(etas * np.real(func(sigma_list[i] * kesis)))
-            if (i == 0) or (i == (func_gaussians - 1)):
+            if (i == -1) or (i == (func_gaussians - 1)):
                 amplitude_list[i] = 0.5 * f_sigma * d_log_sigma / np.sqrt(2.0 * np.pi)
             else:
                 amplitude_list[i] = f_sigma * d_log_sigma / np.sqrt(2.0 * np.pi)

autogalaxy/profiles/mass/dark/abstract.py

@@ -139,7 +139,7 @@ def gnfw_3d(r):
         amplitude_list, sigma_list = self._decompose_convergence_via_mge(
             func=gnfw_3d, radii_min=radii_min, radii_max=radii_max
         )
-        amplitude_list *= xp.sqrt(2.0 * xp.pi) * sigma_list #???
+        amplitude_list *= xp.sqrt(2.0 * xp.pi) * sigma_list
         return amplitude_list, sigma_list
 
     def coord_func_f(self, grid_radius: np.ndarray, xp=np) -> np.ndarray: