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