| /petsc/src/vec/vec/impls/seq/ftn-kernels/ |
| H A D | fmaxpy.F90 | 13 PetscScalar, intent(in) :: y0(*), y1(*), y2(*), y3(*) 17 PETSC_AssertAlignx(16, y0(1)) 22 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n) + a2*y2(1:n) + a3*y3(1:n)) 30 PetscScalar, intent(in) :: y0(*), y1(*), y2(*) 34 PETSC_AssertAlignx(16, y0(1)) 38 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n) + a2*y2(1:n)) 46 PetscScalar, intent(in) :: y0(*), y1(*) 50 PETSC_AssertAlignx(16, y0(1)) 53 x(1:n) = x(1:n) + (a0*y0(1:n) + a1*y1(1:n))
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | pyr_tet.geo | 2 y0 = 0; y1 = 1; 20 Point(1) = {x0, y0, z0}; 21 Point(2) = {x1, y0, z0}; 24 Point(5) = {x0, y0, z1}; 25 Point(6) = {x1, y0, z1};
|
| H A D | gmsh-3d.geo | 6 y0 = 0; y1 = 1; 24 Point(1) = {x0, y0, z0}; 25 Point(2) = {x1, y0, z0}; 28 Point(5) = {x0, y0, z1}; 29 Point(6) = {x1, y0, z1};
|
| H A D | gmsh-hex.geo | 2 y0 = 0; y1 = 1; 20 Point(1) = {x0, y0, z0}; 21 Point(2) = {x1, y0, z0}; 24 Point(5) = {x0, y0, z1}; 25 Point(6) = {x1, y0, z1};
|
| /petsc/src/ts/tutorials/output/ |
| H A D | ex20adj_2.out | 2 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 8 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex20adj_1.out | 2 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 8 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex20adj_imex.out | 2 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 8 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex16fwd_1.out | 7 forward sensitivity: d[y(tf) z(tf)]/d[y0 z0 mu]
|
| H A D | ex20fwd_1.out | 9 forward sensitivity: d[y(tf) z(tf)]/d[y0 z0 mu]
|
| H A D | ex20adj_5.out | 58 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 64 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex20adj_6.out | 57 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 63 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex20adj_3.out | 86 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 92 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| H A D | ex20adj_4.out | 100 sensitivity wrt initial conditions: d[y(tf)]/d[y0] d[y(tf)]/d[z0] 106 sensitivity wrt initial conditions: d[z(tf)]/d[y0] d[z(tf)]/d[z0]
|
| /petsc/src/snes/utils/dm/ |
| H A D | dminterpolatesnes.c | 457 const PetscScalar y0 = vertices[1]; in QuadMap_Private() local 465 const PetscScalar g_1 = y1 - y0; in QuadMap_Private() 467 const PetscScalar g_3 = y3 - y0; in QuadMap_Private() 469 const PetscScalar g_01 = y2 - y1 - y3 + y0; in QuadMap_Private() 481 real[1] = y0 + g_1 * p0 + g_3 * p1 + g_01 * p0 * p1; in QuadMap_Private() 494 const PetscScalar y0 = vertices[1]; in QuadJacobian_Private() local 502 const PetscScalar g_01 = y2 - y1 - y3 + y0; in QuadJacobian_Private() 515 values[2] = (y1 - y0 + g_01 * y) * 0.5; in QuadJacobian_Private() 516 values[3] = (y3 - y0 + g_01 * x) * 0.5; in QuadJacobian_Private() 627 const PetscScalar y0 = vertices[1]; in HexMap_Private() local [all …]
|
| /petsc/src/snes/tutorials/ |
| H A D | ex70.c | 185 Vec y0, y1; in StokesExactSolution() local 189 PetscCall(VecGetSubVector(s->y, s->isg[0], &y0)); in StokesExactSolution() 190 PetscCall(VecGetOwnershipRange(y0, &start, &end)); in StokesExactSolution() 198 PetscCall(VecSetValue(y0, row, val, INSERT_VALUES)); in StokesExactSolution() 200 PetscCall(VecAssemblyBegin(y0)); in StokesExactSolution() 201 PetscCall(VecAssemblyEnd(y0)); in StokesExactSolution() 202 PetscCall(VecRestoreSubVector(s->y, s->isg[0], &y0)); in StokesExactSolution() 686 Vec y0, y1; in StokesCalcError() local 693 PetscCall(VecGetSubVector(s->y, s->isg[0], &y0)); in StokesCalcError() 694 PetscCall(VecNorm(y0, NORM_2, &val)); in StokesCalcError() [all …]
|
| /petsc/src/tao/constrained/impls/admm/ |
| H A D | admm.h | 18 Vec residual, y, yold, y0, yhat, yhatold, constraint; member
|
| H A D | admm.c | 84 PetscCall(VecWAXPY(tempJR2, -1., am->y, am->y0)); in AdaptiveADMMPenaltyUpdate() 446 PetscCall(VecCopy(am->y, am->y0)); in TaoSolve_ADMM() 462 PetscCall(VecCopy(am->y, am->y0)); in TaoSolve_ADMM() 583 if (!am->y0) { in TaoSetUp_ADMM() 584 PetscCall(VecDuplicate(am->Ax, &am->y0)); in TaoSetUp_ADMM() 585 PetscCall(VecSet(am->y0, 0.0)); in TaoSetUp_ADMM() 631 PetscCall(VecDestroy(&am->y0)); in TaoDestroy_ADMM()
|
| /petsc/src/vec/vec/impls/seq/kokkos/ |
| H A D | veckok.kokkos.cxx | 389 ConstPetscScalarKokkosView xv, y0, y1, y2, y3, y4, y5, y6, y7; in VecMultiDot_Verbose() local 398 PetscCall(VecGetKokkosView(yp[0], &y0)); in VecMultiDot_Verbose() 410 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() 415 PetscCall(VecRestoreKokkosView(yp[0], &y0)); in VecMultiDot_Verbose() 428 if (rem > 0) PetscCall(VecGetKokkosView(yp[0], &y0)); in VecMultiDot_Verbose() 441 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() 449 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() 457 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() 465 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() 472 …lsum0 += xv(i) * PetscConj(y0(i)); lsum1 += xv(i) * PetscConj(y1(i)); lsum2 += xv(i) * PetscConj(y… in VecMultiDot_Verbose() [all …]
|
| /petsc/src/dm/impls/moab/ |
| H A D | dmmbfem.cxx | 453 …const PetscReal x0 = coords[/*0 * 3 +*/ 0], y0 = coords[/*0 * 3 +*/ 1], z0 = coords[/*0 * 3 +*/ 2]; in Compute_Lagrange_Basis_3D_Internal() local 462 jacobian[3] = coords[1 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal() 463 jacobian[4] = coords[2 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal() 464 jacobian[5] = coords[3 * 3 + 1] - y0; in Compute_Lagrange_Basis_3D_Internal()
|
| /petsc/src/ts/tutorials/ |
| H A D | ex77.c | 119 const PetscReal y0 = X[1]; in trig_trig_x() local 120 const PetscReal R0 = PetscSqrtReal(x0 * x0 + y0 * y0); in trig_trig_x() 121 const PetscReal theta0 = PetscAtan2Real(y0, x0); in trig_trig_x()
|