Dear developers, I'm writing a sparse matrix into a file by doing
if (dump_mat) {
PetscViewer viewer;
PetscViewerASCIIOpen(PETSC_COMM_WORLD,"mat-par-aux.m",&viewer);
PetscViewerPushFormat(viewer, PETSC_VIEWER_ASCII_MATLAB);
MatView(A,viewer);
}
This works perfectly for small cases.
The program crashes for a case where the matrix A is of order 1 million but
with only 4 million non-zero elements.
Maybe at some point petsc is full-sizing A?
Thank you,
Alfredo
