octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Test failure with gsl package / help needed


From: Julien Bect
Subject: Re: Test failure with gsl package / help needed
Date: Tue, 13 Sep 2016 19:58:35 +0200
User-agent: Mozilla/5.0 (X11; Linux x86_64; rv:45.0) Gecko/20100101 Icedove/45.2.0

Le 13/09/2016 à 09:02, Julien Bect a écrit :
I am seeing "random" test failures with the gsl package (on a Debian testing, x86_64-linux-gnu, Octave 4.0.3, gcc 6.1.1).

The failure occurs after several (might be 2, might be 20...) calls to test_gslsf:

>> pkg load gsl
>> test_gslsf
PASSES 25 out of 25 tests
>> test_gslsf
PASSES 25 out of 25 tests
>> test_gslsf
PASSES 25 out of 25 tests
>> test_gslsf
PASSES 25 out of 25 tests
>> test_gslsf
PASSES 25 out of 25 tests
>> test_gslsf
PASSES 25 out of 25 tests

....  (run many times)

>> test_gslsf
***** test # (int, double, double) to double
 conicalP_cyl_reg_sc=zeros(N,1);
 conicalP_cyl_reg_sc_err=zeros(N,1);
 for i=1:N
[conicalP_cyl_reg_sc(i) conicalP_cyl_reg_sc_err(i)]=conicalP_cyl_reg(xi(i),xd(i),xd2(i));
 end
[conicalP_cyl_reg_vec conicalP_cyl_reg_vec_err]=conicalP_cyl_reg(xi,xd,xd2);
 conicalP_cyl_reg_val=conicalP_cyl_reg(xi,xd,xd2);
 assert(norm(conicalP_cyl_reg_val(:)-conicalP_cyl_reg_vec(:))==0.0)
 assert(norm(conicalP_cyl_reg_sc(:)-conicalP_cyl_reg_vec(:))==0.0)
assert(norm(conicalP_cyl_reg_sc_err(:)-conicalP_cyl_reg_vec_err(:))==0.0)
!!!!! test failed
GSL error 11 at legendre_con.c, line 108: error

Ok, the problem should now be fixed.

It was simply a case of random test data getting too close to zero, where one GSL algorithm apparently struggles to converge...

Thanks again to Tatsuro, Marco and Dmitri.

@++
Julien



reply via email to

[Prev in Thread] Current Thread [Next in Thread]