Problem with hazard Function

[Janes, Robert (Columbus)]

Hello,

I am working with a large Fortran program being ported from Intel Fortran / 
IMSL to gfortran and gsl.  In working with the program I found that 
gsl_sf_hazard does not handle extreme cases the way I believe it should.  I 
have addressed this in my code but suggest that maybe gsl should incorporate 
this change.

To illustrate the problem run the following Linux script code using the 
attached source code files:
gcc -c -o imsl.o imsl_c.c
gfortran -o gsl.exe gsl.f imsl.o -lgsl -lgslcblas
./gsl.exe

The function returns zero for negative values past a certain value until it 
finally throws an underflow error and aborts.  I suggest that returning zero no 
matter how negative the input is a reasonable change.

Thank you.


Rob Janes

SGS Transportation

Senior Software Engineer

SGS – CMX

2860 N. National Road Suite A

US – 47201 – Columbus, Indiana

Phone: +1 - 812 - 378 - 7966

Fax: +1 812 - 378 - 3393

Email: address@hidden



Information in this email and any attachments is confidential and intended 
solely for the use of the individual(s) to whom it is addressed or otherwise 
directed. Please note that any views or opinions presented in this email are 
solely those of the author and do not necessarily represent those of the 
Company. Finally, the recipient should check this email and any attachments for 
the presence of viruses. The Company accepts no liability for any damage caused 
by any virus transmitted by this email. All SGS services are rendered in 
accordance with the applicable SGS conditions of service available on request 
and accessible at https://www.sgs.com/en/terms-and-conditions

gsl.f
Description: gsl.f

imsl_c.c
Description: imsl_c.c