http://gcc.gnu.org/bugzilla/show_bug.cgi?id=52063
--- Comment #4 from Aravind <aravindvijayan224185 at gmail dot com> 2012-01-31 03:38:01 UTC --- Narrowed down to: program x09f95 use plf95demolib integer, parameter :: nptsx=35, xdim=99 real(kind=8) :: xg1(xdim), distort xg1 = foo (arange(0,nptsx) / dble(nptsx-1), distort) contains elemental real(kind=8) function foo( coord, factor ) real(kind=8), intent(in) :: coord, factor foo = 1 end function foo end program x09f95