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

  4: /*@
  5:   PetscViewerExodusIIGetId - Get the file id of the `PETSCVIEWEREXODUSII` file

  7:   Logically Collective

  9:   Input Parameter:
 10: . viewer - the `PetscViewer`

 12:   Output Parameter:
 13: . exoid - The ExodusII file id

 15:   Level: intermediate

 17: .seealso: `PETSCVIEWEREXODUSII`, `PetscViewer`, `PetscViewerFileSetMode()`, `PetscViewerCreate()`, `PetscViewerSetType()`, `PetscViewerBinaryOpen()`
 18: @*/
 19: PetscErrorCode PetscViewerExodusIIGetId(PetscViewer viewer, int *exoid)
 20: {
 21:   PetscFunctionBegin;
 23:   PetscTryMethod(viewer, "PetscViewerGetId_C", (PetscViewer, int *), (viewer, exoid));
 24:   PetscFunctionReturn(PETSC_SUCCESS);
 25: }

 27: /*@
 28:   PetscViewerExodusIISetOrder - Set the elements order in the ExodusII file.

 30:   Collective

 32:   Input Parameters:
 33: + viewer - the `PETSCVIEWEREXODUSII` viewer
 34: - order  - elements order

 36:   Output Parameter:

 38:   Level: beginner

 40: .seealso: `PETSCVIEWEREXODUSII`, `PetscViewer`, `PetscViewerExodusIIGetId()`, `PetscViewerExodusIIGetOrder()`
 41: @*/
 42: PetscErrorCode PetscViewerExodusIISetOrder(PetscViewer viewer, PetscInt order)
 43: {
 44:   PetscFunctionBegin;
 46:   PetscTryMethod(viewer, "PetscViewerSetOrder_C", (PetscViewer, PetscInt), (viewer, order));
 47:   PetscFunctionReturn(PETSC_SUCCESS);
 48: }

 50: /*@
 51:   PetscViewerExodusIIGetOrder - Get the elements order in the ExodusII file.

 53:   Collective

 55:   Input Parameters:
 56: + viewer - the `PETSCVIEWEREXODUSII` viewer
 57: - order  - elements order

 59:   Output Parameter:

 61:   Level: beginner

 63: .seealso: `PETSCVIEWEREXODUSII`, `PetscViewer`, `PetscViewerExodusIIGetId()`, `PetscViewerExodusIISetOrder()`
 64: @*/
 65: PetscErrorCode PetscViewerExodusIIGetOrder(PetscViewer viewer, PetscInt *order)
 66: {
 67:   PetscFunctionBegin;
 69:   PetscTryMethod(viewer, "PetscViewerGetOrder_C", (PetscViewer, PetscInt *), (viewer, order));
 70:   PetscFunctionReturn(PETSC_SUCCESS);
 71: }

 73: /*@
 74:   PetscViewerExodusIIOpen - Opens a file for ExodusII input/output.

 76:   Collective

 78:   Input Parameters:
 79: + comm - MPI communicator
 80: . name - name of file
 81: - type - type of file
 82: .vb
 83:     FILE_MODE_WRITE - create new file for binary output
 84:     FILE_MODE_READ - open existing file for binary input
 85:     FILE_MODE_APPEND - open existing file for binary output
 86: .ve

 88:   Output Parameter:
 89: . exo - `PETSCVIEWEREXODUSII` `PetscViewer` for Exodus II input/output to use with the specified file

 91:   Level: beginner

 93: .seealso: `PETSCVIEWEREXODUSII`, `PetscViewer`, `PetscViewerPushFormat()`, `PetscViewerDestroy()`,
 94:           `DMLoad()`, `PetscFileMode`, `PetscViewerSetType()`, `PetscViewerFileSetMode()`, `PetscViewerFileSetName()`
 95: @*/
 96: PetscErrorCode PetscViewerExodusIIOpen(MPI_Comm comm, const char name[], PetscFileMode type, PetscViewer *exo)
 97: {
 98:   PetscFunctionBegin;
 99:   PetscCall(PetscViewerCreate(comm, exo));
100:   PetscCall(PetscViewerSetType(*exo, PETSCVIEWEREXODUSII));
101:   PetscCall(PetscViewerFileSetMode(*exo, type));
102:   PetscCall(PetscViewerFileSetName(*exo, name));
103:   PetscCall(PetscViewerSetFromOptions(*exo));
104:   PetscFunctionReturn(PETSC_SUCCESS);
105: }

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

110:   Collective

112:   Input Parameters:
113: + comm        - The MPI communicator
114: . filename    - The name of the ExodusII file
115: - interpolate - Create faces and edges in the mesh

117:   Output Parameter:
118: . dm - The `DM` object representing the mesh

120:   Level: beginner

122: .seealso: [](ch_unstructured), `DM`, `PETSCVIEWEREXODUSII`, `DMPLEX`, `DMCreate()`, `DMPlexCreateExodus()`
123: @*/
124: PetscErrorCode DMPlexCreateExodusFromFile(MPI_Comm comm, const char filename[], PetscBool interpolate, DM *dm)
125: {
126:   PetscFunctionBegin;
127: #if defined(PETSC_HAVE_EXODUSII)
128:   PetscMPIInt      rank;
129:   PetscExodusIIInt CPU_word_size = sizeof(PetscReal), IO_word_size = 0, exoid = -1;
130:   float            version;

132:   PetscAssertPointer(filename, 2);
133:   PetscCallMPI(MPI_Comm_rank(comm, &rank));
134:   if (rank == 0) {
135:     exoid = ex_open(filename, EX_READ, &CPU_word_size, &IO_word_size, &version);
136:     PetscCheck(exoid > 0, PETSC_COMM_SELF, PETSC_ERR_LIB, "ex_open(\"%s\",...) did not return a valid file ID", filename);
137:   }
138:   PetscCall(DMPlexCreateExodus(comm, exoid, interpolate, dm));
139:   if (rank == 0) PetscCallExternal(ex_close, exoid);
140:   PetscFunctionReturn(PETSC_SUCCESS);
141: #else
142:   SETERRQ(comm, PETSC_ERR_SUP, "Loading meshes requires EXODUSII support. Reconfigure using --with-exodusii-dir or -download-exodusii");
143: #endif
144: }