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

Reply via email to