Slepc is not thread safe?
Hi everyone.
I'm using OpenMP to parallelize the for-loop in my program, which means one thread should solve an independent eigenproblem at a time. By reading the manual page of Petsc, I realize that Petsc support such "embarrassing" parallelization when configuring with options --with-threadsafety --with-log=0 --with-openmp. I wrote a piece of test code in which several KSP problems are solved in a for-loop and it worked well. But when I replaced the KSP problem with the EPS problem in the test code, I received the seg fault.
I'm wondering whether Slepc is not thread-safe even when it is built based thread-safe version of Petsc?
Thanks.