Lines Matching refs:dm
23 PetscErrorCode DMFieldToDSField(DM dm, DMLabel domain_label, PetscInt dm_field, PetscInt *ds_field)… in DMFieldToDSField() argument
30 PetscCall(DMGetRegionDS(dm, domain_label, &field_is, &ds, NULL)); in DMFieldToDSField()
41 …PetscCheck(*ds_field != -1, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "Could not find dm_fi… in DMFieldToDSField()
80 PetscErrorCode DMPlexCeedElemRestrictionCreate(Ceed ceed, DM dm, DMLabel domain_label, PetscInt lab… in DMPlexCeedElemRestrictionCreate() argument
94 …PetscCall(PetscObjectContainerQuery((PetscObject)dm, container_name, (void **)&container_restricti… in DMPlexCeedElemRestrictionCreate()
105 …PetscCall(DMPlexGetLocalOffsets(dm, domain_label, label_value, height, dm_field, &num_elem, &elem_… in DMPlexCeedElemRestrictionCreate()
114 PetscCall(PetscObjectContainerCompose((PetscObject)dm, container_name, container_restriction, in DMPlexCeedElemRestrictionCreate()
137 PetscErrorCode DMPlexCeedElemRestrictionCoordinateCreate(Ceed ceed, DM dm, DMLabel domain_label, Pe… in DMPlexCeedElemRestrictionCoordinateCreate() argument
142 PetscCall(DMGetCellCoordinateDM(dm, &dm_coord)); in DMPlexCeedElemRestrictionCoordinateCreate()
144 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMPlexCeedElemRestrictionCoordinateCreate()
166 static PetscErrorCode DMPlexCeedElemRestrictionStridedCreate(Ceed ceed, DM dm, DMLabel domain_label… in DMPlexCeedElemRestrictionStridedCreate() argument
176 PetscCall(DMPlexGetDepth(dm, &depth)); in DMPlexCeedElemRestrictionStridedCreate()
177 PetscCall(DMPlexGetDepthLabel(dm, &depth_label)); in DMPlexCeedElemRestrictionStridedCreate()
206 PetscCall(DMGetRegionDS(dm, domain_label, NULL, &ds, NULL)); in DMPlexCeedElemRestrictionStridedCreate()
207 PetscCall(DMFieldToDSField(dm, domain_label, dm_field, &ds_field)); in DMPlexCeedElemRestrictionStridedCreate()
222 PetscCall(DMGetRegionDS(dm, domain_label, NULL, &ds, NULL)); in DMPlexCeedElemRestrictionStridedCreate()
223 PetscCall(DMFieldToDSField(dm, domain_label, dm_field, &ds_field)); in DMPlexCeedElemRestrictionStridedCreate()
252 PetscErrorCode DMPlexCeedElemRestrictionQDataCreate(Ceed ceed, DM dm, DMLabel domain_label, PetscIn… in DMPlexCeedElemRestrictionQDataCreate() argument
255 …PetscCall(DMPlexCeedElemRestrictionStridedCreate(ceed, dm, domain_label, label_value, height, q_da… in DMPlexCeedElemRestrictionQDataCreate()
274 PetscErrorCode DMPlexCeedElemRestrictionCollocatedCreate(Ceed ceed, DM dm, DMLabel domain_label, Pe… in DMPlexCeedElemRestrictionCollocatedCreate() argument
277 …PetscCall(DMPlexCeedElemRestrictionStridedCreate(ceed, dm, domain_label, label_value, height, q_da… in DMPlexCeedElemRestrictionCollocatedCreate()
495 static inline PetscErrorCode GetGlobalDMPlexPolytopeType(DM dm, DMLabel domain_label, PetscInt labe… in GetGlobalDMPlexPolytopeType() argument
506 PetscCall(DMPlexGetDepth(dm, &depth)); in GetGlobalDMPlexPolytopeType()
507 PetscCall(DMPlexGetDepthLabel(dm, &depth_label)); in GetGlobalDMPlexPolytopeType()
512 …PetscCall(DMGetFirstLabeledPoint(dm, dm, domain_label ? domain_label : depth_label, 1, ids, height… in GetGlobalDMPlexPolytopeType()
513 if (first_point != -1) PetscCall(DMPlexGetCellType(dm, first_point, cell_type)); in GetGlobalDMPlexPolytopeType()
519 …duce(&cell_type_local, &cell_type_global, 1, MPIU_INT, MPI_MAX, PetscObjectComm((PetscObject)dm))); in GetGlobalDMPlexPolytopeType()
536 static inline PetscErrorCode GetClosurePermutationAndFieldOffsetAtDepth(DM dm, PetscInt depth, Pets… in GetClosurePermutationAndFieldOffsetAtDepth() argument
545 PetscCall(DMGetDimension(dm, &dim)); in GetClosurePermutationAndFieldOffsetAtDepth()
546 PetscCall(DMGetLocalSection(dm, §ion)); in GetClosurePermutationAndFieldOffsetAtDepth()
557 PetscCall(DMGetField(dm, f, NULL, &obj)); in GetClosurePermutationAndFieldOffsetAtDepth()
572 … PetscCall(PetscSectionGetClosurePermutation(section, (PetscObject)dm, depth, size, permutation)); in GetClosurePermutationAndFieldOffsetAtDepth()
591 static inline PetscErrorCode ComputeFieldTabulationP2C(DM dm, PetscInt field, PetscInt face, PetscT… in ComputeFieldTabulationP2C() argument
607 PetscCall(DMPlexIsSimplex(dm, &is_simplex)); in ComputeFieldTabulationP2C()
609 …PetscCall(GetClosurePermutationAndFieldOffsetAtDepth(dm, dim, field, &permutation, &field_offset)); in ComputeFieldTabulationP2C()
764 PetscErrorCode DMPlexCeedBasisCreate(Ceed ceed, DM dm, DMLabel domain_label, PetscInt label_value, … in DMPlexCeedBasisCreate() argument
766 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in DMPlexCeedBasisCreate()
779 PetscCall(PetscObjectContainerQuery((PetscObject)dm, container_name, (void **)&container_basis)); in DMPlexCeedBasisCreate()
794 PetscCall(DMGetRegionDS(dm, domain_label, NULL, &ds, NULL)); in DMPlexCeedBasisCreate()
795 PetscCall(DMFieldToDSField(dm, domain_label, dm_field, &ds_field)); in DMPlexCeedBasisCreate()
818 PetscCall(GetGlobalDMPlexPolytopeType(dm, domain_label, label_value, height, &cell_type)); in DMPlexCeedBasisCreate()
826 PetscCall(ComputeFieldTabulationP2C(dm, dm_field, face, basis_tabulation, &interp, &grad)); in DMPlexCeedBasisCreate()
847 PetscCall(ComputeFieldTabulationP2C(dm, dm_field, 0, basis_tabulation, &interp_1d, &grad_1d)); in DMPlexCeedBasisCreate()
859 …PetscCall(PetscObjectContainerCompose((PetscObject)dm, container_name, container_basis, (PetscCtxD… in DMPlexCeedBasisCreate()
879 PetscErrorCode DMPlexCeedBasisCoordinateCreate(Ceed ceed, DM dm, DMLabel domain_label, PetscInt lab… in DMPlexCeedBasisCoordinateCreate() argument
883 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMPlexCeedBasisCoordinateCreate()
925 PetscErrorCode DMPlexCeedCoordinateCreateField(Ceed ceed, DM dm, DMLabel domain_label, PetscInt lab… in DMPlexCeedCoordinateCreateField() argument
930 …PetscCall(DMPlexCeedElemRestrictionCoordinateCreate(ceed, dm, domain_label, label_value, height, &… in DMPlexCeedCoordinateCreateField()
935 …if (basis) PetscCall(DMPlexCeedBasisCoordinateCreate(ceed, dm, domain_label, label_value, height, … in DMPlexCeedCoordinateCreateField()
943 PetscCall(DMGetCellCoordinateDM(dm, &cdm)); in DMPlexCeedCoordinateCreateField()
945 PetscCall(DMGetCellCoordinatesLocal(dm, &X_loc)); in DMPlexCeedCoordinateCreateField()
947 PetscCall(DMGetCoordinatesLocal(dm, &X_loc)); in DMPlexCeedCoordinateCreateField()
984 PetscErrorCode DMPlexCeedBasisCellToFaceCreate(Ceed ceed, DM dm, DMLabel domain_label, PetscInt lab… in DMPlexCeedBasisCellToFaceCreate() argument
993 PetscCall(DMGetRegionDS(dm, domain_label, NULL, &ds, NULL)); in DMPlexCeedBasisCellToFaceCreate()
994 PetscCall(DMFieldToDSField(dm, domain_label, dm_field, &ds_field)); in DMPlexCeedBasisCellToFaceCreate()
1007 PetscCall(GetGlobalDMPlexPolytopeType(dm, domain_label, label_value, height, &cell_type)); in DMPlexCeedBasisCellToFaceCreate()
1009 …PetscCheck(elem_topo, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "DMPlex topology not suppor… in DMPlexCeedBasisCellToFaceCreate()
1017 PetscCall(ComputeFieldTabulationP2C(dm, dm_field, face, basis_tabulation, &interp, &grad)); in DMPlexCeedBasisCellToFaceCreate()
1040 PetscErrorCode DMPlexCeedBasisCellToFaceCoordinateCreate(Ceed ceed, DM dm, DMLabel domain_label, Pe… in DMPlexCeedBasisCellToFaceCoordinateCreate() argument
1045 PetscCall(DMGetCellCoordinateDM(dm, &dm_coord)); in DMPlexCeedBasisCellToFaceCoordinateCreate()
1047 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMPlexCeedBasisCellToFaceCoordinateCreate()
1070 PetscInt num_fields, const PetscInt *field_sizes, DM dm) { in DMSetupByOrderBegin_FEM() argument
1074 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in DMSetupByOrderBegin_FEM()
1077 PetscCall(DMPlexIsSimplex(dm, &is_simplex)); in DMSetupByOrderBegin_FEM()
1080 PetscCall(DMGetDimension(dm, &dim)); in DMSetupByOrderBegin_FEM()
1087 PetscCall(DMAddField(dm, NULL, (PetscObject)fe)); in DMSetupByOrderBegin_FEM()
1090 PetscCall(DMCreateDS(dm)); in DMSetupByOrderBegin_FEM()
1100 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMSetupByOrderBegin_FEM()
1101 PetscCall(DMGetCoordinateDim(dm, &num_comp_coord)); in DMSetupByOrderBegin_FEM()
1111 PetscCall(DMSetCoordinateDisc(dm, fe_coord_new, PETSC_FALSE, PETSC_TRUE)); in DMSetupByOrderBegin_FEM()
1127 PetscErrorCode DMSetupByOrderEnd_FEM(PetscBool setup_coords, DM dm) { in DMSetupByOrderEnd_FEM() argument
1131 PetscCall(DMPlexIsSimplex(dm, &is_simplex)); in DMSetupByOrderEnd_FEM()
1134 PetscCall(DMPlexSetClosurePermutationTensor(dm, PETSC_DETERMINE, NULL)); in DMSetupByOrderEnd_FEM()
1138 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMSetupByOrderEnd_FEM()
1142 PetscCall(DMLocalizeCoordinates(dm)); // Must localize after tensor closure setting in DMSetupByOrderEnd_FEM()
1163 PetscInt num_fields, const PetscInt *field_sizes, DM dm) { in DMSetupByOrder_FEM() argument
1165 …erBegin_FEM(setup_faces, setup_coords, degree, coord_order, q_extra, num_fields, field_sizes, dm)); in DMSetupByOrder_FEM()
1166 PetscCall(DMSetupByOrderEnd_FEM(setup_coords, dm)); in DMSetupByOrder_FEM()
1184 static PetscErrorCode DMGetStratumISAtHeight(DM dm, const char name[], PetscInt value, PetscInt hei… in DMGetStratumISAtHeight() argument
1190 PetscCall(DMPlexGetDepth(dm, &depth)); in DMGetStratumISAtHeight()
1191 PetscCall(DMPlexGetDepthLabel(dm, &depth_label)); in DMGetStratumISAtHeight()
1193 PetscCall(DMGetStratumIS(dm, name, value, &label_is)); in DMGetStratumISAtHeight()
1210 PetscErrorCode DMPlexCreateFaceLabel(DM dm, PetscInt dm_face, char **face_label_name) { in DMPlexCreateFaceLabel() argument
1215 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in DMPlexCreateFaceLabel()
1225 PetscCall(DMHasLabel(dm, *face_label_name, &has_label_already)); in DMPlexCreateFaceLabel()
1228 PetscCall(DMCreateLabel(dm, *face_label_name)); in DMPlexCreateFaceLabel()
1229 PetscCall(DMGetLabel(dm, *face_label_name, &face_label)); in DMPlexCreateFaceLabel()
1232 PetscCall(DMGetStratumISAtHeight(dm, "Face Sets", dm_face, face_height, &is_face_points)); in DMPlexCreateFaceLabel()
1244 PetscCall(DMPlexGetSupport(dm, face_point, &face_support)); in DMPlexCreateFaceLabel()
1245 PetscCall(DMPlexGetSupportSize(dm, face_point, &face_support_size)); in DMPlexCreateFaceLabel()
1251 PetscCall(DMPlexGetCone(dm, cell_point, &cell_cone)); in DMPlexCreateFaceLabel()
1252 PetscCall(DMPlexGetConeSize(dm, cell_point, &cell_cone_size)); in DMPlexCreateFaceLabel()
1261 PetscCall(DMPlexLabelAddFaceCells(dm, face_label)); in DMPlexCreateFaceLabel()
1278 PetscErrorCode DMLabelCreateGlobalValueArray(DM dm, DMLabel label, PetscInt *num_values, PetscInt *… in DMLabelCreateGlobalValueArray() argument
1281 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in DMLabelCreateGlobalValueArray()
1339 PetscErrorCode DMGetCoordinateNumComps(DM dm, PetscInt *num_comp) { in DMGetCoordinateNumComps() argument
1345 PetscCall(DMGetCoordinateDM(dm, &dm_coord)); in DMGetCoordinateNumComps()
1358 PetscErrorCode DMGetFieldNumComps(DM dm, PetscInt field, PetscInt *num_comp) { in DMGetFieldNumComps() argument
1362 PetscCall(DMGetLocalSection(dm, §ion)); in DMGetFieldNumComps()