Hi
I am testing the following code with composer_xe_2015
lwork = 256*n
call ZSYTRF( 'U', n, afull, n, ipiv, work, lwork, error )
call ZSYTRI( 'U', n, afull, n, ipiv, work, lwork, error )
using
gfortran -O3 -fopenmp read_blas.f90 -Wl,--start-group ${MKLROOT}/lib/intel64/libmkl_gf_lp64.a ${MKLROOT}/lib/intel64/libmkl_core.a ${MKLROOT}/lib/intel64/libmkl_gnu_thread.a -Wl,--end-group -ldl -lpthread -lm
After setting OMP_NUM_THREADS=16 I can see in "top" that ZSYTRF runs in parallel but ZSYTRI not. The matrix is large enough (n=15120)
Any idea how to solve this problem?
Olaf