Lines Matching refs:dmmoab
97 DM_Moab *dmmoab; in DMMoabCreateMoab() local
103 dmmoab = (DM_Moab *)dmmb->data; in DMMoabCreateMoab()
106 dmmoab->mbiface = new moab::Core(); in DMMoabCreateMoab()
107 dmmoab->icreatedinstance = PETSC_TRUE; in DMMoabCreateMoab()
109 dmmoab->mbiface = mbiface; in DMMoabCreateMoab()
110 dmmoab->icreatedinstance = PETSC_FALSE; in DMMoabCreateMoab()
114 dmmoab->fileset = 0; in DMMoabCreateMoab()
115 dmmoab->hlevel = 0; in DMMoabCreateMoab()
116 dmmoab->nghostrings = 0; in DMMoabCreateMoab()
123 merr = dmmoab->mbiface->create_meshset(moab::MESHSET_SET, partnset); in DMMoabCreateMoab()
127 dmmoab->pcomm = moab::ParallelComm::get_pcomm(dmmoab->mbiface, partnset, &comm); in DMMoabCreateMoab()
131 dmmoab->bs = 1; in DMMoabCreateMoab()
132 dmmoab->numFields = 1; in DMMoabCreateMoab()
133 PetscCall(PetscMalloc(dmmoab->numFields * sizeof(char *), &dmmoab->fieldNames)); in DMMoabCreateMoab()
134 PetscCall(PetscStrallocpy("DEFAULT", (char **)&dmmoab->fieldNames[0])); in DMMoabCreateMoab()
135 dmmoab->rw_dbglevel = 0; in DMMoabCreateMoab()
136 dmmoab->partition_by_rank = PETSC_FALSE; in DMMoabCreateMoab()
137 dmmoab->extra_read_options[0] = '\0'; in DMMoabCreateMoab()
138 dmmoab->extra_write_options[0] = '\0'; in DMMoabCreateMoab()
139 dmmoab->read_mode = READ_PART; in DMMoabCreateMoab()
140 dmmoab->write_mode = WRITE_PART; in DMMoabCreateMoab()
146 merr = dmmoab->mbiface->tag_get_handle(GLOBAL_ID_TAG_NAME, dmmoab->ltog_tag); in DMMoabCreateMoab()
148 if (ltog_tag) *ltog_tag = dmmoab->ltog_tag; in DMMoabCreateMoab()
151 merr = dmmoab->mbiface->tag_get_handle(MATERIAL_SET_TAG_NAME, dmmoab->material_tag); in DMMoabCreateMoab()
202 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMMoabSetInterface() local
208 dmmoab->pcomm = NULL; in DMMoabSetInterface()
210 dmmoab->mbiface = mbiface; in DMMoabSetInterface()
211 dmmoab->icreatedinstance = PETSC_FALSE; in DMMoabSetInterface()
259 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMMoabSetLocalVertices() local
263 dmmoab->vlocal->clear(); in DMMoabSetLocalVertices()
264 dmmoab->vowned->clear(); in DMMoabSetLocalVertices()
266 dmmoab->vlocal->insert(range->begin(), range->end()); in DMMoabSetLocalVertices()
271 …merr = dmmoab->pcomm->filter_pstatus(*dmmoab->vlocal, PSTATUS_NOT_OWNED, PSTATUS_NOT, -1, dmmoab->… in DMMoabSetLocalVertices()
275 tmpvtxs = moab::subtract(*dmmoab->vlocal, *dmmoab->vowned); in DMMoabSetLocalVertices()
276 …merr = dmmoab->pcomm->filter_pstatus(tmpvtxs, PSTATUS_INTERFACE, PSTATUS_OR, -1, dmmoab->vghost… in DMMoabSetLocalVertices()
278 tmpvtxs = moab::subtract(tmpvtxs, *dmmoab->vghost); in DMMoabSetLocalVertices()
279 *dmmoab->vlocal = moab::subtract(*dmmoab->vlocal, tmpvtxs); in DMMoabSetLocalVertices()
281 *dmmoab->vowned = *dmmoab->vlocal; in DMMoabSetLocalVertices()
285 dmmoab->nloc = dmmoab->vowned->size(); in DMMoabSetLocalVertices()
286 dmmoab->nghost = dmmoab->vghost->size(); in DMMoabSetLocalVertices()
288 …PetscCallMPI(MPIU_Allreduce(&dmmoab->nloc, &dmmoab->n, 1, MPI_INTEGER, MPI_SUM, ((PetscObject)dm)-… in DMMoabSetLocalVertices()
290 dmmoab->n = dmmoab->nloc; in DMMoabSetLocalVertices()
381 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMMoabSetLocalElements() local
385 dmmoab->elocal->clear(); in DMMoabSetLocalElements()
386 dmmoab->eghost->clear(); in DMMoabSetLocalElements()
387 dmmoab->elocal->insert(range->begin(), range->end()); in DMMoabSetLocalElements()
390 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->elocal, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMMoabSetLocalElements()
392 *dmmoab->eghost = moab::subtract(*range, *dmmoab->elocal); in DMMoabSetLocalElements()
394 dmmoab->neleloc = dmmoab->elocal->size(); in DMMoabSetLocalElements()
395 dmmoab->neleghost = dmmoab->eghost->size(); in DMMoabSetLocalElements()
397 …PetscCallMPI(MPIU_Allreduce(&dmmoab->neleloc, &dmmoab->nele, 1, MPI_INTEGER, MPI_SUM, ((PetscObjec… in DMMoabSetLocalElements()
398 …%" PetscInt_FMT " local and %" PetscInt_FMT " global elements.\n", dmmoab->neleloc, dmmoab->nele)); in DMMoabSetLocalElements()
400 dmmoab->nele = dmmoab->neleloc; in DMMoabSetLocalElements()
635 DM_Moab *dmmoab; in DMMoabGetMaterialBlock() local
640 dmmoab = (DM_Moab *)dm->data; in DMMoabGetMaterialBlock()
641 *mat = dmmoab->materials[dmmoab->elocal->index(ehandle)]; in DMMoabGetMaterialBlock()
665 DM_Moab *dmmoab; in DMMoabGetVertexCoordinates() local
672 dmmoab = (DM_Moab *)dm->data; in DMMoabGetVertexCoordinates()
675 if (dmmoab->hlevel) { in DMMoabGetVertexCoordinates()
676 …merr = dmmoab->hierarchy->get_coordinates(const_cast<moab::EntityHandle *>(conn), nconn, dmmoab->h… in DMMoabGetVertexCoordinates()
679 merr = dmmoab->mbiface->get_coords(conn, nconn, vpos); in DMMoabGetVertexCoordinates()
704 DM_Moab *dmmoab; in DMMoabGetVertexConnectivity() local
711 dmmoab = (DM_Moab *)dm->data; in DMMoabGetVertexConnectivity()
714 …merr = dmmoab->mbiface->get_adjacencies(&vhandle, 1, 1, true, adj_entities, moab::Interface::UNION… in DMMoabGetVertexConnectivity()
716 merr = dmmoab->mbiface->get_connectivity(&adj_entities[0], adj_entities.size(), connect); in DMMoabGetVertexConnectivity()
772 DM_Moab *dmmoab; in DMMoabGetElementConnectivity() local
781 dmmoab = (DM_Moab *)dm->data; in DMMoabGetElementConnectivity()
784 merr = dmmoab->mbiface->get_connectivity(ehandle, connect, nnodes); in DMMoabGetElementConnectivity()
810 DM_Moab *dmmoab; in DMMoabIsEntityOnBoundary() local
816 dmmoab = (DM_Moab *)dm->data; in DMMoabIsEntityOnBoundary()
819 etype = dmmoab->mbiface->type_from_handle(ent); in DMMoabIsEntityOnBoundary()
823 edim = dmmoab->mbiface->dimension_from_handle(ent); in DMMoabIsEntityOnBoundary()
827 *ent_on_boundary = ((dmmoab->bndyvtx->index(ent) >= 0) ? PETSC_TRUE : PETSC_FALSE); in DMMoabIsEntityOnBoundary()
829 if (edim == dmmoab->dim) { /* check the higher-dimensional elements first */ in DMMoabIsEntityOnBoundary()
830 if (dmmoab->bndyelems->index(ent) >= 0) *ent_on_boundary = PETSC_TRUE; in DMMoabIsEntityOnBoundary()
832 if (dmmoab->bndyfaces->index(ent) >= 0) *ent_on_boundary = PETSC_TRUE; in DMMoabIsEntityOnBoundary()
855 DM_Moab *dmmoab; in DMMoabCheckBoundaryVertices() local
862 dmmoab = (DM_Moab *)dm->data; in DMMoabCheckBoundaryVertices()
864 …for (i = 0; i < nconn; ++i) isbdvtx[i] = (dmmoab->bndyvtx->index(cnt[i]) >= 0 ? PETSC_TRUE : PETSC… in DMMoabCheckBoundaryVertices()
885 DM_Moab *dmmoab; in DMMoabGetBoundaryMarkers() local
889 dmmoab = (DM_Moab *)dm->data; in DMMoabGetBoundaryMarkers()
891 if (bdvtx) *bdvtx = dmmoab->bndyvtx; in DMMoabGetBoundaryMarkers()
892 if (bdfaces) *bdfaces = dmmoab->bndyfaces; in DMMoabGetBoundaryMarkers()
893 if (bdelems) *bdfaces = dmmoab->bndyelems; in DMMoabGetBoundaryMarkers()
901 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMDestroy_Moab() local
906 dmmoab->refct--; in DMDestroy_Moab()
907 if (!dmmoab->refct) { in DMDestroy_Moab()
908 delete dmmoab->vlocal; in DMDestroy_Moab()
909 delete dmmoab->vowned; in DMDestroy_Moab()
910 delete dmmoab->vghost; in DMDestroy_Moab()
911 delete dmmoab->elocal; in DMDestroy_Moab()
912 delete dmmoab->eghost; in DMDestroy_Moab()
913 delete dmmoab->bndyvtx; in DMDestroy_Moab()
914 delete dmmoab->bndyfaces; in DMDestroy_Moab()
915 delete dmmoab->bndyelems; in DMDestroy_Moab()
917 PetscCall(PetscFree(dmmoab->gsindices)); in DMDestroy_Moab()
918 PetscCall(PetscFree2(dmmoab->gidmap, dmmoab->lidmap)); in DMDestroy_Moab()
919 PetscCall(PetscFree(dmmoab->dfill)); in DMDestroy_Moab()
920 PetscCall(PetscFree(dmmoab->ofill)); in DMDestroy_Moab()
921 PetscCall(PetscFree(dmmoab->materials)); in DMDestroy_Moab()
922 if (dmmoab->fieldNames) { in DMDestroy_Moab()
923 for (i = 0; i < dmmoab->numFields; i++) PetscCall(PetscFree(dmmoab->fieldNames[i])); in DMDestroy_Moab()
924 PetscCall(PetscFree(dmmoab->fieldNames)); in DMDestroy_Moab()
927 if (dmmoab->nhlevels) { in DMDestroy_Moab()
928 PetscCall(PetscFree(dmmoab->hsets)); in DMDestroy_Moab()
929 dmmoab->nhlevels = 0; in DMDestroy_Moab()
930 if (!dmmoab->hlevel && dmmoab->icreatedinstance) delete dmmoab->hierarchy; in DMDestroy_Moab()
931 dmmoab->hierarchy = NULL; in DMDestroy_Moab()
934 if (dmmoab->icreatedinstance) { in DMDestroy_Moab()
935 delete dmmoab->pcomm; in DMDestroy_Moab()
936 merr = dmmoab->mbiface->delete_mesh(); in DMDestroy_Moab()
938 delete dmmoab->mbiface; in DMDestroy_Moab()
940 dmmoab->mbiface = NULL; in DMDestroy_Moab()
942 dmmoab->pcomm = NULL; in DMDestroy_Moab()
944 PetscCall(VecScatterDestroy(&dmmoab->ltog_sendrecv)); in DMDestroy_Moab()
945 PetscCall(ISLocalToGlobalMappingDestroy(&dmmoab->ltog_map)); in DMDestroy_Moab()
953 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMSetFromOptions_Moab() local
957 …rbosity level for reading and writing MOAB meshes", "DMView", dmmoab->rw_dbglevel, &dmmoab->rw_dbg… in DMSetFromOptions_Moab()
958 …ion by rank when reading MOAB meshes from file", "DMView", dmmoab->partition_by_rank, &dmmoab->par… in DMSetFromOptions_Moab()
960 …B reader to load DM from file", "DMView", dmmoab->extra_read_options, dmmoab->extra_read_options, … in DMSetFromOptions_Moab()
961 …iter to serialize DM to file", "DMView", dmmoab->extra_write_options, dmmoab->extra_write_options,… in DMSetFromOptions_Moab()
962 …parallel read mode", "DMView", MoabReadModes, (PetscEnum)dmmoab->read_mode, (PetscEnum *)&dmmoab->… in DMSetFromOptions_Moab()
963 …allel write mode", "DMView", MoabWriteModes, (PetscEnum)dmmoab->write_mode, (PetscEnum *)&dmmoab->… in DMSetFromOptions_Moab()
975 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMSetUp_Moab() local
981 …PetscCheck(dmmoab->mbiface != NULL, PETSC_COMM_WORLD, PETSC_ERR_ORDER, "Set the MOAB Interface bef… in DMSetUp_Moab()
983 …PetscCheck(dmmoab->pcomm != NULL, PETSC_COMM_WORLD, PETSC_ERR_ORDER, "Set the MOAB ParallelComm ob… in DMSetUp_Moab()
987 if (dmmoab->vlocal->empty()) { in DMSetUp_Moab()
989 merr = dmmoab->mbiface->get_entities_by_dimension(dmmoab->fileset, 0, *dmmoab->vlocal, false); in DMSetUp_Moab()
994 …merr = dmmoab->pcomm->filter_pstatus(*dmmoab->vlocal, PSTATUS_NOT_OWNED, PSTATUS_NOT, -1, dmmoab->… in DMSetUp_Moab()
999 adjs = moab::subtract(*dmmoab->vlocal, *dmmoab->vowned); in DMSetUp_Moab()
1000 …merr = dmmoab->pcomm->filter_pstatus(adjs, PSTATUS_GHOST | PSTATUS_INTERFACE, PSTATUS_OR, -1, dmmo… in DMSetUp_Moab()
1002 adjs = moab::subtract(adjs, *dmmoab->vghost); in DMSetUp_Moab()
1003 *dmmoab->vlocal = moab::subtract(*dmmoab->vlocal, adjs); in DMSetUp_Moab()
1005 *dmmoab->vowned = *dmmoab->vlocal; in DMSetUp_Moab()
1009 dmmoab->nloc = dmmoab->vowned->size(); in DMSetUp_Moab()
1010 dmmoab->nghost = dmmoab->vghost->size(); in DMSetUp_Moab()
1013 …PetscCallMPI(MPIU_Allreduce(&dmmoab->nloc, &dmmoab->n, 1, MPI_INTEGER, MPI_SUM, ((PetscObject)dm)-… in DMSetUp_Moab()
1014 …scInt_FMT ", ghosted - %" PetscInt_FMT ".\n", dmmoab->fileset, dmmoab->vlocal->size(), dmmoab->nlo… in DMSetUp_Moab()
1016 dmmoab->n = dmmoab->nloc; in DMSetUp_Moab()
1022 dmmoab->eghost->clear(); in DMSetUp_Moab()
1026 dmmoab->elocal->clear(); in DMSetUp_Moab()
1027 merr = dmmoab->mbiface->get_entities_by_dimension(dmmoab->fileset, i, *dmmoab->elocal, false); in DMSetUp_Moab()
1031 if (dmmoab->elocal->size()) { in DMSetUp_Moab()
1032 dmmoab->dim = i; in DMSetUp_Moab()
1037 PetscCall(DMSetDimension(dm, dmmoab->dim)); in DMSetUp_Moab()
1041 *dmmoab->eghost = *dmmoab->elocal; in DMSetUp_Moab()
1042 … merr = dmmoab->pcomm->filter_pstatus(*dmmoab->elocal, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1044 *dmmoab->eghost = moab::subtract(*dmmoab->eghost, *dmmoab->elocal); in DMSetUp_Moab()
1047 dmmoab->neleloc = dmmoab->elocal->size(); in DMSetUp_Moab()
1048 dmmoab->neleghost = dmmoab->eghost->size(); in DMSetUp_Moab()
1051 …PetscCallMPI(MPIU_Allreduce(&dmmoab->neleloc, &dmmoab->nele, 1, MPI_INTEGER, MPI_SUM, ((PetscObjec… in DMSetUp_Moab()
1052 …ned - %" PetscInt_FMT ", ghosted - %" PetscInt_FMT ".\n", dmmoab->dim, dmmoab->neleloc, dmmoab->ne… in DMSetUp_Moab()
1054 dmmoab->nele = dmmoab->neleloc; in DMSetUp_Moab()
1058 bs = dmmoab->bs; in DMSetUp_Moab()
1059 if (!dmmoab->ltog_tag) { in DMSetUp_Moab()
1063 merr = dmmoab->mbiface->tag_get_handle(GLOBAL_ID_TAG_NAME, dmmoab->ltog_tag); in DMSetUp_Moab()
1067 totsize = dmmoab->vlocal->size(); in DMSetUp_Moab()
1068 …dmmoab->nloc + dmmoab->nghost, PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "Mismatch between local … in DMSetUp_Moab()
1069 PetscCall(PetscCalloc1(totsize, &dmmoab->gsindices)); in DMSetUp_Moab()
1072 merr = dmmoab->mbiface->tag_get_data(dmmoab->ltog_tag, *dmmoab->vowned, &dmmoab->gsindices[0]); in DMSetUp_Moab()
1074 if (dmmoab->nghost) { /* next get the ghosted indices */ in DMSetUp_Moab()
1075 …merr = dmmoab->mbiface->tag_get_data(dmmoab->ltog_tag, *dmmoab->vghost, &dmmoab->gsindices[dmmoab-… in DMSetUp_Moab()
1080 dmmoab->lminmax[0] = dmmoab->lminmax[1] = dmmoab->gsindices[0]; in DMSetUp_Moab()
1082 if (dmmoab->lminmax[0] > dmmoab->gsindices[i]) dmmoab->lminmax[0] = dmmoab->gsindices[i]; in DMSetUp_Moab()
1083 if (dmmoab->lminmax[1] < dmmoab->gsindices[i]) dmmoab->lminmax[1] = dmmoab->gsindices[i]; in DMSetUp_Moab()
1086 …PetscCallMPI(MPIU_Allreduce(&dmmoab->lminmax[0], &dmmoab->gminmax[0], 1, MPI_INT, MPI_MIN, ((Petsc… in DMSetUp_Moab()
1087 …PetscCallMPI(MPIU_Allreduce(&dmmoab->lminmax[1], &dmmoab->gminmax[1], 1, MPI_INT, MPI_MAX, ((Petsc… in DMSetUp_Moab()
1090 …for (i = 0; i < totsize; ++i) dmmoab->gsindices[i] -= dmmoab->gminmax[0]; /* zero based index need… in DMSetUp_Moab()
1091 dmmoab->lminmax[0] -= dmmoab->gminmax[0]; in DMSetUp_Moab()
1092 dmmoab->lminmax[1] -= dmmoab->gminmax[0]; in DMSetUp_Moab()
1094 …[%" PetscInt_FMT ", %" PetscInt_FMT "]\n", dmmoab->lminmax[0], dmmoab->lminmax[1], dmmoab->gminmax… in DMSetUp_Moab()
1096 …dmmoab->bs == dmmoab->numFields || dmmoab->bs == 1, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Mismatc… in DMSetUp_Moab()
1097 dmmoab->numFields); in DMSetUp_Moab()
1100 dmmoab->seqstart = dmmoab->mbiface->id_from_handle(dmmoab->vlocal->front()); in DMSetUp_Moab()
1101 dmmoab->seqend = dmmoab->mbiface->id_from_handle(dmmoab->vlocal->back()); in DMSetUp_Moab()
1102 …: Local [min, max] - [%" PetscInt_FMT ", %" PetscInt_FMT "]\n", dmmoab->seqstart, dmmoab->seqend)); in DMSetUp_Moab()
1104 …PetscCall(PetscMalloc2(dmmoab->seqend - dmmoab->seqstart + 1, &dmmoab->gidmap, dmmoab->seqend - dm… in DMSetUp_Moab()
1105 PetscCall(PetscMalloc1(totsize * dmmoab->numFields, &lgmap)); in DMSetUp_Moab()
1109 …for (moab::Range::iterator iter = dmmoab->vowned->begin(); iter != dmmoab->vowned->end(); iter++, … in DMSetUp_Moab()
1110 vent = dmmoab->mbiface->id_from_handle(*iter) - dmmoab->seqstart; in DMSetUp_Moab()
1111 dmmoab->gidmap[vent] = dmmoab->gsindices[i]; in DMSetUp_Moab()
1112 dmmoab->lidmap[vent] = i; in DMSetUp_Moab()
1113 …for (f = 0; f < dmmoab->numFields; f++, j++) lgmap[j] = (bs > 1 ? dmmoab->gsindices[i] * dmmoab->n… in DMSetUp_Moab()
1116 …for (moab::Range::iterator iter = dmmoab->vghost->begin(); iter != dmmoab->vghost->end(); iter++, … in DMSetUp_Moab()
1117 vent = dmmoab->mbiface->id_from_handle(*iter) - dmmoab->seqstart; in DMSetUp_Moab()
1118 dmmoab->gidmap[vent] = dmmoab->gsindices[i]; in DMSetUp_Moab()
1119 dmmoab->lidmap[vent] = i; in DMSetUp_Moab()
1120 …for (f = 0; f < dmmoab->numFields; f++, j++) lgmap[j] = (bs > 1 ? dmmoab->gsindices[i] * dmmoab->n… in DMSetUp_Moab()
1132 PetscCall(VecGetOwnershipRange(global, &dmmoab->vstart, &dmmoab->vend)); in DMSetUp_Moab()
1135 …PetscCall(ISCreateStride(((PetscObject)dm)->comm, dmmoab->nloc * dmmoab->numFields, dmmoab->vstart… in DMSetUp_Moab()
1138 …PetscCall(ISCreateGeneral(((PetscObject)dm)->comm, dmmoab->nloc * dmmoab->numFields, &lgmap[0], PE… in DMSetUp_Moab()
1141 if (!dmmoab->ltog_map) { in DMSetUp_Moab()
1143 …MappingCreate(((PetscObject)dm)->comm, dmmoab->bs, totsize * dmmoab->numFields, lgmap, PETSC_COPY_… in DMSetUp_Moab()
1147 PetscCall(VecScatterCreate(local, from, global, to, &dmmoab->ltog_sendrecv)); in DMSetUp_Moab()
1157 dmmoab->bndyvtx = new moab::Range(); in DMSetUp_Moab()
1158 dmmoab->bndyfaces = new moab::Range(); in DMSetUp_Moab()
1159 dmmoab->bndyelems = new moab::Range(); in DMSetUp_Moab()
1161 if (!dmmoab->hlevel) { in DMSetUp_Moab()
1166 moab::Skinner skinner(dmmoab->mbiface); in DMSetUp_Moab()
1169 …merr = skinner.find_skin(dmmoab->fileset, *dmmoab->elocal, false, *dmmoab->bndyfaces, NULL, true, … in DMSetUp_Moab()
1174 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->bndyfaces, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1176 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->bndyfaces, PSTATUS_INTERFACE, PSTATUS_NOT); in DMSetUp_Moab()
1181 merr = dmmoab->mbiface->get_connectivity(*dmmoab->bndyfaces, *dmmoab->bndyvtx, false); in DMSetUp_Moab()
1183 …merr = dmmoab->mbiface->get_adjacencies(*dmmoab->bndyvtx, dmmoab->dim, false, *dmmoab->bndyelems, … in DMSetUp_Moab()
1187 …for (moab::Range::iterator iter = dmmoab->elocal->begin(); iter != dmmoab->elocal->end(); iter++) { in DMSetUp_Moab()
1189 if (dmmoab->hierarchy->is_entity_on_boundary(elemHandle)) { in DMSetUp_Moab()
1190 dmmoab->bndyelems->insert(elemHandle); in DMSetUp_Moab()
1193 merr = dmmoab->hierarchy->get_connectivity(elemHandle, dmmoab->hlevel, connect); in DMSetUp_Moab()
1196 … if (dmmoab->hierarchy->is_entity_on_boundary(connect[iv])) dmmoab->bndyvtx->insert(connect[iv]); in DMSetUp_Moab()
1199 merr = dmmoab->hierarchy->get_adjacencies(elemHandle, dmmoab->dim - 1, faces); in DMSetUp_Moab()
1202 … if (dmmoab->hierarchy->is_entity_on_boundary(faces[ifa])) dmmoab->bndyfaces->insert(faces[ifa]); in DMSetUp_Moab()
1207 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->bndyvtx, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1209 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->bndyfaces, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1211 merr = dmmoab->pcomm->filter_pstatus(*dmmoab->bndyelems, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1215 … boundary faces and %zu boundary elements.\n", dmmoab->bndyvtx->size(), dmmoab->bndyfaces->size(),… in DMSetUp_Moab()
1219 PetscCall(PetscCalloc1(dmmoab->elocal->size(), &dmmoab->materials)); in DMSetUp_Moab()
1223 …merr = dmmoab->mbiface->get_entities_by_type_and_tag(dmmoab->fileset, moab::MBENTITYSET, &dmmoab->… in DMSetUp_Moab()
1229 merr = dmmoab->mbiface->get_entities_by_dimension(msets[i], dmmoab->dim, msetelems, true); in DMSetUp_Moab()
1233 merr = dmmoab->pcomm->filter_pstatus(msetelems, PSTATUS_NOT_OWNED, PSTATUS_NOT); in DMSetUp_Moab()
1239 … merr = dmmoab->mbiface->tag_get_data(dmmoab->material_tag, &mset, 1, &partID); in DMSetUp_Moab()
1242 …for (unsigned j = 0; j < msetelems.size(); ++j) dmmoab->materials[dmmoab->elocal->index(msetelems[… in DMSetUp_Moab()
1268 DM_Moab *dmmoab; in DMMoabCreateVertices() local
1275 dmmoab = (DM_Moab *)dm->data; in DMMoabCreateVertices()
1278 merr = dmmoab->mbiface->create_vertices(&coords[0], nverts, verts); in DMMoabCreateVertices()
1280 merr = dmmoab->mbiface->add_entities(dmmoab->fileset, verts); in DMMoabCreateVertices()
1308 DM_Moab *dmmoab; in DMMoabCreateElement() local
1315 dmmoab = (DM_Moab *)dm->data; in DMMoabCreateElement()
1318 merr = dmmoab->mbiface->create_element(type, conn, nverts, elem); in DMMoabCreateElement()
1320 merr = dmmoab->mbiface->add_entities(dmmoab->fileset, &elem, 1); in DMMoabCreateElement()
1346 DM_Moab *dmmoab; in DMMoabCreateSubmesh() local
1353 dmmoab = (DM_Moab *)dm->data; in DMMoabCreateSubmesh()
1356 …PetscCall(DMMoabCreateMoab(((PetscObject)dm)->comm, dmmoab->mbiface, &dmmoab->ltog_tag, NULL, newd… in DMMoabCreateSubmesh()
1369 merr = ndmmoab->mbiface->add_entities(ndmmoab->fileset, *dmmoab->vlocal); in DMMoabCreateSubmesh()
1371 merr = ndmmoab->mbiface->add_entities(ndmmoab->fileset, *dmmoab->elocal); in DMMoabCreateSubmesh()
1375 PetscCall(DMMoabSetFieldNames(*newdm, dmmoab->numFields, dmmoab->fieldNames)); in DMMoabCreateSubmesh()
1381 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMMoabView_Ascii() local
1391 …scCall(PetscViewerASCIIPrintf(viewer, "%s in %" PetscInt_FMT " dimensions:\n", name, dmmoab->dim)); in DMMoabView_Ascii()
1392 …else PetscCall(PetscViewerASCIIPrintf(viewer, "Mesh in %" PetscInt_FMT " dimensions:\n", dmmoab->d… in DMMoabView_Ascii()
1396 … ", vertices=%" PetscInt_FMT ", blocks=%" PetscInt_FMT "\n", dmmoab->nele, dmmoab->n, dmmoab->bs)); in DMMoabView_Ascii()
1401 …ewer, "cells=%zu, faces=%zu, vertices=%zu\n", dmmoab->bndyelems->size(), dmmoab->bndyfaces->size()… in DMMoabView_Ascii()
1405 …PetscCall(PetscViewerASCIIPrintf(viewer, "Fields: %" PetscInt_FMT " components\n", dmmoab->numFiel… in DMMoabView_Ascii()
1408 …for (int i = 0; i < dmmoab->numFields; ++i) PetscCall(PetscViewerASCIIPrintf(viewer, "[%" PetscInt… in DMMoabView_Ascii()