@thi.ng/dual-algebra 0.4.45 → 0.4.47

package/CHANGELOG.md

@@ -1,6 +1,6 @@
 # Change Log
 
-- **Last updated**: 2023-12-09T19:12:03Z
+- **Last updated**: 2023-12-18T13:41:19Z
 - **Generator**: [thi.ng/monopub](https://thi.ng/monopub)
 
 All notable changes to this project will be documented in this file.

package/README.md

@@ -138,7 +138,7 @@ For Node.js REPL:
 const dualAlgebra = await import("@thi.ng/dual-algebra");
 ```
 
-Package sizes (brotli'd, pre-treeshake): ESM: 990 bytes
+Package sizes (brotli'd, pre-treeshake): ESM: 991 bytes
 
 ## Dependencies
 

package/api.js

@@ -1 +0,0 @@
-export {};

package/ops.js

@@ -1,185 +1,190 @@
-/**
- * Creates a multivariate dual number.
- *
- * @param real - real value
- * @param n - number of variables (default: 1)
- * @param i - variable index (0 < i <= n)
- */
-export const dual = (real, n = 1, i = 0) => {
-    const out = new Array(n + 1).fill(0, 1);
-    out[0] = real;
-    i > 0 && (out[i] = 1);
-    return out;
+const dual = (real, n = 1, i = 0) => {
+  const out = new Array(n + 1).fill(0, 1);
+  out[0] = real;
+  i > 0 && (out[i] = 1);
+  return out;
 };
-/**
- * Creates a 1-dual number of `r`. Syntax sugar for {@link dual}.
- *
- * @param r -
- * @param i -
- */
-export const $ = (r, i = 0) => [r, i === 1 ? 1 : 0];
-/**
- * Creates a 2-dual number of `r`. Syntax sugar for {@link dual}.
- *
- * @param r -
- * @param i -
- */
-export const $2 = (r, i = 0) => dual(r, 2, i);
-/**
- * Creates a 3-dual number of `r`. Syntax sugar for {@link dual}.
- *
- * @param r -
- * @param i -
- */
-export const $3 = (r, i = 0) => dual(r, 3, i);
-/**
- * Creates a 4-dual number of `r`. Syntax sugar for {@link dual}.
- *
- * @param r -
- * @param i -
- */
-export const $4 = (r, i = 0) => dual(r, 4, i);
-export const defOp = (single, multi, dispatch = 0) => (((...xs) => xs[dispatch].length < 3 ? single(...xs) : multi(...xs)));
-export const add = defOp((a, b) => [a[0] + b[0], a[1] + b[1]], (a, b) => a.map((x, i) => x + b[i]));
-export const sub = defOp((a, b) => [a[0] - b[0], a[1] - b[1]], (a, b) => a.map((x, i) => x - b[i]));
-export const neg = defOp((a) => [-a[0], -a[1]], (a) => a.map((x) => (x !== 0 ? -x : 0)));
-export const mul = defOp(([ar, ad], [br, bd]) => [ar * br, ar * bd + ad * br], (a, b) => {
+const $ = (r, i = 0) => [r, i === 1 ? 1 : 0];
+const $2 = (r, i = 0) => dual(r, 2, i);
+const $3 = (r, i = 0) => dual(r, 3, i);
+const $4 = (r, i = 0) => dual(r, 4, i);
+const defOp = (single, multi, dispatch = 0) => (...xs) => xs[dispatch].length < 3 ? single(...xs) : multi(...xs);
+const add = defOp(
+  (a, b) => [a[0] + b[0], a[1] + b[1]],
+  (a, b) => a.map((x, i) => x + b[i])
+);
+const sub = defOp(
+  (a, b) => [a[0] - b[0], a[1] - b[1]],
+  (a, b) => a.map((x, i) => x - b[i])
+);
+const neg = defOp(
+  (a) => [-a[0], -a[1]],
+  (a) => a.map((x) => x !== 0 ? -x : 0)
+);
+const mul = defOp(
+  ([ar, ad], [br, bd]) => [ar * br, ar * bd + ad * br],
+  (a, b) => {
     const ar = a[0];
     const br = b[0];
     const out = [ar * br];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = ar * b[i] + a[i] * br;
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = ar * b[i] + a[i] * br;
     }
     return out;
-});
-export const div = defOp(([ar, ad], [br, bd]) => [ar / br, (ad * br - ar * bd) / (br * br)], (a, b) => {
+  }
+);
+const div = defOp(
+  ([ar, ad], [br, bd]) => [ar / br, (ad * br - ar * bd) / (br * br)],
+  (a, b) => {
     const ar = a[0];
     const br = b[0];
     const ibr = 1 / (br * br);
     const out = [ar / br];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = (a[i] * br - ar * b[i]) * ibr;
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = (a[i] * br - ar * b[i]) * ibr;
     }
     return out;
-});
-export const abs = defOp(([ar, ad]) => [Math.abs(ar), ad * Math.sign(ar)], (a) => {
+  }
+);
+const abs = defOp(
+  ([ar, ad]) => [Math.abs(ar), ad * Math.sign(ar)],
+  (a) => {
     const s = Math.sign(a[0]);
     const out = [Math.abs(a[0])];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = s * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = s * a[i];
     }
     return out;
-});
-export const sqrt = defOp((a) => {
+  }
+);
+const sqrt = defOp(
+  (a) => {
     const s = Math.sqrt(a[0]);
-    return [s, (0.5 * a[1]) / s];
-}, (a) => {
+    return [s, 0.5 * a[1] / s];
+  },
+  (a) => {
     const s = Math.sqrt(a[0]);
     const si = 0.5 / s;
     const out = [s];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = si * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = si * a[i];
     }
     return out;
-});
-export const exp = defOp(([ar, ad]) => {
+  }
+);
+const exp = defOp(
+  ([ar, ad]) => {
     ar = Math.exp(ar);
     return [ar, ad * ar];
-}, (a) => {
+  },
+  (a) => {
     const ar = Math.exp(a[0]);
     const out = [ar];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = ar * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = ar * a[i];
     }
     return out;
-});
-export const log = defOp(([ar, ad]) => [Math.log(ar), ad / ar], (a) => {
+  }
+);
+const log = defOp(
+  ([ar, ad]) => [Math.log(ar), ad / ar],
+  (a) => {
     const ar = Math.log(a[0]);
     const iar = 1 / ar;
     const out = [ar];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = iar * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = iar * a[i];
     }
     return out;
-});
-export const pow = defOp(([ar, ad], k) => [ar ** k, ad * k * ar ** (k - 1)], (a, k) => {
+  }
+);
+const pow = defOp(
+  ([ar, ad], k) => [ar ** k, ad * k * ar ** (k - 1)],
+  (a, k) => {
     const f = k * a[0] ** (k - 1);
     const out = [a[0] ** k];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = f * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = f * a[i];
     }
     return out;
-});
-export const sin = defOp(([ar, ad]) => [Math.sin(ar), ad * Math.cos(ar)], (a) => {
+  }
+);
+const sin = defOp(
+  ([ar, ad]) => [Math.sin(ar), ad * Math.cos(ar)],
+  (a) => {
     const c = Math.cos(a[0]);
     const out = [Math.sin(a[0])];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = c * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = c * a[i];
     }
     return out;
-});
-export const cos = defOp(([ar, ad]) => [Math.cos(ar), -ad * Math.sin(ar)], (a) => {
+  }
+);
+const cos = defOp(
+  ([ar, ad]) => [Math.cos(ar), -ad * Math.sin(ar)],
+  (a) => {
     const s = -Math.sin(a[0]);
     const out = [Math.cos(a[0])];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = s * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = s * a[i];
     }
     return out;
-});
-export const tan = defOp(([ar, ad]) => {
+  }
+);
+const tan = defOp(
+  ([ar, ad]) => {
     const c = Math.cos(ar);
     return [Math.tan(ar), ad / (c * c)];
-}, (a) => {
+  },
+  (a) => {
     const c = Math.cos(a[0]);
     const ic = 1 / (c * c);
     const out = [Math.tan(a[0])];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = ic * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = ic * a[i];
     }
     return out;
-});
-export const atan = defOp(([ar, ad]) => [Math.atan(ar), ad / (1 + ar * ar)], (a) => {
+  }
+);
+const atan = defOp(
+  ([ar, ad]) => [Math.atan(ar), ad / (1 + ar * ar)],
+  (a) => {
     const ar = a[0];
     const iar = 1 / (1 + ar * ar);
     const out = [Math.atan(ar)];
-    for (let i = a.length; i-- > 1;) {
-        out[i] = iar * a[i];
+    for (let i = a.length; i-- > 1; ) {
+      out[i] = iar * a[i];
     }
     return out;
-});
-/**
- * Linear interpolation for dual numbers: `a + (b - a) * t`
- *
- * @param a -
- * @param b -
- * @param t -
- */
-export const mix = (a, b, t) => add(a, mul(sub(b, a), t));
-/**
- * Higher order function. Takes a 2-multivariate {@link Op2} and returns new
- * function which takes two real numbers `x` and `y` representing variables in
- * the given function. When called, converts `x` and `y` first into dual numbers
- * and then calls `fn` and returns result.
- *
- * @remarks
- * The result tuple can be interpreted the following:
- *
- * - index 0 - real value of `fn`
- * - index 1 - derivative of first var at `x`
- * - index 2 - derivative of second var at `y`
- *
- * @param fn -
- */
-export const evalFn2 = (fn) => (x, y) => fn([x, 1, 0], [y, 0, 1]);
-/**
- * Same as {@link evalFn2}, but 3-multivariate functions.
- *
- * @param fn -
- */
-export const evalFn3 = (fn) => (x, y, z) => fn([x, 1, 0, 0], [y, 0, 1, 0], [z, 0, 0, 1]);
-/**
- * Same as {@link evalFn4}, but 4-multivariate functions.
- *
- * @param fn -
- */
-export const evalFn4 = (fn) => (x, y, z, w) => fn([x, 1, 0, 0, 0], [y, 0, 1, 0, 0], [z, 0, 0, 1, 0], [w, 0, 0, 0, 1]);
+  }
+);
+const mix = (a, b, t) => add(a, mul(sub(b, a), t));
+const evalFn2 = (fn) => (x, y) => fn([x, 1, 0], [y, 0, 1]);
+const evalFn3 = (fn) => (x, y, z) => fn([x, 1, 0, 0], [y, 0, 1, 0], [z, 0, 0, 1]);
+const evalFn4 = (fn) => (x, y, z, w) => fn([x, 1, 0, 0, 0], [y, 0, 1, 0, 0], [z, 0, 0, 1, 0], [w, 0, 0, 0, 1]);
+export {
+  $,
+  $2,
+  $3,
+  $4,
+  abs,
+  add,
+  atan,
+  cos,
+  defOp,
+  div,
+  dual,
+  evalFn2,
+  evalFn3,
+  evalFn4,
+  exp,
+  log,
+  mix,
+  mul,
+  neg,
+  pow,
+  sin,
+  sqrt,
+  sub,
+  tan
+};

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@thi.ng/dual-algebra",
-  "version": "0.4.45",
+  "version": "0.4.47",
   "description": "Multivariate dual number algebra, automatic differentiation",
   "type": "module",
   "module": "./index.js",
@@ -24,7 +24,9 @@
   "author": "Karsten Schmidt (https://thi.ng)",
   "license": "Apache-2.0",
   "scripts": {
-    "build": "yarn clean && tsc --declaration",
+    "build": "yarn build:esbuild && yarn build:decl",
+    "build:decl": "tsc --declaration --emitDeclarationOnly",
+    "build:esbuild": "esbuild --format=esm --platform=neutral --target=es2022 --tsconfig=tsconfig.json --outdir=. src/**/*.ts",
     "clean": "rimraf --glob '*.js' '*.d.ts' '*.map' doc",
     "doc": "typedoc --excludePrivate --excludeInternal --out doc src/index.ts",
     "doc:ae": "mkdir -p .ae/doc .ae/temp && api-extractor run --local --verbose",
@@ -33,10 +35,11 @@
     "test": "bun test"
   },
   "dependencies": {
-    "@thi.ng/api": "^8.9.11"
+    "@thi.ng/api": "^8.9.13"
   },
   "devDependencies": {
     "@microsoft/api-extractor": "^7.38.3",
+    "esbuild": "^0.19.8",
     "rimraf": "^5.0.5",
     "tools": "^0.0.1",
     "typedoc": "^0.25.4",
@@ -59,7 +62,7 @@
     "access": "public"
   },
   "engines": {
-    "node": ">=12.7"
+    "node": ">=18"
   },
   "files": [
     "./*.js",
@@ -89,5 +92,5 @@
     "status": "alpha",
     "year": 2020
   },
-  "gitHead": "25f2ac8ff795a432a930119661b364d4d93b59a0\n"
+  "gitHead": "25a42a81fac8603a1e440a7aa8bc343276211ff4\n"
 }

package/poly.js

@@ -1,71 +1,25 @@
 import { add, mul } from "./ops.js";
-/**
- * Computes: `ax^2 + bx + c`. All args must have same size/arity.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- */
-export const quadratic = (x, a, b, c) => add(add(mul(a, mul(x, x)), mul(b, x)), c);
-/**
- * Same as {@link quadratic}, but for real/scalar inputs. `x` is treated as
- * variable `x+1ε`, the rest as `n+0ε`.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- */
-export const quadraticS = (x, a, b, c) => quadratic([x, 1], [a, 0], [b, 0], [c, 0]);
-/**
- * Computes: `ax^3 + bx^2 + cx + d`. All args must have same size/arity.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- * @param d -
- */
-export const cubic = (x, a, b, c, d) => {
-    const x2 = mul(x, x);
-    return add(add(add(mul(a, mul(x2, x)), mul(b, x2)), mul(c, x)), d);
+const quadratic = (x, a, b, c) => add(add(mul(a, mul(x, x)), mul(b, x)), c);
+const quadraticS = (x, a, b, c) => quadratic([x, 1], [a, 0], [b, 0], [c, 0]);
+const cubic = (x, a, b, c, d) => {
+  const x2 = mul(x, x);
+  return add(add(add(mul(a, mul(x2, x)), mul(b, x2)), mul(c, x)), d);
 };
-/**
- * Same as {@link cubic}, but for real/scalar inputs. `x` is treated as variable
- * `x+1ε`, the rest as `n+0ε`.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- * @param d -
- */
-export const cubicS = (x, a, b, c, d) => cubic([x, 1], [a, 0], [b, 0], [c, 0], [d, 0]);
-/**
- * Computes: `ax^4 + bx^3 + cx^2 + dx + e`. All args must have same size/arity.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- * @param d -
- * @param e -
- */
-export const quartic = (x, a, b, c, d, e) => {
-    const x2 = mul(x, x);
-    const x3 = mul(x2, x);
-    return add(add(add(add(mul(a, mul(x3, x)), mul(b, x3)), mul(c, x2)), mul(d, x)), e);
+const cubicS = (x, a, b, c, d) => cubic([x, 1], [a, 0], [b, 0], [c, 0], [d, 0]);
+const quartic = (x, a, b, c, d, e) => {
+  const x2 = mul(x, x);
+  const x3 = mul(x2, x);
+  return add(
+    add(add(add(mul(a, mul(x3, x)), mul(b, x3)), mul(c, x2)), mul(d, x)),
+    e
+  );
+};
+const quarticS = (x, a, b, c, d, e) => quartic([x, 1], [a, 0], [b, 0], [c, 0], [d, 0], [e, 0]);
+export {
+  cubic,
+  cubicS,
+  quadratic,
+  quadraticS,
+  quartic,
+  quarticS
 };
-/**
- * Same as {@link quartic}, but for real/scalar inputs. `x` is treated as
- * variable `x+1ε`, the rest as `n+0ε`.
- *
- * @param x -
- * @param a -
- * @param b -
- * @param c -
- * @param d -
- * @param e -
- */
-export const quarticS = (x, a, b, c, d, e) => quartic([x, 1], [a, 0], [b, 0], [c, 0], [d, 0], [e, 0]);

package/vector.js

@@ -1,41 +1,14 @@
 import { add, div, dual, mul, sub } from "./ops.js";
-const defVecOp2 = (op) => (a, b) => a.map((a, i) => op(a, b[i]));
-/**
- * Dual vector addition. Applies {@link add} in a component-wise manner. Returns
- * new (dual) vector.
- *
- * @param a -
- * @param b -
- */
-export const vadd = defVecOp2(add);
-/**
- * Dual vector subtraction. Applies {@link sub} in a component-wise manner.
- * Returns new (dual) vector.
- *
- * @param a -
- * @param b -
- */
-export const vsub = defVecOp2(sub);
-/**
- * Dual vector multiplication. Applies {@link mul} in a component-wise manner.
- * Returns new (dual) vector.
- *
- * @param a -
- * @param b -
- */
-export const vmul = defVecOp2(mul);
-/**
- * Dual vector division. Applies {@link div} in a component-wise manner.
- * Returns new (dual) vector.
- *
- * @param a -
- * @param b -
- */
-export const vdiv = defVecOp2(div);
-/**
- * Computes dot product of 2 dual vectors.
- *
- * @param a -
- * @param b -
- */
-export const dot = (a, b) => a.reduce((acc, a, i) => add(acc, mul(a, b[i])), dual(0, a[0].length));
+const defVecOp2 = (op) => (a, b) => a.map((a2, i) => op(a2, b[i]));
+const vadd = defVecOp2(add);
+const vsub = defVecOp2(sub);
+const vmul = defVecOp2(mul);
+const vdiv = defVecOp2(div);
+const dot = (a, b) => a.reduce((acc, a2, i) => add(acc, mul(a2, b[i])), dual(0, a[0].length));
+export {
+  dot,
+  vadd,
+  vdiv,
+  vmul,
+  vsub
+};