Hello,<div><br></div><div>In reply,</div><div><br></div><div>I get a Segment fault even though I only call:</div><div>KSP;</div><div>KSPCreate(MPI_COMM_SELF, &ksp); </div><div><br></div><div>Again... note that this is being done by all cores in the MPI_COMM_WORLD at the same time. <br>
<br><div class="gmail_quote">On Mon, Feb 28, 2011 at 11:11 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov">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 Feb 28, 2011, at 10:07 PM, Joshua Booth wrote:<br>
<br>
> Hello,<br>
><br>
> I was wondering if it is possible to have multiple KSP running, each on their own core but in the same program.<br>
<br>
</div> Yes. If you want a separate linear solve on each MPI process then you use MPI_COMM_SELF for the KSP, note that you also need to use that same MPI_COMM_SELF for creating the Vecs and Mats since each process needs its own for its particular linear system.<br>
<div class="im"><br>
> When I have tried this before using MPI_COMM_SELF, I get an error.<br>
<br>
</div> We would need more information to determine what has gone wrong.<br>
<br>
Barry<br>
<br>
<br>
><br>
> Thank you<br>
><br>
> Joshua Booth<br>
<br>
</blockquote></div><br></div>