[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Parallel use of Prometheus

I have not tried this with any of the examples in the distribution only with my own program.

The steps are relatively straightforward but you must know them already since you have it
working in serial. Nonetheless here is the order of call I have:

call MatCreateMPIBAIJ(...                 ! must be blocked
call MatAssemblyBegin (...                ! Assemble the matrix
call MatAssemblyEnd   (...

call VecAssemblyBegin (...                !Assemble the RHS
call VecAssemblyEnd   (...

call KSPCreate (...
call KSPSetOperators (...
call KSPGetPC (...
call PCSetCoordinates(... ! I'm not fully certain but your problem may have to be defined over R^3 for this to work well
call KSPSetFromOptions(...

call KSPSolve         (...


Nicolas Tardieu wrote: 
Hi Sanjay,

Have you used it on simple examples of the PETSc distribution?

Nicolas

2008/1/16, Sanjay Govindjee <sanjay@xxxxxxxxxxxxxxx <mailto:sanjay@xxxxxxxxxxxxxxx>>:

    We have used it extensively with good success in the past.
    -sanjay

    Nicolas Tardieu wrote:
    > Hi everyone,
    >
    > I have been using Prometheus for scalar and vector problems with
    great
    > success on sequential computers.
    > But when trying to use it in parallel, Prometheus crashes.
    > Then I tried it on some PETSc tests
    (src/ksp/ksp/examples/tests/ex3.c
    > for instance) on parallel and they also crash...
    > Does anyone have experienced Prometheus in parallel?
    >
    > Thanks,
    >
    > Nicolas
    >