Make SparseMatrix::coeff() returns a const reference and add a non const version of SparseMatrix::diagonal()
Document non-const SparseMatrix::diagonal() method.