test/sparse: add test for Accelerate

Author: jcarpent

unittest/CMakeLists.txt

@@ -172,4 +172,10 @@ if(BUILD_TESTING_SCIPY)
       "python")
 
   endif(BUILD_WITH_CHOLMOD_SUPPORT)
+
+  if(BUILD_WITH_ACCELERATE_SUPPORT)
+    add_python_unit_test(
+      "py-Accelerate"
+      "unittest/python/decompositions/sparse/test_Accelerate.py" "python")
+  endif(BUILD_WITH_ACCELERATE_SUPPORT)
 endif()

unittest/python/decompositions/sparse/test_Accelerate.py

@@ -0,0 +1,35 @@
+import eigenpy
+
+import numpy as np
+from scipy.sparse import csc_matrix
+
+
+def test(SolverType: type):
+    dim = 100
+    A = np.random.rand(dim, dim)
+    A = (A + A.T) * 0.5 + np.diag(10.0 + np.random.rand(dim))
+
+    A = csc_matrix(A)
+
+    llt = SolverType(A)
+
+    assert llt.info() == eigenpy.ComputationInfo.Success
+
+    X = np.random.rand(dim, 20)
+    B = A.dot(X)
+    X_est = llt.solve(B)
+    #    import pdb; pdb.set_trace()
+    assert eigenpy.is_approx(X, X_est)
+    assert eigenpy.is_approx(A.dot(X_est), B)
+
+    llt.analyzePattern(A)
+    llt.factorize(A)
+
+
+test(eigenpy.AccelerateLLT)
+test(eigenpy.AccelerateLDLT)
+test(eigenpy.AccelerateLDLTUnpivoted)
+test(eigenpy.AccelerateLDLTSBK)
+test(eigenpy.AccelerateLDLTTPP)
+test(eigenpy.AccelerateQR)
+# test(eigenpy.AccelerateCholeskyAtA)