49   const PetscScalar *y, *ydot;  in IFunctionImplicit()  local
54   PetscCall(VecGetArrayRead(Y, &y));  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 …dot[1] / 1.e6 - 0.0006666766666666667 + PetscExpReal((500 * (y[1] - y[2])) / 13.) / 1.e8 + y[1] / …  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 …/ 1.e6 - 0.0006676566666666666 + (99 * PetscExpReal((500 * (y[1] - y[2])) / 13.)) / 1.e8 + y[3] / …  in IFunctionImplicit()
62   f[4] = (3 * ydot[4]) / 1.e6 - (3 * ydot[3]) / 1.e6 + y[4] / 9000.;  in IFunctionImplicit()
64   PetscCall(VecRestoreArrayRead(Y, &y));  in IFunctionImplicit()
76   const PetscScalar *y, *ydot;  in IJacobianImplicit()  local
80   PetscCall(VecGetArrayRead(Y, &y));  in IJacobianImplicit()
88   J[1][1] = a / 1.e6 + 0.00022222222222222223 + PetscExpReal((500 * (y[1] - y[2])) / 13.) / 2.6e6;  in IJacobianImplicit()
89   J[1][2] = -PetscExpReal((500 * (y[1] - y[2])) / 13.) / 2.6e6;  in IJacobianImplicit()
90   J[2][1] = -PetscExpReal((500 * (y[1] - y[2])) / 13.) / 26000.;  in IJacobianImplicit()
91 …J[2][2] = a / 500000 + 0.00011111111111111112 + PetscExpReal((500 * (y[1] - y[2])) / 13.) / 26000.;  in IJacobianImplicit()
92   J[3][1] = (99 * PetscExpReal((500 * (y[1] - y[2])) / 13.)) / 2.6e6;  in IJacobianImplicit()
93   J[3][2] = (-99 * PetscExpReal((500 * (y[1] - y[2])) / 13.)) / 2.6e6;  in IJacobianImplicit()
101   PetscCall(VecRestoreArrayRead(Y, &y));  in IJacobianImplicit()
120   PetscScalar *y;  in main()  local
140   PetscCall(VecGetArray(Y, &y));  in main()
141   y[0] = 0.0;  in main()
142   y[1] = 3.0;  in main()
143   y[2] = y[1];  in main()
144   y[3] = 6.0;  in main()
145   y[4] = 0.0;  in main()
146   PetscCall(VecRestoreArray(Y, &y));  in main()