34   const PetscScalar *x, *xdot;  in RoberFunction()  local
38   PetscCall(VecGetArrayRead(Xdot, &xdot));  in RoberFunction()
40   f[0] = xdot[0] + 0.04 * x[0] - 1e4 * x[1] * x[2];  in RoberFunction()
41   f[1] = xdot[1] - 0.04 * x[0] + 1e4 * x[1] * x[2] + 3e7 * PetscSqr(x[1]);  in RoberFunction()
42   f[2] = xdot[2] - 3e7 * PetscSqr(x[1]);  in RoberFunction()
44   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in RoberFunction()
53   const PetscScalar *x, *xdot;  in RoberJacobian()  local
57   PetscCall(VecGetArrayRead(Xdot, &xdot));  in RoberJacobian()
69   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in RoberJacobian()
125   const PetscScalar *x, *xdot;  in CEFunction()  local
129   PetscCall(VecGetArrayRead(Xdot, &xdot));  in CEFunction()
131   f[0] = xdot[0] + l * (x[0] - PetscCosReal(t));  in CEFunction()
133 …SC_COMM_WORLD," f(t=%g,x=%g,xdot=%g) = %g\n",(double)t,(double)x[0],(double)xdot[0],(double)f[0]));  in CEFunction()
136   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in CEFunction()
146   const PetscScalar *x, *xdot;  in CEJacobian()  local
150   PetscCall(VecGetArrayRead(Xdot, &xdot));  in CEJacobian()
154   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in CEJacobian()
208   const PetscScalar *x, *xdot;  in OregoFunction()  local
212   PetscCall(VecGetArrayRead(Xdot, &xdot));  in OregoFunction()
214   f[0] = xdot[0] - 77.27 * (x[1] + x[0] * (1. - 8.375e-6 * x[0] - x[1]));  in OregoFunction()
215   f[1] = xdot[1] - 1 / 77.27 * (x[2] - (1. + x[0]) * x[1]);  in OregoFunction()
216   f[2] = xdot[2] - 0.161 * (x[0] - x[2]);  in OregoFunction()
218   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in OregoFunction()
227   const PetscScalar *x, *xdot;  in OregoJacobian()  local
231   PetscCall(VecGetArrayRead(Xdot, &xdot));  in OregoJacobian()
243   PetscCall(VecRestoreArrayRead(Xdot, &xdot));  in OregoJacobian()