Re: lib/50698: hypotf() of small numbers yields infinity

[Andreas Gustafsson <gson%gson.org@localhost>]

The following reply was made to PR lib/50698; it has been noted by GNATS.

From: Andreas Gustafsson <gson%gson.org@localhost>
To: Rin Okuyama <okuyama%flex.phys.tohoku.ac.jp@localhost>
Cc: gnats-bugs%NetBSD.org@localhost
Subject: Re: lib/50698: hypotf() of small numbers yields infinity
Date: Mon, 25 Jan 2016 17:01:25 +0200

 Rin Okuyama wrote:
 >  This is because some magic numbers in e_hypotf.c are wrong. In FreeBSD,
 >  this bug was fixed by commit #23397:
 >  
 >    https://svnweb.freebsd.org/base/head/lib/msun/src/e_hypotf.c#rev23397
 >  
 >  I attached the minimal patch below. With this patch, I obtain (on amd64)
 >  
 >    % ./test-hypotf 1e-18 1e-18
 >    1.41421363268178e-18
 >    %
 
 Thank you for your analysis and patch.  The patch has now been committed.
 
 >  Other changes have been made for libm in FreeBSD. IMO, we should merge
 >  reasonable ones.
 
 Additional patches for the ones you think are reasonable are welcome.
 More test cases would be good, too.
 -- 
 Andreas Gustafsson, gson%gson.org@localhost