[petsc-users] Output data using ExodusIIViewer

[David Andrs]

Hello!

I am trying to store data into an ExodusII file using the ExodusIIViewer,
but running into a segfault inside PETSc. Attached is a minimal example
showing the problem. It can very much be that I am missing something
obvious. However, if I change the code to VTKViewer I get the desired
output file.

Machine: MacBook Pro 2019
OS version/type: Darwin notak.local 21.1.0 Darwin Kernel Version 21.1.0:
Wed Oct 13 17:33:23 PDT 2021; root:xnu-8019.41.5~1/RELEASE_X86_64 x86_64
PETSc: Petsc Release Version 3.16.1, Nov 01, 2021
MPI: MPICH 3.4.2
Compiler: clang-12

Call stack (not sure how relevant that is since it is from opt version):

* thread #1, queue = 'com.apple.main-thread', stop reason = EXC_BAD_ACCESS
(code=1, address=0xc)
    frame #0: 0x0000000102303ba9
libpetsc.3.16.dylib`DMView_PlexExodusII(dm=<unavailable>,
viewer=<unavailable>) at plexexodusii.c:457:45 [opt]
   454        else if (degree == 2) nodes[cs] = nodesHexP2;
   455      }
   456      /* Compute the number of cells not in the connectivity table */
-> 457      cellsNotInConnectivity -= nodes[cs][3]*csSize;
   458
   459      ierr = ISRestoreIndices(stratumIS, &cells);CHKERRQ(ierr);
   460      ierr = ISDestroy(&stratumIS);CHKERRQ(ierr);

With regards,

David Andrs
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211130/c591fe99/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: exo2.c
Type: application/octet-stream
Size: 1219 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20211130/c591fe99/attachment.obj>