Your message dated Tue, 19 May 2020 15:04:19 +0000
with message-id <[email protected]>
and subject line Bug#953437: fixed in flint 2.5.2-22
has caused the Debian Bug report #953437,
regarding libflint-arb2: acb_hypgeom_2f1 segfaults on mips64el for certain
inputs
to be marked as done.
This means that you claim that the problem has been dealt with.
If this is not the case it is now your responsibility to reopen the
Bug report if necessary, and/or fix the problem forthwith.
(NB: If you are a system administrator and have no idea what this
message is talking about, this may indicate a serious mail system
misconfiguration somewhere. Please contact [email protected]
immediately.)
--
953437: https://bugs.debian.org/cgi-bin/bugreport.cgi?bug=953437
Debian Bug Tracking System
Contact [email protected] with problems
--- Begin Message ---
Package: libflint-arb2
Version: 1:2.17.0-1
Severity: normal
Dear Maintainer,
acb_hypgeom_2f1 segfaults on mips64el for certain inputs.
This triggered a FTBFS through the test-suite of Gyoto:
https://buildd.debian.org/status/fetch.php?pkg=gyoto&arch=mips64el&ver=1.4.3-2&stamp=1582934308&raw=0
which was passing previously:
https://buildd.debian.org/status/fetch.php?pkg=gyoto&arch=mips64el&ver=1.4.3-1%2Bb1&stamp=1574106308&raw=0
and still passes on all other architectures.
libflint-arb seems to be the same version but libflint changed from
version 2.5.2-20 to 2.5.2-21. The bug may thus actually be in libflint.
I wrote a small case:
8<----8<----8<----8<---- hypergeom.c 8<----8-----8<----
#include <stdio.h>
#include <stdlib.h>
#include <acb_hypgeom.h>
double hypergeom (double kappaIndex, double thetae) {
// See documentation: http://arblib.org/acb_hypgeom.html#c.acb_hypgeom_2f1
acb_t FF, aa, bb, cc, zed;
acb_init(FF);
acb_init(aa);
acb_init(bb);
acb_init(cc);
acb_init(zed);
acb_set_d_d(aa, kappaIndex-1./3., 0.);
acb_set_d_d(bb, kappaIndex+1., 0.);
acb_set_d_d(cc, kappaIndex+2./3., 0.);
acb_set_d_d(zed, -kappaIndex*thetae, 0.);
slong prec=53; // 53 for double precision
acb_hypgeom_2f1(FF, aa, bb, cc, zed, ACB_HYPGEOM_2F1_AC, prec);
double hypergeom = arf_get_d(&acb_realref(FF)->mid, ARF_RND_NEAR);
// uncertainty
// double rad = mag_get_d(&acb_realref(FF)->rad);
acb_clear(FF);
acb_clear(aa);
acb_clear(bb);
acb_clear(cc);
acb_clear(zed);
return hypergeom;
}
int main(int argc, char **argv) {
double k=5.5, t=2.279935427802561, h=0;
if (argc >= 2) k=atof(argv[1]);
if (argc >= 3) t=atof(argv[2]);
h=hypergeom(k, t);
printf ("hypergeom(%e, %e)=%e\n",k, t, h);
return 0;
}
8<----8<----8<----8<----8<----8<----8<----8<----8<----8<----
Compile with
gcc -o hypergeom main.c -lflint-arb -lflint
The program segfaults for the built-in data as well as for "large"
values of the two arguments:
(sid_mips64el-dchroot)thibaut@eller:~$ ./hypergeom 0.6 0.5
Segmentation fault
(sid_mips64el-dchroot)thibaut@eller:~$ ./hypergeom 0.5 0.6
Segmentation fault
but runs fine for smaller values of each argument:
(sid_mips64el-dchroot)thibaut@eller:~$ ./hypergeom 0.5 0.5
hypergeom(5.000000e-01, 5.000000e-01)=9.539590e-01
Here comes a backtrace:
(gdb) run
Starting program: /home/thibaut/hypergeom
[Thread debugging using libthread_db enabled]
Using host libthread_db library
"/lib/mips64el-linux-gnuabi64/libthread_db.so.1".
Program received signal SIGSEGV, Segmentation fault.
flint_mpz_set_si (r=0xaaaaacfb60, r=0xaaaaacfb60, s=1025435369) at
./gmpcompat.h:66
warning: Source file is more recent than executable.
66 r->_mp_d[0] = s;
(gdb) bt
#0 flint_mpz_set_si (r=0xaaaaacfb60, r=0xaaaaacfb60, s=1025435369) at
./gmpcompat.h:66
#1 fmpz_mul_2exp (f=0xffffffb348, g=<optimized out>, exp=33) at
mul_2exp.c:51
#2 0x000000fff7df2738 in arb_get_interval_fmpz_2exp (a=0xffffffb340,
b=0xffffffb348, exp=0xffffffb350,
x=<optimized out>) at get_interval_fmpz_2exp.c:61
#3 0x000000fff7de8138 in arb_get_unique_fmpz (z=0xaaaaacedf0,
x=0xffffffb3e8) at get_unique_fmpz.c:67
#4 0x000000fff7f5de24 in bernoulli_rev_next (numer=0xaaaaacedf0,
denom=0xaaaaacedf8, iter=0xffffffb4e8)
at rev_next.c:80
#5 0x000000fff7f5d8b8 in bernoulli_cache_compute (n=<optimized out>) at
cache_compute.c:55
#6 0x000000fff7e16f08 in arb_gamma_stirling_coeff (b=0xffffffb678, k=8,
digamma=<optimized out>,
prec=10) at gamma.c:339
#7 0x000000fff7e17348 in arb_gamma_stirling_eval (s=0xffffffb828,
z=0xffffffb7f8, nterms=9,
digamma=<optimized out>, prec=59) at gamma.c:396
#8 0x000000fff7e18ac4 in _arb_gamma (y=y@entry=0xffffffbcf8,
x=0xffffffc370, prec=53,
inverse=inverse@entry=0) at gamma.c:871
#9 0x000000fff7e18de8 in arb_gamma (y=0xaaaaacfb60,
y@entry=0xffffffbcf8, x=0xaaaaacfb60, prec=33)
at gamma.c:898
#10 0x000000fff7e6d558 in acb_gamma (y=0xffffffbcf8, x=<optimized out>,
prec=<optimized out>)
at gamma.c:220
#11 0x000000fff7ee33b8 in acb_hypgeom_2f1_transform_nolimit
(res=0xffffffc250, a=0xffffffc2b0,
b=0xffffffc310, c=0xffffffc370, z=0xffffffc3d0,
regularized=<optimized out>, which=<optimized out>,
prec=53) at 2f1_transform.c:248
#12 0x000000fff7ee3a44 in acb_hypgeom_2f1_transform (res=0xffffffc250,
a=0xffffffc2b0, b=0xffffffc310,
c=0xffffffc370, z=<optimized out>, flags=<optimized out>,
which=<optimized out>, prec=53)
at 2f1_transform.c:458
#13 0x000000fff7ed9f10 in acb_hypgeom_2f1 (res=0xffffffc250,
a=0xffffffc2b0, b=0xffffffc310,
c=0xffffffc370, z=0xffffffc3d0, flags=<optimized out>, prec=53) at
2f1.c:185
#14 0x000000aaaaaab180 in hypergeom ()
#15 0x000000aaaaaab304 in main ()
_mp_d above is initialized through _fmpz_promote which I don't
understand.
Regards, Thibaut.
-- System Information:
Debian Release: bullseye/sid
APT prefers unstable-debug
APT policy: (500, 'unstable-debug'), (500, 'unstable')
Architecture: mips64el (mips64)
Kernel: Linux 4.19.0-8-octeon (SMP w/4 CPU cores)
Locale: LANG=C, LC_CTYPE=C.UTF-8 (charmap=UTF-8), LANGUAGE=C (charmap=UTF-8)
Shell: /bin/sh linked to /bin/dash
Init: unable to detect
Versions of packages libflint-arb2 depends on:
ii libc6 2.29-10
ii libflint-2.5.2 2.5.2-21
ii libgmp10 2:6.2.0+dfsg-4
ii libmpfr6 4.0.2-1
libflint-arb2 recommends no packages.
libflint-arb2 suggests no packages.
-- no debconf information
signature.asc
Description: OpenPGP digital signature
--- End Message ---
--- Begin Message ---
Source: flint
Source-Version: 2.5.2-22
Done: Julien Puydt <[email protected]>
We believe that the bug you reported is fixed in the latest version of
flint, which is due to be installed in the Debian FTP archive.
A summary of the changes between this version and the previous one is
attached.
Thank you for reporting the bug, which will now be closed. If you
have further comments please address them to [email protected],
and the maintainer will reopen the bug report if appropriate.
Debian distribution maintenance software
pp.
Julien Puydt <[email protected]> (supplier of updated flint package)
(This message was generated automatically at their request; if you
believe that there is a problem with it please contact the archive
administrators by mailing [email protected])
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA512
Format: 1.8
Date: Tue, 19 May 2020 16:29:46 +0200
Source: flint
Architecture: source
Version: 2.5.2-22
Distribution: unstable
Urgency: medium
Maintainer: Debian Science Maintainers
<[email protected]>
Changed-By: Julien Puydt <[email protected]>
Closes: 953437
Changes:
flint (2.5.2-22) unstable; urgency=medium
.
* Add patch fixing segfaults on mips64el (Closes: #953437).
* Bump std-ver to 4.5.0.
Checksums-Sha1:
d47fc4ed71a9bb19bc92f3c00eb1acd09260993f 2143 flint_2.5.2-22.dsc
cf146e9c6a72cf0a0750325770df061d8a0d4ce2 9040 flint_2.5.2-22.debian.tar.xz
e1576935cd937e0ec543f020323a788fb19537c2 7800 flint_2.5.2-22_source.buildinfo
Checksums-Sha256:
1fa7c6ee180b24a9ec2417bee69659a9b1257f6530ffd29570dc92e2cf1f8f88 2143
flint_2.5.2-22.dsc
70cd8611a9ed36e1a089be3ee7808f9a3a5a8ddc2c2c4a338504c32fbfd8b822 9040
flint_2.5.2-22.debian.tar.xz
e134edb5aeed1c1650319a1a525c9c78d25073948f35b1f2cde8fd56490cf78c 7800
flint_2.5.2-22_source.buildinfo
Files:
d1ce4d418866701de062b73d5615aaf4 2143 math optional flint_2.5.2-22.dsc
4fe6f5bf3b85517aa4ef2239bf627de1 9040 math optional
flint_2.5.2-22.debian.tar.xz
4ff41211b1e8f8b0512f85c0bf4bb330 7800 math optional
flint_2.5.2-22_source.buildinfo
-----BEGIN PGP SIGNATURE-----
iQJGBAEBCgAwFiEEgS7v2KP7pKzk3xFLBMU71/4DBVEFAl7D7X8SHGpwdXlkdEBk
ZWJpYW4ub3JnAAoJEATFO9f+AwVR0dsQAKZaTT/XHSOJIbA/kjKbo6C/hLky1Yrd
gv62JMqyJWvVx1+0vSRPxHnuf2QebsEwcV8U30T4SjCf2MQzuqCQhPJ9Z+medqa7
S3cDXvdQ15gwlx4JC8FYO4Bdz0WyyTCNuKvOuHuKwJeT6lS8Yuvn+P/SC963XbfN
isH3OSQJtsuaxG9nxGgfjULUOack4ATUAGnHb40fVK5xHch/IiJgBC+wBxJ++jJ/
7ZXbAKOOwEzycgoRqvEvLk6UKlIVGTyzfQktv3+tA+NKvy0auwOisVLLmLvF/Lh6
I8UX9WQqaImqji7CvKzKdw3YdoSgj8gQ/IIMMNRSQaskpTGsJAJouvpunfqFGYHG
aN1BJY05afIwG0wQAGBV7LvOBWkrj0yLEq0n4pB8OfIq2UVx2N3QtS9NC7DAZSnJ
Fiyya1pNJSX+B97ntaErjDp7YsiJ5UbrUgR48Lvsvp+Bk+P4eYVKFAfxCS5xGKEe
97D3ihMNcAAId6twIORTA2WDtxXS1995VJmWrp7vhbK98vM0aQGrT9G0xrcMm20v
axbOO/HIEYSZFc5Oh0U2C30WRD3r1Kf/4HMy6NDFXfmkLZIrhD41N/JlwXe937Op
pOptWkFMAlEYVjgAkAiLcwM23IhR9VLjVFgCrSGtcaTvTAUyJI7PKpAp7iYumefw
rJP6pU7rFTbo
=v3kQ
-----END PGP SIGNATURE-----
--- End Message ---
--
debian-science-maintainers mailing list
[email protected]
https://alioth-lists.debian.net/cgi-bin/mailman/listinfo/debian-science-maintainers
