Merge pull request #2905 from juj/musl_mathfuncs

Musl mathfuncs

Author: kripken

emscripten-version.txt

@@ -1,2 +1,2 @@
-1.28.0
+1.28.1
 

src/library.js

@@ -4518,202 +4518,11 @@ LibraryManager.library = {
   llvm_log_f64: 'Math_log',
   llvm_exp_f32: 'Math_exp',
   llvm_exp_f64: 'Math_exp',
-  cbrt: function(x) {
-    return Math.pow(x, 1/3);
-  },
-  cbrtf: 'cbrt',
-  cbrtl: 'cbrt',
-
-  modf: function(x, intpart) {
-    {{{ makeSetValue('intpart', 0, '(x < 0) ? Math.ceil(x) : Math.floor(x)', 'double') }}};
-    return x - {{{ makeGetValue('intpart', 0, 'double') }}};
-  },
-  modff: function(x, intpart) {
-    {{{ makeSetValue('intpart', 0, '(x < 0) ? Math.ceil(x) : Math.floor(x)', 'float') }}};
-    return x - {{{ makeGetValue('intpart', 0, 'float') }}};
-  },
-  finite: function(x) {
-    return isFinite(x);
-  },
-  __finite: 'finite',
-  isinf: function(x) {
-    return !isNaN(x) && !isFinite(x);
-  },
-  __isinf: 'isinf',
-  isnan: function(x) {
-    return isNaN(x);
-  },
-  __isnan: 'isnan',
 
   _reallyNegative: function(x) {
     return x < 0 || (x === 0 && (1/x) === -Infinity);
   },
 
-  copysign__deps: ['_reallyNegative'],
-  copysign: function(a, b) {
-    return __reallyNegative(a) === __reallyNegative(b) ? a : -a;
-  },
-  copysignf: 'copysign',
-  copysignl: 'copysign',
-  __signbit__deps: ['copysign'],
-  __signbit: function(x) {
-    // We implement using copysign so that we get support
-    // for negative zero (once copysign supports that).
-    return _copysign(1.0, x) < 0;
-  },
-  __signbitd: '__signbit',
-  __signbitf: '__signbit',
-  __signbitl: '__signbit',
-  hypot: function(a, b) {
-     return Math.sqrt(a*a + b*b);
-  },
-  hypotf: 'hypot',
-  hypotl: 'hypot',
-  sinh: function(x) {
-    var p = Math.pow(Math.E, x);
-    return (p - (1 / p)) / 2;
-  },
-  sinhf: 'sinh',
-  sinhl: 'sinh',
-  cosh: function(x) {
-    var p = Math.pow(Math.E, x);
-    return (p + (1 / p)) / 2;
-  },
-  coshf: 'cosh',
-  coshl: 'cosh',
-  tanh__deps: ['sinh', 'cosh'],
-  tanh: function(x) {
-    return _sinh(x) / _cosh(x);
-  },
-  tanhf: 'tanh',
-  tanhl: 'tanh',
-  asinh: function(x) {
-    return Math.log(x + Math.sqrt(x * x + 1));
-  },
-  asinhf: 'asinh',
-  asinhl: 'asinh',
-  acosh: function(x) {
-    return Math.log(x * 1 + Math.sqrt(x * x - 1));
-  },
-  acoshf: 'acosh',
-  acoshl: 'acosh',
-  atanh: function(x) {
-    return Math.log((1 + x) / (1 - x)) / 2;
-  },
-  atanhf: 'atanh',
-  atanhl: 'atanh',
-  exp2: function(x) {
-    return Math.pow(2, x);
-  },
-  exp2f: 'exp2',
-  exp2l: 'exp2',
-  expm1: function(x) {
-    return Math.exp(x) - 1;
-  },
-  expm1f: 'expm1',
-  expm1l: 'expm1',
-  round: function(x) {
-    return (x < 0) ? -Math.round(-x) : Math.round(x);
-  },
-  roundf: 'round',
-  roundl: 'round',
-  lround: 'round',
-  lroundf: 'round',
-  lroundl: 'round',
-  llround: 'round',
-  llroundf: 'round',
-  llroundl: 'round',
-  rint: function(x) {
-    if (Math.abs(x % 1) !== 0.5) return Math.round(x);
-    return x + x % 2 + ((x < 0) ? 1 : -1);
-  },
-  rintf: 'rint',
-  rintl: 'rint',
-  lrint: 'rint',
-  lrintf: 'rint',
-  lrintl: 'rint',
-#if USE_TYPED_ARRAYS == 2
-  llrint: function(x) {
-    x = (x < 0) ? -Math.round(-x) : Math.round(x);
-    {{{ makeStructuralReturn(splitI64('x')) }}};
-  },
-#else
-  llrint: 'rint',
-#endif
-  llrintf: 'llrint',
-  llrintl: 'llrint',
-  nearbyint: 'rint',
-  nearbyintf: 'rint',
-  nearbyintl: 'rint',
-  trunc: function(x) {
-    return (x < 0) ? Math.ceil(x) : Math.floor(x);
-  },
-  truncf: 'trunc',
-  truncl: 'trunc',
-  fdim: function(x, y) {
-    return (x > y) ? x - y : 0;
-  },
-  fdimf: 'fdim',
-  fdiml: 'fdim',
-  fmax: function(x, y) {
-    return isNaN(x) ? y : isNaN(y) ? x : Math.max(x, y);
-  },
-  fmaxf: 'fmax',
-  fmaxl: 'fmax',
-  fmin: function(x, y) {
-    return isNaN(x) ? y : isNaN(y) ? x : Math.min(x, y);
-  },
-  fminf: 'fmin',
-  fminl: 'fmin',
-  fma: function(x, y, z) {
-    return x * y + z;
-  },
-  fmaf: 'fma',
-  fmal: 'fma',
-  fmod: function(x, y) {
-    return x % y;
-  },
-  fmodf: 'fmod',
-  fmodl: 'fmod',
-  remainder: 'fmod',
-  remainderf: 'fmod',
-  remainderl: 'fmod',
-  log10: function(x) {
-    return Math.log(x) / Math.LN10;
-  },
-  log10f: 'log10',
-  log10l: 'log10',
-  log1p: function(x) {
-    return Math.log(1 + x);
-  },
-  log1pf: 'log1p',
-  log1pl: 'log1p',
-  log2: function(x) {
-    return Math.log(x) / Math.LN2;
-  },
-  log2f: 'log2',
-  log2l: 'log2',
-  nan: function(x) {
-    return NaN;
-  },
-  nanf: 'nan',
-  nanl: 'nan',
-
-  sincos: function(x, sine, cosine) {
-    var sineVal = Math.sin(x),
-        cosineVal = Math.cos(x);
-    {{{ makeSetValue('sine', '0', 'sineVal', 'double') }}};
-    {{{ makeSetValue('cosine', '0', 'cosineVal', 'double') }}};
-  },
-  sincosl: 'sincos',
-
-  sincosf: function(x, sine, cosine) {
-    var sineVal = Math.sin(x),
-        cosineVal = Math.cos(x);
-    {{{ makeSetValue('sine', '0', 'sineVal', 'float') }}};
-    {{{ makeSetValue('cosine', '0', 'cosineVal', 'float') }}};
-  },
-
   div: function(divt, numer, denom) {
     var quot = (numer / denom) | 0;
     var rem = numer - quot * denom;
@@ -4722,17 +4531,6 @@ LibraryManager.library = {
     return divt;
   },
 
-  __fpclassify: function(x) {
-    if (isNaN(x)) return {{{ cDefine('FP_NAN') }}};
-    if (!isFinite(x)) return {{{ cDefine('FP_INFINITE') }}};
-    if (x == 0) return {{{ cDefine('FP_ZERO') }}};
-    // FP_SUBNORMAL..?
-    return {{{ cDefine('FP_NORMAL') }}};
-  },
-  __fpclassifyd: '__fpclassify', // Needed by tests/python/python.le32.bc
-  __fpclassifyf: '__fpclassify',
-  __fpclassifyl: '__fpclassify',
-
   // ==========================================================================
   // sys/utsname.h
   // ==========================================================================

system/lib/libc.symbols

@@ -1,26 +1,83 @@
          T __floatscan
+         T __fpclassify
+         T __fpclassifyf
+         T __fpclassifyl
          T __intscan
          T __overflow
          T __shgetc
          T __shlim
+         T __signbit
+         T __signbitf
+         T __signbitl
          W __strtod_l
          W __strtof_l
          W __strtold_l
          T __toread
          T __towrite
          T __uflow
+         T acosh
+         T acoshf
+         T acoshl
+         T asinh
+         T asinhf
+         T asinhl
+         T atanh
+         T atanhf
+         T atanhl
          T atof
          T atoi
          T atol
          W bulk_free
          W calloc
+         T cbrt
+         T cbrtf
+         T cbrtl
+         T copysign
+         T copysignf
+         T copysignl
+         T cosh
+         T coshf
+         T coshl
          T drand48
          T erand48
+         T exp2
+         T exp2f
+         T exp2l
+         T expm1
+         T expm1f
+         T expm1l
+         T feclearexcept
+         T feraiseexcept
+         T fetestexcept
+         T fegetround
+         T __fesetround
+         T fegetenv
+         T fesetenv
+         T fdim
+         T fdimf
+         T fdiml
          T ffs
+         T finite
+         T finitef
+         T fma
+         T fmaf
+         T fmal
+         T fmax
+         T fmaxf
+         T fmaxl
+         T fmin
+         T fminf
+         T fminl
+         T fmod
+         T fmodf
+         T fmodl
          W free
          T frexp
          T frexpf
          T frexpl
+         T hypot
+         T hypotf
+         T hypotl
          W independent_calloc
          W independent_comalloc
          T initstate
@@ -31,8 +88,29 @@
          T lcg31
          T lcg64
          T lcong48
+         T llrint
+         T llrintf
+         T llrintl
+         T llround
+         T llroundf
+         T llroundl
          T loadstate
+         T log10
+         T log10f
+         T log10l
+         T log1p
+         T log1pf
+         T log1pl
+         T log2
+         T log2f
+         T log2l
          T lrand48
+         T lrint
+         T lrintf
+         T lrintl
+         T lround
+         T lroundf
+         T lroundl
          W mallinfo
          W malloc
          W malloc_footprint
@@ -44,8 +122,17 @@
          W malloc_usable_size
          W mallopt
          W memalign
+         T modf
+         T modff
+         T modfl
          T MUSL_vfprintf
          T mrand48
+         T nan
+         T nanf
+         T nanl
+         T nearbyint
+         T nearbyintf
+         T nearbyintl
          T nrand48
          W posix_memalign
          W pvalloc
@@ -54,9 +141,30 @@
          T random
          W realloc
          W realloc_in_place
+         T remainder
+         T remainderf
+         T remainderl
+         T rint
+         T rintf
+         T rintl
+         T round
+         T roundf
+         T roundl
          T scalbn
          T scalbnl
+         T sincos
+         T sincosf
+         T sincosl
+         T sinh
+         T sinhf
+         T sinhl
          T srand
+         T tanh
+         T tanhf
+         T tanhl
+         T trunc
+         T truncf
+         T truncl
          T memchr
          T memcmp
          T memcpy