Made sure that preferred matrix type for biphasic problems in unsymme…

…tric.
febiosoftware · Oct 29, 2024 · 602986e · 602986e
1 parent e6a1885
commit 602986e
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/FEBioMix/FEBiphasicSolver.h b/FEBioMix/FEBiphasicSolver.h
@@ -103,6 +103,8 @@ class FEBIOMIX_API FEBiphasicSolver : public FENewtonSolver
 	void GetDisplacementData(vector<double>& di, vector<double>& ui);
 	void GetPressureData(vector<double>& pi, vector<double>& ui);
 
+	Matrix_Type PreferredMatrixType() const override { return REAL_UNSYMMETRIC; };
+
 public:
 	// additional convergence norms
 	double	m_Dtol;			//!< displacement tolerance

diff --git a/FECore/FESolver.h b/FECore/FESolver.h
@@ -183,8 +183,8 @@ class FECORE_API FESolver : public FECoreBase
 	// return the solution vector
 	virtual std::vector<double> GetSolutionVector() const;
 
-private:
-	Matrix_Type PreferredMatrixType() const;
+protected:
+	virtual Matrix_Type PreferredMatrixType() const;
 
 public: //TODO Move these parameters elsewhere
 	bool				m_bwopt;	    //!< bandwidth optimization flag