------- Comment #4 from jb at gcc dot gnu dot org 2007-04-27 11:29 ------- (In reply to comment #3) > 1. Convert to single precision > 2. Calculate rcp(s|p)s or rsqrt(p|s)s > 3. Refine with newton iteration > > vs. just using div(p|s)d or sqrt(p|s)d?
This should be 1. Convert to single precision 2. Calculate rcp(s|p)s or rsqrt(p|s)s 3. Convert back to double precision 4. Refine with newton iteration -- http://gcc.gnu.org/bugzilla/show_bug.cgi?id=31723
