88   const PetscScalar *vl, *vr, *u;  in RHSHessianProductUU()  local
97   PetscCall(VecGetArrayRead(Vr, &vr));  in RHSHessianProductUU()
106       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in RHSHessianProductUU()
111   PetscCall(VecRestoreArrayRead(Vr, &vr));  in RHSHessianProductUU()
118   const PetscScalar *vl, *vr, *u;  in RHSHessianProductUP()  local
126   PetscCall(VecGetArrayRead(Vr, &vr));  in RHSHessianProductUP()
134       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdP[i][j][k] * vr[k];  in RHSHessianProductUP()
139   PetscCall(VecRestoreArrayRead(Vr, &vr));  in RHSHessianProductUP()
146   const PetscScalar *vl, *vr, *u;  in RHSHessianProductPU()  local
154   PetscCall(VecGetArrayRead(Vr, &vr));  in RHSHessianProductPU()
162       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in RHSHessianProductPU()
167   PetscCall(VecRestoreArrayRead(Vr, &vr));  in RHSHessianProductPU()
249   const PetscScalar *vl, *vr, *u;  in IHessianProductUU()  local
258   PetscCall(VecGetArrayRead(Vr, &vr));  in IHessianProductUU()
267       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in IHessianProductUU()
272   PetscCall(VecRestoreArrayRead(Vr, &vr));  in IHessianProductUU()
279   const PetscScalar *vl, *vr, *u;  in IHessianProductUP()  local
287   PetscCall(VecGetArrayRead(Vr, &vr));  in IHessianProductUP()
295       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdP[i][j][k] * vr[k];  in IHessianProductUP()
300   PetscCall(VecRestoreArrayRead(Vr, &vr));  in IHessianProductUP()
307   const PetscScalar *vl, *vr, *u;  in IHessianProductPU()  local
315   PetscCall(VecGetArrayRead(Vr, &vr));  in IHessianProductPU()
323       for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k];  in IHessianProductPU()
328   PetscCall(VecRestoreArrayRead(Vr, &vr));  in IHessianProductPU()