Function PetscMonitor

src/solver/impls/petsc/petsc.cxx:931–980

}

PetscErrorCode PetscMonitor(TS ts, PetscInt UNUSED(step), PetscReal t, Vec X, void* ctx) {
  PetscErrorCode ierr;
  auto* s = static_cast<PetscSolver*>(ctx);
  PetscReal tfinal, dt;
  Vec interpolatedX;
  const PetscScalar* x;
  static int i = 0;

  PetscFunctionBegin;
  ierr = TSGetTimeStep(ts, &dt);
  CHKERRQ(ierr);

#if PETSC_VERSION_GE(3, 8, 0)
  ierr = TSGetMaxTime(ts, &tfinal);
  CHKERRQ(ierr);
#else
  ierr = TSGetDuration(ts, nullptr, &tfinal);
  CHKERRQ(ierr);
#endif

  /* Duplicate the solution vector X into a work vector */
  ierr = VecDuplicate(X, &interpolatedX);
  CHKERRQ(ierr);
  while (s->next_output <= t && s->next_output <= tfinal) {
    if (s->interpolate) {
      ierr = TSInterpolate(ts, s->next_output, interpolatedX);
      CHKERRQ(ierr);
    }

    /* Place the interpolated values into the global variables */
    ierr = VecGetArrayRead(interpolatedX, &x);
    CHKERRQ(ierr);
    s->load_vars(const_cast<BoutReal*>(x));
    ierr = VecRestoreArrayRead(interpolatedX, &x);
    CHKERRQ(ierr);

    if (s->call_monitors(simtime, i++, s->getNumberOutputSteps())) {
      PetscFunctionReturn(1);
    }

    s->next_output += s->getOutputTimestep();
    simtime = s->next_output;
  }

  /* Done with vector, so destroy it */
  ierr = VecDestroy(&interpolatedX);
  CHKERRQ(ierr);

  PetscFunctionReturn(0);
}

PetscErrorCode PetscSNESMonitor(SNES snes, PetscInt its, PetscReal norm, void* ctx) {
  PetscErrorCode ierr;