http://gcc.gnu.org/bugzilla/show_bug.cgi?id=50721
--- Comment #8 from Bernd Schmidt <bernds at gcc dot gnu.org> 2011-10-14 13:23:58 UTC --- It looks like gsl_poly_dd_taylor overwrites part of the caller's stack frame. Specifically, this statement w[size - 1] = 1.0; appears to go out of bounds. Can someone who understands Fortran check the source whether it's actually valid?