<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=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">Begin forwarded message:</div><br class="Apple-interchange-newline"><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">Date: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">December 19, 2014 at 9:44:31 AM CST<br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">From: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">George Bosilca <<a href="mailto:bosilca@icl.utk.edu" class="">bosilca@icl.utk.edu</a>><br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">To: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">Open MPI Users <<a href="mailto:users@open-mpi.org" class="">users@open-mpi.org</a>><br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">Cc: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class="">"<a href="mailto:petsc-maint@mcs.anl.gov" class="">petsc-maint@mcs.anl.gov</a>" <<a href="mailto:petsc-maint@mcs.anl.gov" class="">petsc-maint@mcs.anl.gov</a>><br class=""></span></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px;" class=""><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif; color:rgba(0, 0, 0, 1.0);" class=""><b class="">Subject: </b></span><span style="font-family: -webkit-system-font, Helvetica Neue, Helvetica, sans-serif;" class=""><b class="">Re: [petsc-maint] [OMPI users] Deadlock in OpenMPI 1.8.3 and PETSc 3.4.5</b><br class=""></span></div><br class=""><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class=""><div class="gmail_extra"><div class="gmail_quote">On Fri, Dec 19, 2014 at 8:58 AM, Jeff Squyres (jsquyres) <span dir="ltr" class=""><<a href="mailto:jsquyres@cisco.com" target="_blank" class="">jsquyres@cisco.com</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">George:<br class="">
<br class="">
(I'm not a member of petsc-maint; I have no idea whether my mail will actually go through to that list)<br class="">
<br class="">
TL;DR: I do not think that George's change was correct. PETSC is relying on undefined behavior in the MPI standard and should probably update to use a different scheme.<br class=""></blockquote><div class=""><br class=""></div>I will detail below why the change is correct or at least not worst than what we have today.</div><div class="gmail_quote"><br class=""></div><div class="gmail_quote">Regarding your second point, while I do tend to agree that such issue is better addressed in the MPI Forum, the last attempt to fix this was certainly not a resounding success.</div><div class="gmail_quote"><br class=""></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><br class="">
More detail:<br class="">
<br class="">
More specifically, George's change can lead to inconsistency/incorrectness in the presence of multiple threads simultaneously executing attribute actions on a single entity.<br class=""></blockquote><div class="gmail_quote"><br class=""></div>Indeed, there is a slight window of opportunity for inconsistencies in the recursive behavior. But the inconsistencies were already in the code, especially in the single threaded case. As we never received any complaints related to this topic I did not deemed interesting to address them with my last commit. Moreover, the specific behavior needed by PETSc is available in Open MPI when compiled without thread support, as the only thing that "protects" the attributes is that global mutex.<div class=""><br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">For example, in ompi_attr_delete_all(), it gets the count of all attributes and then loops <count> times to delete each attribute.  But each attribute callback can now insert or delete attributes on that entity.  This can mean that the loop can either fail to delete an attribute (because some attribute callback already deleted it) or fail to delete *all* attributes (because some attribute callback added more).<br class=""></blockquote><div class=""><br class=""></div><div class="">To be extremely precise the deletion part is always correct as it copies the values to be deleted into a temporary array before calling any callbacks (and before releasing the mutex), so we only remove what was in the object attribute hash when the function was called. Don't misunderstand we have an extremely good reason to do it this way, we need to call the callbacks in the order in which they were created (mandated by the MPI standard).</div><div class=""> <br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">ompi_attr_copy_all() has similar problems -- in general, the hash that it is looping over can change underneath it.<br class=""></blockquote><div class=""><br class=""></div><div class=""><div class="">For the copy it is a little bit more tricky, as the calling order is not imposed. Our peculiar implementation of the hash table (with array) makes the code work, with a single (possible minor) exception when the hash table itself is grown between 2 calls. However, as stated before this issue was already present in the code in single threaded cases for years. Addressing it is another 2 line patch, but I leave this exercise to an interested reader.</div></div><div class=""><br class=""></div><div class="">   George.</div><div class=""><br class=""></div><div class=""><br class=""></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">This is, unfortunately, an undefined area of the MPI specification.  I do believe that our previous behavior was *correct* -- it just deadlocks with PETSC because PETSC is relying on undefined behavior.<br class="">
<br class="">
For those who care, Microsoft/Cisco proposed a new attribute system to the Forum a while ago that removes all these kinds of ambiguities (see <a href="http://meetings.mpi-forum.org/secretary/2013/09/slides/jsquyres-attributes-revamp.pdf" target="_blank" class="">http://meetings.mpi-forum.org/secretary/2013/09/slides/jsquyres-attributes-revamp.pdf</a>).  However, we didn't get a huge amount of interest, and therefore lost our window of availability opportunity to be able to advance the proposal.  I'd be more than happy to talk anyone through the proposal if they have interest/cycles in taking it over and advancing it with the Forum.<br class="">
<br class="">
Two additional points from the PDF listed above:<br class="">
<br class="">
- on slide 21, it was decided to no allow the recursive behavior (i.e., you can ignore the "This is under debate" bullet.<br class="">
- the "destroy" callback was not judged to be useful; you can ignore slides 22 and 23.<br class="">
<div class=""><div class="h5"><br class="">
<br class="">
On Dec 17, 2014, at 11:16 PM, George Bosilca <<a href="mailto:bosilca@icl.utk.edu" class="">bosilca@icl.utk.edu</a>> wrote:<br class="">
<br class="">
> Ben,<br class="">
><br class="">
> I can't find anything in the MPI standard suggesting that a recursive behavior of the attribute deletion is enforced/supported by the MPI standard. Thus, the current behavior of Open MPI (a single lock for all attributes), while maybe a little strict, is standard compliant (and thus correct). I guess this is one of these peculiar differences between different MPI libraries that makes other libraries (taking advantage of particular implementation of MPI functions) less portable.<br class="">
><br class="">
> That being said, changing Open MPI to a less restrictive behavior is a 4 lines patch [1]. Unfortunately, this will not make your life easier, you will still have to cope with all the existing versions of Open MPI that do not support this recursive behavior.<br class="">
><br class="">
> I leave the answer to the question of including this patch into the next release to our release managers.<br class="">
><br class="">
>   George.<br class="">
><br class="">
> [1] <a href="https://github.com/open-mpi/ompi/commit/4d55ae838d5eff2818707da7ba60ffd640144360" target="_blank" class="">https://github.com/open-mpi/ompi/commit/4d55ae838d5eff2818707da7ba60ffd640144360</a><br class="">
><br class="">
><br class="">
> On Wed, Dec 17, 2014 at 10:11 PM, Howard Pritchard <<a href="mailto:hppritcha@gmail.com" class="">hppritcha@gmail.com</a>> wrote:<br class="">
> Hi Ben,<br class="">
><br class="">
> Would you mind checking if you still observe this deadlock condition if you use<br class="">
> the 1.8.4 rc4 candidate?<br class="">
><br class="">
> openmpi-1.8.4rc4.tar.gz<br class="">
><br class="">
> I realize the behavior will likely be the same, but this is just to double check.<br class="">
><br class="">
> The Open MPI man page for MPI_Attr_get (hmm.. no MPI_Comm_attr_get man page,<br class="">
> needs to be fixed) says nothing about issues with recursion with respect to invoking<br class="">
> this function within an attribute delete callback, so I would treat this as a bug.<br class="">
><br class="">
> Thanks for your patience,<br class="">
><br class="">
> Howard<br class="">
><br class="">
><br class="">
> 2014-12-17 17:07 GMT-07:00 Ben Menadue <<a href="mailto:ben.menadue@nci.org.au" class="">ben.menadue@nci.org.au</a>>:<br class="">
> Hi PETSc and OpenMPI teams,<br class="">
><br class="">
> I'm running into a deadlock in PETSc 3.4.5 with OpenMPI 1.8.3:<br class="">
><br class="">
>  1. PetscCommDestroy calls MPI_Attr_delete<br class="">
>  2. MPI_Attr_delete acquires a lock<br class="">
>  3. MPI_Attr_delete calls Petsc_DelComm_Outer (through a callback)<br class="">
>  4. Petsc_DelComm_Outer calls MPI_Attr_get<br class="">
>  5. MPI_Attr_get wants to also acquire the lock from step 2.<br class="">
><br class="">
> Looking at the OpenMPI source code, it looks like you can't call an<br class="">
> MPI_Attr_* function from inside the registered deletion callback. The<br class="">
> OpenMPI source code notes that all of the functions acquire a global lock,<br class="">
> which is where the problem is coming from - here are the comments and the<br class="">
> lock definition, in ompi/attribute/attribute.c of OpenMPI 1.8.3:<br class="">
><br class="">
>     404 /*<br class="">
>     405  * We used to have multiple locks for semi-fine-grained locking.<br class="">
> But<br class="">
>     406  * the code got complex, and we had to spend time looking for subtle<br class="">
>     407  * bugs.  Craziness -- MPI attributes are *not* high performance, so<br class="">
>     408  * just use a One Big Lock approach: there is *no* concurrent<br class="">
> access.<br class="">
>     409  * If you have the lock, you can do whatever you want and no data<br class="">
> will<br class="">
>     410  * change/disapear from underneath you.<br class="">
>     411  */<br class="">
>     412 static opal_mutex_t attribute_lock;<br class="">
><br class="">
> To get it to work, I had to modify the definition of this lock to use a<br class="">
> recursive mutex:<br class="">
><br class="">
>     412 static opal_mutex_t attribute_lock = { .m_lock_pthread =<br class="">
> PTHREAD_RECURSIVE_MUTEX_INITIALIZER_NP };<br class="">
><br class="">
> but this is non-portable.<br class="">
><br class="">
> Is the behaviour expected from new versions OpenMPI? In which case a new<br class="">
> approach might be needed in PETSc. Otherwise, maybe a per-attribute lock is<br class="">
> needed in OpenMPI - but I'm not sure whether the get in the callback is on<br class="">
> the same attribute as is being deleted.<br class="">
><br class="">
> Thanks,<br class="">
> Ben<br class="">
><br class="">
> #0  0x00007fd7d5de4264 in __lll_lock_wait () from /lib64/libpthread.so.0<br class="">
> #1  0x00007fd7d5ddf508 in _L_lock_854 () from /lib64/libpthread.so.0<br class="">
> #2  0x00007fd7d5ddf3d7 in pthread_mutex_lock () from /lib64/libpthread.so.0<br class="">
> #3  0x00007fd7d27d91bc in ompi_attr_get_c () from<br class="">
> /apps/openmpi/1.8.3/lib/libmpi.so.1<br class="">
> #4  0x00007fd7d2803f03 in PMPI_Attr_get () from<br class="">
> /apps/openmpi/1.8.3/lib/libmpi.so.1<br class="">
> #5  0x00007fd7d7716006 in Petsc_DelComm_Outer (comm=0x7fd7d2a83b30,<br class="">
> keyval=128, attr_val=0x7fff00a20f00, extra_state=0xffffffffffffffff) at<br class="">
> pinit.c:406<br class="">
> #6  0x00007fd7d27d8cad in ompi_attr_delete_impl () from<br class="">
> /apps/openmpi/1.8.3/lib/libmpi.so.1<br class="">
> #7  0x00007fd7d27d8f2f in ompi_attr_delete () from<br class="">
> /apps/openmpi/1.8.3/lib/libmpi.so.1<br class="">
> #8  0x00007fd7d2803dfc in PMPI_Attr_delete () from<br class="">
> /apps/openmpi/1.8.3/lib/libmpi.so.1<br class="">
> #9  0x00007fd7d78bf5c5 in PetscCommDestroy (comm=0x7fd7d2a83b30) at<br class="">
> tagm.c:256<br class="">
> #10 0x00007fd7d7506f58 in PetscHeaderDestroy_Private (h=0x7fd7d2a83b30) at<br class="">
> inherit.c:114<br class="">
> #11 0x00007fd7d75038a0 in ISDestroy (is=0x7fd7d2a83b30) at index.c:225<br class="">
> #12 0x00007fd7d75029b7 in PCReset_ILU (pc=0x7fd7d2a83b30) at ilu.c:42<br class="">
> #13 0x00007fd7d77a9baa in PCReset (pc=0x7fd7d2a83b30) at precon.c:81<br class="">
> #14 0x00007fd7d77a99ae in PCDestroy (pc=0x7fd7d2a83b30) at precon.c:117<br class="">
> #15 0x00007fd7d7557c1a in KSPDestroy (ksp=0x7fd7d2a83b30) at itfunc.c:788<br class="">
> #16 0x00007fd7d91cdcca in linearSystemPETSc<double>::~linearSystemPETSc<br class="">
> (this=0x7fd7d2a83b30) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Solver/li<br class="">
> nearSystemPETSc.hpp:73<br class="">
> #17 0x00007fd7d8ddb63b in GFaceCompound::parametrize (this=0x7fd7d2a83b30,<br class="">
> step=128, tom=10620672) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Geo/GFace<br class="">
> Compound.cpp:1672<br class="">
> #18 0x00007fd7d8dda0fe in GFaceCompound::parametrize (this=0x7fd7d2a83b30)<br class="">
> at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Geo/GFace<br class="">
> Compound.cpp:916<br class="">
> #19 0x00007fd7d8f98b0e in checkMeshCompound (gf=0x7fd7d2a83b30, edges=...)<br class="">
> at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Mesh/mesh<br class="">
> GFace.cpp:2588<br class="">
> #20 0x00007fd7d8f95c7e in meshGenerator (gf=0xd13020, RECUR_ITER=0,<br class="">
> repairSelfIntersecting1dMesh=true, onlyInitialMesh=false, debug=false,<br class="">
> replacement_edges=0x0)<br class="">
>     at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Mesh/mesh<br class="">
> GFace.cpp:1075<br class="">
> #21 0x00007fd7d8f9a41e in meshGFace::operator() (this=0x7fd7d2a83b30,<br class="">
> gf=0x80, print=false) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Mesh/mesh<br class="">
> GFace.cpp:2562<br class="">
> #22 0x00007fd7d8f8c327 in Mesh2D (m=0x7fd7d2a83b30) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Mesh/Gene<br class="">
> rator.cpp:407<br class="">
> #23 0x00007fd7d8f8ad0b in GenerateMesh (m=0x7fd7d2a83b30, ask=128) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Mesh/Gene<br class="">
> rator.cpp:641<br class="">
> #24 0x00007fd7d8e43126 in GModel::mesh (this=0x7fd7d2a83b30, dimension=128)<br class="">
> at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Geo/GMode<br class="">
> l.cpp:535<br class="">
> #25 0x00007fd7d8c1acd2 in GmshBatch () at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Common/Gm<br class="">
> sh.cpp:240<br class="">
> #26 0x000000000040187a in main (argc=-760726736, argv=0x80) at<br class="">
> /short/z00/bjm900/build/fluidity/intel15-ompi183/gmsh-2.8.5-source/Common/Ma<br class="">
> in.cpp:27<br class="">
><br class="">
><br class="">
> _______________________________________________<br class="">
> users mailing list<br class="">
> <a href="mailto:users@open-mpi.org" class="">users@open-mpi.org</a><br class="">
> Subscription: <a href="http://www.open-mpi.org/mailman/listinfo.cgi/users" target="_blank" class="">http://www.open-mpi.org/mailman/listinfo.cgi/users</a><br class="">
> Link to this post: <a href="http://www.open-mpi.org/community/lists/users/2014/12/26018.php" target="_blank" class="">http://www.open-mpi.org/community/lists/users/2014/12/26018.php</a><br class="">
><br class="">
> _______________________________________________<br class="">
> users mailing list<br class="">
> <a href="mailto:users@open-mpi.org" class="">users@open-mpi.org</a><br class="">
> Subscription: <a href="http://www.open-mpi.org/mailman/listinfo.cgi/users" target="_blank" class="">http://www.open-mpi.org/mailman/listinfo.cgi/users</a><br class="">
> Link to this post: <a href="http://www.open-mpi.org/community/lists/users/2014/12/26025.php" target="_blank" class="">http://www.open-mpi.org/community/lists/users/2014/12/26025.php</a><br class="">
><br class="">
> _______________________________________________<br class="">
> users mailing list<br class="">
> <a href="mailto:users@open-mpi.org" class="">users@open-mpi.org</a><br class="">
> Subscription: <a href="http://www.open-mpi.org/mailman/listinfo.cgi/users" target="_blank" class="">http://www.open-mpi.org/mailman/listinfo.cgi/users</a><br class="">
</div></div>> Link to this post: <a href="http://www.open-mpi.org/community/lists/users/2014/12/26028.php" target="_blank" class="">http://www.open-mpi.org/community/lists/users/2014/12/26028.php</a><br class="">
<span class=""><br class="">
<br class="">
--<br class="">
Jeff Squyres<br class="">
<a href="mailto:jsquyres@cisco.com" class="">jsquyres@cisco.com</a><br class="">
For corporate legal information go to: <a href="http://www.cisco.com/web/about/doing_business/legal/cri/" target="_blank" class="">http://www.cisco.com/web/about/doing_business/legal/cri/</a><br class="">
<br class="">
</span><span class="">_______________________________________________<br class="">
users mailing list<br class="">
<a href="mailto:users@open-mpi.org" class="">users@open-mpi.org</a><br class="">
Subscription: <a href="http://www.open-mpi.org/mailman/listinfo.cgi/users" target="_blank" class="">http://www.open-mpi.org/mailman/listinfo.cgi/users</a><br class="">
</span>Link to this post: <a href="http://www.open-mpi.org/community/lists/users/2014/12/26041.php" target="_blank" class="">http://www.open-mpi.org/community/lists/users/2014/12/26041.php</a><br class="">
</blockquote></div><br class=""></div></div>
</div></blockquote></div><br class=""></body></html>