Home
last modified time | relevance | path

Searched refs:vhv (Results 1 – 3 of 3) sorted by relevance

/petsc/src/ts/tutorials/
H A Dex20opt_p.c89 PetscScalar *vhv; in RHSHessianProductUU() local
98 PetscCall(VecGetArray(VHV[0], &vhv)); in RHSHessianProductUU()
104 vhv[j] = 0; in RHSHessianProductUU()
106 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k]; in RHSHessianProductUU()
112 PetscCall(VecRestoreArray(VHV[0], &vhv)); in RHSHessianProductUU()
119 PetscScalar *vhv; in RHSHessianProductUP() local
127 PetscCall(VecGetArray(VHV[0], &vhv)); in RHSHessianProductUP()
132 vhv[j] = 0; in RHSHessianProductUP()
134 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdP[i][j][k] * vr[k]; in RHSHessianProductUP()
140 PetscCall(VecRestoreArray(VHV[0], &vhv)); in RHSHessianProductUP()
[all …]
H A Dex20opt_ic.c81 PetscScalar *vhv; in RHSHessianProductUU() local
90 PetscCall(VecGetArray(VHV[0], &vhv)); in RHSHessianProductUU()
96 vhv[j] = 0; in RHSHessianProductUU()
98 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k]; in RHSHessianProductUU()
103 PetscCall(VecRestoreArray(VHV[0], &vhv)); in RHSHessianProductUU()
178 PetscScalar *vhv; in IHessianProductUU() local
187 PetscCall(VecGetArray(VHV[0], &vhv)); in IHessianProductUU()
192 vhv[j] = 0; in IHessianProductUU()
194 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k]; in IHessianProductUU()
199 PetscCall(VecRestoreArray(VHV[0], &vhv)); in IHessianProductUU()
/petsc/src/ts/tutorials/optimal_control/
H A Dex1.c100 PetscScalar *vhv; in RHSHessianProductPP() local
112 PetscCall(VecGetArray(VHV[0], &vhv)); in RHSHessianProductPP()
122 vhv[2 * step + j] = 0; in RHSHessianProductPP()
124 for (i = 0; i < 2; i++) vhv[2 * step + j] += vl[i] * dJpdP[i][j][k] * vr[2 * step + k]; in RHSHessianProductPP()
129 PetscCall(VecRestoreArray(VHV[0], &vhv)); in RHSHessianProductPP()
138 PetscScalar *vhv; in IntegrandHessianProductUU() local
148 PetscCall(VecGetArray(VHV[0], &vhv)); in IntegrandHessianProductUU()
154 vhv[j] = 0; in IntegrandHessianProductUU()
155 for (k = 0; k < 2; k++) vhv[j] += dRudU[j][k] * vr[k]; in IntegrandHessianProductUU()
159 PetscCall(VecRestoreArray(VHV[0], &vhv)); in IntegrandHessianProductUU()