complex system GETRF+GETRS

diedro · ‎12-09-2014

Dear all,

I would like to solve a linear complex system with MKL libraries. As I have done with real system I use GETRF with GETRS. The MKL reference says that I can also use getrs also for comple system. Here my example code:

program testmkl
use LAPACK95
implicit none
complex    ,allocatable,dimension(:,:)::AA
complex    ,allocatable,dimension(:)  ::BB
integer    ,allocatable,dimension(:)::IPV
integer :: info,n

n=10
allocate(AA(n,n))
allocate(BB(n))
allocate(IPV(n))

 call GETRF(AA,IPV,info)
 call GETRS(AA,IPV,BB,info)
 
endprogram

However, I am not able to compile it. This is my error:

There is no matching specific subroutine for this generic subroutine call.   [GETRS]
 call GETRS(AA,IPV,BB,info)

Where that am i wrong?

Thanks

diedro · ‎12-09-2014

Dear all,

I am sorry. I get the error. Here, the correct code

program testmkl
use LAPACK95
implicit none
complex    ,allocatable,dimension(:,:)::AA
complex    ,allocatable,dimension(:)  ::BB
integer    ,allocatable,dimension(:)::IPV
integer :: info,n

n=10
allocate(AA(n,n))
allocate(BB(n))
allocate(IPV(n))

 call GETRF(AA,IPV)
 call GETRS(AA,IPV,BB,'N',info)
 
endprogram