<br><br><div class="gmail_quote">On Wed, Oct 10, 2012 at 11:38 AM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

<div class="im"><br>
On Oct 10, 2012, at 8:50 AM, Jed Brown <<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>> wrote:<br>
<br>
> How did you get the MatNest?<br>
><br>
> What you describe is the right approach, but proper preallocation is messy in the general case. There is a partial implementation in petsc-dev, but it needs to be finished/debugged.<br>
<br>
</div>   Thomas,<br>
<br>
      We should have a reasonable MatConvert_MatNest_MATAIJ. Perhaps you can look at what Jed has started in petsc-dev and improve it enough to suit your purpose.<br></blockquote><div>There is this commit from last summer: <a href="http://petsc.cs.iit.edu/petsc/petsc-dev/rev/5289fd118d48">http://petsc.cs.iit.edu/petsc/petsc-dev/rev/5289fd118d48</a></div>

<div>It's untested, but I pushed it back then hoping someone could beat it into shape.</div><div>I also have a patch that eliminates ISAllGather() in favor of PetscSF operations, but that's even</div><div>less tested (if such a thing is possible).  I can clean it up and push, if somebody wants to play with it.</div>

<div><br></div><div>Dmitry.</div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<span class="HOEnZb"><font color="#888888"><br>
    Barry<br>
</font></span><div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
><br>
> On Oct 10, 2012 8:45 AM, "Thomas Witkowski" <<a href="mailto:thomas.witkowski@tu-dresden.de">thomas.witkowski@tu-dresden.de</a>> wrote:<br>
> You are right. But what if a have a MATNEST and i must convert it to MATAIJ. Any better way than what I wrote in my first mail?<br>
><br>
> Thomas<br>
><br>
> Am 10.10.2012 15:33, schrieb Jed Brown:<br>
>> The best way, by far, is to assemble the matrix as AIJ if that is what you want to use. See MatGetLocalSubMatrix and SNES ex28. Using MATNEST should he a run-time choice.<br>
>><br>
>> On Oct 10, 2012 8:30 AM, "Thomas Witkowski" <<a href="mailto:thomas.witkowski@tu-dresden.de">thomas.witkowski@tu-dresden.de</a>> wrote:<br>
>> Is there a fast way to convert a MatNest to a MatAIJ? I tried to use MatConvert, but it seems that this convert is not yet implemented. As I need this functionality to factorize the nested matrix with MUMPS, is there any better way than to extract the nested matrices, iterate over all rows of these matrices and use MatSetValues to add them again to the final MatAIJ?<br>


>><br>
>> Thomas<br>
><br>
<br>
</div></div></blockquote></div><br>