1: #define PETSCDM_DLL
  2: #include <petsc/private/dmpleximpl.h>

  4: /*@
  5:   DMPlexCreateExodusFromFile - Create a `DMPLEX` mesh from an ExodusII file.

  7:   Collective

  9:   Input Parameters:
 10: + comm        - The MPI communicator
 11: . filename    - The name of the ExodusII file
 12: - interpolate - Create faces and edges in the mesh

 14:   Output Parameter:
 15: . dm - The `DM` object representing the mesh

 17:   Level: beginner

 19: .seealso: [](ch_unstructured), `DM`, `PETSCVIEWEREXODUSII`, `DMPLEX`, `DMCreate()`, `DMPlexCreateExodus()`
 20: @*/
 21: PetscErrorCode DMPlexCreateExodusFromFile(MPI_Comm comm, const char filename[], PetscBool interpolate, DM *dm)
 22: {
 23:   PetscFunctionBegin;
 24: #if defined(PETSC_HAVE_EXODUSII)
 25:   PetscMPIInt      rank;
 26:   PetscExodusIIInt CPU_word_size = sizeof(PetscReal), IO_word_size = 0, exoid = -1;
 27:   float            version;

 29:   PetscAssertPointer(filename, 2);
 30:   PetscCallMPI(MPI_Comm_rank(comm, &rank));
 31:   if (rank == 0) {
 32:     exoid = ex_open(filename, EX_READ, &CPU_word_size, &IO_word_size, &version);
 33:     PetscCheck(exoid > 0, PETSC_COMM_SELF, PETSC_ERR_LIB, "ex_open(\"%s\",...) did not return a valid file ID", filename);
 34:   }
 35:   PetscCall(DMPlexCreateExodus(comm, exoid, interpolate, dm));
 36:   if (rank == 0) PetscCallExternal(ex_close, exoid);
 37:   PetscFunctionReturn(PETSC_SUCCESS);
 38: #else
 39:   SETERRQ(comm, PETSC_ERR_SUP, "Loading meshes requires EXODUSII support. Reconfigure using --with-exodusii-dir or -download-exodusii");
 40: #endif
 41: }