fix(pyabacus): call static base wrappers directly

Author: zxy.monado

python/pyabacus/CONTRIBUTING.md

@@ -375,7 +375,7 @@ To contribute to the `pyabacus` project, follow these steps:
 
        @staticmethod
        def sphbes_zeros(l: int, n: int, zeros: NDArray[np.float64]) -> None: 
-           super().sphbes_zeros(l, n, zeros)
+           _Sphbes.sphbes_zeros(l, n, zeros)
    ```
 
 ## Conclusion

python/pyabacus/src/pyabacus/ModuleBase/_module_base.py

@@ -48,7 +48,7 @@ def dsphbesj(self, *args, **kwargs):
 
     @staticmethod
     def sphbes_zeros(l: int, n: int, zeros: NDArray[np.float64]) -> None: 
-        super().sphbes_zeros(l, n, zeros)
+        _Sphbes.sphbes_zeros(l, n, zeros)
 
 class Integral(_Integral):
     def __init__(self) -> None: 
@@ -81,7 +81,7 @@ def Simpson_Integral_0toall(
         rab: NDArray[np.float64], 
         asum: NDArray[np.float64]
     ) -> None: 
-        super().Simpson_Integral_0toall(mesh, func, rab, asum)
+        _Integral.Simpson_Integral_0toall(mesh, func, rab, asum)
 
     @staticmethod
     def Simpson_Integral_alltoinf(
@@ -90,7 +90,7 @@ def Simpson_Integral_alltoinf(
         rab: NDArray[np.float64], 
         asum: NDArray[np.float64]
     ) -> None: 
-        super().Simpson_Integral_alltoinf(mesh, func, rab, asum)
+        _Integral.Simpson_Integral_alltoinf(mesh, func, rab, asum)
 
     @overload
     @staticmethod

python/pyabacus/tests/test_base_math.py

@@ -13,6 +13,23 @@ def test_sphbes():
     assert s.sphbesj(0, 0.0) == 1.0
     assert s.sphbesj(1, np.array([0.0]), 1, 1, np.zeros(1)) == None
 
+def test_static_array_wrappers():
+    zeros = np.zeros(2)
+    base.Sphbes.sphbes_zeros(0, 2, zeros)
+    np.testing.assert_allclose(zeros, [np.pi, 2 * np.pi])
+
+    mesh = 3
+    func = np.ones(mesh)
+    rab = np.ones(mesh)
+    integral_from_zero = np.zeros(mesh)
+    integral_to_infinity = np.zeros(mesh)
+
+    base.Integral.Simpson_Integral_0toall(mesh, func, rab, integral_from_zero)
+    base.Integral.Simpson_Integral_alltoinf(mesh, func, rab, integral_to_infinity)
+
+    np.testing.assert_allclose(integral_from_zero, [0.0, 1.0, 2.0])
+    np.testing.assert_allclose(integral_to_infinity, [2.0, 1.0, 0.0])
+
 def test_sbt():
     sbt = base.SphericalBesselTransformer()