33   PetscScalar *u;  in MyMonitor()  local
41     PetscCall(MatDenseGetColumn(sp, 2, &u));  in MyMonitor()
43 …5lf %20.15lf %20.15lf\n", (double)time, (double)PetscRealPart(u[0]), (double)PetscRealPart(u[1])));  in MyMonitor()
44     PetscCall(MatDenseRestoreColumn(sp, &u));  in MyMonitor()
53   const PetscScalar *u;  in EventFunction()  local
56   PetscCall(VecGetArrayRead(U, &u));  in EventFunction()
58     fvalue[0] = PetscRealPart(u[1] - actx->lambda1 * u[0]);  in EventFunction()
60     fvalue[0] = PetscRealPart(u[1] - actx->lambda2 * u[0]);  in EventFunction()
62   PetscCall(VecRestoreArrayRead(U, &u));  in EventFunction()
70   PetscScalar *u;  in ShiftGradients()  local
76   PetscCall(VecGetArray(U, &u));  in ShiftGradients()
79     denorm   = -actx->lambda1 * (u[0] - 100. * u[1]) + 1. * (10. * u[0] + u[1]);  in ShiftGradients()
80     A1[0][0] = 110. * u[1] * (-actx->lambda1) / denorm + 1.;  in ShiftGradients()
81     A1[1][0] = -110. * u[0] * (-actx->lambda1) / denorm;  in ShiftGradients()
82     A1[0][1] = 110. * u[1] * 1. / denorm;  in ShiftGradients()
83     A1[1][1] = -110. * u[0] * 1. / denorm + 1.;  in ShiftGradients()
85     A2[0] = 110. * u[1] * (-u[0]) / denorm;  in ShiftGradients()
86     A2[1] = -110. * u[0] * (-u[0]) / denorm;  in ShiftGradients()
88     denorm   = -actx->lambda2 * (u[0] + 10. * u[1]) + 1. * (-100. * u[0] + u[1]);  in ShiftGradients()
89     A1[0][0] = 110. * u[1] * (actx->lambda2) / denorm + 1;  in ShiftGradients()
90     A1[1][0] = -110. * u[0] * (actx->lambda2) / denorm;  in ShiftGradients()
91     A1[0][1] = -110. * u[1] * 1. / denorm;  in ShiftGradients()
92     A1[1][1] = 110. * u[0] * 1. / denorm + 1.;  in ShiftGradients()
98   PetscCall(VecRestoreArray(U, &u));  in ShiftGradients()
147   const PetscScalar *u, *udot;  in IFunction()  local
151   PetscCall(VecGetArrayRead(U, &u));  in IFunction()
156     f[0] = udot[0] - u[0] + 100 * u[1];  in IFunction()
157     f[1] = udot[1] - 10 * u[0] - u[1];  in IFunction()
159     f[0] = udot[0] - u[0] - 10 * u[1];  in IFunction()
160     f[1] = udot[1] + 100 * u[0] - u[1];  in IFunction()
163   PetscCall(VecRestoreArrayRead(U, &u));  in IFunction()
177   const PetscScalar *u, *udot;  in IJacobian()  local
180   PetscCall(VecGetArrayRead(U, &u));  in IJacobian()
196   PetscCall(VecRestoreArrayRead(U, &u));  in IJacobian()
223   PetscScalar *u;  in main()  local
271   PetscCall(VecGetArray(U, &u));  in main()
272   u[0] = 0;  in main()
273   u[1] = 1;  in main()
274   PetscCall(VecRestoreArray(U, &u));  in main()