MINRES: Backward compatibility to call from eigenpy and eigenpy.solvers

Author: Lucas-Haubert

src/decompositions/decompositions.cpp

@@ -3,6 +3,7 @@
  */
 
 #include "eigenpy/decompositions/decompositions.hpp"
+#include "eigenpy/solvers/MINRES.hpp"
 
 #include "eigenpy/fwd.hpp"
 
@@ -31,6 +32,11 @@ void exposePermutationMatrix();
 void exposeBDCSVDSolver();
 void exposeJacobiSVDSolver();
 
+void exposeBackwardCompatibilityAliases() {
+  typedef Eigen::MatrixXd MatrixXd;
+  MINRESSolverVisitor<MatrixXd>::expose("MINRES");
+}
+
 void exposeDecompositions() {
   using namespace Eigen;
 
@@ -80,5 +86,7 @@ void exposeDecompositions() {
 #ifdef EIGENPY_WITH_ACCELERATE_SUPPORT
   exposeAccelerate();
 #endif
+
+  exposeBackwardCompatibilityAliases();
 }
 }  // namespace eigenpy

unittest/python/test_MINRES.py

@@ -11,6 +11,11 @@
 X = rng.random((dim, 20))
 B = A.dot(X)
 X_est = minres.solve(B)
-print("A.dot(X_est):", A.dot(X_est))
-print("B:", B)
+assert eigenpy.is_approx(A.dot(X_est), B, 1e-6)
+
+minres_back = eigenpy.MINRES(A)
+
+X = rng.random((dim, 20))
+B = A.dot(X)
+X_est = minres_back.solve(B)
 assert eigenpy.is_approx(A.dot(X_est), B, 1e-6)