bugGNU 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

Submitted by:  Patrick Alken <psa>
Submitted on:  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:

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.

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.

https://github.com/conda-forge/gsl-feedstock/pull/24#issuecomment-370585026

Dougal Sutherland <dougals>
Thu 18 Jan 2018 11:13:07 PM UTC, original submission:

from eshell =at= ucsc =dot= edu

subject: GSL 2.4 make check failing on Rocks cluster running CentOS 6.9

note: it looks like the error term (r) is causing the failure, not the computed j2 value itself


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:

> Hello,
>
> It looks like its calculating the value correctly, but the test
> tolerance needs to be relaxed slightly. Can you locate this line in
> specfun/test_bessel.c:
>
> 186 TEST_SF(s, gsl_sf_bessel_j2_e, (1048576.0, &r),
> -3.1518539455252413111e-07, TEST_TOL3, GSL_SUCCESS);
>
> and change the TEST_TOL3 to TEST_TOL4, and let me know if the test passes?
> If not try TEST_TOL5 and then TEST_TOL6, and tell me which one allows the
> test to pass.
>
> Thanks,
> Patrick
>
>
> On 01/18/2018 03:56 PM, Eric Shell wrote:
>
>> Hello,
>>
>> I'm trying to install GSL 2.4 on a Rocks cluster running CentOS 6.9, gcc
>> version 4.4.7. I am running configure with just a --prefix argument.
>> make
>> succeeds without errors, but make check is failing on the specfunc test.
>> Here are the contents of the log file:
>>
>> ------------------------------------------------------------
>> --------------------
>>
>> ======================================
>> 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: 4.5474735088646412e-13
>> value/expected not consistent within reported error
>> -3.151853945525253879e-07 2.799408656462224591e-22
>> FAIL: Bessel Functions [407]
>>
>> ------------------------------------------------------------
>> --------------------
>>
>> Is this a known issue? How can I address the underlying issue?
>>
>> Thanks!
>>
>> - Eric
>>

Patrick Alken <psa>
Project Administrator

 

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by vladk (Posted a comment)
  • -email is unavailable- added by dougals (Posted a comment)
  • -email is unavailable- added by psa (Submitted the item)
  •  

     

     

    No changes have been made to this item

    Back to the top


    Powered by Savane 3.3