Lines Matching refs:ct
464 DMPolytopeType ct; in DMPlexTransformCreateOffset_Internal() local
474 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformCreateOffset_Internal()
481 if (DMPolytopeTypeGetDim(ct) < 0 || DMPolytopeTypeGetDim(ctNew) < 0) { in DMPlexTransformCreateOffset_Internal()
523 const DMPolytopeType ct = (DMPolytopeType)c; in DMPlexTransformCreateOffset_Internal() local
530 …if (DMPolytopeTypeGetDim(ct) < 0 || ct == DM_POLYTOPE_UNKNOWN_CELL || ct == DM_POLYTOPE_UNKNOWN_FA… in DMPlexTransformCreateOffset_Internal()
531 off[ct * DM_NUM_POLYTOPES + ctNew] = -1; in DMPlexTransformCreateOffset_Internal()
534 off[ct * DM_NUM_POLYTOPES + ctNew] = 0; in DMPlexTransformCreateOffset_Internal()
540 if (ict == ct) { in DMPlexTransformCreateOffset_Internal()
543 if (n == Nct) off[ct * DM_NUM_POLYTOPES + ctNew] = -1; in DMPlexTransformCreateOffset_Internal()
547 …if (rct[n] == ctNew) off[ct * DM_NUM_POLYTOPES + ctNew] += (ctStart[ictn] - ctStart[ict]) * rsize[… in DMPlexTransformCreateOffset_Internal()
610 DMPolytopeType ct; in DMPlexTransformSetUp() local
615 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetUp()
616 …PetscCheck(ct != DM_POLYTOPE_UNKNOWN && ct != DM_POLYTOPE_UNKNOWN_CELL && ct != DM_POLYTOPE_UNKNOW… in DMPlexTransformSetUp()
617 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &cone, &ornt)); in DMPlexTransformSetUp()
628 DMPolytopeType ct; in DMPlexTransformSetUp() local
633 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetUp()
634 …PetscCheck(ct != DM_POLYTOPE_UNKNOWN && ct != DM_POLYTOPE_UNKNOWN_CELL && ct != DM_POLYTOPE_UNKNOW… in DMPlexTransformSetUp()
635 ++ctC[ct]; in DMPlexTransformSetUp()
636 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &cone, &ornt)); in DMPlexTransformSetUp()
812 static PetscErrorCode DMPlexTransformGetCoordinateFE(DMPlexTransform tr, DMPolytopeType ct, PetscFE… in DMPlexTransformGetCoordinateFE() argument
815 if (!tr->coordFE[ct]) { in DMPlexTransformGetCoordinateFE()
818 dim = DMPolytopeTypeGetDim(ct); in DMPlexTransformGetCoordinateFE()
820 …(PetscFECreateLagrangeByCell(PETSC_COMM_SELF, dim, cdim, ct, 1, PETSC_DETERMINE, &tr->coordFE[ct])… in DMPlexTransformGetCoordinateFE()
830 PetscCall(DMPlexTransformGetCellVertices(tr, ct, &Nq, &Xq)); in DMPlexTransformGetCoordinateFE()
837 PetscCall(PetscFESetQuadrature(tr->coordFE[ct], quad)); in DMPlexTransformGetCoordinateFE()
839 PetscCall(PetscFEGetDualSpace(tr->coordFE[ct], &dsp)); in DMPlexTransformGetCoordinateFE()
841 PetscCall(PetscFEGeomCreate(quad, 1, cdim, PETSC_FEGEOM_BASIC, &tr->refGeom[ct])); in DMPlexTransformGetCoordinateFE()
842 cg = tr->refGeom[ct]; in DMPlexTransformGetCoordinateFE()
847 *fe = tr->coordFE[ct]; in DMPlexTransformGetCoordinateFE()
1000 PetscErrorCode DMPlexTransformGetTargetPoint(DMPlexTransform tr, DMPolytopeType ct, DMPolytopeType … in DMPlexTransformGetTargetPoint() argument
1006 …PetscInt ctS = tr->ctStart[ct], ctE = tr->ctStart[tr->ctOrderOld[tr->ctOrderInvOld[ct] + 1]… in DMPlexTransformGetTargetPoint()
1011 …scInt_FMT " is not a %s [%" PetscInt_FMT ", %" PetscInt_FMT ")", p, DMPolytopeTypes[ct], ctS, ctE); in DMPlexTransformGetTargetPoint()
1012 PetscCall(DMPlexTransformCellTransform(tr, ct, p, &rt, &Nct, &rct, &rsize, &cone, &ornt)); in DMPlexTransformGetTargetPoint()
1016 … %s point %" PetscInt_FMT " does not have refine type %" PetscInt_FMT, DMPolytopeTypes[ct], p, rt); in DMPlexTransformGetTargetPoint()
1018 cind = ct; in DMPlexTransformGetTargetPoint()
1022 …PetscInt_FMT " does not produce type %s for transform %s", DMPolytopeTypes[ct], rt, p, DMPolytopeT… in DMPlexTransformGetTargetPoint()
1027 …) for subcell type %s in cell type %s", r, rsize[n], DMPolytopeTypes[rct[n]], DMPolytopeTypes[ct]); in DMPlexTransformGetTargetPoint()
1057 PetscErrorCode DMPlexTransformGetSourcePoint(DMPlexTransform tr, PetscInt pNew, DMPolytopeType *ct,… in DMPlexTransformGetSourcePoint() argument
1153 if (ct) *ct = (DMPolytopeType)ctO; in DMPlexTransformGetSourcePoint()
1400 DMPolytopeType ct; in DMPlexTransformSetConeSizes() local
1405 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetConeSizes()
1406 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformSetConeSizes()
1409 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &pNew)); in DMPlexTransformSetConeSizes()
1440 …etCone_Internal(DMPlexTransform tr, PetscInt p, PetscInt o, DMPolytopeType ct, DMPolytopeType ctNe… in DMPlexTransformGetCone_Internal() argument
1454 if (DMPolytopeTypeGetDim(ctNew) == dim && DMPolytopeTypeGetDim(ct) == dim - 1) { in DMPlexTransformGetCone_Internal()
1455 PetscCall(DMPlexTransformGetSubcellOrientation(tr, ct, p, o, ctNew, cr, co, &nr, &no)); in DMPlexTransformGetCone_Internal()
1462 …DMPolytopeType pct = ct; /* Parent type: Cell type for parent o… in DMPlexTransformGetCone_Internal()
1529 DMPolytopeType ct; in DMPlexTransformSetCones() local
1546 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetCones()
1547 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformSetCones()
1552 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &pNew)); in DMPlexTransformSetCones()
1553 …PetscCall(DMPlexTransformGetCone_Internal(tr, p, 0, ct, ctNew, rcone, &coff, rornt, &ooff, coneNew… in DMPlexTransformSetCones()
1572 DMPolytopeType ct, qct; in DMPlexTransformGetConeOriented() local
1585 PetscCall(DMPlexTransformGetSourcePoint(tr, q, &ct, &qct, &p, &r)); in DMPlexTransformGetConeOriented()
1586 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformGetConeOriented()
1604 …PetscCall(DMPlexTransformGetCone_Internal(tr, p, po, ct, qct, rcone, &coff, rornt, &ooff, qcone, q… in DMPlexTransformGetConeOriented()
1613 DMPolytopeType ct, qct; in DMPlexTransformGetCone() local
1626 PetscCall(DMPlexTransformGetSourcePoint(tr, q, &ct, &qct, &p, &r)); in DMPlexTransformGetCone()
1627 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformGetCone()
1645 …PetscCall(DMPlexTransformGetCone_Internal(tr, p, 0, ct, qct, rcone, &coff, rornt, &ooff, qcone, qo… in DMPlexTransformGetCone()
1672 const DMPolytopeType ct = (DMPolytopeType)ict; in DMPlexTransformCreateCellVertices_Internal() local
1685 …if (DMPolytopeTypeGetDim(ct) <= 0 || ct == DM_POLYTOPE_UNKNOWN_CELL || ct == DM_POLYTOPE_UNKNOWN_F… in DMPlexTransformCreateCellVertices_Internal()
1686 PetscCall(DMPlexCreateReferenceCell(PETSC_COMM_SELF, ct, &refdm)); in DMPlexTransformCreateCellVertices_Internal()
1696 tr->trNv[ct] = vEnd - vStart; in DMPlexTransformCreateCellVertices_Internal()
1699 …ct] * trdim == Nc, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Cell type %s, transformed coordinate size %" … in DMPlexTransformCreateCellVertices_Internal()
1700 PetscCall(PetscCalloc1(Nc, &tr->trVerts[ct])); in DMPlexTransformCreateCellVertices_Internal()
1702 PetscCall(PetscArraycpy(tr->trVerts[ct], coords, Nc)); in DMPlexTransformCreateCellVertices_Internal()
1705 PetscCall(PetscCalloc1(DM_NUM_POLYTOPES, &tr->trSubVerts[ct])); in DMPlexTransformCreateCellVertices_Internal()
1706 PetscCall(DMPlexTransformCellTransform(reftr, ct, 0, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateCellVertices_Internal()
1710 PetscCall(PetscCalloc1(rsize[n], &tr->trSubVerts[ct][rct[n]])); in DMPlexTransformCreateCellVertices_Internal()
1715 … PetscCall(PetscCalloc1(DMPolytopeTypeGetNumVertices(rct[n]), &tr->trSubVerts[ct][rct[n]][r])); in DMPlexTransformCreateCellVertices_Internal()
1716 PetscCall(DMPlexTransformGetTargetPoint(reftr, ct, rct[n], 0, r, &pNew)); in DMPlexTransformCreateCellVertices_Internal()
1721 if ((sv >= vStart) && (sv < vEnd)) tr->trSubVerts[ct][rct[n]][r][Nv++] = sv - vStart; in DMPlexTransformCreateCellVertices_Internal()
1724 …ll %s", Nv, DMPolytopeTypeGetNumVertices(rct[n]), DMPolytopeTypes[rct[n]], r, DMPolytopeTypes[ct]); in DMPlexTransformCreateCellVertices_Internal()
1732 …tscPrintf(PETSC_COMM_SELF, "%s: %" PetscInt_FMT " vertices\n", DMPolytopeTypes[ct], tr->trNv[ct])); in DMPlexTransformCreateCellVertices_Internal()
1733 for (v = 0; v < tr->trNv[ct]; ++v) { in DMPlexTransformCreateCellVertices_Internal()
1735 …++d) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%g ", (double)PetscRealPart(tr->trVerts[ct][off++]))); in DMPlexTransformCreateCellVertices_Internal()
1739 … PetscCall(DMPlexTransformCellTransform(reftr, ct, 0, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateCellVertices_Internal()
1742 …s: %s subvertices %" PetscInt_FMT "\n", DMPolytopeTypes[ct], DMPolytopeTypes[rct[n]], tr->trNv[ct]… in DMPlexTransformCreateCellVertices_Internal()
1745 …v) PetscCall(PetscPrintf(PETSC_COMM_SELF, "%" PetscInt_FMT " ", tr->trSubVerts[ct][rct[n]][r][v])); in DMPlexTransformCreateCellVertices_Internal()
1772 PetscErrorCode DMPlexTransformGetCellVertices(DMPlexTransform tr, DMPolytopeType ct, PetscInt *Nv, … in DMPlexTransformGetCellVertices() argument
1776 if (Nv) *Nv = tr->trNv[ct]; in DMPlexTransformGetCellVertices()
1777 if (trVerts) *trVerts = tr->trVerts[ct]; in DMPlexTransformGetCellVertices()
1797 PetscErrorCode DMPlexTransformGetSubcellVertices(DMPlexTransform tr, DMPolytopeType ct, DMPolytopeT… in DMPlexTransformGetSubcellVertices() argument
1801 …k(tr->trSubVerts[ct][rct], PetscObjectComm((PetscObject)tr), PETSC_ERR_ARG_WRONG, "Cell type %s do… in DMPlexTransformGetSubcellVertices()
1802 if (subVerts) *subVerts = tr->trSubVerts[ct][rct][r]; in DMPlexTransformGetSubcellVertices()
1807 …Barycenter_Internal(DMPlexTransform tr, DMPolytopeType pct, DMPolytopeType ct, PetscInt p, PetscIn… in DMPlexTransformMapCoordinatesBarycenter_Internal() argument
1812 …PetscCheck(ct == DM_POLYTOPE_POINT, PETSC_COMM_SELF, PETSC_ERR_SUP, "Not for refined point type %s… in DMPlexTransformMapCoordinatesBarycenter_Internal()
1842 …sformMapCoordinates(DMPlexTransform tr, DMPolytopeType pct, DMPolytopeType ct, PetscInt p, PetscIn… in DMPlexTransformMapCoordinates() argument
1845 if (Nv) PetscUseTypeMethod(tr, mapcoordinates, pct, ct, p, r, Nv, dE, in, out); in DMPlexTransformMapCoordinates()
1868 …DMLabel label, PetscInt pp, DMPolytopeType pct, PetscInt p, DMPolytopeType ct, PetscInt r, PetscIn… in DMPlexTransformLabelProducedPoint_Private() argument
1871 if (tr->labelMatchStrata && pct != ct) PetscFunctionReturn(PETSC_SUCCESS); in DMPlexTransformLabelProducedPoint_Private()
1903 DMPolytopeType ct; in RefineLabel_Internal() local
1908 PetscCall(DMPlexGetCellType(dm, point, &ct)); in RefineLabel_Internal()
1909 … PetscCall(DMPlexTransformCellTransform(tr, ct, point, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in RefineLabel_Internal()
1912 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], point, r, &pNew)); in RefineLabel_Internal()
1913 …PetscCall(DMPlexTransformLabelProducedPoint_Private(tr, labelNew, point, ct, pNew, rct[n], r, valu… in RefineLabel_Internal()
2022 DMPolytopeType ct; in DMPlexTransformCreateSF() local
2027 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformCreateSF()
2028 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateSF()
2038 DMPolytopeType ct; in DMPlexTransformCreateSF() local
2043 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformCreateSF()
2044 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateSF()
2057 DMPolytopeType ct; in DMPlexTransformCreateSF() local
2064 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformCreateSF()
2065 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateSF()
2068 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &pNew)); in DMPlexTransformCreateSF()
2080 DMPolytopeType ct; in DMPlexTransformCreateSF() local
2086 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformCreateSF()
2087 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformCreateSF()
2090 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &pNew)); in DMPlexTransformCreateSF()
2148 static PetscErrorCode DMPlexTransformMapLocalizedCoordinates(DMPlexTransform tr, DMPolytopeType ct,… in DMPlexTransformMapLocalizedCoordinates() argument
2154 PetscCall(DMPlexTransformGetCoordinateFE(tr, ct, &fe)); in DMPlexTransformMapLocalizedCoordinates()
2155 PetscCall(DMPlexTransformGetSubcellVertices(tr, ct, rct, r, &subcellV)); in DMPlexTransformMapLocalizedCoordinates()
2157 …Vertices(rct); ++v) PetscCall(PetscFEInterpolate_Static(fe, x, tr->refGeom[ct], subcellV[v], &xr[v… in DMPlexTransformMapLocalizedCoordinates()
2230 DMPolytopeType ct; in DMPlexTransformSetCoordinates() local
2235 PetscCall(DMPlexGetCellType(dm, c, &ct)); in DMPlexTransformSetCoordinates()
2236 dim = DMPolytopeTypeGetDim(ct); in DMPlexTransformSetCoordinates()
2237 … PetscCall(DMPlexTransformCellTransform(tr, ct, c, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformSetCoordinates()
2245 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], c, r, &cNew)); in DMPlexTransformSetCoordinates()
2282 DMPolytopeType ct; in DMPlexTransformSetCoordinates() local
2288 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetCoordinates()
2289 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformSetCoordinates()
2307 if (ct != DM_POLYTOPE_POINT) { in DMPlexTransformSetCoordinates()
2321 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &vNew)); in DMPlexTransformSetCoordinates()
2323 PetscCall(DMPlexTransformMapCoordinates(tr, ct, rct[n], p, r, Nv, dEo, icoords, vcoords)); in DMPlexTransformSetCoordinates()
2356 DMPolytopeType ct; in DMPlexTransformSetCoordinates() local
2361 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetCoordinates()
2362 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexTransformSetCoordinates()
2371 if (DMPolytopeTypeGetDim(ct) != DMPolytopeTypeGetDim(rct[n])) continue; in DMPlexTransformSetCoordinates()
2378 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, r, &pNew)); in DMPlexTransformSetCoordinates()
2380 … PetscCall(DMPlexTransformMapLocalizedCoordinates(tr, ct, rct[n], r, pcoords, &coordsNew[offNew])); in DMPlexTransformSetCoordinates()