jyapayne/musl
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

math: simplify hypot and hypotf using scalbn

Browse source 
this also fixes overflow/underflow raising and excess
precision issues (as those are handled well in scalbn)
This commit is contained in:
Szabolcs Nagy 2012-11-13 21:54:32 +01:00
commit 529a8de46c
2 changed files with 4 additions and 11 deletions
Download patch file Download diff file Expand all files Collapse all files

9
src/math/hypot.c
View file

@ -117,12 +117,7 @@ double hypot(double x, double y)
t2 = a - t1; t2 = a - t1;
w = sqrt(t1*y1-(w*(-w)-(t1*y2+t2*b))); w = sqrt(t1*y1-(w*(-w)-(t1*y2+t2*b)));
} }
if (k != 0) { if (k)
uint32_t high; w = scalbn(w, k);
t1 = 1.0;
GET_HIGH_WORD(high, t1);
SET_HIGH_WORD(t1, high+(k<<20));
return t1*w;
}
return w; return w;
} }

6
src/math/hypotf.c
View file

@ -80,9 +80,7 @@ float hypotf(float x, float y)
t2 = a - t1; t2 = a - t1;
w = sqrtf(t1*y1-(w*(-w)-(t1*y2+t2*b))); w = sqrtf(t1*y1-(w*(-w)-(t1*y2+t2*b)));
} }
if (k != 0) { if (k)
SET_FLOAT_WORD(t1, 0x3f800000+(k<<23)); w = scalbnf(w, k);
return t1*w;
}
return w; return w;
} }