[Bug fortran/41408] Segmentation fault calling DGETRF from gfortran

[dominiq at lps dot ens dot fr]

------- 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