@thi.ng/math 5.7.6 → 5.7.7

package/int.js

@@ -1,64 +1,128 @@
-const M8 = 0xff;
-const M16 = 0xffff;
-export const signExtend8 = (a) => ((a &= M8), a & 0x80 ? a | ~M8 : a);
-export const signExtend16 = (a) => ((a &= M16), a & 0x8000 ? a | ~M16 : a);
-export const addi8 = (a, b) => signExtend8((a | 0) + (b | 0));
-export const divi8 = (a, b) => signExtend8((a | 0) / (b | 0));
-export const muli8 = (a, b) => signExtend8((a | 0) * (b | 0));
-export const subi8 = (a, b) => signExtend8((a | 0) - (b | 0));
-export const andi8 = (a, b) => signExtend8((a | 0) & (b | 0));
-export const ori8 = (a, b) => signExtend8(a | 0 | (b | 0));
-export const xori8 = (a, b) => signExtend8((a | 0) ^ (b | 0));
-export const noti8 = (a) => signExtend8(~a);
-export const lshifti8 = (a, b) => signExtend8((a | 0) << (b | 0));
-export const rshifti8 = (a, b) => signExtend8((a | 0) >> (b | 0));
-export const addi16 = (a, b) => signExtend16((a | 0) + (b | 0));
-export const divi16 = (a, b) => signExtend16((a | 0) / (b | 0));
-export const muli16 = (a, b) => signExtend16((a | 0) * (b | 0));
-export const subi16 = (a, b) => signExtend16((a | 0) - (b | 0));
-export const andi16 = (a, b) => signExtend16((a | 0) & (b | 0));
-export const ori16 = (a, b) => signExtend16(a | 0 | (b | 0));
-export const xori16 = (a, b) => signExtend16((a | 0) ^ (b | 0));
-export const noti16 = (a) => signExtend16(~a);
-export const lshifti16 = (a, b) => signExtend16((a | 0) << (b | 0));
-export const rshifti16 = (a, b) => signExtend16((a | 0) >> (b | 0));
-export const addi32 = (a, b) => ((a | 0) + (b | 0)) | 0;
-export const divi32 = (a, b) => ((a | 0) / (b | 0)) | 0;
-export const muli32 = (a, b) => ((a | 0) * (b | 0)) | 0;
-export const subi32 = (a, b) => ((a | 0) - (b | 0)) | 0;
-export const andi32 = (a, b) => (a | 0) & (b | 0);
-export const ori32 = (a, b) => a | 0 | (b | 0);
-export const xori32 = (a, b) => (a | 0) ^ (b | 0);
-export const lshifti32 = (a, b) => (a | 0) << (b | 0);
-export const rshifti32 = (a, b) => (a | 0) >> (b | 0);
-export const noti32 = (a) => ~a;
-export const addu8 = (a, b) => ((a & M8) + (b & M8)) & M8;
-export const divu8 = (a, b) => ((a & M8) / (b & M8)) & M8;
-export const mulu8 = (a, b) => ((a & M8) * (b & M8)) & M8;
-export const subu8 = (a, b) => ((a & M8) - (b & M8)) & M8;
-export const andu8 = (a, b) => a & M8 & (b & M8) & M8;
-export const oru8 = (a, b) => ((a & M8) | (b & M8)) & M8;
-export const xoru8 = (a, b) => ((a & M8) ^ (b & M8)) & M8;
-export const notu8 = (a) => ~a & M8;
-export const lshiftu8 = (a, b) => ((a & M8) << (b & M8)) & M8;
-export const rshiftu8 = (a, b) => ((a & M8) >>> (b & M8)) & M8;
-export const addu16 = (a, b) => ((a & M16) + (b & M16)) & M16;
-export const divu16 = (a, b) => ((a & M16) / (b & M16)) & M16;
-export const mulu16 = (a, b) => ((a & M16) * (b & M16)) & M16;
-export const subu16 = (a, b) => ((a & M16) - (b & M16)) & M16;
-export const andu16 = (a, b) => a & M16 & (b & M16) & M16;
-export const oru16 = (a, b) => ((a & M16) | (b & M16)) & M16;
-export const xoru16 = (a, b) => ((a & M16) ^ (b & M16)) & M16;
-export const notu16 = (a) => ~a & M16;
-export const lshiftu16 = (a, b) => ((a & M16) << (b & M16)) & M16;
-export const rshiftu16 = (a, b) => ((a & M16) >>> (b & M16)) & M16;
-export const addu32 = (a, b) => ((a >>> 0) + (b >>> 0)) >>> 0;
-export const divu32 = (a, b) => ((a >>> 0) / (b >>> 0)) >>> 0;
-export const mulu32 = (a, b) => ((a >>> 0) * (b >>> 0)) >>> 0;
-export const subu32 = (a, b) => ((a >>> 0) - (b >>> 0)) >>> 0;
-export const andu32 = (a, b) => ((a >>> 0) & (b >>> 0)) >>> 0;
-export const oru32 = (a, b) => ((a >>> 0) | (b >>> 0)) >>> 0;
-export const xoru32 = (a, b) => ((a >>> 0) ^ (b >>> 0)) >>> 0;
-export const notu32 = (a) => ~a >>> 0;
-export const lshiftu32 = (a, b) => ((a >>> 0) << (b >>> 0)) >>> 0;
-export const rshiftu32 = (a, b) => ((a >>> 0) >>> (b >>> 0)) >>> 0;
+const M8 = 255;
+const M16 = 65535;
+const signExtend8 = (a) => (a &= M8, a & 128 ? a | ~M8 : a);
+const signExtend16 = (a) => (a &= M16, a & 32768 ? a | ~M16 : a);
+const addi8 = (a, b) => signExtend8((a | 0) + (b | 0));
+const divi8 = (a, b) => signExtend8((a | 0) / (b | 0));
+const muli8 = (a, b) => signExtend8((a | 0) * (b | 0));
+const subi8 = (a, b) => signExtend8((a | 0) - (b | 0));
+const andi8 = (a, b) => signExtend8((a | 0) & (b | 0));
+const ori8 = (a, b) => signExtend8(a | 0 | (b | 0));
+const xori8 = (a, b) => signExtend8((a | 0) ^ (b | 0));
+const noti8 = (a) => signExtend8(~a);
+const lshifti8 = (a, b) => signExtend8((a | 0) << (b | 0));
+const rshifti8 = (a, b) => signExtend8((a | 0) >> (b | 0));
+const addi16 = (a, b) => signExtend16((a | 0) + (b | 0));
+const divi16 = (a, b) => signExtend16((a | 0) / (b | 0));
+const muli16 = (a, b) => signExtend16((a | 0) * (b | 0));
+const subi16 = (a, b) => signExtend16((a | 0) - (b | 0));
+const andi16 = (a, b) => signExtend16((a | 0) & (b | 0));
+const ori16 = (a, b) => signExtend16(a | 0 | (b | 0));
+const xori16 = (a, b) => signExtend16((a | 0) ^ (b | 0));
+const noti16 = (a) => signExtend16(~a);
+const lshifti16 = (a, b) => signExtend16((a | 0) << (b | 0));
+const rshifti16 = (a, b) => signExtend16((a | 0) >> (b | 0));
+const addi32 = (a, b) => (a | 0) + (b | 0) | 0;
+const divi32 = (a, b) => (a | 0) / (b | 0) | 0;
+const muli32 = (a, b) => (a | 0) * (b | 0) | 0;
+const subi32 = (a, b) => (a | 0) - (b | 0) | 0;
+const andi32 = (a, b) => (a | 0) & (b | 0);
+const ori32 = (a, b) => a | 0 | (b | 0);
+const xori32 = (a, b) => (a | 0) ^ (b | 0);
+const lshifti32 = (a, b) => (a | 0) << (b | 0);
+const rshifti32 = (a, b) => (a | 0) >> (b | 0);
+const noti32 = (a) => ~a;
+const addu8 = (a, b) => (a & M8) + (b & M8) & M8;
+const divu8 = (a, b) => (a & M8) / (b & M8) & M8;
+const mulu8 = (a, b) => (a & M8) * (b & M8) & M8;
+const subu8 = (a, b) => (a & M8) - (b & M8) & M8;
+const andu8 = (a, b) => a & M8 & (b & M8) & M8;
+const oru8 = (a, b) => (a & M8 | b & M8) & M8;
+const xoru8 = (a, b) => (a & M8 ^ b & M8) & M8;
+const notu8 = (a) => ~a & M8;
+const lshiftu8 = (a, b) => (a & M8) << (b & M8) & M8;
+const rshiftu8 = (a, b) => (a & M8) >>> (b & M8) & M8;
+const addu16 = (a, b) => (a & M16) + (b & M16) & M16;
+const divu16 = (a, b) => (a & M16) / (b & M16) & M16;
+const mulu16 = (a, b) => (a & M16) * (b & M16) & M16;
+const subu16 = (a, b) => (a & M16) - (b & M16) & M16;
+const andu16 = (a, b) => a & M16 & (b & M16) & M16;
+const oru16 = (a, b) => (a & M16 | b & M16) & M16;
+const xoru16 = (a, b) => (a & M16 ^ b & M16) & M16;
+const notu16 = (a) => ~a & M16;
+const lshiftu16 = (a, b) => (a & M16) << (b & M16) & M16;
+const rshiftu16 = (a, b) => (a & M16) >>> (b & M16) & M16;
+const addu32 = (a, b) => (a >>> 0) + (b >>> 0) >>> 0;
+const divu32 = (a, b) => (a >>> 0) / (b >>> 0) >>> 0;
+const mulu32 = (a, b) => (a >>> 0) * (b >>> 0) >>> 0;
+const subu32 = (a, b) => (a >>> 0) - (b >>> 0) >>> 0;
+const andu32 = (a, b) => (a >>> 0 & b >>> 0) >>> 0;
+const oru32 = (a, b) => (a >>> 0 | b >>> 0) >>> 0;
+const xoru32 = (a, b) => (a >>> 0 ^ b >>> 0) >>> 0;
+const notu32 = (a) => ~a >>> 0;
+const lshiftu32 = (a, b) => a >>> 0 << (b >>> 0) >>> 0;
+const rshiftu32 = (a, b) => a >>> 0 >>> (b >>> 0) >>> 0;
+export {
+  addi16,
+  addi32,
+  addi8,
+  addu16,
+  addu32,
+  addu8,
+  andi16,
+  andi32,
+  andi8,
+  andu16,
+  andu32,
+  andu8,
+  divi16,
+  divi32,
+  divi8,
+  divu16,
+  divu32,
+  divu8,
+  lshifti16,
+  lshifti32,
+  lshifti8,
+  lshiftu16,
+  lshiftu32,
+  lshiftu8,
+  muli16,
+  muli32,
+  muli8,
+  mulu16,
+  mulu32,
+  mulu8,
+  noti16,
+  noti32,
+  noti8,
+  notu16,
+  notu32,
+  notu8,
+  ori16,
+  ori32,
+  ori8,
+  oru16,
+  oru32,
+  oru8,
+  rshifti16,
+  rshifti32,
+  rshifti8,
+  rshiftu16,
+  rshiftu32,
+  rshiftu8,
+  signExtend16,
+  signExtend8,
+  subi16,
+  subi32,
+  subi8,
+  subu16,
+  subu32,
+  subu8,
+  xori16,
+  xori32,
+  xori8,
+  xoru16,
+  xoru32,
+  xoru8
+};

