Problems using CGEQRF

Problems using CGEQRF

I am trying to use Lapack's cgeqrf routine but there is something wrong:

!---------------------------------------------------------------

program main

implicit none

integer,parameter :: p = 5
complex :: tau(p), u(p,p), work(p)

u = cmplx(1.0,0.0)
tau = cmplx(0.0,0.0)
work = cmplx(0.0,0.0)

call cgeqrf(p, p, u, p, tau, work, p)

end program main

!---------------------------------------------------------------

Then, when trying to run:

† perjanus † Area51 $ ifort -L/opt/intel/composerxe/mkl/lib/intel64/ -traceback -mkl ./main.f90 -o main.exe
† perjanus † Area51 $ ./main.exe
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
libmkl_intel_lp64 00007F92E6F46A64 Unknown Unknown Unknown
main.exe 0000000000403880 MAIN__ 13 main.f90
main.exe 000000000040371C Unknown Unknown Unknown
libc.so.6 00007F92E3F334BD Unknown Unknown Unknown
main.exe 0000000000403619 Unknown Unknown Unknown
† perjanus † Area51 $

Am I doing something wrong? I am using ifort 12.1.5 20120612 (don't know how to find mkl version, but it came in the same bundle with ifort and icc).

Thanks in advance

2 Beiträge / 0 neu
Letzter Beitrag
Nähere Informationen zur Compiler-Optimierung finden Sie in unserem Optimierungshinweis.

How did you link this example? Please check with the Linker Adviser how to do that. see the link to the adviser follow this link: http://software.intel.com/sites/products/mkl/

Kommentar hinterlassen

Bitte anmelden, um einen Kommentar hinzuzufügen. Sie sind noch nicht Mitglied? Jetzt teilnehmen