| /petsc/src/dm/impls/plex/transform/impls/refine/bl/ |
| H A D | plexrefbl.c | 18 DMPolytopeType ct; in DMPlexTransformSetUp_BL() local 46 PetscCall(DMPlexGetCellType(dm, point, &ct)); in DMPlexTransformSetUp_BL() 47 val = (PetscInt)ct; in DMPlexTransformSetUp_BL() 48 if (ct > DM_POLYTOPE_POINT_PRISM_TENSOR) ++val; in DMPlexTransformSetUp_BL() 49 if (ct > DM_POLYTOPE_SEG_PRISM_TENSOR) ++val; in DMPlexTransformSetUp_BL() 50 if (ct > DM_POLYTOPE_TRI_PRISM_TENSOR) ++val; in DMPlexTransformSetUp_BL() 51 if (ct > DM_POLYTOPE_QUAD_PRISM_TENSOR) ++val; in DMPlexTransformSetUp_BL() 52 switch (ct) { in DMPlexTransformSetUp_BL() 73 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetUp_BL() 74 val = (PetscInt)ct; in DMPlexTransformSetUp_BL() [all …]
|
| /petsc/src/dm/impls/plex/transform/impls/extrude/ |
| H A D | plextrextrude.c | 100 DMPolytopeType ct; in DMPlexTransformExtrudeComputeExtrusionDim() local 112 PetscCall(DMPlexGetCellType(dm, points[p], &ct)); in DMPlexTransformExtrudeComputeExtrusionDim() 113 dimExtPoint = PetscMax(dimExtPoint, DMPolytopeTypeGetDim(ct)); in DMPlexTransformExtrudeComputeExtrusionDim() 147 const DMPolytopeType ct = DM_POLYTOPE_POINT; in DMPlexTransformExtrudeSetUp_Point() local 152 ex->Nt[ct] = 2; in DMPlexTransformExtrudeSetUp_Point() 155 …PetscCall(PetscMalloc4(ex->Nt[ct], &ex->target[ct], ex->Nt[ct], &ex->size[ct], Nc, &ex->cone[ct], … in DMPlexTransformExtrudeSetUp_Point() 156 ex->target[ct][0] = DM_POLYTOPE_POINT; in DMPlexTransformExtrudeSetUp_Point() 157 ex->target[ct][1] = ex->useTensor ? DM_POLYTOPE_POINT_PRISM_TENSOR : DM_POLYTOPE_SEGMENT; in DMPlexTransformExtrudeSetUp_Point() 158 ex->size[ct][0] = Np; in DMPlexTransformExtrudeSetUp_Point() 159 ex->size[ct][1] = Nl; in DMPlexTransformExtrudeSetUp_Point() [all …]
|
| H A D | plextrcohesive.c | 779 PetscInt ct, val; in DMPlexTransformSetUp_Cohesive() local 781 PetscCall(DMLabelGetValue(celltype, p, &ct)); in DMPlexTransformSetUp_Cohesive() 786 PetscCall(DMLabelSetValue(tr->trType, p, ct * 2 + 0)); in DMPlexTransformSetUp_Cohesive() 794 …ct == DM_POLYTOPE_POINT_PRISM_TENSOR || ct == DM_POLYTOPE_SEG_PRISM_TENSOR || ct == DM_POLYTOPE_TR… in DMPlexTransformSetUp_Cohesive() 795 else PetscCall(DMLabelSetValue(tr->trType, p, ct * 2 + 1)); in DMPlexTransformSetUp_Cohesive() 799 if (!fsplit) PetscCall(DMLabelSetValue(tr->trType, p, ct * 2 + 0)); in DMPlexTransformSetUp_Cohesive() 800 else PetscCall(DMLabelSetValue(tr->trType, p, (ct * 2 + 0) * 100 + fsplit)); in DMPlexTransformSetUp_Cohesive() 804 PetscCall(DMLabelSetValue(tr->trType, p, (ct * 2 + 1) * 100 + funsplit)); in DMPlexTransformSetUp_Cohesive() 833 PetscInt ct; in DMPlexTransformDestroy_Cohesive() local 837 …for (ct = 0; ct < DM_NUM_POLYTOPES * 2 * 100; ++ct) PetscCall(PetscFree4(ex->target[ct], ex->size[… in DMPlexTransformDestroy_Cohesive() [all …]
|
| /petsc/src/dm/impls/plex/tutorials/ |
| H A D | ex11.c | 69 DMPolytopeType ct; in CheckCellVertices() local 77 PetscCall(DMPlexGetCellType(dm, cell, &ct)); in CheckCellVertices() 86 …ct), comm, PETSC_ERR_ARG_WRONG, "Cell %" PetscInt_FMT " has %" PetscInt_FMT " vertices != %" Petsc… in CheckCellVertices() 87 arrVerts = DMPolytopeTypeGetVertexArrangement(ct, o); in CheckCellVertices() 136 DMPolytopeType ct; in GenerateArrangements() local 143 PetscCall(DMPlexGetCellType(dm, 0, &ct)); in GenerateArrangements() 144 No = DMPolytopeTypeGetNumArrangements(ct) / 2; in GenerateArrangements() 160 DMPolytopeType ct; in VerifyCayleyTable() local 169 PetscCall(DMPlexGetCellType(dm, 0, &ct)); in VerifyCayleyTable() 171 No = DMPolytopeTypeGetNumArrangements(ct) / 2; in VerifyCayleyTable() [all …]
|
| /petsc/src/dm/impls/plex/transform/interface/ |
| H A D | plextransform.c | 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() [all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | mattransposematmult.c | 17 PetscCall(VecDestroy(&atb->ct)); in MatProductCtxDestroy_SeqDense_MatTransMatMult() 50 PetscCall(MatCreateVecs(atb->mA, &atb->ct, &atb->bt)); in MatTMatTMultSymbolic_SeqAIJ_SeqDense() 69 Vec bt, ct; in MatTMatTMultNumeric_SeqAIJ_SeqDense() local 77 ct = atb->ct; in MatTMatTMultNumeric_SeqAIJ_SeqDense() 94 PetscCall(MatMultTranspose(atb->mA, bt, ct)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() 97 PetscCall(VecGetArrayRead(ct, &ctarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() 101 PetscCall(VecRestoreArrayRead(ct, &ctarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() 107 PetscCall(VecPlaceArray(ct, Barray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() 112 PetscCall(VecGetArrayWrite(ct, &ctarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() 116 PetscCall(VecRestoreArrayWrite(ct, &ctarray)); in MatTMatTMultNumeric_SeqAIJ_SeqDense() [all …]
|
| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | xxt.c | 465 PetscInt i, ct, id; in det_separators() local 542 for (ct = i = 0; i < n; i++) { in det_separators() 545 ct++; in det_separators() 564 (id < mask) ? (sum[2] = ct) : (sum[3] = ct); in det_separators() 565 (id < mask) ? (rsum[2] = ct) : (rsum[3] = ct); in det_separators() 575 for (ct = i = 0; i < n; i++) { in det_separators() 577 ct++; in det_separators() 586 if (ct > 1) PetscCall(PCTFS_ivec_sort(iptr, ct)); in det_separators() 588 lnsep[edge] = ct; in det_separators() 595 for (ct = i = 0; i < n; i++) { in det_separators() [all …]
|
| H A D | xyt.c | 575 PetscInt i, ct, id; in det_separators() local 645 for (ct = i = 0; i < n; i++) { in det_separators() 648 ct++; in det_separators() 658 (id < mask) ? (sum[2] = ct) : (sum[3] = ct); in det_separators() 661 ct = 0; in det_separators() 666 ct++; in det_separators() 677 ct++; in det_separators() 686 if (ct > 1) PetscCall(PCTFS_ivec_sort(iptr, ct)); in det_separators() 687 lnsep[edge] = ct; in det_separators() 699 for (ct = i = 0; i < n; i++) { in det_separators() [all …]
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpimattransposematmult.c | 17 PetscCall(VecDestroy(&atb->ct)); in MatProductCtxDestroy_MPIDense_MatTransMatMult() 47 PetscCall(MatCreateVecs(atb->mA, &atb->ct, &atb->bt)); in MatTransposeMatMultSymbolic_MPIAIJ_MPIDense() 62 Vec bt, ct; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() local 74 ct = atb->ct; in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() 86 PetscCall(MatMultTranspose(atb->mA, bt, ct)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() 91 PetscCall(VecGetArrayRead(ct, &ctarray)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense() 94 PetscCall(VecRestoreArrayRead(ct, &ctarray)); in MatTransposeMatMultNumeric_MPIAIJ_MPIDense()
|
| /petsc/src/dm/impls/plex/ |
| H A D | plexinterpolate.c | 46 PetscErrorCode DMPlexGetRawFaces_Internal(DM dm, DMPolytopeType ct, const PetscInt cone[], PetscInt… in DMPlexGetRawFaces_Internal() argument 72 switch (ct) { in DMPlexGetRawFaces_Internal() 475 …Object)dm), PETSC_ERR_ARG_OUTOFRANGE, "No face description for cell type %s", DMPolytopeTypes[ct]); in DMPlexGetRawFaces_Internal() 480 PetscErrorCode DMPlexRestoreRawFaces_Internal(DM dm, DMPolytopeType ct, const PetscInt cone[], Pets… in DMPlexRestoreRawFaces_Internal() argument 519 DMPolytopeType ct; in DMPlexInterpolateFaces_Internal() local 522 PetscCall(DMPlexGetCellType(dm, c, &ct)); in DMPlexInterpolateFaces_Internal() 527 …if (ct != DM_POLYTOPE_SEGMENT && ct != DM_POLYTOPE_POINT_PRISM_TENSOR) PetscCall(DMPlexGetRawFaces… in DMPlexInterpolateFaces_Internal() 539 DMPolytopeType ct; in DMPlexInterpolateFaces_Internal() local 542 PetscCall(DMPlexGetCellType(dm, c, &ct)); in DMPlexInterpolateFaces_Internal() 545 if (ct != DM_POLYTOPE_SEGMENT && ct != DM_POLYTOPE_POINT_PRISM_TENSOR) { in DMPlexInterpolateFaces_Internal() [all …]
|
| H A D | plexextrude.c | 97 DMPolytopeType ct; in DMPlexExtrude() local 102 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexExtrude() 103 PetscCall(DMPlexTransformCellTransform(tr, ct, p, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMPlexExtrude() 107 if (ct == rct[n] || pdim > dim) { in DMPlexExtrude() 108 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], p, 0, &pNew)); in DMPlexExtrude()
|
| H A D | plexfluent.c | 509 DMPolytopeType ct; in ReorderWedge() local 511 PetscCall(DMPlexGetCellType(dm, cone[c], &ct)); in ReorderWedge() 512 if (ct == DM_POLYTOPE_TRIANGLE) { in ReorderWedge() 531 DMPolytopeType ct; in ReorderWedge() local 535 PetscCall(DMPlexGetCellType(dm, cone[c], &ct)); in ReorderWedge() 536 if (ct != DM_POLYTOPE_QUADRILATERAL) continue; in ReorderWedge() 602 static PetscErrorCode ReorderCell(PetscViewer viewer, DM dm, PetscInt cell, DMPolytopeType ct) in ReorderCell() argument 605 switch (ct) { in ReorderCell() 620 … PetscCheck(0, PETSC_COMM_SELF, PETSC_ERR_SUP, "Celltype %s is unsupported", DMPolytopeTypes[ct]); in ReorderCell() 626 static PetscErrorCode GetNumCellFaces(int nd, PetscInt *numCellFaces, DMPolytopeType *ct) in GetNumCellFaces() argument [all …]
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 66 (*q)->ct = DM_POLYTOPE_UNKNOWN; in PetscQuadratureCreate() 93 DMPolytopeType ct; in PetscQuadratureDuplicate() local 101 PetscCall(PetscQuadratureGetCellType(q, &ct)); in PetscQuadratureDuplicate() 102 PetscCall(PetscQuadratureSetCellType(*r, ct)); in PetscQuadratureDuplicate() 156 PetscErrorCode PetscQuadratureGetCellType(PetscQuadrature q, DMPolytopeType *ct) in PetscQuadratureGetCellType() argument 160 PetscAssertPointer(ct, 2); in PetscQuadratureGetCellType() 161 *ct = q->ct; in PetscQuadratureGetCellType() 178 PetscErrorCode PetscQuadratureSetCellType(PetscQuadrature q, DMPolytopeType ct) in PetscQuadratureSetCellType() argument 182 q->ct = ct; in PetscQuadratureSetCellType() 355 …if (A->ct != B->ct || A->dim != B->dim || A->Nc != B->Nc || A->order != B->order || A->numPoints !… in PetscQuadratureEqual() [all …]
|
| /petsc/include/ |
| H A D | petscdm.h | 468 static inline PetscBool DMPolytopeTypeIsHybrid(DMPolytopeType ct) in DMPolytopeTypeIsHybrid() argument 470 switch (ct) { in DMPolytopeTypeIsHybrid() 481 static inline PetscInt DMPolytopeTypeGetDim(DMPolytopeType ct) in DMPolytopeTypeGetDim() argument 483 switch (ct) { in DMPolytopeTypeGetDim() 507 static inline PetscInt DMPolytopeTypeGetConeSize(DMPolytopeType ct) in DMPolytopeTypeGetConeSize() argument 509 switch (ct) { in DMPolytopeTypeGetConeSize() 539 static inline PetscInt DMPolytopeTypeGetNumVertices(DMPolytopeType ct) in DMPolytopeTypeGetNumVertices() argument 541 switch (ct) { in DMPolytopeTypeGetNumVertices() 576 static inline PetscInt DMPolytopeTypeGetNumArrangements(DMPolytopeType ct) in DMPolytopeTypeGetNumArrangements() argument 578 switch (ct) { in DMPolytopeTypeGetNumArrangements() [all …]
|
| /petsc/src/dm/impls/plex/tests/ |
| H A D | ex41.c | 63 DMPolytopeType ct; in ConstructRefineTree() local 69 PetscCall(DMPlexGetCellType(odm, c, &ct)); in ConstructRefineTree() 70 dim = DMPolytopeTypeGetDim(ct); in ConstructRefineTree() 71 PetscCall(DMPlexTransformCellTransform(tr, ct, c, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in ConstructRefineTree() 75 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], c, r, &pNew)); in ConstructRefineTree()
|
| H A D | ex72.c | 43 DMPolytopeType ct; in CreateDiscretization() local 49 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in CreateDiscretization() 50 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, 1, ct, NULL, -1, &fe)); in CreateDiscretization()
|
| H A D | ex33.c | 43 DMPolytopeType ct; in CreateDiscretization() local 49 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in CreateDiscretization() 50 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, 1, ct, NULL, PETSC_DETERMINE, &fe)); in CreateDiscretization()
|
| H A D | ex102.c | 63 DMPolytopeType ct; in SetupFE() local 70 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in SetupFE() 72 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, dim, ct, name ? prefix : NULL, -1, &fe)); in SetupFE()
|
| /petsc/src/dm/label/impls/ephemeral/ |
| H A D | dmlabeleph.c | 32 DMPolytopeType ct; in DMLabelEphemeralComputeStratumSize_Private() local 37 PetscCall(DMPlexGetCellType(dm, point, &ct)); in DMLabelEphemeralComputeStratumSize_Private() 38 … PetscCall(DMPlexTransformCellTransform(tr, ct, point, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMLabelEphemeralComputeStratumSize_Private() 69 DMPolytopeType ct; in DMLabelGetStratumIS_Ephemeral() local 74 PetscCall(DMPlexGetCellType(dm, point, &ct)); in DMLabelGetStratumIS_Ephemeral() 75 … PetscCall(DMPlexTransformCellTransform(tr, ct, point, NULL, &Nct, &rct, &rsize, &rcone, &rornt)); in DMLabelGetStratumIS_Ephemeral() 78 PetscCall(DMPlexTransformGetTargetPoint(tr, ct, rct[n], point, r, &pNew)); in DMLabelGetStratumIS_Ephemeral()
|
| /petsc/src/mat/tests/ |
| H A D | ex47.c | 12 PetscInt i, j, row, m, n, ncols1, ncols2, ct, m2, n2; in main() local 81 for (ct = 0; ct < 100; ct++) { in main()
|
| /petsc/src/dm/impls/plex/transform/impls/refine/regular/ |
| H A D | plexrefregular.c | 45 PetscErrorCode DMPlexRefineRegularGetAffineFaceTransforms(DMPlexTransform tr, DMPolytopeType ct, Pe… in DMPlexRefineRegularGetAffineFaceTransforms() argument 92 switch (ct) { in DMPlexRefineRegularGetAffineFaceTransforms() 108 …SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unsupported polytope type %s", DMPolytopeTypes[ct]); in DMPlexRefineRegularGetAffineFaceTransforms() 133 PetscErrorCode DMPlexRefineRegularGetAffineTransforms(DMPlexTransform tr, DMPolytopeType ct, PetscI… in DMPlexRefineRegularGetAffineTransforms() argument 271 switch (ct) { in DMPlexRefineRegularGetAffineTransforms() 303 …SETERRQ(PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Unsupported polytope type %s", DMPolytopeTypes[ct]); in DMPlexRefineRegularGetAffineTransforms() 309 …lexCellRefinerGetCellVertices_Regular(DMPlexCellRefiner cr, DMPolytopeType ct, PetscInt *Nv, Petsc… 328 switch (ct) { 334 …TERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "No subcell vertices for cell type %s", DMPolytopeTypes[ct]); 339 …CellRefinerGetSubcellVertices_Regular(DMPlexCellRefiner cr, DMPolytopeType ct, DMPolytopeType rct,… [all …]
|
| /petsc/src/dm/impls/plex/transform/impls/refine/1d/ |
| H A D | plexref1d.c | 18 DMPolytopeType ct; in DMPlexTransformSetUp_1D() local 21 PetscCall(DMPlexGetCellType(dm, p, &ct)); in DMPlexTransformSetUp_1D() 22 switch (ct) { in DMPlexTransformSetUp_1D() 33 … SETERRQ(PETSC_COMM_SELF, PETSC_ERR_SUP, "Cannot handle points of type %s", DMPolytopeTypes[ct]); in DMPlexTransformSetUp_1D()
|
| /petsc/src/dm/impls/plex/transform/impls/filter/ |
| H A D | plextrfilter.c | 43 DMPolytopeType ct; in DMPlexTransformSetUp_Filter() local 48 PetscCall(DMPlexGetCellType(dm, closure[cl], &ct)); in DMPlexTransformSetUp_Filter() 49 PetscCall(DMLabelSetValue(tr->trType, closure[cl], ct)); in DMPlexTransformSetUp_Filter()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex20.c | 101 DMPolytopeType ct; in SetupDiscretization() local 109 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in SetupDiscretization() 110 …simplex = DMPolytopeTypeGetNumVertices(ct) == DMPolytopeTypeGetDim(ct) + 1 ? PETSC_TRUE : PETSC_FA… in SetupDiscretization()
|
| /petsc/src/dm/impls/swarm/tests/ |
| H A D | ex11.c | 114 DMPolytopeType ct; in CreateDiscretization() local 120 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in CreateDiscretization() 126 PetscCall(PetscFVCreateDualSpace(fv, ct)); in CreateDiscretization() 134 PetscCall(PetscFVCreateDualSpace(fv, ct)); in CreateDiscretization() 139 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, 1, ct, NULL, -1, &fe)); in CreateDiscretization() 143 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, dim, ct, NULL, -1, &fe)); in CreateDiscretization() 164 DMPolytopeType ct; in CreateGradDiscretization() local 170 PetscCall(DMPlexGetCellType(dm, cStart, &ct)); in CreateGradDiscretization() 171 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, dim, ct, NULL, -1, &fe)); in CreateGradDiscretization() 175 PetscCall(PetscFECreateByCell(PETSC_COMM_SELF, dim, 2 * dim, ct, NULL, -1, &fe)); in CreateGradDiscretization()
|