------- Comment #11 from sgk at troutmask dot apl dot washington dot edu 2009-04-30 16:09 ------- Subject: Re: gfortran support for non-standard sind,cosd and friends intrinsics
On Thu, Apr 30, 2009 at 07:38:31AM -0000, ruben at tapir dot caltech dot edu wrote: > > > ------- Comment #9 from ruben at tapir dot caltech dot edu 2009-04-30 07:38 > ------- > Someone suggested multiplying or dividing the argument by a factor. > That does not do - the floating point precision achieved would be > insufficient. > Yep. I even wrote the necessary code. You're more than welcomed to improve it. http://gcc.gnu.org/ml/fortran/2006-12/msg00329.html You need to (1) Treat special cases. (2) Add argument reduction in degrees. (3) Treat special cases of the reduced argument. (4) Convert to radians and use the hardware FPU, or implement a Remes polynomial approximation to evaluation the function in degrees. -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=27452
