<div dir="ltr">On Tue, Feb 5, 2013 at 12:01 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
  Ok, one decision made. Dynamic Library loads are collective.<br>
<br>
   Now we need to make several more decisions.<br>
<br>
1) are PETSc package registrations collective? sys, vec, mat, dm, ksp, snes, ts<br>
     currently as Jed noted they are not except with dynamic loading of PETSc libs<br></blockquote><div><br></div><div style>I vote yes. Its much easier, and I see no real upside to being more flexible.</div><div> </div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
2) are all PETSc packages registered upfront during PetscInitialize()?<br>
    currently only with dynamic loading of petsc libs<br></blockquote><div><br></div><div style>Yes again. I think mirroring the dynamic way makes things easier too.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

2a) If all PETSc packages are registered up front what is the mechanism to turn off registering some? Thought Matt disagrees there is always some asshole who says, I don't use TS so I don't want it registered.<br>
</blockquote><div><br></div><div style>Tell them its not loaded and see if they can figure out that it is.</div><div style><br></div><div style>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

…..<br>
<br>
Including where the comms are passed around for all the methods.<br>
<span class="HOEnZb"><font color="#888888"><br>
  Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
On Feb 4, 2013, at 10:53 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
<br>
> On Mon, Feb 4, 2013 at 11:40 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
> On Mon, Feb 4, 2013 at 10:05 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
>    This is currently a mess.<br>
><br>
>    Say one process calls PetscFunctionListAdd() with a function pointer, but another calls it with the string name of the function. Now both processes call PetscFunctionListFind() with a common comm. The process with the function pointer will return immediately with the answer. The one without the function pointer will start mucking around with dynamic libraries which "sometimes" could be collective on the comm so it would block?<br>

><br>
>   These sets of routines evolved organically overtime. We need to refactor the whole hierarchy of these routines and figure out what collectivity is needed.  There are too many potential comms since they were kind of shoved in over time.<br>

><br>
>   It may be simplest if we treat accessing the dynamic libraries as completely non-collective, this means removing things like PetscDLLibraryRetrieve() which, while a way cool concept has never proven to be practical during its 15 years of existence.<br>

><br>
>    So are we able to treat accessing dynamic libraries as completely non-collective? Will this lose a valuable feature?<br>
><br>
> Sort of. The problem is that independent access to the file system is already so slow on current hardware that shared libraries bring those expensive machines to their knees. When we worked on this problem for Python (which is _heavily_ dependent on dynamical loading), we patched glibc-rtld so we could get hooks into a library I called "collfs" (collective file system) that would do a collective open implemented using MPI_Bcast. Most of this circus could go away if libc provided "dlopenfd()", in which case we could use shm_open() and avoid touching the file system at all.<br>

><br>
> From the glibc implementation, I don't think anyone was trying to make adding dlopenfd() easy to implement, so we probably have to deal with paths. Still, if we have a working shm_open and a communicator, we can avoid the libc-rtld hocus pocus with a fast collective load implemented as:<br>

><br>
>   rank 0 mmaps the file<br>
>   everyone else does shm_open and mmap<br>
>   MPI_Bcast<br>
>   dlopen("/dev/shm/thelib.so",)<br>
><br>
> In summary, I think collective loads are useful even without the "retrieval" stuff.<br>
><br>
> I agree with Jed here. We should keep the collective semantics for loading.<br>
><br>
> Now it's cleaner for modularity to load the entire plugin library up-front and let PetscDLLibraryRegister_thelib call MatRegister for everything that it provides. It's easy to manage collectivity this way, but unfortunately, it eats up startup time and memory. (PETSc's current dynamic registration is like Emacs "autoload".)<br>

><br>
> I think this was a hack to begin with (and I did it), so calling the library loads up front does not bother me. I don't<br>
> think anyone in the world is interested an a lightweight, partial PETSc.<br>
><br>
>    Matt<br>
><br>
> At a cost of at least one reduction per library per communicator, we could keep track of the scope on which each library has been loaded so that all loads are safe. Of course performance would go way down if many callers brought in the library on a small object, but that may be unavoidable.<br>

><br>
><br>
>    Barry<br>
><br>
><br>
> On Feb 4, 2013, at 9:22 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
><br>
> > On Sat, Feb 2, 2013 at 3:30 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
> ><br>
> >    Yeah I noticed this problem but didn't want to deal with it when I changed the code.<br>
> ><br>
> > So if we believe the documentation of PetscFunctionListAdd, XXInitializePackage() is effectively collective on COMM_WORLD (though not documented as such). This means that if !defined(PETSC_USE_DYNAMIC_LIBRARIES), the following could deadlock:<br>

> ><br>
> > if (!rank) {<br>
> >   VecCreate(PETSC_COMM_SELF,....);<br>
> > }<br>
> ><br>
> > which would be awfully bad behavior. In reality, PetscFunctionListAdd() does not reference comm at all. Why did you add the comm argument? "Consistency"?<br>
> ><br>
> > Whatever the "next" documentation system is, it should be taught to trace the "collective" attribute and complain if a "Not Collective" function calls a Collective function with an argument other than COMM_SELF.<br>

> ><br>
> ><br>
> >     Yes we should remove the "Formally Collective", I was drinking that week :-)<br>
> ><br>
> >    Barry<br>
> ><br>
> > On Feb 2, 2013, at 2:54 PM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
> ><br>
> > > In [1], PetscFunctionListAdd became implicitly collective on COMM_WORLD, but the all the XXRegisterDynamic() say "Not collective". These all have to be updated if this is the case, but I'm not sure it's even a good thing. What if we have a big multi-domain simulation in which we initialize each of the components on their own subcomm. Those sub-components would not be allowed to register methods (or load plugins) that they might use because registration was implicitly more global.<br>

> > ><br>
> > > The comm is used by PetscLs and others. This is important because file systems are terrible at independent access. (Same for loading shared libraries; it's potentially much easier to do it by broadcasting the library, though portability is tricky.)<br>

> > ><br>
> > > Anyway, it would be really bad to PetscDLLibraryAppend() on a subcomm and have the registration function in the shared lib call PCRegisterDynamic() that promotes itself to COMM_WORLD.<br>
> > ><br>
> > > Maybe we need to pass an explicit comm to all the registration functions.<br>
> > ><br>
> > > [1] <a href="https://bitbucket.org/petsc/petsc-dev/commits/07f9e01e040feeb4162253a60ca63556436f4135" target="_blank">https://bitbucket.org/petsc/petsc-dev/commits/07f9e01e040feeb4162253a60ca63556436f4135</a><br>

> > ><br>
> > > What does "Formally collective" mean anyway? Either it's always safe to call independently, it's "Logically collective" so that there is no performance impact, but it still needs to be collective to have consistent state, or it's Not Collective. This falls under Not Collective because it can deadlock if you call it independently.<br>

> ><br>
> ><br>
><br>
><br>
><br>
><br>
><br>
> --<br>
> What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
> -- Norbert Wiener<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>