Actual source code: dmmbio.cxx

  1: #include <petsc/private/dmmbimpl.h>
  2: #include <petscdmmoab.h>

  4: static PetscErrorCode DMMoab_GetWriteOptions_Private(PetscInt fsetid, PetscInt numproc, PetscInt dim, MoabWriteMode mode, PetscInt dbglevel, const char *dm_opts, const char *extra_opts, const char **write_opts)
  5: {
  6:   char *wopts;
  7:   char  wopts_par[PETSC_MAX_PATH_LEN];
  8:   char  wopts_parid[PETSC_MAX_PATH_LEN];
  9:   char  wopts_dbg[PETSC_MAX_PATH_LEN];

 11:   PetscMalloc1(PETSC_MAX_PATH_LEN, &wopts);
 12:   PetscMemzero(&wopts_par, PETSC_MAX_PATH_LEN);
 13:   PetscMemzero(&wopts_parid, PETSC_MAX_PATH_LEN);
 14:   PetscMemzero(&wopts_dbg, PETSC_MAX_PATH_LEN);

 16:   // do parallel read unless only one processor
 17:   if (numproc > 1) {
 18:     PetscSNPrintf(wopts_par, PETSC_MAX_PATH_LEN, "PARALLEL=%s;", MoabWriteModes[mode]);
 19:     if (fsetid >= 0) PetscSNPrintf(wopts_parid, PETSC_MAX_PATH_LEN, "PARALLEL_COMM=%d;", fsetid);
 20:   }

 22:   if (dbglevel) PetscSNPrintf(wopts_dbg, PETSC_MAX_PATH_LEN, "CPUTIME;DEBUG_IO=%d;", dbglevel);

 24:   PetscSNPrintf(wopts, PETSC_MAX_PATH_LEN, "%s%s%s%s%s", wopts_par, wopts_parid, wopts_dbg, (extra_opts ? extra_opts : ""), (dm_opts ? dm_opts : ""));
 25:   *write_opts = wopts;
 26:   return 0;
 27: }

 29: /*@C
 30:   DMMoabOutput - Output the solution vectors that are stored in the DMMoab object as tags
 31:   along with the complete mesh data structure in the native H5M or VTK format. The H5M output file
 32:   can be visualized directly with Paraview (if compiled with appropriate plugin) or converted
 33:   with MOAB/tools/mbconvert to a VTK or Exodus file.

 35:   This routine can also be used for check-pointing purposes to store a complete history of
 36:   the solution along with any other necessary data to restart computations.

 38:   Collective

 40:   Input Parameters:
 41: + dm     - the discretization manager object containing solution in MOAB tags.
 42: .  filename - the name of the output file: e.g., poisson.h5m
 43: -  usrwriteopts - the parallel write options needed for serializing a MOAB mesh database. Can be NULL.
 44:    Reference (Parallel Mesh Initialization: http://ftp.mcs.anl.gov/pub/fathom/moab-docs/contents.html#fivetwo)

 46:   Level: intermediate

 48: .seealso: `DMMoabLoadFromFile()`, `DMMoabSetGlobalFieldVector()`
 49: @*/
 50: PetscErrorCode DMMoabOutput(DM dm, const char *filename, const char *usrwriteopts)
 51: {
 52:   DM_Moab        *dmmoab;
 53:   const char     *writeopts;
 54:   PetscBool       isftype;
 55:   moab::ErrorCode merr;

 58:   dmmoab = (DM_Moab *)(dm)->data;

 60:   PetscStrendswith(filename, "h5m", &isftype);

 62:   /* add mesh loading options specific to the DM */
 63:   if (isftype) {
 64: #ifdef MOAB_HAVE_MPI
 65:     DMMoab_GetWriteOptions_Private(dmmoab->pcomm->get_id(), dmmoab->pcomm->size(), dmmoab->dim, dmmoab->write_mode, dmmoab->rw_dbglevel, dmmoab->extra_write_options, usrwriteopts, &writeopts);
 66: #else
 67:     DMMoab_GetWriteOptions_Private(0, 1, dmmoab->dim, dmmoab->write_mode, dmmoab->rw_dbglevel, dmmoab->extra_write_options, usrwriteopts, &writeopts);
 68: #endif
 69:     PetscInfo(dm, "Writing file %s with options: %s\n", filename, writeopts);
 70:   } else {
 71:     writeopts = NULL;
 72:   }

 74:   /* output file, using parallel write */
 75:   merr = dmmoab->mbiface->write_file(filename, NULL, writeopts, &dmmoab->fileset, 1);
 76:   MBERRVM(dmmoab->mbiface, "Writing output of DMMoab failed.", merr);
 77:   PetscFree(writeopts);
 78:   return 0;
 79: }