Fix Trilinos nullspace

polyfem · Aug 12, 2024 · f920ea4 · f920ea4
1 parent 8ab068d
commit f920ea4
Show file tree

Hide file tree

Showing 3 changed files with 0 additions and 15 deletions.
diff --git a/src/polysolve/linear/Solver.hpp b/src/polysolve/linear/Solver.hpp
@@ -4,9 +4,6 @@
 
 #include <memory>
 
-#include <Eigen/Dense>
-#include <Eigen/Sparse>
-#include <vector>
 
 #define POLYSOLVE_DELETE_MOVE_COPY(Base) \
     Base(Base &&) = delete;              \

diff --git a/src/polysolve/linear/TrilinosSolver.cpp b/src/polysolve/linear/TrilinosSolver.cpp
@@ -12,7 +12,6 @@ namespace polysolve::linear
     TrilinosSolver::TrilinosSolver()
     {
         precond_num_ = 0;
-// #ifdef HAVE_MPI
         int done_already;
         MPI_Initialized(&done_already);
         if (!done_already)
@@ -25,9 +24,6 @@ namespace polysolve::linear
             MPI_Init(&argc, &argvv);
             CommPtr = new Epetra_MpiComm(MPI_COMM_WORLD);
         }
-// #else
-    //  CommPtr=new Epetra_SerialComm;
-// #endif
     }
 
     ////////////////////////////////////////////////////////////////
@@ -50,10 +46,6 @@ namespace polysolve::linear
             {
                 conv_tol_ = params["Trilinos"]["tolerance"];
             }
-            if (params["Trilinos"].contains("is_nullspace"))
-            {
-                is_nullspace_ = params["Trilinos"]["is_nullspace"];
-            }
         }
     }
 
@@ -289,9 +281,7 @@ namespace polysolve::linear
         delete A;
         delete rowMap;
         delete MLPrec;   
-// #ifdef HAVE_MPI
         MPI_Finalize() ;
-// #endif
     }
 }
 

diff --git a/src/polysolve/linear/TrilinosSolver.hpp b/src/polysolve/linear/TrilinosSolver.hpp
@@ -84,8 +84,6 @@ namespace polysolve::linear
         double conv_tol_ = 1e-8;
         size_t iterations_;
         double residual_error_;
-        bool is_nullspace_ = true;
-        Eigen::MatrixXd reduced_vertices;
         ML_Epetra::MultiLevelPreconditioner* MLPrec=NULL;
         Epetra_Map *rowMap=NULL;