Home
last modified time | relevance | path

Searched refs:ww (Results 1 – 15 of 15) sorted by relevance

/petsc/src/vec/vec/impls/seq/
H A Dbvec2.c14 PetscScalar *ww, *xx, *yy; /* cannot make xx or yy const since might be ww */ in VecPointwiseApply_Seq() local
19 PetscCall(VecGetArray(win, &ww)); in VecPointwiseApply_Seq()
20 for (PetscInt i = 0; i < n; ++i) ww[i] = func(xx[i], yy[i]); in VecPointwiseApply_Seq()
23 PetscCall(VecRestoreArray(win, &ww)); in VecPointwiseApply_Seq()
75 PetscScalar *ww, *xx, *yy; /* cannot make xx or yy const since might be ww */ in VecPointwiseMult_Seq() local
80 PetscCall(VecGetArray(win, &ww)); in VecPointwiseMult_Seq()
81 if (ww == xx) { in VecPointwiseMult_Seq()
82 for (i = 0; i < n; i++) ww[i] *= yy[i]; in VecPointwiseMult_Seq()
83 } else if (ww == yy) { in VecPointwiseMult_Seq()
84 for (i = 0; i < n; i++) ww[i] *= xx[i]; in VecPointwiseMult_Seq()
[all …]
H A Ddvec2.c801 PetscScalar *ww; in VecWAXPY_Seq() local
806 PetscCall(VecGetArray(win, &ww)); in VecWAXPY_Seq()
810 for (PetscInt i = 0; i < n; i++) ww[i] = yy[i] + xx[i]; in VecWAXPY_Seq()
813 for (PetscInt i = 0; i < n; i++) ww[i] = yy[i] - xx[i]; in VecWAXPY_Seq()
815 PetscCall(PetscArraycpy(ww, yy, n)); in VecWAXPY_Seq()
819 fortranwaxpy_(&n, &alpha, xx, yy, ww); in VecWAXPY_Seq()
821 for (PetscInt i = 0; i < n; i++) ww[i] = yy[i] + alpha * xx[i]; in VecWAXPY_Seq()
826 PetscCall(VecRestoreArray(win, &ww)); in VecWAXPY_Seq()
/petsc/src/mat/impls/mffd/
H A Dmffd.c388 PetscScalar h, *aa, *ww, v; in MatGetDiagonal_MFFD() local
406 PetscCall(VecGetArray(w, &ww)); in MatGetDiagonal_MFFD()
407 h = ww[i - rstart]; in MatGetDiagonal_MFFD()
413 ww[i - rstart] += h; in MatGetDiagonal_MFFD()
414 PetscCall(VecRestoreArray(w, &ww)); in MatGetDiagonal_MFFD()
418 PetscCall(VecGetArray(w, &ww)); in MatGetDiagonal_MFFD()
419 ww[i - rstart] -= h; in MatGetDiagonal_MFFD()
420 PetscCall(VecRestoreArray(w, &ww)); in MatGetDiagonal_MFFD()
/petsc/src/sys/classes/draw/impls/x/
H A Dxinit.c274 unsigned int ww = 0, hh = 0, dummy; in PetscDrawXiGetGeometry() local
278 XGetGeometry(XiWin->disp, XiWin->win, &parent, &xx, &yy, &ww, &hh, &dummy, &dummy); in PetscDrawXiGetGeometry()
286 XGetGeometry(XiWin->disp, XiWin->drw, &root, &xx, &yy, &ww, &hh, &dummy, &dummy); in PetscDrawXiGetGeometry()
290 if (w) *w = (int)ww; in PetscDrawXiGetGeometry()
/petsc/src/ts/tests/
H A Dex30.c333 …dim + 0] : 1, w = fact * wp[p] * ctx->n_0 * ctx->masses[ctx->species_offset[grid]], ww = w * N_inv; in gridToParticles_private() local
335 moments[3] -= ww * PetscLogReal(ww); in gridToParticles_private()
336 …PetscCheck(ww < 1 - PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "ww (%g) > 1", (doubl… in gridToParticles_private()
429 …im + 0] : 1, w = fact * wp[pp] * ctx->n_0 * ctx->masses[ctx->species_offset[grid]], ww = w * N_inv; in PostStep() local
431 moments[3] -= ww * PetscLogReal(ww); in PostStep()
675 …w = fact * ctx->n_0 * ctx->masses[ctx->species_offset[grid]] * wp_t[grid][tid][pp], ww = w * N_inv; in go() local
677 moments_0[3] -= ww * PetscLogReal(ww); in go()
678 …PetscCheck(ww < 1 - PETSC_MACHINE_EPSILON, PETSC_COMM_WORLD, PETSC_ERR_ARG_OUTOFRANGE, "ww (%g) > … in go()
812 …im + 0] : 1, w = fact * wp[pp] * ctx->n_0 * ctx->masses[ctx->species_offset[grid]], ww = w * N_inv; in go() local
814 moments_1a[3] -= ww * PetscLogReal(ww); in go()
[all …]
/petsc/src/ts/tutorials/
H A Dex10.c514 static const PetscReal ww[1] = {1.}; in RDGetQuadrature() local
522 refweight = ww; in RDGetQuadrature()
535 static const PetscReal ww[2] = {0.5, 0.5}; in RDGetQuadrature() local
537 refweight = ww; in RDGetQuadrature()
552 static const PetscReal ww[3] = {5. / 18, 8. / 18, 5. / 18}; in RDGetQuadrature() local
554 refweight = ww; in RDGetQuadrature()
571 …static const PetscReal ww[] = {0.17392742256872692, 0.3260725774312731, 0.3260725774312731, 0.1739… in RDGetQuadrature() local
574 refweight = ww; in RDGetQuadrature()
587 static const PetscReal ww[2] = {0.5, 0.5}; in RDGetQuadrature() local
589 refweight = ww; in RDGetQuadrature()
[all …]
/petsc/src/dm/impls/da/
H A Dda2.c155 PetscScalar h, *aa, *ww, v; in DMDAGetDiagonal_MFFD() local
165 PetscCall(VecGetArray(U, &ww)); in DMDAGetDiagonal_MFFD()
169 h = ww[gI]; in DMDAGetDiagonal_MFFD()
175 ww[gI] += h; in DMDAGetDiagonal_MFFD()
178 ww[gI] -= h; in DMDAGetDiagonal_MFFD()
181 PetscCall(VecRestoreArray(U, &ww)); in DMDAGetDiagonal_MFFD()
/petsc/src/dm/dt/interface/
H A Ddtaltv.c389 const PetscReal *ww = NULL; in PetscDTAltVPullback() local
400 ww = walloc; in PetscDTAltVPullback()
402 ww = w; in PetscDTAltVPullback()
415 prod = isOdd ? -ww[i] : ww[i]; in PetscDTAltVPullback()
H A Ddt.c1915 PetscReal *x, *w, *xw, *ww; in PetscDTGaussTensorQuadrature() local
1934 PetscCall(PetscMalloc1(npoints, &ww)); in PetscDTGaussTensorQuadrature()
1935 PetscCall(PetscDTGaussQuadrature(npoints, a, b, x, ww)); in PetscDTGaussTensorQuadrature()
1937 for (PetscInt c = 0; c < Nc; ++c) w[i * Nc + c] = ww[i]; in PetscDTGaussTensorQuadrature()
1938 PetscCall(PetscFree(ww)); in PetscDTGaussTensorQuadrature()
1942 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
1943 PetscCall(PetscDTGaussQuadrature(npoints, a, b, xw, ww)); in PetscDTGaussTensorQuadrature()
1948 for (PetscInt c = 0; c < Nc; ++c) w[(i * npoints + j) * Nc + c] = ww[i] * ww[j]; in PetscDTGaussTensorQuadrature()
1951 PetscCall(PetscFree2(xw, ww)); in PetscDTGaussTensorQuadrature()
1955 PetscCall(PetscMalloc2(npoints, &xw, npoints, &ww)); in PetscDTGaussTensorQuadrature()
[all …]
/petsc/src/mat/impls/aij/seq/aijperm/
H A Daijperm.c458 static PetscErrorCode MatMultAdd_SeqAIJPERM(Mat A, Vec xx, Vec ww, Vec yy) in MatMultAdd_SeqAIJPERM() argument
494 PetscCall(VecGetArrayPair(yy, ww, &y, &w)); in MatMultAdd_SeqAIJPERM()
591 PetscCall(VecRestoreArrayPair(yy, ww, &y, &w)); in MatMultAdd_SeqAIJPERM()
/petsc/src/vec/vec/impls/nest/
H A Dvecnest.c607 Vec *ww; in VecCreateLocalVector_Nest() local
619 PetscCall(PetscMalloc2(bv->nb, &ww, bv->nb, &wis)); in VecCreateLocalVector_Nest()
620 for (i = 0; i < bv->nb; i++) PetscCall(VecCreateLocalVector(bv->v[i], &ww[i])); in VecCreateLocalVector_Nest()
625 PetscCall(ISOnComm(bv->is[i], PetscObjectComm((PetscObject)ww[i]), PETSC_COPY_VALUES, &wis[i])); in VecCreateLocalVector_Nest()
628 PetscCall(VecCreateNest(PETSC_COMM_SELF, bv->nb, wis, ww, w)); in VecCreateLocalVector_Nest()
630 PetscCall(VecDestroy(&ww[i])); in VecCreateLocalVector_Nest()
633 PetscCall(PetscFree2(ww, wis)); in VecCreateLocalVector_Nest()
/petsc/src/dm/dt/tests/
H A Dex7.c12 const PetscReal *ww = NULL; in CheckPullback() local
22 ww = walloc; in CheckPullback()
24 ww = w; in CheckPullback()
75 PetscCall(PetscDTAltVApply(M, k, ww, Lx, &wLx)); in CheckPullback()
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c211 ….x, *yy = (PetscReal *)ctx->SData_d.y, *zz = (PetscReal *)ctx->SData_d.z, *ww = (PetscReal *)ctx->… in LandauFormJacobian_Internal() local
341 …etscReal vj[3] = {xx[jpidx_glb], yy[jpidx_glb], zz ? zz[jpidx_glb] : 0}, wj = ww[jpidx_glb]; in LandauFormJacobian_Internal()
353 const PetscReal wi = ww[ipidx], x = xx[ipidx], y = yy[ipidx]; in LandauFormJacobian_Internal()
431 PetscReal wj = ww[jpidx_glb]; in LandauFormJacobian_Internal()
1699 PetscReal *invJe, *ww, *xx, *yy, *zz = NULL, *invJ_a; in CreateStaticData() local
1715 PetscCall(PetscMalloc4(nip_glb, &ww, nip_glb, &xx, nip_glb, &yy, nip_glb * dim * dim, &invJ_a)); in CreateStaticData()
1763 ww[gidx] = detJj[qj] * quadWeights[qj]; in CreateStaticData()
1764 if (dim == 2) ww[gidx] *= vj[qj * dim + 0]; /* cylindrical coordinate, w/o 2pi */ in CreateStaticData()
1816 …et, nu_alpha, nu_beta, invMass, (PetscReal *)ctx->lambdas, invJ_a, xx, yy, zz, ww, &ctx->SData_d)); in CreateStaticData()
1818 PetscCall(PetscFree4(ww, xx, yy, invJ_a)); in CreateStaticData()
[all …]
/petsc/share/petsc/datafiles/meshes/
H A Dmesh-3d-box-innersphere_bin.msh17 …`�?N����Y�yO�?+�>'�?0�4?����O����0d��?�? �g�If������[��P�������r��o��ww��?Z�9�����Q���������…
H A Dtestcase3D.cas2203 (ww-stress/interpolate #t)
2785 (ww-stress/extrapolate/cells #t)
3549 (ww-stress/patch 1)
8211 (les-ww-wall? #f)