[Octave-bug-tracker] [bug #45507] acos returns different results on big input values

[Marco Caliari]

Follow-up Comment #14, bug #45507 (project octave):

Hi,

I tried 


write(6,*) acos((1.0d8,0.0d0))


in fortran and got (0.0000000000000000,-19.113827924512311). A part from the
sign of the imaginary part (which can be easly changed), no overflow due to
catastrophic cancellation. Why is it not possible to use the same standard
math library used by fortran (gfortran-4.8 in my case) and fix the sign later?

    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?45507>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/