It was wrongly assumed that the matrix is a MatDense but it can be a MatSparse in which case the casting to MatDense gives a nullptr and a segfault.