diff --git a/FEBioMix/FEBiphasicSolver.h b/FEBioMix/FEBiphasicSolver.h
index 0aa7c8c35..ea828f84f 100644
--- 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
index 5c0acf162..226776c6a 100644
--- 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