<div class="gmail_quote">On Thu, Jan 26, 2012 at 04:04, Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com">stefano.zampini@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div id=":4pk">I'm partially refactoring BDDC code I contributed to petsc-dev, trying to limit the extent of hand-written code and to exploit already existing PETSc code. I want to know if there are DM objects in PETSc suitable to store informations related to the local to global map of MATIS objects which I can use to manage communications between different levels (through DMCoarsen, DMGetInterpolation and DMGetMatrix). DMMesh can do the job?<br>
</div></blockquote></div><br><div>Sorry for the slow reply. The question, I think, is how to produce a Galerkin coarse grid operator of type MATIS. I think we can write MatPtAP_IS without too much trouble, although I don't know how it would perform.</div>
<div><br></div><div>I don't see how there is enough information in a MATIS to define coarsening and interpolation algebraically. But maybe you are thinking of making PCBDDC algebraically construct a grid hierarchy? I would not be inclined to put that into DM, instead, I would call PCMG functions directly.</div>