You also want:
PetscCall(PetscViewerPopFormat(viewer));
PetscCall(PetscViewerDestroy(&viewer));
This should not be a problem.
If this is a segv and you configure it with '--with-debugging=1', you
should get a stack trace, which would help immensely.
Or run in a debugger to get a stack trace.
Thanks,
Mark
On Fri, Aug 12, 2022 at 11:26 AM Alfredo Jaramillo <
[email protected]> wrote:
> 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
>
