summaryrefslogtreecommitdiff
path: root/src/MSL_C.PPCEABI.bare.H/s_ldexp.c
diff options
context:
space:
mode:
authormrb0nk500 <b0nk@b0nk.xyz>2023-02-02 17:29:19 -0400
committermrb0nk500 <b0nk@b0nk.xyz>2023-02-03 13:27:06 -0400
commitf2eabdb6257c09cf2890dac5e9737912728542af (patch)
tree3b46f6787185d65605651a0f48776dc9779ce648 /src/MSL_C.PPCEABI.bare.H/s_ldexp.c
parenteef1dd840b7cecac28c2e6b0574707b90a37d4e7 (diff)
global: Add rest of Dolphin SDK proper, add MSL, and MetroTRK
Finally, it links properly.
Diffstat (limited to 'src/MSL_C.PPCEABI.bare.H/s_ldexp.c')
-rw-r--r--src/MSL_C.PPCEABI.bare.H/s_ldexp.c54
1 files changed, 54 insertions, 0 deletions
diff --git a/src/MSL_C.PPCEABI.bare.H/s_ldexp.c b/src/MSL_C.PPCEABI.bare.H/s_ldexp.c
new file mode 100644
index 0000000..aeeda31
--- /dev/null
+++ b/src/MSL_C.PPCEABI.bare.H/s_ldexp.c
@@ -0,0 +1,54 @@
+/* @(#)s_ldexp.c 1.2 95/01/04 */
+/*
+ * ====================================================
+ * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
+ *
+ * Developed at SunPro, a Sun Microsystems, Inc. business.
+ * Permission to use, copy, modify, and distribute this
+ * software is freely granted, provided that this notice
+ * is preserved.
+ * ====================================================
+ */
+
+#include "fdlibm.h"
+
+static const double two54 = 1.80143985094819840000e+16, /* 0x43500000, 0x00000000 */
+ twom54 = 5.55111512312578270212e-17, /* 0x3C900000, 0x00000000 */
+ big = 1.0e+300, tiny = 1.0e-300;
+
+double ldexp(double x, int n) {
+ _INT32 k, hx, lx;
+ if (!isfinite(x) || x == 0.0)
+ return x;
+
+ hx = __HI(x);
+ lx = __LO(x);
+ k = (hx & 0x7ff00000) >> 20; /* extract exponent */
+ if (k == 0) { /* 0 or subnormal x */
+ if ((lx | (hx & 0x7fffffff)) == 0)
+ return x; /* +-0 */
+ x *= two54;
+ hx = __HI(x);
+ k = ((hx & 0x7ff00000) >> 20) - 54;
+ if (n < -50000)
+ return tiny * x; /*underflow*/
+ }
+ if (k == 0x7ff)
+ return x + x; /* NaN or Inf */
+ k = k + n;
+ if (k > 0x7fe)
+ return big * copysign(big, x); /* overflow */
+ if (k > 0) /* normal result */
+ {
+ __HI(x) = (hx & 0x800fffff) | (k << 20);
+ return x;
+ }
+ if (k <= -54)
+ if (n > 50000) /* in case integer overflow in n+k */
+ return big * copysign(big, x); /*overflow*/
+ else
+ return tiny * copysign(tiny, x); /*underflow*/
+ k += 54; /* subnormal result */
+ __HI(x) = (hx & 0x800fffff) | (k << 20);
+ return x * twom54;
+}