[petsc-dev] OpenMP in PETSc when calling from Fortran?

[Barry Smith]

   1) Work with petsc-dev http://www.mcs.anl.gov/petsc/developers/index.html  the threaded stuff is more mature.

    2) You should not need petscthreadcomm.h in Fortran. Simply using OpenMP progmas in your portion of the code.

    3) See http://www.mcs.anl.gov/petsc/features/threads.html for a little bit of detail and feel free to post additional questions on this mailing list.

   Barry


On Feb 28, 2013, at 4:35 AM, Åsmund Ervik <asmund.ervik at ntnu.no> wrote:

> Hi,
> 
> Having read that petsc-dev has the possibility of using OpenMP to run PC and KSP routines in parallel, and that people are getting nice speedup with this, I wanted to make use of this feature in our code. It is an incompressible two-phase flow solver programmed in Fortran which used PETSc for solving the pressure Poisson equation.
> 
> I am, however, unable to use PETSC with OpenMP from Fortran. In particular, I can't seem to find a Fortran header file corresponding to petscthreadcomm.h (which is for C/C++). Is this not implemented yet? If so, is only the header file missing, or is more work required?
> 
> Regards,
> Åsmund Ervik