On Wed, Oct 10, 2012 at 9:45 AM, Thomas Witkowski <span dir="ltr"><<a href="mailto:thomas.witkowski@tu-dresden.de" target="_blank">thomas.witkowski@tu-dresden.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">

  
    
  
  <div text="#000000" bgcolor="#FFFFFF">
    <div>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></div></div></blockquote><div><br></div><div>I think you miss Jed's point. You should NEVER have a MATNEST. You should have a Mat. Its type</div><div>should be determined by MatSetFromOptions(). If you want MATAIJ, you give -mat_type aij, and if you</div>
<div>want MATNEST, you give -mat_type nest.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF"><div>
      Thomas<br>
      <br>
      Am 10.10.2012 15:33, schrieb Jed Brown:<br>
    </div>
    <blockquote type="cite">
      <p>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.</p>
      <div class="gmail_quote">On Oct 10, 2012 8:30 AM, "Thomas
        Witkowski" <<a href="mailto:thomas.witkowski@tu-dresden.de" target="_blank">thomas.witkowski@tu-dresden.de</a>>
        wrote:<br type="attribution">
        <blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
          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>
        </blockquote>
      </div>
    </blockquote>
    <br>
  </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<br>