<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 11, 2014 at 11:47 AM, Eric Chamberland <span dir="ltr"><<a href="mailto:Eric.Chamberland@giref.ulaval.ca" target="_blank">Eric.Chamberland@giref.ulaval.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">On 09/11/2014 12:26 PM, Matthew Knepley wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I am trying to understand the problem. This matrix is created by by the<br>
MatGetSubMatrix() call.<br>
</blockquote>
<br>
If this is right, then I misunderstood the effect of MatGetSubMatrix on a MatNest.  We created a MatNest which holds the a_11 matrix...<br></blockquote><div><br></div><div>MatNest is absolutely the worst thing in the PETSc interface. You should never ever ever ever</div><div>be calling MatNest directly. You should be assembling into one matrix from views. Then MatNest</div><div>can be used for optimization in the background.</div><div><br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
   1) What options are you trying to pass to it?<br>
</blockquote>
<br>
mat_type sbaij<br>
mat_block_size  3 (not in this example however)<br></blockquote><div><br></div><div>These are reasonable, but they really apply at creation time (since you would not want</div><div>to convert after values have been set), and it sounds like that is what you are doing.</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>
   2) Why would they interfere with the KSP options?<br>
</blockquote>
<br>
I didn't mentioned that they could interfere, just that the name was changed and that it was "curious" from my point of view.<br></blockquote><div><br></div><div>Okay, so the name change is strange, and happens because MatNest returns a reference</div><div>to the inner matrix rather than some view which gets created and destroyed. Let me talk</div><div>to Jed.</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">
But if the sub-matrix extracted from the matnest is not the a_11 matrix we gave to the matnest, then maybe we did something wrong?  Or you can't verify that the extracted matrix did exists before the extraction?<br>
<br>
Again, the options we gave are treated correctly, since we create the matrix and give it to the nest.  The "curious" thing is when petsc retreive (the same?) matrix and changes its name...<br>
<br>
Thanks again!<span class="HOEnZb"><font color="#888888"><br>
<br>
Eric<br>
<br>
</font></span></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>