openblas openmp fortran multi- treaded
===

I have a code using fortran f95 and the gfortran compiler. I am also using openMp and i have to handle very big arrays. In my code i also have to solve a system of linear equations using the solver DGTSV from openblas. I want to parallelize this solver as well using openblas which should be capable of that. But I have trouble with the syntax. Using the attached pseudo code all 4 cpus are used to almost 100% but I am not sure if each kernel solves the linear equations separatly or if they split it into parts and calculating it parallel.
The whole stuff is compiled using gfortran -fopenmp -lblas a.f95 -o a.out

So my pseudo code looks like

    program a
    implicit none
    integer, parameter	::	N	=		200
    real*8, dimension(numx)	::	D	=		0.0
    real*8, dimension(numx-1)::	DL	=		0.0
    real*8, dimension(numx-1)::	DU	=		0.0
    real*8, dimension(numx)	::	b	=		0.0
    integer			::	info	=		0
    integer :: numthread=4
    ...
    !$OMP PARALLEL NUM_THREADS(numthread)
    ...
    !$OMP DO
    ...
    !$OMP END DO
    CALL DGTSV(N,1,DL,D,DU,b,N,info)
    !$OMP DO
    ...
    !$OMP END DO
    ...
    !$OMP END PARALLEL
    end program a

What does i have to do to make the solver parallelized, so each kernel calculates parts of the solver ? Thank you very much for your help!