#include #include struct _n_PetscViewers { MPI_Comm comm; PetscViewer *viewer; int n; }; /*@C PetscViewersDestroy - Destroys a set of `PetscViewer`s created with `PetscViewersCreate()`. Collective on v Input Parameters: . v - the `PetscViewer`s to be destroyed. Level: intermediate .seealso: `PetscViewer`, `PetscViewerSocketOpen()`, `PetscViewerASCIIOpen()`, `PetscViewerCreate()`, `PetscViewerDrawOpen()`, `PetscViewersCreate()` @*/ PetscErrorCode PetscViewersDestroy(PetscViewers *v) { int i; PetscFunctionBegin; if (!*v) PetscFunctionReturn(0); for (i = 0; i < (*v)->n; i++) PetscCall(PetscViewerDestroy(&(*v)->viewer[i])); PetscCall(PetscFree((*v)->viewer)); PetscCall(PetscFree(*v)); PetscFunctionReturn(0); } /*@C PetscViewersCreate - Creates a container to hold a set of `PetscViewer`s. The container is essentially a sparse, growable in length array of `PetscViewer`s Collective Input Parameter: . comm - the MPI communicator Output Parameter: . v - the collection of `PetscViewer`s Level: intermediate .seealso: `PetscViewer`, `PetscViewerCreate()`, `PetscViewersDestroy()` @*/ PetscErrorCode PetscViewersCreate(MPI_Comm comm, PetscViewers *v) { PetscFunctionBegin; PetscValidPointer(v, 2); PetscCall(PetscNew(v)); (*v)->n = 64; (*v)->comm = comm; PetscCall(PetscCalloc1(64, &(*v)->viewer)); PetscFunctionReturn(0); } /*@C PetscViewersGetViewer - Gets a `PetscViewer` from a `PetscViewer` collection Not Collective, but the resulting `PetscViewer` will be collective object on viewers Input Parameters: + viewers - object created with `PetscViewersCreate()` - n - number of `PetscViewer `you want Output Parameter: . viewer - the `PetscViewer` Level: intermediate .seealso: `PetscViewer`, `PetscViewersCreate()`, `PetscViewersDestroy()` @*/ PetscErrorCode PetscViewersGetViewer(PetscViewers viewers, PetscInt n, PetscViewer *viewer) { PetscFunctionBegin; PetscValidPointer(viewers, 1); PetscValidPointer(viewer, 3); PetscCheck(n >= 0, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Cannot access using a negative index - %" PetscInt_FMT, n); if (n >= viewers->n) { PetscViewer *v; int newn = n + 64; /* add 64 new ones at a time */ PetscCall(PetscCalloc1(newn, &v)); PetscCall(PetscArraycpy(v, viewers->viewer, viewers->n)); PetscCall(PetscFree(viewers->viewer)); viewers->viewer = v; } if (!viewers->viewer[n]) PetscCall(PetscViewerCreate(viewers->comm, &viewers->viewer[n])); *viewer = viewers->viewer[n]; PetscFunctionReturn(0); } /* PetscMonitorCompare - Checks if two monitors are identical; if they are then it destroys the new one Not collective Input Parameters: + nmon - The new monitor . nmctx - The new monitor context, or NULL . nmdestroy - The new monitor destroy function, or NULL . mon - The old monitor . mctx - The old monitor context, or NULL - mdestroy - The old monitor destroy function, or NULL Output Parameter: . identical - `PETSC_TRUE` if the monitors are the same Level: developer .seealso: `DMMonitorSetFromOptions()`, `KSPMonitorSetFromOptions()`, `SNESMonitorSetFromOptions()` */ PetscErrorCode PetscMonitorCompare(PetscErrorCode (*nmon)(void), void *nmctx, PetscErrorCode (*nmdestroy)(void **), PetscErrorCode (*mon)(void), void *mctx, PetscErrorCode (*mdestroy)(void **), PetscBool *identical) { PetscFunctionBegin; PetscValidBoolPointer(identical, 7); *identical = PETSC_FALSE; if (nmon == mon && nmdestroy == mdestroy) { if (nmctx == mctx) *identical = PETSC_TRUE; else if (nmdestroy == (PetscErrorCode(*)(void **))PetscViewerAndFormatDestroy) { PetscViewerAndFormat *old = (PetscViewerAndFormat *)mctx, *newo = (PetscViewerAndFormat *)nmctx; if (old->viewer == newo->viewer && old->format == newo->format) *identical = PETSC_TRUE; } if (*identical) { if (mdestroy) PetscCall((*mdestroy)(&nmctx)); } } PetscFunctionReturn(0); }