31   const PetscScalar *u;  in EventFunction()  local
34   PetscCall(VecGetArrayRead(U, &u));  in EventFunction()
36     fvalue[0] = PetscRealPart(u[1] - actx->lambda1 * u[0]);  in EventFunction()
38     fvalue[0] = PetscRealPart(u[1] - actx->lambda2 * u[0]);  in EventFunction()
40   PetscCall(VecRestoreArrayRead(U, &u));  in EventFunction()
48   const PetscScalar *u;  in ShiftGradients()  local
54   PetscCall(VecGetArrayRead(U, &u));  in ShiftGradients()
57     denorm1  = -actx->lambda1 * (u[0] - 100. * u[1]) + 1. * (10. * u[0] + u[1]);  in ShiftGradients()
58     denorm2  = -actx->lambda1 * (u[0] + 10. * u[1]) + 1. * (-100. * u[0] + u[1]);  in ShiftGradients()
59     A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm1 + 1.;  in ShiftGradients()
60     A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm1;  in ShiftGradients()
61     A1[1][0] = 110. * u[1] * 1. / denorm1;  in ShiftGradients()
62     A1[1][1] = -110. * u[0] * 1. / denorm1 + 1.;  in ShiftGradients()
64     A2[0] = 110. * u[1] * (-u[0]) / denorm2;  in ShiftGradients()
65     A2[1] = -110. * u[0] * (-u[0]) / denorm2;  in ShiftGradients()
67     denorm2  = -actx->lambda2 * (u[0] + 10. * u[1]) + 1. * (-100. * u[0] + u[1]);  in ShiftGradients()
68     A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm2 + 1;  in ShiftGradients()
69     A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm2;  in ShiftGradients()
70     A1[1][0] = 110. * u[1] * 1. / denorm2;  in ShiftGradients()
71     A1[1][1] = -110. * u[0] * 1. / denorm2 + 1.;  in ShiftGradients()
77   PetscCall(VecRestoreArrayRead(U, &u));  in ShiftGradients()
124   const PetscScalar *u, *udot;  in IFunction()  local
128   PetscCall(VecGetArrayRead(U, &u));  in IFunction()
133     f[0] = udot[0] - u[0] + 100 * u[1];  in IFunction()
134     f[1] = udot[1] - 10 * u[0] - u[1];  in IFunction()
136     f[0] = udot[0] - u[0] - 10 * u[1];  in IFunction()
137     f[1] = udot[1] + 100 * u[0] - u[1];  in IFunction()
140   PetscCall(VecRestoreArrayRead(U, &u));  in IFunction()
154   const PetscScalar *u, *udot;  in IJacobian()  local
157   PetscCall(VecGetArrayRead(U, &u));  in IJacobian()
173   PetscCall(VecRestoreArrayRead(U, &u));  in IJacobian()
193   PetscScalar *u;  in main()  local
237   PetscCall(VecGetArray(U, &u));  in main()
238   u[0] = 0;  in main()
239   u[1] = 1;  in main()
240   PetscCall(VecRestoreArray(U, &u));  in main()
273   PetscCall(VecGetArray(U, &u));  in main()
274   tmp[0] = u[0];  in main()
275   tmp[1] = u[1];  in main()
277   u[0] = 0 + delta;  in main()
278   u[1] = 1;  in main()
279   PetscCall(VecRestoreArray(U, &u));  in main()
283   PetscCall(VecGetArray(U, &u));  in main()
284   sensi[0] = (u[0] - tmp[0]) / delta;  in main()
285   sensi[1] = (u[1] - tmp[1]) / delta;  in main()
287   u[0] = 0;  in main()
288   u[1] = 1 + delta;  in main()
289   PetscCall(VecRestoreArray(U, &u));  in main()
293   PetscCall(VecGetArray(U, &u));  in main()
294   sensi[0] = (u[0] - tmp[0]) / delta;  in main()
295   sensi[1] = (u[1] - tmp[1]) / delta;  in main()
297   u[0]        = 0;  in main()
298   u[1]        = 1;  in main()
300   PetscCall(VecRestoreArray(U, &u));  in main()
304   PetscCall(VecGetArray(U, &u));  in main()
305   sensi[0] = (u[0] - tmp[0]) / delta;  in main()
306   sensi[1] = (u[1] - tmp[1]) / delta;  in main()
308   PetscCall(VecRestoreArray(U, &u));  in main()