| /petsc/src/ksp/ksp/tutorials/amrex/ |
| H A D | MyTest.cxx | 96 ParmParse pp; in readParameters() local 97 pp.query("max_level", max_level); in readParameters() 98 pp.query("n_cell", n_cell); in readParameters() 99 pp.query("max_grid_size", max_grid_size); in readParameters() 100 pp.query("is_periodic", is_periodic); in readParameters() 101 pp.query("eb_is_dirichlet", eb_is_dirichlet); in readParameters() 103 pp.query("plot_file", plot_file_name); in readParameters() 113 pp.queryarr("scalars", scalars); in readParameters() 115 pp.query("verbose", verbose); in readParameters() 116 pp.query("bottom_verbose", bottom_verbose); in readParameters() [all …]
|
| H A D | initEB.cxx | 16 ParmParse pp("eb2"); in initializeEB() local 18 pp.get("geom_type", geom_type); in initializeEB()
|
| /petsc/src/ts/tests/ |
| H A D | ex30.c | 407 for (PetscInt pp = 0; pp < NN; pp++) { in PostStep() local 408 …scReal v2 = 0, fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp] * … in PostStep() 409 for (PetscInt i = 0; i < dim; ++i) v2 += PetscSqr(coords[pp * dim + i]); in PostStep() 411 moments[1] += w * ctx->v_0 * coords[pp * dim + 1]; // z-momentum in PostStep() 428 for (PetscInt pp = 0; pp < NN; pp++) { in PostStep() local 429 …PetscReal fact = (dim == 2) ? 2.0 * PETSC_PI * coords[pp * dim + 0] : 1, w = fact * wp[pp] * ctx->… in PostStep() 577 for (PetscInt pj = 0, pp = 0; pj < Npj; pj++) { in go() local 579 for (PetscInt pi = 0; pi < Npi; pi++, pp++) { in go() 581 wp_t[grid][tid][pp] = 0; in go() 583 xx_t[grid][tid][pp] = lo[0] + hp[0] / 2.0 + pi * hp[0]; in go() [all …]
|
| /petsc/src/sys/classes/viewer/impls/socket/mex-scripts/ |
| H A D | bread.c | 82 char *pp = (char *)p; in PetscBinaryRead() local 95 err = (int)read(fd, pp, (int)wsize); in PetscBinaryRead() 100 pp += err; in PetscBinaryRead() 126 char *pp = (char *)p; in PetscBinaryWrite() local 147 err = (int)write(fd, pp, (int)wsize); in PetscBinaryWrite() 155 pp += err; in PetscBinaryWrite()
|
| /petsc/src/ksp/ksp/impls/cg/pipecg2/ |
| H A D | pipecg2.c | 85 PetscScalar *PETSC_RESTRICT pp, *PETSC_RESTRICT pq; in VecMergedOpsShort_Private() local 94 PetscCall(VecGetArray(vp, (PetscScalar **)&pp)); in VecMergedOpsShort_Private() 120 pp[j] = pu[j]; in VecMergedOpsShort_Private() 126 px[j] = px[j] + alpha0 * pp[j]; in VecMergedOpsShort_Private() 141 pp[j] = pu[j] + beta1 * pp[j]; in VecMergedOpsShort_Private() 145 px[j] = px[j] + alpha1 * pp[j]; in VecMergedOpsShort_Private() 175 pp[j] = pu[j]; in VecMergedOpsShort_Private() 181 px[j] = px[j] + alpha0 * pp[j]; in VecMergedOpsShort_Private() 196 pp[j] = pu[j] + beta1 * pp[j]; in VecMergedOpsShort_Private() 200 px[j] = px[j] + alpha1 * pp[j]; in VecMergedOpsShort_Private() [all …]
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex7.c | 286 for (int pi = 0, pp = 0; pi < Np2[0]; pi++) { in go() local 287 for (int pj = 0; pj < Np2[1]; pj++, pp++) { in go() 288 xx_t[tid][pp] = lo[0] + hp[0] / 2. + pi * hp[0]; in go() 289 yy_t[tid][pp] = lo[1] + hp[1] / 2. + pj * hp[1]; in go() 291 PetscReal x[] = {xx_t[tid][pp], yy_t[tid][pp]}; in go() 292 PetscCall(maxwellian(2, x, 1.0, vol / (PetscReal)Np_t[tid], &wp_t[tid][pp])); in go() 295 moments_0[0] += wp_t[tid][pp]; in go() 296 moments_0[1] += wp_t[tid][pp] * xx_t[tid][pp]; // x-momentum in go() 297 moments_0[2] += wp_t[tid][pp] * (PetscSqr(xx_t[tid][pp]) + PetscSqr(yy_t[tid][pp])); in go()
|
| /petsc/src/ksp/pc/impls/wb/ |
| H A D | wb.c | 20 …PetscInt mwidth, nwidth, pwidth, cnt, mp, np, pp, Ntotal, gl[26], *globals, Ng, *IIi… in DMDAGetWireBasketInterpolation() local 35 …PetscCall(DMDAGetInfo(da, &dim, NULL, NULL, NULL, &mp, &np, &pp, &dof, NULL, NULL, NULL, NULL, NUL… in DMDAGetWireBasketInterpolation() 230 …(np + 1) * (pp + 1) + mp * np * (pp + 1) + mp * pp * (np + 1) + np * pp * (mp + 1) + mp * (np + 1)… in DMDAGetWireBasketInterpolation() 292 PetscCall(PetscMalloc1(26 * mp * np * pp, &globals)); in DMDAGetWireBasketInterpolation() 297 for (i = 0, cnt = 0; i < 26 * mp * np * pp; i++) { in DMDAGetWireBasketInterpolation() 368 …PetscInt mwidth, nwidth, pwidth, cnt, mp, np, pp, Ntotal, gl[6], *globals, Ng, *IIin… in DMDAGetFaceInterpolation() local 383 …PetscCall(DMDAGetInfo(da, &dim, NULL, NULL, NULL, &mp, &np, &pp, &dof, NULL, NULL, NULL, NULL, NUL… in DMDAGetFaceInterpolation() 552 Ntotal = mp * np * (pp + 1) + mp * pp * (np + 1) + np * pp * (mp + 1); in DMDAGetFaceInterpolation() 581 PetscCall(PetscMalloc1(6 * mp * np * pp, &globals)); in DMDAGetFaceInterpolation() 586 for (i = 0, cnt = 0; i < 6 * mp * np * pp; i++) { in DMDAGetFaceInterpolation()
|
| /petsc/doc/community/meetings/2025/ |
| H A D | index.md | 41 …tube.com/watch?v=vl5WOooGg_M&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=14&pp=iAQB) | 42 …tube.com/watch?v=efOjxXJeEXU&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=10&pp=iAQB) | 43 …outube.com/watch?v=-TYu8x8hkcU&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=9&pp=iAQB) | 44 …tube.com/watch?v=XwEZ6-VqEt8&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=18&pp=iAQB) | 46 …utube.com/watch?v=eKks1IXHQ2E&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=17&pp=iAQB) | 47 ….youtube.com/watch?v=IrQ0ptBbzSU&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=11&pp=iAQB) | 48 …tube.com/watch?v=dHF2D2NIdgA&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=15&pp=iAQB) | 50 …www.youtube.com/watch?v=MDYmajrhLsw&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=3&pp=iAQB) | 51 …ube.com/watch?v=4T4sU4pmfEc&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=1&t=26s&pp=iAQB) | 52 …utube.com/watch?v=l-z2Hj9pGtY&list=PLgFMPm27S9JrnzQmEZq54TGe5Iv8P6V_V&index=4&pp=iAQB) | [all …]
|
| /petsc/src/ts/tutorials/ |
| H A D | ex14.c | 803 const PetscReal pp = phi[l], *dp = dphi[l]; in THIFunctionLocal_3D() local 804 …* dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * du[2] + pp * jw * thi->rhog * d… in THIFunctionLocal_3D() 805 …* dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * dv[2] + pp * jw * thi->rhog * d… in THIFunctionLocal_3D() 806 fn[l]->u += pp * jw * udot * thi->inertia * pp; in THIFunctionLocal_3D() 807 fn[l]->v += pp * jw * vdot * thi->inertia * pp; in THIFunctionLocal_3D() 842 const PetscReal pp = phi[l]; in THIFunctionLocal_3D() local 843 fn[ls + l]->u += pp * jw * beta2 * u; in THIFunctionLocal_3D() 844 fn[ls + l]->v += pp * jw * beta2 * v; in THIFunctionLocal_3D() 1136 const PetscReal pp = phi[l], *restrict dp = dphi[l]; in THIJacobianLocal_Momentum() local 1153 Ke[l * 2 + 0][ll * 2 + 0] += pp * jw * thi->inertia * pp; in THIJacobianLocal_Momentum() [all …]
|
| /petsc/src/sys/fileio/ |
| H A D | sysio.c | 402 const char *pp = (char *)p; in PetscBinaryWrite() local 433 pp = fname; in PetscBinaryWrite() 443 pv = (PetscReal *)pp; in PetscBinaryWrite() 445 pp = (char *)ppp; in PetscBinaryWrite() 472 err = (size_t)write(fd, pp, wsize); in PetscBinaryWrite() 476 pp += wsize; in PetscBinaryWrite()
|
| /petsc/share/petsc/matlab/ |
| H A D | laplacian.m | 327 pp = p1(mod(p1,2)==0); variable 328 if ~isempty(pp) 359 pp = p2(mod(p2,2)==0); variable 360 if ~isempty(pp) 390 pp = p1(mod(p3,2)==0); variable 391 if ~isempty(pp)
|
| /petsc/src/snes/tutorials/network/power/ |
| H A D | case9.m | 9 % Transactions on, vol. 26, no. 1, pp. 12-19, Feb. 2011
|
| /petsc/src/snes/tutorials/ |
| H A D | ex48.c | 801 const PetscReal pp = phi[l], *dp = dphi[l]; in THIFunctionLocal() local 802 …* dv[1]) + dp[1] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * du[2] + pp * jw * thi->rhog * d… in THIFunctionLocal() 803 …* dv[1]) + dp[0] * jw * eta * (du[1] + dv[0]) + dp[2] * jw * eta * dv[2] + pp * jw * thi->rhog * d… in THIFunctionLocal() 838 const PetscReal pp = phi[l]; in THIFunctionLocal() local 839 fn[ls + l]->u += pp * jw * beta2 * u; in THIFunctionLocal() 840 fn[ls + l]->v += pp * jw * beta2 * v; in THIFunctionLocal() 1004 const PetscReal pp = phi[l], *dp = dphi[l]; in THIJacobianLocal_2D() local 1011 …l * 2 + 0] += dp[0] * jw * eta * 4. * dpl[0] + dp[1] * jw * eta * dpl[1] + pp * jw * (beta2 / h) *… in THIJacobianLocal_2D() 1014 …l * 2 + 1] += dp[1] * jw * eta * 4. * dpl[1] + dp[0] * jw * eta * dpl[0] + pp * jw * (beta2 / h) *… in THIJacobianLocal_2D() 1016 …* (4. * du[0] + 2. * dv[1]) + dp[1] * jw * deta * dgdu * (du[1] + dv[0]) + pp * jw * (dbeta2 / h) … in THIJacobianLocal_2D() [all …]
|
| H A D | ex77.c | 149 PetscReal cofu_x[9 /* Ncomp*dim */], detu_x, pp, pm, p = PetscRealPart(u[Ncomp]); in g3_uu_3d() local 155 pp = p / detu_x + kappa; in g3_uu_3d() 167 …g3[((compI * Ncomp + compJ) * dim + d1) * dim + d2] = g * G * mu + pp * cofu_x[compI * dim + d1] *… in g3_uu_3d()
|
| /petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/ |
| H A D | jacobian.m | 7 % SIAM J. Optimization, Vol. 20 (1), pp.172-191, 2009.
|
| H A D | dfoxs.m | 8 % SIAM J. Optimization, Vol. 20 (1), pp.172-191, 2009.
|
| H A D | perf_profile.m | 6 % SIAM J. Optimization, Vol. 20 (1), pp.172-191, 2009.
|
| H A D | dfovec.m | 7 % SIAM J. Optimization, Vol. 20 (1), pp.172-191, 2009.
|
| /petsc/src/mat/impls/maij/ |
| H A D | maij.c | 586 Mat_SeqMAIJ *pp = (Mat_SeqMAIJ *)PP->data; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() local 587 Mat P = pp->AIJ; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() 593 const PetscInt am = A->rmap->N, cn = C->cmap->N, cm = C->rmap->N, ppdof = pp->dof; in MatPtAPNumeric_SeqAIJ_SeqMAIJ() 671 Mat_SeqMAIJ *pp = (Mat_SeqMAIJ *)PP->data; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() local 672 Mat P = pp->AIJ; in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() 676 …const PetscInt an = A->cmap->N, am = A->rmap->N, pn = P->cmap->N, pm = P->rmap->N, ppdof = pp-… in MatPtAPSymbolic_SeqAIJ_SeqMAIJ() 777 PetscCall(MatSetBlockSize(C, pp->dof)); in MatPtAPSymbolic_SeqAIJ_SeqMAIJ()
|
| /petsc/doc/miscellaneous/ |
| H A D | prizes.md | 9 Software Tools in Scientific Computing, pp. 163-202 (1997)
|
| /petsc/src/dm/impls/plex/transform/interface/ |
| H A D | plextransform.c | 1460 …PetscInt pp = p; /* Parent point: Point in the original… in DMPlexTransformGetCone_Internal() local 1478 ppp = pp; in DMPlexTransformGetCone_Internal() 1479 pp = pcone[pcp]; in DMPlexTransformGetCone_Internal() 1480 PetscCall(DMPlexGetCellType(dm, pp, &pct)); in DMPlexTransformGetCone_Internal() 1483 PetscCall(DMPlexGetOrientedCone(dm, pp, &pcone, NULL)); in DMPlexTransformGetCone_Internal() 1488 if (lc) PetscCall(DMPlexRestoreOrientedCone(dm, pp, &pcone, NULL)); in DMPlexTransformGetCone_Internal() 1491 PetscCall(DMPlexTransformGetSubcellOrientation(tr, pct, pp, fn ? po : o, ft, pr, fo, &pr, &fo)); in DMPlexTransformGetCone_Internal() 1492 PetscCall(DMPlexTransformGetTargetPoint(tr, pct, ft, pp, pr, &coneNew[c])); in DMPlexTransformGetCone_Internal() 1868 …formLabelProducedPoint_Private(DMPlexTransform tr, DMLabel label, PetscInt pp, DMPolytopeType pct,… in DMPlexTransformLabelProducedPoint_Private() argument
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexglvis.c | 1002 const PetscInt pp = closure[2 * i]; in DMPlexView_GLVis_ASCII() local 1004 if (pp >= vStart && pp < vEnd) PetscCall(PetscBTSet(vown, pp - vStart)); in DMPlexView_GLVis_ASCII()
|
| H A D | plexgeometry.c | 611 PetscReal pp[3]; in DMPlexLocatePoint_Hex_3D_Linear_Internal() local 623 pp[0] = PetscRealPart(coords[faces[f * 4 + 0] * 3 + 0] - point[0]); in DMPlexLocatePoint_Hex_3D_Linear_Internal() 624 pp[1] = PetscRealPart(coords[faces[f * 4 + 0] * 3 + 1] - point[1]); in DMPlexLocatePoint_Hex_3D_Linear_Internal() 625 pp[2] = PetscRealPart(coords[faces[f * 4 + 0] * 3 + 2] - point[2]); in DMPlexLocatePoint_Hex_3D_Linear_Internal() 626 dot = normal[0] * pp[0] + normal[1] * pp[1] + normal[2] * pp[2]; in DMPlexLocatePoint_Hex_3D_Linear_Internal()
|
| /petsc/src/ksp/pc/impls/bddc/ |
| H A D | bddcprivate.c | 2654 PetscInt s, pp, p = closure[k], off, dof, cdof; in PCBDDCDetectDisconnectedComponents() local 2665 PetscCall(DMPlexGetTreeParent(dm, p, &pp, NULL)); in PCBDDCDetectDisconnectedComponents() 2666 if (pp != p) { in PCBDDCDetectDisconnectedComponents() 2667 PetscCall(PetscSectionGetConstraintDof(subSection, pp, &cdof)); in PCBDDCDetectDisconnectedComponents() 2668 PetscCall(PetscSectionGetOffset(subSection, pp, &off)); in PCBDDCDetectDisconnectedComponents() 2669 PetscCall(PetscSectionGetDof(subSection, pp, &dof)); in PCBDDCDetectDisconnectedComponents() 2672 if (PetscBTLookup(btvc, pp - pStart)) pids[cump++] = off + s; /* subdomain corner */ in PCBDDCDetectDisconnectedComponents() 2815 PetscInt p, *pp; in PCBDDCBenignDetectSaddlePoint() local 2818 PetscCall(PetscMalloc1(pcbddc->n_ISForDofsLocal, &pp)); in PCBDDCBenignDetectSaddlePoint() 2821 …PetscCall(PetscOptionsIntArray("-pc_bddc_pressure_field", "Field id for pressures", NULL, pp, &n, … in PCBDDCBenignDetectSaddlePoint() [all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 823 PetscReal pp; in PetscDTJacobiEval_Internal() local 826 pp = (cnm1 + cnm1x * x) * pm1 - cnm2 * pm2; in PetscDTJacobiEval_Internal() 828 pm1 = pp; in PetscDTJacobiEval_Internal() 830 p[l] = pp; in PetscDTJacobiEval_Internal()
|