[MOAB-dev] Polyhedra elements in parallell

[WARNER, JAMES E. (LARC-D309)]

Hi Iulian,

The face that the two polyhedron elements have in common was given a different ID in each element's connectivity… so that is probably the issue. I'm going to look into fixing that and I'll let you know if it's resolved afterwards. Sorry about that!

-Jim


From: <Grindeanu>, "Iulian R." <iulian at mcs.anl.gov<mailto:iulian at mcs.anl.gov>>
Date: Monday, April 27, 2015 4:42 PM
To: "WARNER, JAMES E. (LARC-D309)" <james.e.warner at nasa.gov<mailto:james.e.warner at nasa.gov>>, "moab-dev at mcs.anl.gov<mailto:moab-dev at mcs.anl.gov>" <moab-dev at mcs.anl.gov<mailto:moab-dev at mcs.anl.gov>>
Subject: RE: [MOAB-dev] Polyhedra elements in parallell


my guess is that mbpart (with default options) distributed all faces to different processors, so maybe meshes are not "complete" in the partition; I think that mbpart will need to change for polyhedral meshes.
As Rajeev said, can you send your mesh files (before and after partitioning) ?
Thanks,
Iulian
________________________________
From: moab-dev-bounces at mcs.anl.gov<mailto:moab-dev-bounces at mcs.anl.gov> [moab-dev-bounces at mcs.anl.gov<mailto:moab-dev-bounces at mcs.anl.gov>] on behalf of WARNER, JAMES E. (LARC-D309) [james.e.warner at nasa.gov<mailto:james.e.warner at nasa.gov>]
Sent: Monday, April 27, 2015 3:12 PM
To: moab-dev at mcs.anl.gov<mailto:moab-dev at mcs.anl.gov>
Subject: [MOAB-dev] Polyhedra elements in parallell

Hello,

I generated a MOAB finite element mesh with polyhedra elements that works fine on one processor. However, when I use mbpart to generate two partitions and run my finite element code in parallel, I get the following error when calling DMMoabLoadFromFile() :

"
ParallelComm.cpp:2928: moab::ErrorCode moab::ParallelComm::get_local_handles(moab::EntityHandle*, int, const std::vector<long unsigned int>&): Assertion `ID_FROM_HANDLE(from_vec[i]) < (int) new_ents.size()' failed.
"

Ultimately, I'd like to use MOAB's polyhedra elements to implement arbitrary user-defined finite elements in the code I work on, but as a simple initial test I am trying to create a mesh with just two 3D bricks using MBPOLYHEDRON rather than MBHEX. I first generate 2D quad elements for the six facets of each 3D element and use an array of these facets as the connectivity for each polyhedron element, according to the documentation. I don't get any warning or error messages when generating the input mesh file or when calling mbpart to partition the mesh, only when trying to read it into the finite element code in parallel. And again, it appears to be working correctly in serial.

It's possible I am overlooking something in my code for generating the mesh, so if you can provide any info about the error above and what may be causing it I'd appreciate it!

Thanks a lot,
-Jim




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20150427/b32c79f9/attachment.html>