25   const PetscScalar *u;  in RHSFunctionslow()  local
29   PetscCall(VecGetArrayRead(U, &u));  in RHSFunctionslow()
31 …+ u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) + 0.05 * (-2.0 + u[1] * u[1] - PetscCosScalar(5.…  in RHSFunctionslow()
32   PetscCall(VecRestoreArrayRead(U, &u));  in RHSFunctionslow()
39   const PetscScalar *u;  in RHSFunctionfast()  local
43   PetscCall(VecGetArrayRead(U, &u));  in RHSFunctionfast()
45 …f[0] = 0.05 * (-1.0 + u[0] * u[0] - PetscCosScalar(t)) / (2.0 * u[0]) - 5.0 * PetscSinScalar(5.0 *…  in RHSFunctionfast()
46   PetscCall(VecRestoreArrayRead(U, &u));  in RHSFunctionfast()
53   const PetscScalar *u, *udot;  in IFunctionfast()  local
57   PetscCall(VecGetArrayRead(U, &u));  in IFunctionfast()
60   f[0] = udot[0] + (-2.0 + u[1] * u[1] - PetscCosScalar(5.0 * t)) / (2.0 * u[1]);  in IFunctionfast()
62   PetscCall(VecRestoreArrayRead(U, &u));  in IFunctionfast()
70   const PetscScalar *u;  in IJacobianfast()  local
74   PetscCall(VecGetArrayRead(U, &u));  in IJacobianfast()
75   J[0][0] = a + (2.0 + PetscCosScalar(5.0 * t)) / (2.0 * u[1] * u[1]) + 0.5;  in IJacobianfast()
77   PetscCall(VecRestoreArrayRead(U, &u));  in IJacobianfast()
93   PetscScalar *u;  in sol_true()  local
96   PetscCall(VecGetArray(U, &u));  in sol_true()
97   u[0] = PetscSqrtScalar(1.0 + PetscCosScalar(t));  in sol_true()
98   u[1] = PetscSqrtScalar(2.0 + PetscCosScalar(5.0 * t));  in sol_true()
99   PetscCall(VecRestoreArray(U, &u));  in sol_true()
111   PetscScalar *u;  in main()  local
159   PetscCall(VecGetArray(U, &u));  in main()
160   u[0] = PetscSqrtScalar(2.0);  in main()
161   u[1] = PetscSqrtScalar(3.0);  in main()
162   PetscCall(VecRestoreArray(U, &u));  in main()