GNU Scientific Library - Bugs: bug #52927, make check fails on Bessel j2 test
You are not allowed to post comments on this tracker with your current authentication level.
bug #52927: make check fails on Bessel j2 test
Submitter: | Patrick Alken <psa> | ||
Submitted: | Thu 18 Jan 2018 11:13:07 PM UTC | ||
Category: | Runtime error | Severity: | 3 - Normal |
Operating System: | Status: | None | |
Assigned to: | None | Open/Closed: | Open |
Release: |
Sat 21 Jul 2018 12:41:59 PM UTC, comment #2: |
Vladimir Khodygo <vladk> |
Mon 05 Mar 2018 11:01:38 PM UTC, comment #1: Just a note that we're getting the same error on CentOS 6, specifically the condaforge/linux-anvil Docker build ( https://hub.docker.com/r/condaforge/linux-anvil/ ); exact same output as below. TEST_TOL6 also isn't enough for the test to pass; I guess we're just commenting out the test for now.
|
Dougal Sutherland <dougals> |
Thu 18 Jan 2018 11:13:07 PM UTC, original submission:
from eshell =at= ucsc =dot= edu
Thanks much for the quick reply! I incremented the tolerance value, but it is still failing even at TEST_TOL6: ====================================== gsl 2.4: specfunc/test-suite.log ====================================== # TOTAL: 1 # PASS: 0 # SKIP: 0 # XFAIL: 0 # FAIL: 1 # XPASS: 0 # ERROR: 0 .. contents:: :depth: 2 FAIL: test ========== FAIL: gsl_sf_bessel_j2_e(1048576.0, &r) [168] expected: -3.1518539455252412e-07 obtained: -3.1518539455252539e-07 ± 2.7994086564622246e-22 (rel=8.88178e-16) fracdiff: 2.0155588470164931e-15 tolerance: 2.3283064365386963e-10 value/expected not consistent within reported error -3.151853945525253879e-07 2.799408656462224591e-22 FAIL: Bessel Functions [407] On Thu, Jan 18, 2018 at 3:01 PM, Patrick Alken <alken@colorado.edu> wrote:
|
Patrick Alken <psa>![]() |
No files currently attached
Depends on the following items: None found
Items that depend on this one: None found
I do confirm this error. Our cluster uses CentOS 6 as well and gives almost _exactly_ the same values:
Updating gsl to version 2.5 fixes this problem.
At the same time I have no problems with my desktop:
I also would like to report a bug here since it's related to Bessel functions. According to this an example from the documentation gives incorrect answers.
What you get: JO(5) = -1.917848549326277019e-01.
What you expect (see documentation): JO(5) = -1.775967713143382920e-01.
When I replace gsl_sf_bessel_j0 with gsl_sf_bessel_J0 I get JO(5) = -1.775967713143382642e-01 which is close to the desired output.