Implementing an ADI method in parallel

Satish Balay balay at mcs.anl.gov
Fri Aug 24 11:40:27 CDT 2007


Or for  uniproc solves - you can create as many solver objects
as you need with MPI_COMM_SELF or PETSC_COMM_SELF, on each proc

Satish

On Fri, 24 Aug 2007, Matthew Knepley wrote:

> You can easily setup concurrent solves by splitting the communicator. For
> example, if you divide the communicator into two groups using MPI_Comm_split()
> and then create two KSP objects, each with a subcommunicator obtained from
> the split, Then these two solves can run concurrently.
> 
>   Thanks,
> 
>     Matt
> 
> On 8/24/07, David Knezevic <david.knezevic at balliol.ox.ac.uk> wrote:
> > Hello,
> >
> > I was hoping to get some advice on whether the following is possible in
> > PETSc: I want to implement an ADI-like method, in which I want to do
> > multiple uniprocessor linear solves concurrently (each linear solve is
> > quite small, but there are a lot of them).
> >
> > This isn't what PETSc is designed for I guess, since one would typically
> > do one large solve with multiple processors, but I was wondering if
> > there is a good way of implementing this kind of ADI-like method in
> > PETSc? Or should I instead implement this kind of thing using a serial
> > linear solver, and control the concurrent solves with MPI? I'd prefer to
> > stick to PETSc if possible, since it has so much nice functionality
> > bundled in it, and has a lot of the MPI stuff already taken care of...
> >
> > Thanks very much for the help.
> >
> > Regards,
> > David Knezevic
> >
> >
> 
> 
> 




More information about the petsc-users mailing list