<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Thanks Vijay,</div>
<div class=""><br class="">
</div>
<div class="">Only to report, deleting and creating new pcomm, solved the problem.</div>
<div class=""><br class="">
</div>
Lukasz
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On 27 Nov 2017, at 22:15, Lukasz Kaczmarczyk <<a href="mailto:Lukasz.Kaczmarczyk@glasgow.ac.uk" class="">Lukasz.Kaczmarczyk@glasgow.ac.uk</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Hello,</div>
<div class=""><br class="">
See attached figure. You can see red surface with is imprinted on the mesh. Surface is imprinted by edge splitting. With this, I create good quality mesh with node merging and flips. The surface is evolving and I solve PDE on it.</div>
<div class=""><br class="">
What would you suggest, to create new pcomm on cloned mpi comm? <br class="">
<br class="">
Regards,<br class="">
Lukasz</div>
<div class=""><br class="">
</div>
<br class="">
<blockquote type="cite" class="">On 27 Nov 2017, at 22:06, Vijay S. Mahadevan <<a href="mailto:vijay.m@gmail.com" class="">vijay.m@gmail.com</a>> wrote:<br class="">
<br class="">
You could try calling the remove_pcomm to delete it from the global<br class="">
tag/internal set. Now calling add_pcomm should return a ParallelComm<br class="">
object with the same ID, potentially. This should theoretically help<br class="">
what you want to do, but I still don't entirely see the use of<br class="">
recreating the PComm object on the same MPI_Comm.<br class="">
<br class="">
Vijay<br class="">
<br class="">
<br class="">
On Mon, Nov 27, 2017 at 4:58 PM, Lukasz Kaczmarczyk<br class="">
<<a href="mailto:Lukasz.Kaczmarczyk@glasgow.ac.uk" class="">Lukasz.Kaczmarczyk@glasgow.ac.uk</a>> wrote:<br class="">
<blockquote type="cite" class="">Hello,<br class="">
<br class="">
Thanks for fast response. Yes, changes are in the original file set. On some subset of entities and new create new entities. I do mesh cutting with the arbitrary surface.<br class="">
<br class="">
Currently I resolve ents on that meshset, but before I reset tags, pstatus_tag(), sharedp_tag(), sharedps_tag(), sharedh_tag(), sharedhs_tag(), and in principle it works. I do that on the same pcomm.<br class="">
<br class="">
I create pcomm like this<br class="">
<br class="">
ParallelComm *pcomm =<br class="">
ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX, &moab_comm_world);<br class="">
<br class="">
Is it a way to delete this and create a new pcomm on the same index?<br class="">
<br class="">
Regards,<br class="">
Lukasz<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<br class="">
<blockquote type="cite" class="">On 27 Nov 2017, at 21:38, Vijay S. Mahadevan <<a href="mailto:vijay.m@gmail.com" class="">vijay.m@gmail.com</a>> wrote:<br class="">
<br class="">
Lukasz,<br class="">
<br class="">
Are the topological changes still part of the original file_set ? If<br class="">
yes, have you tried to just call resolve_shared_ents on that file set<br class="">
again ? Does this not solve the issue ? There is no global way to<br class="">
reset all the ParallelComm tags currently I think, but you could<br class="">
potentially create a new ParallelComm object spanning the same<br class="">
communicator as before and redo your operations on the older or a new<br class="">
file set containing the entity changes.<br class="">
<br class="">
Also, when you say repartition, do you mean dynamically load balance<br class="">
and redistribute entities ? If not, please explain.<br class="">
<br class="">
If neither of these suggestions are helpful, it would be useful if you<br class="">
can provide an example use case for us to better understand your<br class="">
workflow.<br class="">
<br class="">
Thanks,<br class="">
Vijay<br class="">
<br class="">
On Mon, Nov 27, 2017 at 4:02 PM, Lukasz Kaczmarczyk<br class="">
<<a href="mailto:Lukasz.Kaczmarczyk@glasgow.ac.uk" class="">Lukasz.Kaczmarczyk@glasgow.ac.uk</a>> wrote:<br class="">
<blockquote type="cite" class="">Hello,<br class="">
<br class="">
I do intense topological changes and need to keep precisely the same mesh an all procs, at least is the simplest solution. I partition mesh and resolve shared entities. All calculations are distributed, and I build approx. fields only on partitions. It works
great.<br class="">
<br class="">
However, I am not clear what is the best way to reset whole ParallelComm and ParallelComm tags, etc. and once again resolve shared entities. I repartition problem to balance mesh after some topological changes on the mesh. Can you give me any guidance and advice?<br class="">
<br class="">
Kind regards,<br class="">
Lukasz<br class="">
<br class="">
<br class="">
[University of Glasgow: The Times Scottish University of the Year 2018]<br class="">
</blockquote>
</blockquote>
<br class="">
</blockquote>
</blockquote>
<span id="cid:B87BD9EF-5557-49C2-806C-C0BFC0B63522"><cut_sylinedr_with_holes.png></span></div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>