Hi, does anybody know for which x huge+x>one does _not_hold? Where huge = 1.000e+300, one = 1.00000000000000000000e+00 |x| < 2**-27 I've no idea, but I need that to build a test case.
Anyways we think t=0 is fine for |x| < 2**-27 as in that range x = asin(x) wrt. to the preciseness of double. If t = 0, p will evalutate to p=0, q to q=1 and thus x will be returned. Best regards, Goetz. > -----Original Message----- > From: core-libs-dev [mailto:[email protected]] On Behalf > Of Andrew Haley > Sent: Montag, 5. Dezember 2016 11:10 > To: [email protected] > Subject: Re: potential error in fdlibm asin > > On 02/12/16 14:24, Lindenmaier, Goetz wrote: > > > > I found one copy of this code in the internet where the "else" is removed, > > but it's not completely the same: > > https://github.com/jerryscript-project/jerryscript/blob/master/jerry- > libm/asin.c > > > > Does anybody have the knowledge to tell me what's wrong? > > Does anybody know how to contribute a fix to netlib fdlibm? > > I think the only way to proceed with this is to find an input argument > which is affected by this very odd code. If there is one. does it > give a correct or an incorrect result? If not, we don't care and it > can be tidied up later. > > Andrew.
