| /petsc/src/ts/tutorials/ |
| H A D | ex20opt_p.c | 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() [all …]
|
| H A D | ex20opt_ic.c | 80 const PetscScalar *vl, *vr, *u; in RHSHessianProductUU() local 89 PetscCall(VecGetArrayRead(Vr, &vr)); in RHSHessianProductUU() 98 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k]; in RHSHessianProductUU() 102 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductUU() 177 const PetscScalar *vl, *vr, *u; in IHessianProductUU() local 186 PetscCall(VecGetArrayRead(Vr, &vr)); in IHessianProductUU() 194 for (i = 0; i < 2; i++) vhv[j] += vl[i] * dJdU[i][j][k] * vr[k]; in IHessianProductUU() 198 PetscCall(VecRestoreArrayRead(Vr, &vr)); in IHessianProductUU()
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | eptorsion3.c | 261 PetscReal zero = 0.0, vb, vl, vr, vt, dvdx, dvdy, flin = 0.0, fquad = 0.0; in FormFunction() local 275 vr = zero; in FormFunction() 278 if (i < nx - 1 && j > -1) vr = x[k + 1]; in FormFunction() 280 dvdx = (vr - v) / hx; in FormFunction() 283 flin -= cdiv3 * (v + vr + vt); in FormFunction() 333 PetscReal vb, vl, vr, vt, dvdx, dvdy; in FormGradient() local 349 vr = zero; in FormGradient() 352 if (i < nx - 1 && j > -1) vr = x[k + 1]; in FormGradient() 354 dvdx = (vr - v) / hx; in FormGradient() 537 PetscReal v, vb, vl, vr, vt, hxhx, hyhy; in HessianProduct() local [all …]
|
| H A D | eptorsion1.c | 245 PetscReal zero = 0.0, vb, vl, vr, vt, dvdx, dvdy, flin = 0.0, fquad = 0.0; in FormFunction() local 259 vr = zero; in FormFunction() 262 if (i < nx - 1 && j > -1) vr = x[k + 1]; in FormFunction() 264 dvdx = (vr - v) / hx; in FormFunction() 267 flin -= cdiv3 * (v + vr + vt); in FormFunction() 317 PetscReal vb, vl, vr, vt, dvdx, dvdy; in FormGradient() local 333 vr = zero; in FormGradient() 336 if (i < nx - 1 && j > -1) vr = x[k + 1]; in FormGradient() 338 dvdx = (vr - v) / hx; in FormGradient() 521 PetscReal v, vb, vl, vr, vt, hxhx, hyhy; in HessianProduct() local [all …]
|
| H A D | eptorsion2.c | 205 PetscReal v, vb, vl, vr, vt, dvdx, dvdy; in FormFunctionGradient() local 248 vr = zero; in FormFunctionGradient() 251 if (i < mx - 1 && j > -1) vr = x[k + 1]; in FormFunctionGradient() 253 dvdx = (vr - v) / hx; in FormFunctionGradient() 271 flin -= cdiv3 * (v + vr + vt); in FormFunctionGradient()
|
| H A D | eptorsion2f.F90 | 119 PetscReal v, vb, vl, vr, vt, dvdx 181 vr = zero 184 if (i < mx - 1 .and. j > -1) vr = lx_v(k + 2) 186 dvdx = (vr - v)/hx 204 flin = flin - cdiv3*(v + vr + vt)
|
| /petsc/src/mat/tests/ |
| H A D | ex221.c | 289 Vec vr, vl; in main() local 291 PetscCall(MatCreateVecs(A, &vr, &vl)); in main() 293 PetscCall(VecSetRandom(vr, NULL)); in main() 296 PetscCall(VecSet(vr, test % 2 ? 0.15 : 1.0 / 0.15)); in main() 299 PetscCall(MatDiagonalScale(A, vl, vr)); in main() 300 PetscCall(MatDiagonalScale(S, vl, vr)); in main() 301 PetscCall(VecDestroy(&vr)); in main()
|
| /petsc/src/ts/tutorials/optimal_control/ |
| H A D | ex1.c | 99 const PetscScalar *v, *w, *vl, *vr, *u; in RHSHessianProductPP() local 110 PetscCall(VecGetArrayRead(Vr, &vr)); 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() 128 PetscCall(VecRestoreArrayRead(Vr, &vr)); in RHSHessianProductPP() 137 const PetscScalar *v, *w, *vr, *u; in IntegrandHessianProductUU() local 147 PetscCall(VecGetArrayRead(Vr, &vr)); in IntegrandHessianProductUU() 155 for (k = 0; k < 2; k++) vhv[j] += dRudU[j][k] * vr[k]; in IntegrandHessianProductUU() 158 PetscCall(VecRestoreArrayRead(Vr, &vr)); in IntegrandHessianProductUU()
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busdmnetwork.c | 60 PetscScalar vr; /* Real component of bus voltage */ member 194 bus[i * 9 + j].vr = varr[2 * j]; in read_data() 397 xarr[offset] = bus->vr; in SetInitialGuess() 400 Vr = bus->vr; in SetInitialGuess() 553 Vr0 = bus->vr; in FormIFunction() 803 Vr0 = bus->vr; in AlgFunction()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 82 static PetscErrorCode VecMergedOpsShort_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec… in VecMergedOpsShort_Private() argument 91 PetscCall(VecGetArray(vr, (PetscScalar **)&pr)); in VecMergedOpsShort_Private() 281 PetscCall(VecRestoreArray(vr, (PetscScalar **)&pr)); in VecMergedOpsShort_Private() 304 static PetscErrorCode VecMergedOps_Private(Vec vx, Vec vr, Vec vz, Vec vw, Vec vp, Vec vq, Vec vc, … in VecMergedOps_Private() argument 313 PetscCall(VecGetArray(vr, (PetscScalar **)&pr)); in VecMergedOps_Private() 510 PetscCall(VecRestoreArray(vr, (PetscScalar **)&pr)); in VecMergedOps_Private()
|
| /petsc/src/binding/petsc4py/test/ |
| H A D | test_mat_py.py | 285 def diagonalScale(self, mat, vl, vr): argument 288 if vr: 289 self.D.pointwiseMult(self.D, vr)
|
| /petsc/src/binding/petsc4py/src/petsc4py/PETSc/ |
| H A D | KSP.pyx | 2033 cdef PetscVec *vr=NULL, *vl=NULL 2038 CHKERR(KSPCreateVecs(self.ksp, nr, &vr, nl, &vr)) 2041 vecsr.append(ref_Vec(vr[i])) 2045 if nr > 0 and vr != NULL: 2046 VecDestroyVecs(nr, &vr) # XXX errors?
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aij.c | 2275 PetscScalar vc, vr; in MatIsTranspose_SeqAIJ() local 2280 vr = vb[bptr[idc]]; in MatIsTranspose_SeqAIJ() 2281 if (i != idr || PetscAbsScalar(vc - vr) > tol) { in MatIsTranspose_SeqAIJ() 2327 PetscScalar vc, vr; in MatIsHermitianTranspose_SeqAIJ() local 2332 vr = vb[bptr[idc]]; in MatIsHermitianTranspose_SeqAIJ() 2333 if (i != idr || PetscAbsScalar(vc - PetscConj(vr)) > tol) { in MatIsHermitianTranspose_SeqAIJ()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 4612 (hyb-init/turb-vr 10.) 6217 (udf/vr-rate-fcn "none")
|