Switch MINRES in eigenpy.solvers scope

Author: Lucas-Haubert

CMakeLists.txt

@@ -191,6 +191,7 @@ set(${PROJECT_NAME}_SOLVERS_HEADERS
     include/eigenpy/solvers/IterativeSolverBase.hpp
     include/eigenpy/solvers/LeastSquaresConjugateGradient.hpp
     include/eigenpy/solvers/BiCGSTAB.hpp
+    include/eigenpy/solvers/MINRES.hpp
     include/eigenpy/solvers/ConjugateGradient.hpp
     include/eigenpy/solvers/SparseSolverBase.hpp
     include/eigenpy/solvers/BasicPreconditioners.hpp
@@ -248,7 +249,6 @@ set(${PROJECT_NAME}_DECOMPOSITIONS_HEADERS
     include/eigenpy/decompositions/CompleteOrthogonalDecomposition.hpp
     include/eigenpy/decompositions/FullPivHouseholderQR.hpp
     include/eigenpy/decompositions/SelfAdjointEigenSolver.hpp
-    include/eigenpy/decompositions/MINRES.hpp
     include/eigenpy/decompositions/SVDBase.hpp
     include/eigenpy/decompositions/BDCSVD.hpp
     include/eigenpy/decompositions/JacobiSVD.hpp)
@@ -335,7 +335,6 @@ set(${PROJECT_NAME}_DECOMPOSITIONS_SOURCES
     src/decompositions/tridiagonalization.cpp
     src/decompositions/real-schur.cpp
     src/decompositions/partialpivlu-solver.cpp
-    src/decompositions/minres-solver.cpp
     src/decompositions/sparse-lu-solver.cpp
     src/decompositions/sparse-qr-solver.cpp
     src/decompositions/qr-solvers.cpp

include/eigenpy/decompositions/MINRES.hpp renamed to include/eigenpy/solvers/MINRES.hpp

@@ -23,7 +23,6 @@ void exposeLDLTSolver();
 void exposeFullPivLUSolver();
 void exposePartialPivLUSolver();
 void exposeQRSolvers();
-void exposeMINRESSolver();
 void exposeSimplicialLLTSolver();
 void exposeSimplicialLDLTSolver();
 void exposeIncompleteCholesky();
@@ -52,7 +51,6 @@ void exposeDecompositions() {
   exposeFullPivLUSolver();
   exposePartialPivLUSolver();
   exposeQRSolvers();
-  exposeMINRESSolver();
   exposeBDCSVDSolver();
   exposeJacobiSVDSolver();
 

src/decompositions/decompositions.cpp

@@ -14,6 +14,7 @@
 #endif
 
 #include "eigenpy/solvers/BiCGSTAB.hpp"
+#include "eigenpy/solvers/MINRES.hpp"
 
 namespace eigenpy {
 void exposeSolvers() {
@@ -54,6 +55,8 @@ void exposeSolvers() {
 
   BiCGSTABVisitor<BiCGSTAB<MatrixXd>>::expose("BiCGSTAB");
   BiCGSTABVisitor<IdentityBiCGSTAB>::expose("IdentityBiCGSTAB");
+
+  MINRESSolverVisitor<MatrixXd>::expose("MINRES");
 }
 }  // namespace eigenpy
 

src/decompositions/minres-solver.cpp

@@ -6,7 +6,7 @@
 rng = np.random.default_rng()
 A = np.eye(dim)
 
-minres = eigenpy.MINRES(A)
+minres = eigenpy.solvers.MINRES(A)
 
 X = rng.random((dim, 20))
 B = A.dot(X)