[Bug 1811798] Re: gcc 8 miscompiles scipy/optimize/minpack/qrsolv.f

[Matthias Klose]

** Bug watch added: GCC Bugzilla #88856
   https://gcc.gnu.org/bugzilla/show_bug.cgi?id=88856

** Also affects: scipy via
   https://gcc.gnu.org/bugzilla/show_bug.cgi?id=88856
   Importance: Unknown
       Status: Unknown

** Project changed: scipy => gcc

-- 
You received this bug notification because you are a member of Ubuntu
Foundations Bugs, which is subscribed to gcc-8 in Ubuntu.
https://bugs.launchpad.net/bugs/1811798

Title:
  gcc 8 miscompiles scipy/optimize/minpack/qrsolv.f

Status in gcc:
  Unknown
Status in gcc-8 package in Ubuntu:
  New
Status in imexam package in Ubuntu:
  New
Status in python-scipy package in Ubuntu:
  New
Status in python-scipy package in Debian:
  New

Bug description:
  Originally reported in https://bugs.debian.org/906198

  The following simple test fails on s390x:

  ---------------------------8<------------------------------
  import numpy as np
  from scipy.optimize import leastsq

  y = np.array([0., 1., 1., 2., 1., 1., 0.])
  x = np.arange(len(y))

  def func(par):
      return par[2] * np.exp(-(x - par[0])**2/par[1]) - y

  print(leastsq(func, [0,1,0]))
  ---------------------------8<------------------------------

  The expected result is

  (array([3.        , 4.42280548, 1.67210345]), 1)

  while on S390x, one gets

  (array([0.        , 1.        , 0.34027645]), 3)

  
  Further information from https://bugs.debian.org/915738

  The bug seems to be a miscompilation of scipy/optimize/minpack/qrsolv.f
  with -funroll-loops. Removing -funroll-loops is enough to get the
  leastsq testcase to pass.

To manage notifications about this bug go to:
https://bugs.launchpad.net/gcc/+bug/1811798/+subscriptions