88   const PetscScalar *vl, *vr, *u;  in RHSHessianProductUU()  local
96   PetscCall(VecGetArrayRead(Vl[0], &vl));  in RHSHessianProductUU()
106       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in RHSHessianProductUU()
110   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in RHSHessianProductUU()
118   const PetscScalar *vl, *vr, *u;  in RHSHessianProductUP()  local
125   PetscCall(VecGetArrayRead(Vl[0], &vl));  in RHSHessianProductUP()
134       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdP[i][j][k] * vr[k];  in RHSHessianProductUP()
138   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in RHSHessianProductUP()
146   const PetscScalar *vl, *vr, *u;  in RHSHessianProductPU()  local
153   PetscCall(VecGetArrayRead(Vl[0], &vl));  in RHSHessianProductPU()
162       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in RHSHessianProductPU()
166   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in RHSHessianProductPU()
249   const PetscScalar *vl, *vr, *u;  in IHessianProductUU()  local
257   PetscCall(VecGetArrayRead(Vl[0], &vl));  in IHessianProductUU()
267       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in IHessianProductUU()
271   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in IHessianProductUU()
279   const PetscScalar *vl, *vr, *u;  in IHessianProductUP()  local
286   PetscCall(VecGetArrayRead(Vl[0], &vl));  in IHessianProductUP()
295       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdP[i][j][k] * vr[k];  in IHessianProductUP()
299   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in IHessianProductUP()
307   const PetscScalar *vl, *vr, *u;  in IHessianProductPU()  local
314   PetscCall(VecGetArrayRead(Vl[0], &vl));  in IHessianProductPU()
323       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in IHessianProductPU()
327   PetscCall(VecRestoreArrayRead(Vl[0], &vl));  in IHessianProductPU()