215   EmaxCtx      emaxCtx;      // Information for fit to decay profile  member
288   options->emaxCtx.s              = 50;  in ProcessOptions()
289   options->emaxCtx.per            = 100;  in ProcessOptions()
314 …", "First time step to use for Emax fits", __FILE__, options->emaxCtx.s, &options->emaxCtx.s, NULL…  in ProcessOptions()
315 …"Number of time steps between Emax fits", __FILE__, options->emaxCtx.per, &options->emaxCtx.per, N…  in ProcessOptions()
812 …PetscCall(PetscDrawLGGetData(ctx->drawlgE, NULL, &ctx->emaxCtx.e, &ctx->emaxCtx.t, &ctx->emaxCtx.E…  in MonitorEField()
813   if (!rank && !(ctx->emaxCtx.e % ctx->emaxCtx.per)) {  in MonitorEField()
824     PetscCall(VecCreateSeq(PETSC_COMM_SELF, ctx->emaxCtx.e - ctx->emaxCtx.s, &r));  in MonitorEField()
825     if (fitLog) PetscCall(TaoSetResidualRoutine(tao, r, ComputeLogEmaxResidual, &ctx->emaxCtx));  in MonitorEField()
826     else PetscCall(TaoSetResidualRoutine(tao, r, ComputeEmaxResidual, &ctx->emaxCtx));  in MonitorEField()
828     PetscCall(MatCreateSeqDense(PETSC_COMM_SELF, ctx->emaxCtx.e - ctx->emaxCtx.s, 4, NULL, &J));  in MonitorEField()
829 …fitLog) PetscCall(TaoSetJacobianResidualRoutine(tao, J, J, ComputeLogEmaxJacobian, &ctx->emaxCtx));  in MonitorEField()
830     else PetscCall(TaoSetJacobianResidualRoutine(tao, J, J, ComputeEmaxJacobian, &ctx->emaxCtx));  in MonitorEField()
842       for (PetscInt i = 0; i < ctx->emaxCtx.e; ++i) {  in MonitorEField()
844         PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", ctx->emaxCtx.t[i]));  in MonitorEField()
848       for (PetscInt i = 0; i < ctx->emaxCtx.e; ++i) {  in MonitorEField()
850         PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g", ctx->emaxCtx.Emax[i]));  in MonitorEField()