[petsc-users] PETSC with OPENMP

Enrico degregori at dkrz.de
Fri Nov 6 03:41:02 CST 2020


Hello,

I need to call PETSC functions within an OPENMP parallel region. The 
reason is not related to PETSC itself but to parallel IO. To finalize 
the parallel IO I need to launch an OPENMP task. Thus, all the PETSC 
calls are in a OPENMP single region and executed by one thread while 
another thread finalizes the IO.

This seems to be pretty standard even if I'm using OPENMP with PETSC but 
when I try to run, the code crashes the first time I call a PETSC 
function while the other thread is working on the IO.

I've installed PETSC with threadsave and openmp mode and with log and 
debugging set to 0 so this should be fine.

Do you know about a branch that support this kind of issue or a way to 
solve it?

Enrico


More information about the petsc-users mailing list