------- Comment #1 from dominiq at lps dot ens dot fr 2009-09-19 14:01 ------- Have a look at http://linux.die.net/man/l/dgetrf for the calling convention. The following works for me:
program det_test implicit none integer, parameter :: p15=selected_real_kind(15) integer :: ipiv(2) real(p15) :: A(2,2) integer :: info A(1,1) = 1.0 A(1,2) = 0.0 A(2,1) = 3.0 A(2,2) = 0.5 info = 1 call DGETRF(2,2,A,2,ipiv,info) write(*,*) A(1,1) write(*,*) A(1,2) write(*,*) A(2,1) write(*,*) A(2,2) end program det_test [ibook-dhum] bug/timing% gfc pr41408_db.f90 -Wl,-framework -Wl,vecLib [ibook-dhum] bug/timing% a.out 3.0000000000000000 0.50000000000000000 0.33333333333333331 -0.16666666666666666 -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=41408