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     denorm   = -actx->lambda1 * (u[0] - 100. * u[1]) + 1. * (10. * u[0] + u[1]);  in ShiftGradients()
58     A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm + 1.;  in ShiftGradients()
59     A1[0][1] = -110. * u[0] * (-actx->lambda1) / denorm;  in ShiftGradients()
60     A1[1][0] = 110. * u[1] * 1. / denorm;  in ShiftGradients()
61     A1[1][1] = -110. * u[0] * 1. / denorm + 1.;  in ShiftGradients()
63     A2[0] = 110. * u[1] * (-u[0]) / denorm;  in ShiftGradients()
64     A2[1] = -110. * u[0] * (-u[0]) / denorm;  in ShiftGradients()
66     denorm   = -actx->lambda2 * (u[0] + 10. * u[1]) + 1. * (-100. * u[0] + u[1]);  in ShiftGradients()
67     A1[0][0] = 110. * u[1] * (actx->lambda2) / denorm + 1;  in ShiftGradients()
68     A1[0][1] = -110. * u[0] * (actx->lambda2) / denorm;  in ShiftGradients()
69     A1[1][0] = -110. * u[1] * 1. / denorm;  in ShiftGradients()
70     A1[1][1] = 110. * u[0] * 1. / denorm + 1.;  in ShiftGradients()
76   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()
200   PetscScalar *u, *v;  in main()  local
245   PetscCall(VecGetArray(U, &u));  in main()
246   u[0] = 0;  in main()
247   u[1] = 1;  in main()
248   PetscCall(VecRestoreArray(U, &u));  in main()
295   PetscCall(VecGetArray(lambda[0], &u));  in main()
296   u[0] = 1.;  in main()
297   PetscCall(VecRestoreArray(lambda[0], &u));  in main()
298   PetscCall(VecGetArray(lambda[1], &u));  in main()
299   u[1] = 1.;  in main()
300   PetscCall(VecRestoreArray(lambda[1], &u));  in main()
316   PetscCall(VecGetArray(mu[0], &u));  in main()
319 …RLD, f, "%20.15lf %20.15lf %20.15lf\n", (double)tend, (double)PetscRealPart(u[0]), (double)PetscRe…  in main()
320   PetscCall(VecRestoreArray(mu[0], &u));  in main()