<div dir="ltr">Thanks for reporting this.  <div><div>There appears to be a bug in MatGetSubMatricesParallel().</div><div>I'll look into it and add a test for the case above.</div><div><br></div><div><br></div><div>In the meantime, maybe you could tell me more about what you are trying</div>

<div>to do? PCGASM: it is rather experimental and doesn't get much use ordinarily. </div><div>Depending on your particular use case, we can either try improving PCGASM, </div></div><div>or find a different way to accomplish the same thing.</div>

<div><br></div><div>Thanks.</div><div>Dmitry.</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Nov 25, 2013 at 9:07 AM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@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">Laurent Berenguer <<a href="mailto:laurent.berenguer@gmail.com">laurent.berenguer@gmail.com</a>> writes:<br>


<br>
>>How exactly did you set this subksp?  Did you modify the code or change<br>
>>the command line arguments?  Did the test run without errors when you<br>
>>don't modify anything?<br>
><br>
> I tried both methods, the error is the same in both cases :<br>
> - I added the option -user_set_subdomain_solvers and modified the code<br>
> from<br>
>   ierr = PCSetType(subpc,PCILU);CHKERRQ(ierr);<br>
> to<br>
> ierr = PCSetType(subpc,PCJACOBI);CHKERRQ(ierr);<br>
> - Or I only added  -sub_ksp_type gmres without changing the code.<br>
><br>
> The example runs correctly when I change nothing. But I'm not sure<br>
> what is the default subksp. When I change nothing in the code but it<br>
> seems to be KSPPREONLY.<br>
<br>
</div>It should say if you run with -ksp_view.<br>
<br>
PCGASM is Dmitry's project and I don't think it has gotten much use yet.<br>
I'll leave it to him to explain what's going on here.  If your cases are<br>
supposed to work, they should be added as tests.<br>
</blockquote></div><br><br clear="all"><div><br></div>
</div></div>