package/interval.js

@@ -1,207 +1,75 @@
-/**
- * Clamps value `x` to given closed interval.
- *
- * @param x - value to clamp
- * @param min - lower bound
- * @param max - upper bound
- */
-export const clamp = (x, min, max) => (x < min ? min : x > max ? max : x);
-/**
- * Clamps value `x` to closed [0 .. ∞] interval.
- *
- * @param x -
- */
-export const clamp0 = (x) => (x > 0 ? x : 0);
-/**
- * Clamps value `x` to closed [0 .. 1] interval.
- *
- * @param x -
- */
-export const clamp01 = (x) => (x < 0 ? 0 : x > 1 ? 1 : x);
-/**
- * Clamps value `x` to closed [-1 .. 1] interval.
- *
- * @param x -
- */
-export const clamp11 = (x) => (x < -1 ? -1 : x > 1 ? 1 : x);
-/**
- * Clamps value `x` to closed [0 .. 0.5] interval.
- *
- * @param x -
- */
-export const clamp05 = (x) => (x < 0 ? 0 : x > 0.5 ? 0.5 : x);
-/**
- * Returns 2-tuple of [min(x,y), max(x,y)].
- *
- * @param x
- * @param y
- */
-export const minMax = (x, y) => x < y ? [x, y] : [y, x];
-/**
- * Folds `x` back inside closed [min..max] interval. Also see
- * {@link wrapOnce}.
- *
- * @param x -
- * @param min -
- * @param max -
- */
-export const wrap = (x, min, max) => {
-    if (min === max)
-        return min;
-    if (x > max) {
-        const d = max - min;
-        x -= d;
-        if (x > max)
-            x -= d * (((x - min) / d) | 0);
-    }
-    else if (x < min) {
-        const d = max - min;
-        x += d;
-        if (x < min)
-            x += d * (((min - x) / d + 1) | 0);
-    }
-    return x;
+const clamp = (x, min, max) => x < min ? min : x > max ? max : x;
+const clamp0 = (x) => x > 0 ? x : 0;
+const clamp01 = (x) => x < 0 ? 0 : x > 1 ? 1 : x;
+const clamp11 = (x) => x < -1 ? -1 : x > 1 ? 1 : x;
+const clamp05 = (x) => x < 0 ? 0 : x > 0.5 ? 0.5 : x;
+const minMax = (x, y) => x < y ? [x, y] : [y, x];
+const wrap = (x, min, max) => {
+  if (min === max)
+    return min;
+  if (x > max) {
+    const d = max - min;
+    x -= d;
+    if (x > max)
+      x -= d * ((x - min) / d | 0);
+  } else if (x < min) {
+    const d = max - min;
+    x += d;
+    if (x < min)
+      x += d * ((min - x) / d + 1 | 0);
+  }
+  return x;
 };
-/**
- * Like {@link wrap}, but optimized for cases where `x` is guaranteed to
- * be in `[min - d, max + d]` interval, where `d = max - min`. Result
- * will be in closed `[min..max]` interval.
- *
- * @param x -
- * @param min -
- * @param max -
- */
-export const wrapOnce = (x, min, max) => x < min ? x - min + max : x > max ? x - max + min : x;
-/**
- * Similar to {@link wrapOnce} for [0..1] interval.
- *
- * @param x -
- */
-export const wrap01 = (x) => (x < 0 ? x + 1 : x > 1 ? x - 1 : x);
-/**
- * Similar to {@link wrapOnce} for [-1..1] interval.
- *
- * @param x -
- */
-export const wrap11 = (x) => (x < -1 ? x + 2 : x > 1 ? x - 2 : x);
-export const min2id = (a, b) => (a <= b ? 0 : 1);
-export const min3id = (a, b, c) => a <= b ? (a <= c ? 0 : 2) : b <= c ? 1 : 2;
-export const min4id = (a, b, c, d) => a <= b
-    ? a <= c
-        ? a <= d
-            ? 0
-            : 3
-        : c <= d
-            ? 2
-            : 3
-    : b <= c
-        ? b <= d
-            ? 1
-            : 3
-        : c <= d
-            ? 2
-            : 3;
-export const max2id = (a, b) => (a >= b ? 0 : 1);
-export const max3id = (a, b, c) => a >= b ? (a >= c ? 0 : 2) : b >= c ? 1 : 2;
-export const max4id = (a, b, c, d) => a >= b
-    ? a >= c
-        ? a >= d
-            ? 0
-            : 3
-        : c >= d
-            ? 2
-            : 3
-    : b >= c
-        ? b >= d
-            ? 1
-            : 3
-        : c >= d
-            ? 2
-            : 3;
-/**
- * Returns the non-zero minimum value of the given `a`, `b` args.
- *
- * @param a -
- * @param b -
- */
-export const minNonZero2 = (a, b) => a !== 0 ? (b !== 0 ? Math.min(a, b) : a) : b;
-/**
- * Returns the non-zero minimum value of the given `a`, `b`, `c` args.
- *
- * @param a -
- * @param b -
- * @param c -
- */
-export const minNonZero3 = (a, b, c) => minNonZero2(minNonZero2(a, b), c);
-/**
- * See `smax()`.
- *
- * @param a -
- * @param b -
- * @param k - smooth exponent (MUST be > 0)
- */
-export const smin = (a, b, k) => smax(a, b, -k);
-/**
- * Smooth maximum. Note: Result values will be slightly larger than max value
- * near max(a,b) + eps due to exponential decay. Higher `k` values reduce the
- * error, but also reduce the smoothing. Recommended k=16.
- *
- * https://en.wikipedia.org/wiki/Smooth_maximum
- *
- * @param a -
- * @param b -
- * @param k - smooth exponent (MUST be > 0)
- */
-export const smax = (a, b, k) => {
-    const ea = Math.exp(a * k);
-    const eb = Math.exp(b * k);
-    return (a * ea + b * eb) / (ea + eb);
+const wrapOnce = (x, min, max) => x < min ? x - min + max : x > max ? x - max + min : x;
+const wrap01 = (x) => x < 0 ? x + 1 : x > 1 ? x - 1 : x;
+const wrap11 = (x) => x < -1 ? x + 2 : x > 1 ? x - 2 : x;
+const min2id = (a, b) => a <= b ? 0 : 1;
+const min3id = (a, b, c) => a <= b ? a <= c ? 0 : 2 : b <= c ? 1 : 2;
+const min4id = (a, b, c, d) => a <= b ? a <= c ? a <= d ? 0 : 3 : c <= d ? 2 : 3 : b <= c ? b <= d ? 1 : 3 : c <= d ? 2 : 3;
+const max2id = (a, b) => a >= b ? 0 : 1;
+const max3id = (a, b, c) => a >= b ? a >= c ? 0 : 2 : b >= c ? 1 : 2;
+const max4id = (a, b, c, d) => a >= b ? a >= c ? a >= d ? 0 : 3 : c >= d ? 2 : 3 : b >= c ? b >= d ? 1 : 3 : c >= d ? 2 : 3;
+const minNonZero2 = (a, b) => a !== 0 ? b !== 0 ? Math.min(a, b) : a : b;
+const minNonZero3 = (a, b, c) => minNonZero2(minNonZero2(a, b), c);
+const smin = (a, b, k) => smax(a, b, -k);
+const smax = (a, b, k) => {
+  const ea = Math.exp(a * k);
+  const eb = Math.exp(b * k);
+  return (a * ea + b * eb) / (ea + eb);
+};
+const sclamp = (x, min, max, k) => smin(smax(x, min, k), max, k);
+const absMin = (a, b) => Math.abs(a) < Math.abs(b) ? a : b;
+const absMax = (a, b) => Math.abs(a) > Math.abs(b) ? a : b;
+const foldback = (e, x) => x < -e || x > e ? Math.abs(Math.abs((x - e) % (4 * e)) - 2 * e) - e : x;
+const foldback01 = (x) => (x = Math.abs(x) % 2) > 1 ? 2 - x : x;
+const inRange = (x, min, max) => x >= min && x <= max;
+const inOpenRange = (x, min, max) => x > min && x < max;
+export {
+  absMax,
+  absMin,
+  clamp,
+  clamp0,
+  clamp01,
+  clamp05,
+  clamp11,
+  foldback,
+  foldback01,
+  inOpenRange,
+  inRange,
+  max2id,
+  max3id,
+  max4id,
+  min2id,
+  min3id,
+  min4id,
+  minMax,
+  minNonZero2,
+  minNonZero3,
+  sclamp,
+  smax,
+  smin,
+  wrap,
+  wrap01,
+  wrap11,
+  wrapOnce
 };
-/**
- * Same as `smin(smax(x, min, k), max, k)`.
- *
- * @param x -
- * @param min -
- * @param max -
- * @param k -
- */
-export const sclamp = (x, min, max, k) => smin(smax(x, min, k), max, k);
-export const absMin = (a, b) => (Math.abs(a) < Math.abs(b) ? a : b);
-export const absMax = (a, b) => (Math.abs(a) > Math.abs(b) ? a : b);
-/**
- * If `abs(x) > abs(e)`, recursively mirrors `x` back into `[-e .. +e]`
- * interval at respective positive/negative boundary.
- *
- * @remarks
- * References:
- * - https://www.desmos.com/calculator/lkyf2ag3ta
- * - https://www.musicdsp.org/en/latest/Effects/203-fold-back-distortion.html
- *
- * @param e - threshold (> 0)
- * @param x - input value
- */
-export const foldback = (e, x) => x < -e || x > e ? Math.abs(Math.abs((x - e) % (4 * e)) - 2 * e) - e : x;
-/**
- * Similar to {@link foldback}, but with fixed target range: Folds `x` into the
- * closed [0..1] interval, using infinite internal reflection on either side of
- * the interval.
- *
- * @param x
- */
-export const foldback01 = (x) => ((x = Math.abs(x) % 2) > 1 ? 2 - x : x);
-/**
- * Returns true iff `x` is in closed interval `[min .. max]`
- *
- * @param x -
- * @param min -
- * @param max -
- */
-export const inRange = (x, min, max) => x >= min && x <= max;
-/**
- * Returns true iff `x` is in open interval `(min .. max)`
- *
- * @param x -
- * @param min -
- * @param max -
- */
-export const inOpenRange = (x, min, max) => x > min && x < max;

package/libc.js

@@ -1,109 +1,40 @@
-/**
- * Returns a value with the magnitude of `x` and the sign of `y`.
- *
- * @param x -
- * @param y -
- */
-export const copysign = (x, y) => Math.sign(y) * Math.abs(x);
-/**
- * Returns `2^x`.
- *
- * @param x -
- */
-export const exp2 = (x) => 2 ** x;
-/**
- * Returns the positive difference between `x` and `y`, i.e. `x - y` iff `x > y`,
- * otherwise zero.
- *
- * @param x -
- * @param y -
- */
-export const fdim = (x, y) => Math.max(x - y, 0);
-/**
- * Returns `x * y + z`.
- *
- * @param x -
- * @param y -
- * @param z -
- */
-export const fma = (x, y, z) => x * y + z;
-/**
- * Similar to {@link mod}, {@link remainder}. Returns `x - y * trunc(x / y)`,
- * i.e. essentially the same as JS `%` operator. Result will always have the
- * sign of `x`.
- *
- * @remarks
- * **Caution:** Due to the introduction of libc math functions in v4.0.0 and the
- * resulting name/behavior clashes between the modulo logic in JS, C & GLSL, the
- * previous `fmod` function has been renamed to {@link mod} to align w/ its GLSL
- * version and now exhibits a different behavior to this current {@link fmod}
- * function.
- *
- * Reference: https://www.cplusplus.com/reference/cmath/fmod/
- *
- * @param x -
- * @param y -
- */
-export const fmod = (x, y) => x % y;
-//export const fmod: FnN2 = (x, y) => x - y * Math.trunc(x / y);
-/**
- * Inverse op of {@link ldexp}. Breaks the number `x` into its binary
- * significand (a floating point with an abs value in `[0.5,1.0)` interval and
- * an integral exponent for 2, such that: `x = significand * 2^exp`. Returns
- * tuple of `[sig, exp]`.
- *
- * @remarks
- * - If `x` is zero, both parts (significand and exponent) are zero.
- * - If `x` is negative, the significand returned by this function is negative.
- *
- * Based on:
- * https://github.com/locutusjs/locutus/blob/master/src/c/math/frexp.js
- *
- * @param x -
- */
-export const frexp = (x) => {
-    if (x === 0 || !isFinite(x))
-        return [x, 0];
-    const abs = Math.abs(x);
-    let exp = Math.max(-1023, Math.floor(Math.log2(abs)) + 1);
-    let y = abs * 2 ** -exp;
-    while (y < 0.5) {
-        y *= 2;
-        exp--;
-    }
-    while (y >= 1) {
-        y *= 0.5;
-        exp++;
-    }
-    return [x < 0 ? -y : y, exp];
+const copysign = (x, y) => Math.sign(y) * Math.abs(x);
+const exp2 = (x) => 2 ** x;
+const fdim = (x, y) => Math.max(x - y, 0);
+const fma = (x, y, z) => x * y + z;
+const fmod = (x, y) => x % y;
+const frexp = (x) => {
+  if (x === 0 || !isFinite(x))
+    return [x, 0];
+  const abs = Math.abs(x);
+  let exp = Math.max(-1023, Math.floor(Math.log2(abs)) + 1);
+  let y = abs * 2 ** -exp;
+  while (y < 0.5) {
+    y *= 2;
+    exp--;
+  }
+  while (y >= 1) {
+    y *= 0.5;
+    exp++;
+  }
+  return [x < 0 ? -y : y, exp];
 };
-/**
- * Inverse op of {@link frexp}. Returns `x * 2^exp`
- *
- * @param x -
- * @param exp -
- */
-export const ldexp = (x, exp) => x * 2 ** exp;
-/**
- * Similar to {@link fmod}, {@link mod}. Returns `x - y * round(x / y)`.
- *
- * @remarks
- * https://www.cplusplus.com/reference/cmath/remainder/
- *
- * @param x -
- * @param y -
- */
-export const remainder = (x, y) => x - y * Math.round(x / y);
-/**
- * Computes both the quotient and remainder of the integer division of the
- * numerator `x` by the denominator `y`.
- *
- * @param x -
- * @param y -
- */
-export const ldiv = (x, y) => {
-    x |= 0;
-    y |= 0;
-    const q = (x / y) | 0;
-    return [q, x - q * y];
+const ldexp = (x, exp) => x * 2 ** exp;
+const remainder = (x, y) => x - y * Math.round(x / y);
+const ldiv = (x, y) => {
+  x |= 0;
+  y |= 0;
+  const q = x / y | 0;
+  return [q, x - q * y];
+};
+export {
+  copysign,
+  exp2,
+  fdim,
+  fma,
+  fmod,
+  frexp,
+  ldexp,
+  ldiv,
+  remainder
 };