<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:13px">Can you send your mesh file?<br><div id="yui_3_16_0_1_1429984829891_153878"><span></span></div><div id="yui_3_16_0_1_1429984829891_153916" class="signature">Rajeev Jain <br>630-252-3176 / 630-252-5986 (fax)<br>jain@mcs.anl.gov</div><br>  <div id="yui_3_16_0_1_1429984829891_153881" style="font-family: HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 13px;"> <div id="yui_3_16_0_1_1429984829891_153880" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div id="yui_3_16_0_1_1429984829891_153879" dir="ltr"> <hr size="1">  <font id="yui_3_16_0_1_1429984829891_153917" face="Arial" size="2"> <b><span style="font-weight:bold;">From:</span></b> "WARNER, JAMES E. (LARC-D309)" <james.e.warner@nasa.gov><br> <b><span style="font-weight: bold;">To:</span></b> "moab-dev@mcs.anl.gov" <moab-dev@mcs.anl.gov> <br> <b><span style="font-weight: bold;">Sent:</span></b> Monday, April 27, 2015 3:12 PM<br> <b><span style="font-weight: bold;">Subject:</span></b> [MOAB-dev] Polyhedra elements in parallell<br> </font> </div> <div id="yui_3_16_0_1_1429984829891_153882" class="y_msg_container"><br><div id="yiv2220346356">

 

<div id="yui_3_16_0_1_1429984829891_153884">
<div id="yui_3_16_0_1_1429984829891_153913">
<div id="yui_3_16_0_1_1429984829891_153912">
<div id="yui_3_16_0_1_1429984829891_153911">Hello,</div>
</div>
</div>
<div><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153883">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()
 :</div>
<div><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153910">"</div>
<div id="yui_3_16_0_1_1429984829891_153909">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.</div>
<div id="yui_3_16_0_1_1429984829891_153885">"</div>
<div id="yui_3_16_0_1_1429984829891_153908"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153886">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. </div>
<div id="yui_3_16_0_1_1429984829891_153890"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153889">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!</div>
<div id="yui_3_16_0_1_1429984829891_153887"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153888">Thanks a lot,</div>
<div id="yui_3_16_0_1_1429984829891_153901">-Jim </div>
<div id="yui_3_16_0_1_1429984829891_153918"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153919"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153920"><br>
</div>
<div id="yui_3_16_0_1_1429984829891_153921"><br>
</div>
</div>
</div><br><br></div> </div> </div>  </div></body></html>