| /petsc/src/dm/impls/moab/ |
| H A D | dmmoab.cxx | 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()
[all …]
|
| H A D | dmmbfield.cxx | 22 DM_Moab *dmmoab; in DMMoabSetFieldVector() local
31 dmmoab = (DM_Moab *)dm->data; in DMMoabSetFieldVector()
33 …(ifield >= dmmoab->numFields), PETSC_COMM_SELF, PETSC_ERR_ARG_OUTOFRANGE, "The field %d should be … in DMMoabSetFieldVector()
36 …merr = dmmoab->mbiface->tag_get_handle(dmmoab->fieldNames[ifield], 1, moab::MB_TYPE_DOUBLE, ntag, … in DMMoabSetFieldVector()
41 merr = dmmoab->mbiface->tag_get_name(vtag, tag_name); in DMMoabSetFieldVector()
45 merr = dmmoab->mbiface->tag_set_data(ntag, *dmmoab->vowned, (const void *)varray); in DMMoabSetFieldVector()
49 PetscCall(PetscMalloc1(dmmoab->nloc, &farray)); in DMMoabSetFieldVector()
51 merr = dmmoab->mbiface->tag_get_data(vtag, *dmmoab->vowned, (void *)farray); in DMMoabSetFieldVector()
53 merr = dmmoab->mbiface->tag_set_data(ntag, *dmmoab->vowned, (const void *)farray); in DMMoabSetFieldVector()
59 merr = dmmoab->pcomm->exchange_tags(ntag, *dmmoab->vowned); in DMMoabSetFieldVector()
[all …]
|
| H A D | dmmbvec.cxx | 122 DM_Moab *dmmoab; in DMMoabVecGetArray() local
135 dmmoab = (DM_Moab *)dm->data; in DMMoabVecGetArray()
166 merr = dmmoab->pcomm->exchange_tags(vtag, *dmmoab->vlocal); in DMMoabVecGetArray()
170 PetscCall(PetscMalloc1((dmmoab->nloc + dmmoab->nghost) * dmmoab->numFields, varray)); in DMMoabVecGetArray()
173 …merr = dmmoab->mbiface->tag_iterate(vtag, dmmoab->vlocal->begin(), dmmoab->vlocal->end(), count, r… in DMMoabVecGetArray()
175 …dmmoab->nloc + dmmoab->nghost, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Mismatch between local verti… in DMMoabVecGetArray()
178 …for (moab::Range::iterator iter = dmmoab->vlocal->begin(); iter != dmmoab->vlocal->end(); iter++) { in DMMoabVecGetArray()
179 …for (f = 0; f < dmmoab->numFields; f++, i++) (*varray)[dmmoab->lidmap[(PetscInt)*iter - dmmoab->se… in DMMoabVecGetArray()
202 DM_Moab *dmmoab; in DMMoabVecRestoreArray() local
215 dmmoab = (DM_Moab *)dm->data; in DMMoabVecRestoreArray()
[all …]
|
| H A D | dmmbmat.cxx | 13 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMCreateMatrix_Moab() local
26 nlsiz = (tmp ? dmmoab->nloc : dmmoab->nloc * dmmoab->numFields); in DMCreateMatrix_Moab()
36 …PetscCall(MatSetSizes(A, dmmoab->nloc * dmmoab->numFields, dmmoab->nloc * dmmoab->numFields, PETSC… in DMCreateMatrix_Moab()
38 PetscCall(MatSetBlockSize(A, dmmoab->bs)); in DMCreateMatrix_Moab()
42 …PetscCheck(dmmoab->ltog_map, ((PetscObject)dm)->comm, PETSC_ERR_ORDER, "Cannot create a DMMoab Mat… in DMCreateMatrix_Moab()
43 PetscCall(MatSetLocalToGlobalMapping(A, dmmoab->ltog_map, dmmoab->ltog_map)); in DMCreateMatrix_Moab()
48 PetscCall(MatSeqBAIJSetPreallocation(A, dmmoab->bs, innz, nnz)); in DMCreateMatrix_Moab()
49 PetscCall(MatMPIBAIJSetPreallocation(A, dmmoab->bs, innz, nnz, ionz, onz)); in DMCreateMatrix_Moab()
67 DM_Moab *dmmoab = (DM_Moab *)dm->data; in DMMoab_Compute_NNZ_From_Connectivity() local
75 bs = dmmoab->bs; in DMMoab_Compute_NNZ_From_Connectivity()
[all …]
|
| H A D | dmmbio.cxx | 55 DM_Moab *dmmoab; in DMMoabOutput() local
62 dmmoab = (DM_Moab *)dm->data; in DMMoabOutput()
69 …tWriteOptions_Private(dmmoab->pcomm->get_id(), dmmoab->pcomm->size(), dmmoab->dim, dmmoab->write_m… in DMMoabOutput()
71 …etscCall(DMMoab_GetWriteOptions_Private(0, 1, dmmoab->dim, dmmoab->write_mode, dmmoab->rw_dbglevel… in DMMoabOutput()
79 merr = dmmoab->mbiface->write_file(filename, NULL, writeopts, &dmmoab->fileset, 1); in DMMoabOutput()
80 MBERRVM(dmmoab->mbiface, "Writing output of DMMoab failed.", merr); in DMMoabOutput()
|
| H A D | dmmbutil.cxx | 503 DM_Moab *dmmoab; in DMMoabCreateBoxMesh() local
539 dmmoab = (DM_Moab *)(*dm)->data; in DMMoabCreateBoxMesh()
540 mbImpl = dmmoab->mbiface; in DMMoabCreateBoxMesh()
542 pcomm = dmmoab->pcomm; in DMMoabCreateBoxMesh()
548 global_id_tag = dmmoab->ltog_tag; in DMMoabCreateBoxMesh()
549 dmmoab->dim = dim; in DMMoabCreateBoxMesh()
550 dmmoab->nghostrings = nghost; in DMMoabCreateBoxMesh()
551 dmmoab->refct = 1; in DMMoabCreateBoxMesh()
554 merr = mbImpl->create_meshset(moab::MESHSET_SET, dmmoab->fileset); in DMMoabCreateBoxMesh()
693 if (dmmoab->fileset) { in DMMoabCreateBoxMesh()
[all …]
|
| H A D | dmmbmg.cxx | 34 DM_Moab *dmmoab; in DMMoabGenerateHierarchy() local
41 dmmoab = (DM_Moab *)dm->data; in DMMoabGenerateHierarchy()
49 dmmoab->nhlevels = nlevels; in DMMoabGenerateHierarchy()
53 …dmmoab->hierarchy = new moab::NestedRefine(dynamic_cast<moab::Core *>(dmmoab->mbiface), dmmoab->pc… in DMMoabGenerateHierarchy()
55 …dmmoab->hierarchy = new moab::NestedRefine(dynamic_cast<moab::Core *>(dmmoab->mbiface), NULL, dmmo… in DMMoabGenerateHierarchy()
58 PetscCall(PetscMalloc1(nlevels + 1, &dmmoab->hsets)); in DMMoabGenerateHierarchy()
61 merr = dmmoab->hierarchy->generate_mesh_hierarchy(nlevels, pdegrees, hsets, false); in DMMoabGenerateHierarchy()
65 if (dmmoab->pcomm->size() > 1) { in DMMoabGenerateHierarchy()
66 merr = dmmoab->hierarchy->exchange_ghosts(hsets, dmmoab->nghostrings); in DMMoabGenerateHierarchy()
72 dmmoab->hsets[0] = hsets[0]; in DMMoabGenerateHierarchy()
[all …]
|