43   const PetscScalar *u, *udot;  in IFunctionPassive()  local
47   PetscCall(VecGetArrayRead(U, &u));  in IFunctionPassive()
50   f[0] = udot[0] + ctx->k * u[0] * u[1];  in IFunctionPassive()
51   f[1] = udot[1] + ctx->k * u[0] * u[1];  in IFunctionPassive()
52   f[2] = udot[2] - ctx->k * u[0] * u[1];  in IFunctionPassive()
55   PetscCall(VecRestoreArrayRead(U, &u));  in IFunctionPassive()
65   const PetscScalar *u, *udot;  in IFunctionActive1()  local
72   PetscCall(VecGetArrayRead(U, &u));  in IFunctionActive1()
78   u_a[0] <<= u[0];  in IFunctionActive1()
79   u_a[1] <<= u[1];  in IFunctionActive1()
80   u_a[2] <<= u[2]; /* Mark independence */  in IFunctionActive1()
92   PetscCall(VecRestoreArrayRead(U, &u));  in IFunctionActive1()
102   const PetscScalar *u, *udot;  in IFunctionActive2()  local
109   PetscCall(VecGetArrayRead(U, &u));  in IFunctionActive2()
118   f_a[0] = udot_a[0] + ctx->k * u[0] * u[1];  in IFunctionActive2()
119   f_a[1] = udot_a[1] + ctx->k * u[0] * u[1];  in IFunctionActive2()
120   f_a[2] = udot_a[2] - ctx->k * u[0] * u[1];  in IFunctionActive2()
129   PetscCall(VecRestoreArrayRead(U, &u));  in IFunctionActive2()
140   const PetscScalar *u;  in IJacobian()  local
143   PetscCall(VecGetArrayRead(U, &u));  in IJacobian()
144   PetscCall(PetscAdolcComputeIJacobian(1, 2, A, u, a, appctx->adctx));  in IJacobian()
145   PetscCall(VecRestoreArrayRead(U, &u));  in IJacobian()
155   PetscScalar       *u, d0, q;  in Solution()  local
159   PetscCall(VecGetArray(U, &u));  in Solution()
163   u[0] = uinit[0] / (1.0 + uinit[1] * q);  in Solution()
164   u[1] = u[0] - d0;  in Solution()
165   u[2] = uinit[1] + uinit[2] - u[1];  in Solution()
166   PetscCall(VecRestoreArray(U, &u));  in Solution()
180   PetscScalar      *u;  in main()  local
212   PetscCall(VecGetArray(ctx.initialsolution, &u));  in main()
213   u[0] = 1;  in main()
214   u[1] = .7;  in main()
215   u[2] = 0;  in main()
216   PetscCall(VecRestoreArray(ctx.initialsolution, &u));  in main()