diff --git a/MathLib/LinAlg/Dense/GlobalDenseMatrix.h b/MathLib/LinAlg/Dense/GlobalDenseMatrix.h
index ddf1939c87caba0c72f0a30047672d0b34349cb4..0c9bc24f243d8b2b6a82e50e17e35b91549b3f3b 100644
--- a/MathLib/LinAlg/Dense/GlobalDenseMatrix.h
+++ b/MathLib/LinAlg/Dense/GlobalDenseMatrix.h
@@ -70,6 +70,16 @@ public:
 	virtual bool addValue(IDX_TYPE row, IDX_TYPE col, FP_TYPE val);
 
 
+	/// Add sub-matrix at positions \c row_pos and same column positions as the
+	/// given row positions.
+	template<class T_DENSE_MATRIX>
+	void addSubMatrix(std::vector<IDX_TYPE> const& row_pos,
+			const T_DENSE_MATRIX &sub_matrix,
+			FP_TYPE fkt = static_cast<FP_TYPE>(1.0))
+	{
+		this->addSubMatrix(row_pos, row_pos, sub_matrix, fkt);
+	}
+
 	template<class T_DENSE_MATRIX>
 	void addSubMatrix(std::vector<IDX_TYPE> const& row_pos,
 			std::vector<IDX_TYPE> const& col_pos, const T_DENSE_MATRIX &sub_matrix,