Implementing an ADI method in parallel

Matthew Knepley knepley at gmail.com
Fri Aug 24 11:21:33 CDT 2007


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
>
>


-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener




More information about the petsc-users mailing list