void InLDLTPreconditioner::resetPreconditionerSystem(const BlockedSymSpMatrix& mat){
		int columnCount = mat.getNumColumns();
		values.reserve(columnCount);
		rows.reserve(columnCount);
		pcol.reserve(columnCount);
		invDiagonal.resize(columnCount);
		matFullIndices = mat.getFullIndices();
	}
	InLDLTPreconditioner::InLDLTPreconditioner(const BlockedSymSpMatrix& mat, Scalar sainvEps, Scalar ldltEps)
	: sainvEpsilon(sainvEpsilon), ldltEpsilon(ldltEpsilon), matFullIndices(mat.getFullIndices()) {

		int columnCount = mat.getNumColumns();
		values.reserve(columnCount);
		rows.reserve(columnCount);
		pcol.reserve(columnCount);
		invDiagonal.resize(columnCount);
	}