49   const PetscScalar *y, *ydot;  in IFunctionImplicit()  local
55   PetscCall(VecGetArrayRead(Ydot, &ydot));  in IFunctionImplicit()
58   f[0] = ydot[0] / 1.e6 - ydot[1] / 1.e6 - PetscSinReal(200 * PETSC_PI * t) / 2500. + y[0] / 1000.;  in IFunctionImplicit()
59 …f[1] = -ydot[0] / 1.e6 + ydot[1] / 1.e6 - 0.0006666766666666667 + PetscExpReal((500 * (y[1] - y[2]…  in IFunctionImplicit()
60 …f[2] = ydot[2] / 500000. + 1.e-6 - PetscExpReal((500 * (y[1] - y[2])) / 13.) / 1.e6 + y[2] / 9000.;  in IFunctionImplicit()
61 …f[3] = (3 * ydot[3]) / 1.e6 - (3 * ydot[4]) / 1.e6 - 0.0006676566666666666 + (99 * PetscExpReal((5…  in IFunctionImplicit()
62   f[4] = (3 * ydot[4]) / 1.e6 - (3 * ydot[3]) / 1.e6 + y[4] / 9000.;  in IFunctionImplicit()
65   PetscCall(VecRestoreArrayRead(Ydot, &ydot));  in IFunctionImplicit()
76   const PetscScalar *y, *ydot;  in IJacobianImplicit()  local
81   PetscCall(VecGetArrayRead(Ydot, &ydot));  in IJacobianImplicit()
102   PetscCall(VecRestoreArrayRead(Ydot, &ydot));  in IJacobianImplicit()