bug #59710 (Undefined behavior error) marked fixed but remains in the recent GSL

[Zhoulai Fu]

Dear GSL developers,

The following undefined behavior bug (which I reported) was marked "fixed".
But as I tested today, the bug remained there in the latest GSL. Can you
please double check?
*bug #59710 <https://savannah.gnu.org/bugs/?59710>*: Undefined behavior
error in Nonlinear Least-Squares Fitting
(https://savannah.gnu.org/bugs/?59710)

To test if the bug is well fixed, you can use the nlfit_debug.c provided in
my previous bug report (attached for your convenience in this email). Below
is the error message I got as of today.

% ./a.out 1>/dev/null
iter  0: A = 1.0000, lambda = 1.0000, b = 0.0000, cond(J) =      inf,
|f(x)| = 101.0200
iter  1: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =  92.8216,
|f(x)| = -nan
iter  2: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  3: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  4: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
*nielsen.c:98:7: runtime error: left shift of 4611686018427387904 by 1
places cannot be represented in type 'long int'*
iter  5: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  6: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  7: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  8: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter  9: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
iter 10: A = 3.5110, lambda = -12.8820, b = 1.2364, cond(J) =      nan,
|f(x)| = -nan
.....


BR,
Zhoulai

nlfit_bug.c
Description: Text Data