<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">If you need edges, create them before ghosting; I modified your example, and it now works better<br>
<br>
If edges get created after ghosting, they are not shared correctly between processors, some are duplicated in the file.
<br>
<br>
Or you can try merging them.<br>
<br>
Still, I would delete them if you do not need them.<br>
I agree it is a bug that you cannot read that file with those duplicated edges. <br>
<br>
Iulian<br>
<br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF397506"><font size="2" color="#000000" face="Tahoma"><b>From:</b> kanaev@ibrae.ac.ru [kanaev@ibrae.ac.ru]<br>
<b>Sent:</b> Friday, April 18, 2014 2:42 PM<br>
<b>To:</b> Grindeanu, Iulian R.<br>
<b>Subject:</b> RE: exchanging both quads and edges while exchanging ghosts cells<br>
</font><br>
</div>
<div></div>
<div>
<p>I actually don't need them so writing hexas+quads should work for me </p>
<p>Here's a small program that reads .h5m file in parallel (you should give the name in command line) then gets edges adjacent to every cell and writes the mesh back to ptest.h5m
</p>
<p>You can see you cannot read ptest.h5m then with the same program because of the odd edges
</p>
<p>Anton </p>
<p><br>
On Thu, 17 Apr 2014 22:10:22 +0000, "Grindeanu, Iulian R." wrote: </p>
<blockquote style="padding-left:5px; border-left-color:#1010ff; border-left-width:2px; border-left-style:solid; margin-left:5px; width:100%">
P {margin-top:0;margin-bottom:0;}
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">do you need edges or not?<br>
what is the bug? Are edges duplicated when you ghost? can you replicate a problem in a small moab test?<br>
<br>
Can you delete edges before writing the file?<br>
Or when you write, write only the hexas + quads (put in a set what you need, and write in parallel just that set)<br>
<br>
Why can't you read in parallel after writing?<br>
<br>
<div style="font-family:'Times New Roman'; color:#000000; font-size:16px">
<hr>
<div id="divRpF885029" style="direction:ltr"><font size="2" color="#000000" face="Tahoma"><strong>From:</strong> kanaev@ibrae.ac.ru [kanaev@ibrae.ac.ru]<br>
<strong>Sent:</strong> Thursday, April 17, 2014 4:18 PM<br>
<strong>To:</strong> Grindeanu, Iulian R.<br>
<strong>Cc:</strong> MOAB dev<br>
<strong>Subject:</strong> exchanging both quads and edges while exchanging ghosts cells<br>
</font><br>
</div>
<div>
<p>Hello Iulian, </p>
<p>I exchange one layer of ghost dim 3 </p>
<p> result = pcomm->exchange_ghost_cells(3, // int ghost_dim, </p>
<p> 0, // int bridge_dim, </p>
<p> 1, //int num_layers, </p>
<p> 2, //int addl_ents, </p>
<p> true); // bool store_remote_handles); </p>
<p> </p>
<p>then if i do: </p>
<p> </p>
<p> result = mb->get_entities_by_dimension(0, 3, ents); </p>
<p> PRINT_LAST_ERROR; </p>
<p> </p>
<p> for (moab::Range::iterator it = ents.begin(); it != ents.end(); it++) { </p>
<p> EntityHandle ent = *it; </p>
<p> HandleVec edges; </p>
<p> result = mb->get_adjacencies(&(ent), 1, 1, true, edges); </p>
<p> } </p>
<p> </p>
<p>i'll have additional edges created which i don't need (because after i save file writing everything i cannot load from it in parallel afterwards)
</p>
<p>I think i have to echange edges besides of quads during exchange_ghost_cells but it's not an option.
</p>
<p>How can i avoid unnecessary doubling of edges in this case? </p>
<p> </p>
<p>Thanks </p>
<p>Anton </p>
</div>
</div>
</div>
</blockquote>
<p> </p>
</div>
</div>
</div>
</body>
</html>