[petsc-users] Multiple Copies of KSP

[Satish Balay]

It should work. you'll have to use gdb to determine exact location and
reason for this crash in your code.

satish

On Mon, 28 Feb 2011, Joshua Booth wrote:

> Hello,
> 
> In reply,
> 
> I get a Segment fault even though I only call:
> KSP;
> KSPCreate(MPI_COMM_SELF, &ksp);
> 
> Again... note that this is being done by all cores in the MPI_COMM_WORLD at
> the same time.
> 
> On Mon, Feb 28, 2011 at 11:11 PM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
> >
> > On Feb 28, 2011, at 10:07 PM, Joshua Booth wrote:
> >
> > > Hello,
> > >
> > > I was wondering if it is possible to have multiple KSP running, each on
> > their own core but in the same program.
> >
> >   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.
> >
> > > When I have tried this before using MPI_COMM_SELF,  I get an error.
> >
> >    We would need more information to determine what has gone wrong.
> >
> >   Barry
> >
> >
> > >
> > > Thank you
> > >
> > > Joshua Booth
> >
> >
>