@danielsimonjr/mathts-functions 0.2.12 → 0.2.13

package/dist/index.js

@@ -459,6 +459,44 @@ __export(typed_exports, {
 // src/typed/arithmetic.ts
 import { mathTyped, Complex, Fraction, BigNumber } from "@danielsimonjr/mathts-core";
 import { computePool } from "@danielsimonjr/mathts-parallel";
+
+// src/wasm/elementwise/wasm-bridge.ts
+var WASM_ELEMENTWISE_THRESHOLD = 1024;
+function getWasm() {
+  try {
+    return wasmLoader.getModule();
+  } catch {
+    return null;
+  }
+}
+var scratchBase = -1;
+function elementwiseUnaryDispatch(op, xs) {
+  const n = xs.length;
+  if (n < WASM_ELEMENTWISE_THRESHOLD) return null;
+  const wasm = getWasm();
+  if (!wasm) return null;
+  const fn = wasm[`simd_${op}_array`];
+  const mem = wasm.memory;
+  if (typeof fn !== "function" || !(mem instanceof WebAssembly.Memory)) return null;
+  try {
+    const need = 2 * n * 8;
+    if (scratchBase < 0) scratchBase = mem.buffer.byteLength;
+    if (scratchBase + need > mem.buffer.byteLength) {
+      mem.grow(Math.ceil((scratchBase + need - mem.buffer.byteLength) / 65536));
+    }
+    const inOff = scratchBase;
+    const outOff = scratchBase + n * 8;
+    new Float64Array(mem.buffer, inOff, n).set(xs);
+    fn(inOff, outOff, n);
+    const result = new Float64Array(n);
+    result.set(new Float64Array(mem.buffer, outOff, n));
+    return result;
+  } catch {
+    return null;
+  }
+}
+
+// src/typed/arithmetic.ts
 var add = mathTyped("add", {
   "number, number": (a, b) => a + b,
   "bigint, bigint": (a, b) => a + b,
@@ -594,6 +632,8 @@ var abs = mathTyped("abs", {
   BigNumber: (a) => a.abs(),
   // Parallel array abs
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("abs", a);
+    if (wasm) return wasm;
     const result = await computePool.abs(a);
     return result.result;
   }
@@ -693,6 +733,8 @@ var exp = mathTyped("exp", {
   BigNumber: (a) => a.exp(),
   // Parallel array exp
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("exp", a);
+    if (wasm) return wasm;
     const result = await computePool.exp(a);
     return result.result;
   }
@@ -704,6 +746,8 @@ var log = mathTyped("log", {
   "number, number": (a, base) => Math.log(a) / Math.log(base),
   // Parallel array log
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("log", a);
+    if (wasm) return wasm;
     const result = await computePool.log(a);
     return result.result;
   }
@@ -1214,6 +1258,8 @@ var sin = mathTyped2("sin", {
   BigNumber: (a) => a.sin(),
   // Parallel Float64Array sin
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("sin", a);
+    if (wasm) return wasm;
     const result = await computePool2.sin(a);
     return result.result;
   }
@@ -1224,6 +1270,8 @@ var cos = mathTyped2("cos", {
   BigNumber: (a) => a.cos(),
   // Parallel Float64Array cos
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("cos", a);
+    if (wasm) return wasm;
     const result = await computePool2.cos(a);
     return result.result;
   }
@@ -1234,6 +1282,8 @@ var tan = mathTyped2("tan", {
   BigNumber: (a) => a.tan(),
   // Parallel Float64Array tan
   Float64Array: async (a) => {
+    const wasm = elementwiseUnaryDispatch("tan", a);
+    if (wasm) return wasm;
     const result = await computePool2.tan(a);
     return result.result;
   }
@@ -1974,7 +2024,7 @@ import { computePool as computePool4 } from "@danielsimonjr/mathts-parallel";
 
 // src/wasm/signal/wasm-bridge.ts
 var WASM_SIGNAL_THRESHOLD = 4096;
-function getWasm() {
+function getWasm2() {
   try {
     return wasmLoader.getModule();
   } catch {
@@ -2184,7 +2234,7 @@ function chirpZTransformJS(samples, m, phiStartRe, phiStartIm, phiStepRe, phiSte
 }
 function applyWindowDispatch(samples, windowName) {
   const wt = windowTypeCode(windowName);
-  const wasm = getWasm();
+  const wasm = getWasm2();
   if (wasm) {
     const rustFn = wasm["apply_window_f64"];
     if (typeof rustFn === "function") {
@@ -2215,7 +2265,7 @@ function welchPSDDispatch(samples, frameLength, overlap, windowName) {
   const n = samples.length;
   const wt = windowTypeCode(windowName);
   if (n >= WASM_SIGNAL_THRESHOLD && n >= frameLength) {
-    const wasm = getWasm();
+    const wasm = getWasm2();
     if (wasm) {
       const rustFn = wasm["welch_psd_f64"];
       if (typeof rustFn === "function") {
@@ -2254,7 +2304,7 @@ function welchPSDDispatch(samples, frameLength, overlap, windowName) {
 function bartlettPSDDispatch(samples, frameLength) {
   const n = samples.length;
   if (n >= WASM_SIGNAL_THRESHOLD && n >= frameLength) {
-    const wasm = getWasm();
+    const wasm = getWasm2();
     if (wasm) {
       const rustFn = wasm["bartlett_psd_f64"];
       if (typeof rustFn === "function") {
@@ -2288,7 +2338,7 @@ function bartlettPSDDispatch(samples, frameLength) {
 function goertzelDispatch(samples, targetFreq, sampleRate) {
   const n = samples.length;
   if (n >= WASM_SIGNAL_THRESHOLD) {
-    const wasm = getWasm();
+    const wasm = getWasm2();
     if (wasm) {
       const rustFn = wasm["goertzel_f64"];
       if (typeof rustFn === "function") {
@@ -2320,7 +2370,7 @@ function goertzelDispatch(samples, targetFreq, sampleRate) {
 function chirpZTransformDispatch(samples, m, phiStartRe, phiStartIm, phiStepRe, phiStepIm) {
   const n = samples.length;
   if (Math.max(n, m) >= WASM_SIGNAL_THRESHOLD) {
-    const wasm = getWasm();
+    const wasm = getWasm2();
     if (wasm) {
       const rustFn = wasm["chirp_z_transform_f64"];
       if (typeof rustFn === "function") {
@@ -3409,7 +3459,7 @@ var BINARY_OPS = {
 var UNARY_OPS = {
   bitNot: { rust: "bitNotArray", as: "bitNot_i32_array" }
 };
-function getWasm2() {
+function getWasm3() {
   try {
     return wasmLoader.getModule();
   } catch {
@@ -3417,7 +3467,7 @@ function getWasm2() {
   }
 }
 function getBinaryKernel(op) {
-  const wasm = getWasm2();
+  const wasm = getWasm3();
   if (!wasm) return null;
   const names = BINARY_OPS[op];
   const rust = wasm[names.rust];
@@ -3431,7 +3481,7 @@ function getBinaryKernel(op) {
   return null;
 }
 function getUnaryKernel(op) {
-  const wasm = getWasm2();
+  const wasm = getWasm3();
   if (!wasm) return null;
   const names = UNARY_OPS[op];
   const rust = wasm[names.rust];
@@ -4120,7 +4170,7 @@ import { computePool as computePool6 } from "@danielsimonjr/mathts-parallel";
 
 // src/wasm/special/wasm-bridge.ts
 var WASM_SPECIAL_THRESHOLD = 1024;
-function getWasm3() {
+function getWasm4() {
   try {
     return wasmLoader.getModule();
   } catch {
@@ -4186,7 +4236,7 @@ function airyBiJS(xs) {
 function besselJ0Dispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_j0_f64"];
@@ -4217,7 +4267,7 @@ function besselJ0Dispatch(xs) {
 function besselJ1Dispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_j1_f64"];
@@ -4250,7 +4300,7 @@ function besselJDispatch(order, xs) {
   if (order === 1) return besselJ1Dispatch(xs);
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_j_f64"];
@@ -4281,7 +4331,7 @@ function besselJDispatch(order, xs) {
 function besselY0Dispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_y0_f64"];
@@ -4312,7 +4362,7 @@ function besselY0Dispatch(xs) {
 function besselY1Dispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_y1_f64"];
@@ -4345,7 +4395,7 @@ function besselYDispatch(order, xs) {
   if (order === 1) return besselY1Dispatch(xs);
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["bessel_y_f64"];
@@ -4376,7 +4426,7 @@ function besselYDispatch(order, xs) {
 function airyAiDispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["airy_ai_f64"];
@@ -4409,7 +4459,7 @@ function airyAiDispatch(xs) {
 function airyBiDispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["airy_bi_f64"];
@@ -4472,7 +4522,7 @@ function lgammaJS(xs) {
 function lgammaDispatch(xs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["lgamma_f64"];
@@ -4709,7 +4759,7 @@ function ellipticPiIncompleteScalar(n, phi, m) {
 function carlsonRCDispatch(xs, ys) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["carlson_rc_f64"];
@@ -4740,7 +4790,7 @@ function carlsonRCDispatch(xs, ys) {
 function carlsonRFDispatch(xs, ys, zs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["carlson_rf_f64"];
@@ -4773,7 +4823,7 @@ function carlsonRFDispatch(xs, ys, zs) {
 function carlsonRDDispatch(xs, ys, zs) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["carlson_rd_f64"];
@@ -4806,7 +4856,7 @@ function carlsonRDDispatch(xs, ys, zs) {
 function carlsonRJDispatch(xs, ys, zs, ps) {
   const n = xs.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["carlson_rj_f64"];
@@ -4841,7 +4891,7 @@ function carlsonRJDispatch(xs, ys, zs, ps) {
 function ellipticFIncompleteDispatch(phis, ms) {
   const n = phis.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["elliptic_f_incomplete_f64"];
@@ -4872,7 +4922,7 @@ function ellipticFIncompleteDispatch(phis, ms) {
 function ellipticEIncompleteDispatch(phis, ms) {
   const n = phis.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["elliptic_e_incomplete_f64"];
@@ -4903,7 +4953,7 @@ function ellipticEIncompleteDispatch(phis, ms) {
 function ellipticPiIncompleteDispatch(ns, phis, ms) {
   const n = ns.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["elliptic_pi_incomplete_f64"];
@@ -4994,7 +5044,7 @@ function ellipticEJS(ms) {
 function ellipticKDispatch(ms) {
   const n = ms.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["elliptic_k_f64"];
@@ -5027,7 +5077,7 @@ function ellipticKDispatch(ms) {
 function ellipticEDispatch(ms) {
   const n = ms.length;
   if (n >= WASM_SPECIAL_THRESHOLD) {
-    const wasm = getWasm3();
+    const wasm = getWasm4();
     if (wasm) {
       try {
         const rustFn = wasm["elliptic_e_f64"];
@@ -7529,7 +7579,7 @@ async function distanceMatrix(points) {
 
 // src/wasm/poly/wasm-bridge.ts
 var WASM_POLY_THRESHOLD = 256;
-function getWasm4() {
+function getWasm5() {
   try {
     return wasmLoader.getModule();
   } catch {
@@ -7574,7 +7624,7 @@ function polyDivModJS(num, den) {
 function polyMulDispatch(a, b) {
   const bigEnough = a.length >= WASM_POLY_THRESHOLD || b.length >= WASM_POLY_THRESHOLD;
   if (bigEnough) {
-    const wasm = getWasm4();
+    const wasm = getWasm5();
     if (wasm) {
       try {
         const rustFn = wasm["poly_mul_f64"];
@@ -7606,7 +7656,7 @@ function polyDivModDispatch(num, den) {
   if (den.length === 0) throw new Error("Division by zero polynomial");
   const bigEnough = num.length >= WASM_POLY_THRESHOLD;
   if (bigEnough) {
-    const wasm = getWasm4();
+    const wasm = getWasm5();
     if (wasm) {
       try {
         const rustFn = wasm["poly_div_mod_f64"];
@@ -7734,7 +7784,7 @@ function discriminantJS(p) {
 function resultantDispatch(p, q) {
   const bigEnough = p.length >= WASM_POLY_THRESHOLD || q.length >= WASM_POLY_THRESHOLD;
   if (bigEnough) {
-    const wasm = getWasm4();
+    const wasm = getWasm5();
     if (wasm) {
       try {
         const rustFn = wasm["poly_resultant_f64"];
@@ -7761,7 +7811,7 @@ function resultantDispatch(p, q) {
 function discriminantDispatch(p) {
   const bigEnough = p.length >= WASM_POLY_THRESHOLD;
   if (bigEnough) {
-    const wasm = getWasm4();
+    const wasm = getWasm5();
     if (wasm) {
       try {
         const rustFn = wasm["poly_discriminant_f64"];
@@ -7882,7 +7932,7 @@ function fitDispatch(xs, ys, degree2, rustName, asName, jsFallback) {
   if (xs.length !== ys.length) throw new Error("xs and ys must have the same length");
   const bigEnough = xs.length >= WASM_POLY_FIT_THRESHOLD;
   if (bigEnough) {
-    const wasm = getWasm4();
+    const wasm = getWasm5();
     if (wasm) {
       try {
         const rustFn = wasm[rustName];
@@ -8881,7 +8931,7 @@ async function romberg(f, a, b, tol = 1e-12, options = {}) {
 // src/wasm/interpolation/wasm-bridge.ts
 var WASM_TRIDIAG_THRESHOLD = 1024;
 var WASM_INTERP_THRESHOLD = 256;
-function getWasm5() {
+function getWasm6() {
   try {
     return wasmLoader.getModule();
   } catch {
@@ -8917,7 +8967,7 @@ function tridiagSolveJS(diag2, lower, upper, rhs) {
 function tridiagSolveDispatch(diag2, lower, upper, rhs) {
   const n = diag2.length;
   if (n >= WASM_TRIDIAG_THRESHOLD) {
-    const wasm = getWasm5();
+    const wasm = getWasm6();
     if (wasm) {
       try {
         const rustFn = wasm["tridiag_solve_f64"];
@@ -8978,7 +9028,7 @@ function dividedDifferenceJS(xs, ys) {
 function dividedDifferenceDispatch(xs, ys) {
   const n = xs.length;
   if (n >= WASM_INTERP_THRESHOLD) {
-    const wasm = getWasm5();
+    const wasm = getWasm6();
     if (wasm) {
       try {
         const rustFn = wasm["divided_difference_f64"];

package/dist/typed/arithmetic.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAIH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAS1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AASlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCAwCd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAyBnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAgDnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAyBjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAQd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAYjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAUd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAWd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAahB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAahB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAmBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAqBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAWd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAMhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAKlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAKjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAKpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAKnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAKlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAmBnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAmBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}
+{"version":3,"file":"arithmetic.d.ts","sourceRoot":"","sources":["../../src/typed/arithmetic.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAIH,OAAO,EAAe,WAAW,EAAE,MAAM,gCAAgC,CAAC;AAU1E,4BAA4B;AAC5B,KAAK,GAAG,GAAG,MAAM,CAAC;AASlB;;;;;;;;;;;;;GAaG;AACH,eAAO,MAAM,GAAG,wCAwCd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAyBnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAgDnB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAyBjB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,UAAU,wCAYrB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAOpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAgBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAQd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAYjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAiBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAIlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAad,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAahB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAahB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAmBhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAKd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAqBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAWd,CAAC;AAEH;;;GAGG;AACH,eAAO,MAAM,IAAI,wCA6Bf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAMhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAKlB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,MAAM,wCAKjB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAKpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAKnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,OAAO,wCAKlB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAcd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,QAAQ,wCAmBnB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAmBd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAiBd,CAAC;AAMH;;;;;;;GAOG;AACH,wBAAsB,MAAM,CAC1B,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,EACV,KAAK,EAAE,GAAG,EACV,CAAC,EAAE,YAAY,EACf,KAAK,EAAE,GAAG,GACT,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,SAAS,CAAC,IAAI,EAAE,YAAY,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,GAAG,OAAO,CAAC,YAAY,CAAC,CAG/F;AAED;;GAEG;AACH,wBAAsB,MAAM,CAC1B,MAAM,EAAE,YAAY,EACpB,IAAI,EAAE,GAAG,EACT,IAAI,EAAE,GAAG,EACT,MAAM,EAAE,YAAY,GACnB,OAAO,CAAC,YAAY,CAAC,CAGvB;AAED;;GAEG;AACH,wBAAsB,KAAK,CAAC,CAAC,EAAE,YAAY,EAAE,CAAC,EAAE,YAAY,GAAG,OAAO,CAAC,YAAY,CAAC,CAGnF;AAMD;;GAEG;AACH,wBAAsB,cAAc,IAAI,OAAO,CAAC,IAAI,CAAC,CAEpD;AAED;;GAEG;AACH,wBAAsB,aAAa,IAAI,OAAO,CAAC,IAAI,CAAC,CAEnD;AAED;;GAEG;AACH,wBAAgB,iBAAiB,CAAC,YAAY,EAAE,GAAG,GAAG,OAAO,CAE5D;AAED;;GAEG;AACH,wBAAgB,cAAc,IAAI,WAAW,CAE5C;AAMD,eAAO,MAAM,eAAe;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;CA+E3B,CAAC"}

package/dist/typed/trigonometry.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"trigonometry.d.ts","sourceRoot":"","sources":["../../src/typed/trigonometry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAiBH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAUd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAUd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAUd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAEhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAKhB,CAAC;AAMH,eAAO,MAAM,iBAAiB;;;;;;;;;;;;;;;;;;;;CA2B7B,CAAC"}
+{"version":3,"file":"trigonometry.d.ts","sourceRoot":"","sources":["../../src/typed/trigonometry.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAkBH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAYd,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,GAAG,wCAed,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAqBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAoBf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAef,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAEhB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,IAAI,wCAGf,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAehB,CAAC;AAMH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,SAAS,wCAEpB,CAAC;AAEH;;GAEG;AACH,eAAO,MAAM,KAAK,wCAKhB,CAAC;AAMH,eAAO,MAAM,iBAAiB;;;;;;;;;;;;;;;;;;;;CA2B7B,CAAC"}

package/dist/wasm/elementwise/wasm-bridge.d.ts

@@ -0,0 +1,13 @@
+/** Element-count threshold above which the WASM kernel beats JS. */
+export declare const WASM_ELEMENTWISE_THRESHOLD = 1024;
+/** Unary ops with a Rust `simd_<op>_array` kernel that is net-faster than JS. */
+export declare const WASM_ELEMENTWISE_OPS: readonly ["abs", "sin", "cos", "tan", "exp", "log"];
+export type WasmElementwiseOp = (typeof WASM_ELEMENTWISE_OPS)[number];
+/**
+ * Apply a unary elementwise op over `xs` via the Rust `simd_<op>_array` kernel.
+ * Returns the result, or `null` (below threshold, wasm unavailable, kernel
+ * missing, or any failure) — in which case the caller uses its JS/parallel path.
+ * Never throws.
+ */
+export declare function elementwiseUnaryDispatch(op: WasmElementwiseOp, xs: Float64Array): Float64Array | null;
+//# sourceMappingURL=wasm-bridge.d.ts.map

package/dist/wasm/elementwise/wasm-bridge.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"wasm-bridge.d.ts","sourceRoot":"","sources":["../../../src/wasm/elementwise/wasm-bridge.ts"],"names":[],"mappings":"AAoBA,oEAAoE;AACpE,eAAO,MAAM,0BAA0B,OAAO,CAAC;AAE/C,iFAAiF;AACjF,eAAO,MAAM,oBAAoB,qDAAsD,CAAC;AACxF,MAAM,MAAM,iBAAiB,GAAG,CAAC,OAAO,oBAAoB,CAAC,CAAC,MAAM,CAAC,CAAC;AAyBtE;;;;;GAKG;AACH,wBAAgB,wBAAwB,CAAC,EAAE,EAAE,iBAAiB,EAAE,EAAE,EAAE,YAAY,GAAG,YAAY,GAAG,IAAI,CA0BrG"}

package/dist/wasm/wasm-manifest.json

@@ -1,3 +1,3 @@
-{
-  "mathts.wasm": "sha384-oONMJGNzrFf5H4dPkE7U6ceOD2fBdFtJbrZitiNZOb/CcMvZVXY3g6ZTB/DtVDKy"
-}
+{
+  "mathts.wasm": "sha384-oONMJGNzrFf5H4dPkE7U6ceOD2fBdFtJbrZitiNZOb/CcMvZVXY3g6ZTB/DtVDKy"
+}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-functions",
-  "version": "0.2.12",
+  "version": "0.2.13",
   "description": "Mathematical functions for MathTS - arithmetic, algebra, trigonometry, statistics, and more",
   "author": "Daniel Simon Jr.",
   "license": "MIT",
@@ -22,7 +22,7 @@
     "build": "tsup src/index.ts --format esm --clean && tsc -p tsconfig.dts.json && node scripts/copy-wasm.mjs",
     "dev": "tsup src/index.ts --format esm --dts --watch",
     "test": "vitest run",
-    "test:diff": "node tests/diff-special.test.mjs",
+    "test:diff": "node tests/diff-special.test.mjs && node tests/diff-elementwise.test.mjs",
     "golden:gen": "python -X utf8 tests/golden/gen_special_goldens.py",
     "test:watch": "vitest",
     "test:coverage": "vitest run --coverage",