@danielsimonjr/mathts-matrix 0.6.2 → 0.7.0

package/dist/index.d.ts

@@ -2278,6 +2278,23 @@ interface LUResult {
  * @throws {Error} if `A` is singular (zero pivot encountered).
  */
 declare function lu(A: DenseMatrix): LUResult;
+/**
+ * Solve the dense linear system `A · x = b` using a precomputed LU factorisation
+ * (from {@link lu}), where `P · A = L · U`.
+ *
+ * Applies the row permutation to `b` (`(P·b)[i] = b[P[i]]`), forward-substitutes
+ * the unit-lower-triangular `L`, then back-substitutes the upper-triangular `U`.
+ * Each call is O(n²), so factor **once** with `lu()` and call `luSolve` repeatedly
+ * to solve many right-hand sides against the same matrix — the pattern the stiff
+ * Rosenbrock/RODAS ODE steps use (the iteration matrix is solved against 3–6
+ * successive RHS per step).
+ *
+ * @param fac An `LUResult` produced by {@link lu}.
+ * @param b   Right-hand side vector of length `n` (`number[]` or `Float64Array`).
+ * @returns   The solution vector `x` as a `number[]` of length `n`.
+ * @throws {Error} if `b.length` does not equal the factor dimension `n`.
+ */
+declare function luSolve(fac: LUResult, b: ArrayLike<number>): number[];
 
 /**
  * Cholesky Decomposition (DenseMatrix primitive)
@@ -3024,4 +3041,4 @@ declare function initializeParallelMatrix(): Promise<void>;
  */
 declare function terminateParallelMatrix(): Promise<void>;
 
-export { BUILTIN_SHADERS, type BackendHints, BackendManager, BackendRegistry, type BackendType, BatchExecutor, type CholeskyResult, DEFAULT_BACKEND_HINTS, DEFAULT_EXTENDED_HINTS, DenseMatrix, type EigOptions, type EigResult, type ExpmOptions, type ExtendedBackendHints, GPUBackend, type GPUBackendOptions, type GPUBackendStatus, GPUMatrixBackend, type GPUMatrixBackendConfig, JSBackend, type LQResult, type LUResult, type LogmOptions, Matrix, type MatrixBackend, type MatrixEntry, type MatrixIndex, type MatrixType, type OperationType, ParallelBackend, type ParallelBackendConfig, type PinvOptions, type QLResult, type QROptions, type QRPivotedOptions, type QRPivotedResult, type QRResult, type RQResult, type SVDOptions, type SVDResult, type SchurOptions, type SchurResult, type SliceSpec, SparseMatrix, type SqrtmOptions, type SyncConfig, SyncManager, type SyncStrategy, WASMBackend, type WASMBackendConfig, type WasmFeatures, abs, add, backendManager, backendRegistry, cholesky, clearFeatureCache, column, cond, condest, createBackendManager, createGPUMatrixBackend, createParallelBackend, createSyncManager, createWASMBackend, destroyGlobalGPUBackend, detectWasmFeatures, diag, diagonal, divide, dotMultiply, eig, eigWasm, eigvals, eigvalsWasm, exp, getCachedFeatures, getGlobalGPUBackend, gpuMatrixBackend, identity, initializeGlobalGPUBackend, initializeParallelMatrix, isAtomicsAvailable, isDenseMatrix, isMatrix, isSharedMemoryAvailable, isSparseMatrix, isWasmAvailable, jsBackend, log, lowRankApprox, lq, lu, matrix, matrixExpm, matrixLogm, pinv as matrixPinv, matrixSchur, matrixSqrtm, max, mean, min, multiply, norm, norm2, normFro, ones, parallelBackend, parallelDiag, parallelDotMultiply, parallelIdentity, parallelMatrix, parallelMatrixAbs, parallelMatrixAdd, parallelMatrixColumn, parallelMatrixCos, parallelMatrixDiagonal, parallelMatrixDistance, parallelMatrixDivide, parallelMatrixDot, parallelMatrixExp, parallelMatrixHistogram, parallelMatrixLog, parallelMatrixMatvec, parallelMatrixMax, parallelMatrixMean, parallelMatrixMin, parallelMatrixMultiply, parallelMatrixNorm, parallelMatrixOperations, parallelMatrixOuter, parallelMatrixRow, parallelMatrixSin, parallelMatrixSize, parallelMatrixSqrt, parallelMatrixSquare, parallelMatrixStd, parallelMatrixSubset, parallelMatrixSubtract, parallelMatrixSum, parallelMatrixTan, parallelMatrixTrace, parallelMatrixTranspose, parallelMatrixVariance, parallelOnes, parallelRandom, parallelUnaryMinus, parallelZeros, pinv$1 as pinv, pow, powerIteration, ql, qr, qrPivoted, random, row, rq, singularValues, size, spectralRadiusWasm, sqrt, square, subset, subtract, sum, svd, svdWasm, terminateParallelMatrix, trace, transpose, typedMatrixOperations, unaryMinus, wasmBackend, zeros };
+export { BUILTIN_SHADERS, type BackendHints, BackendManager, BackendRegistry, type BackendType, BatchExecutor, type CholeskyResult, DEFAULT_BACKEND_HINTS, DEFAULT_EXTENDED_HINTS, DenseMatrix, type EigOptions, type EigResult, type ExpmOptions, type ExtendedBackendHints, GPUBackend, type GPUBackendOptions, type GPUBackendStatus, GPUMatrixBackend, type GPUMatrixBackendConfig, JSBackend, type LQResult, type LUResult, type LogmOptions, Matrix, type MatrixBackend, type MatrixEntry, type MatrixIndex, type MatrixType, type OperationType, ParallelBackend, type ParallelBackendConfig, type PinvOptions, type QLResult, type QROptions, type QRPivotedOptions, type QRPivotedResult, type QRResult, type RQResult, type SVDOptions, type SVDResult, type SchurOptions, type SchurResult, type SliceSpec, SparseMatrix, type SqrtmOptions, type SyncConfig, SyncManager, type SyncStrategy, WASMBackend, type WASMBackendConfig, type WasmFeatures, abs, add, backendManager, backendRegistry, cholesky, clearFeatureCache, column, cond, condest, createBackendManager, createGPUMatrixBackend, createParallelBackend, createSyncManager, createWASMBackend, destroyGlobalGPUBackend, detectWasmFeatures, diag, diagonal, divide, dotMultiply, eig, eigWasm, eigvals, eigvalsWasm, exp, getCachedFeatures, getGlobalGPUBackend, gpuMatrixBackend, identity, initializeGlobalGPUBackend, initializeParallelMatrix, isAtomicsAvailable, isDenseMatrix, isMatrix, isSharedMemoryAvailable, isSparseMatrix, isWasmAvailable, jsBackend, log, lowRankApprox, lq, lu, luSolve, matrix, matrixExpm, matrixLogm, pinv as matrixPinv, matrixSchur, matrixSqrtm, max, mean, min, multiply, norm, norm2, normFro, ones, parallelBackend, parallelDiag, parallelDotMultiply, parallelIdentity, parallelMatrix, parallelMatrixAbs, parallelMatrixAdd, parallelMatrixColumn, parallelMatrixCos, parallelMatrixDiagonal, parallelMatrixDistance, parallelMatrixDivide, parallelMatrixDot, parallelMatrixExp, parallelMatrixHistogram, parallelMatrixLog, parallelMatrixMatvec, parallelMatrixMax, parallelMatrixMean, parallelMatrixMin, parallelMatrixMultiply, parallelMatrixNorm, parallelMatrixOperations, parallelMatrixOuter, parallelMatrixRow, parallelMatrixSin, parallelMatrixSize, parallelMatrixSqrt, parallelMatrixSquare, parallelMatrixStd, parallelMatrixSubset, parallelMatrixSubtract, parallelMatrixSum, parallelMatrixTan, parallelMatrixTrace, parallelMatrixTranspose, parallelMatrixVariance, parallelOnes, parallelRandom, parallelUnaryMinus, parallelZeros, pinv$1 as pinv, pow, powerIteration, ql, qr, qrPivoted, random, row, rq, singularValues, size, spectralRadiusWasm, sqrt, square, subset, subtract, sum, svd, svdWasm, terminateParallelMatrix, trace, transpose, typedMatrixOperations, unaryMinus, wasmBackend, zeros };

package/dist/index.js

@@ -2499,7 +2499,7 @@ var WASMBackend = class {
    */
   async resolveAsWasmPath() {
     const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
-    const { resolvePackagedWasm, defaultWasmLocation } = await import("./resolve-VJX3YTS6.js");
+    const { resolvePackagedWasm, defaultWasmLocation } = await import("./resolve-DDCEWHWX.js");
     if (isNode) {
       const found = await resolvePackagedWasm(import.meta.url, "mathts-as.wasm");
       if (found) return found;
@@ -2525,7 +2525,7 @@ var WASMBackend = class {
       Date
     };
     const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
-    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-X3CLNVGW.js");
+    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-LIWCCDIJ.js");
     let instance;
     if (isNode) {
       const fs = await import("fs");
@@ -4289,7 +4289,7 @@ var WasmLoader = class _WasmLoader {
     if (this.compiledModule) return;
     const path = wasmPath || await this.getDefaultWasmPath();
     const startTime = performance.now();
-    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-X3CLNVGW.js");
+    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-LIWCCDIJ.js");
     if (this.isNode) {
       const fs = await import("fs");
       const { promisify } = await import("util");
@@ -4350,20 +4350,20 @@ var WasmLoader = class _WasmLoader {
   async getDefaultWasmPath() {
     const wasmFile = "mathts-as.wasm";
     if (this.isNode) {
-      const { resolvePackagedWasm } = await import("./resolve-VJX3YTS6.js");
+      const { resolvePackagedWasm } = await import("./resolve-DDCEWHWX.js");
       const found = await resolvePackagedWasm(import.meta.url, wasmFile);
       if (found) return found;
-      const { defaultWasmLocation: defaultWasmLocation2 } = await import("./resolve-VJX3YTS6.js");
+      const { defaultWasmLocation: defaultWasmLocation2 } = await import("./resolve-DDCEWHWX.js");
       return defaultWasmLocation2(import.meta.url, wasmFile);
     }
-    const { defaultWasmLocation } = await import("./resolve-VJX3YTS6.js");
+    const { defaultWasmLocation } = await import("./resolve-DDCEWHWX.js");
     return defaultWasmLocation(import.meta.url, wasmFile, { browser: true });
   }
   async loadNodeWasm(path, totalStart) {
     const fs = await import("fs");
     const { promisify } = await import("util");
     const readFile = promisify(fs.readFile);
-    const { verifyWasmIntegrity } = await import("./integrity-X3CLNVGW.js");
+    const { verifyWasmIntegrity } = await import("./integrity-LIWCCDIJ.js");
     const readStart = performance.now();
     const buffer = await readFile(path);
     const readEnd = performance.now();
@@ -4384,7 +4384,7 @@ var WasmLoader = class _WasmLoader {
     return instance.exports;
   }
   async loadBrowserWasm(path, totalStart) {
-    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-X3CLNVGW.js");
+    const { loadWasmManifest, verifyWasmIntegrity } = await import("./integrity-LIWCCDIJ.js");
     const manifest = await loadWasmManifest(path);
     if (!manifest && typeof WebAssembly.instantiateStreaming === "function") {
       const instStart2 = performance.now();
@@ -6492,6 +6492,30 @@ function lu(A) {
     P: perm
   };
 }
+function luSolve(fac, b) {
+  const n = fac.P.length;
+  if (b.length !== n) {
+    throw new Error(
+      `luSolve: right-hand side length ${b.length} does not match factor dimension ${n}`
+    );
+  }
+  const L = fac.L.toFloat64Array();
+  const U = fac.U.toFloat64Array();
+  const P = fac.P;
+  const y = new Array(n);
+  for (let i = 0; i < n; i++) {
+    let s = b[P[i]];
+    for (let k = 0; k < i; k++) s -= L[i * n + k] * y[k];
+    y[i] = s;
+  }
+  const x = new Array(n);
+  for (let i = n - 1; i >= 0; i--) {
+    let s = y[i];
+    for (let k = i + 1; k < n; k++) s -= U[i * n + k] * x[k];
+    x[i] = s / U[i * n + i];
+  }
+  return x;
+}
 
 // src/operations/cholesky.ts
 init_DenseMatrix();
@@ -7556,7 +7580,7 @@ function ql(A) {
 }
 
 // src/operations/condest.ts
-function luSolve(L, U, perm, n, b) {
+function luSolve2(L, U, perm, n, b) {
   const pb = new Array(n);
   for (let i = 0; i < n; i++) pb[i] = b[perm[i]];
   const y = new Array(n);
@@ -7646,7 +7670,7 @@ function condest(A, p = 1) {
     return Infinity;
   }
   const norm1Inv = hagerNormEstimate(
-    (v) => luSolve(L, U, perm, n, v),
+    (v) => luSolve2(L, U, perm, n, v),
     (v) => luSolveTranspose(L, U, perm, n, v),
     n
   );
@@ -8544,6 +8568,7 @@ export {
   lowRankApprox,
   lq,
   lu,
+  luSolve,
   matrix,
   matrixExpm,
   matrixLogm,

package/dist/integrity-LIWCCDIJ.js

@@ -0,0 +1,13 @@
+import "./chunk-4VMDO6W2.js";
+
+// src/backends/wasm/integrity.ts
+import {
+  sha384OfBuffer,
+  loadWasmManifest,
+  verifyWasmIntegrity
+} from "@danielsimonjr/mathts-core/internal";
+export {
+  loadWasmManifest,
+  sha384OfBuffer,
+  verifyWasmIntegrity
+};

package/dist/resolve-DDCEWHWX.js

@@ -0,0 +1,8 @@
+import "./chunk-4VMDO6W2.js";
+
+// src/backends/wasm/resolve.ts
+import { resolvePackagedWasm, defaultWasmLocation } from "@danielsimonjr/mathts-core/internal";
+export {
+  defaultWasmLocation,
+  resolvePackagedWasm
+};

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-matrix",
-  "version": "0.6.2",
+  "version": "0.7.0",
   "description": "Matrix operations for MathTS with WASM/WebGPU backend support",
   "author": "Daniel Simon Jr.",
   "license": "MIT",
@@ -32,9 +32,9 @@
     "build:prod": "tsup src/index.ts --format esm --dts --clean --minify --treeshake && node scripts/copy-wasm.mjs"
   },
   "dependencies": {
-    "@danielsimonjr/mathts-core": "^0.12.0",
+    "@danielsimonjr/mathts-core": "^0.13.0",
     "@danielsimonjr/mathts-gpu": "^0.2.0",
-    "@danielsimonjr/mathts-parallel": "^0.6.2"
+    "@danielsimonjr/mathts-parallel": "^0.6.3"
   },
   "devDependencies": {
     "@types/node": "^25.5.2",

package/dist/integrity-X3CLNVGW.js

@@ -1,78 +0,0 @@
-import "./chunk-4VMDO6W2.js";
-
-// src/backends/wasm/integrity.ts
-var ALGO = "sha384";
-async function sha384OfBuffer(buffer) {
-  const bytes = buffer instanceof Uint8Array ? buffer : new Uint8Array(buffer);
-  const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
-  let digest;
-  if (isNode) {
-    const { createHash } = await import("crypto");
-    const hash = createHash(ALGO);
-    hash.update(bytes);
-    digest = new Uint8Array(hash.digest());
-  } else {
-    const copy = new ArrayBuffer(bytes.byteLength);
-    new Uint8Array(copy).set(bytes);
-    const out = await crypto.subtle.digest("SHA-384", copy);
-    digest = new Uint8Array(out);
-  }
-  let binary = "";
-  for (let i = 0; i < digest.length; i++) binary += String.fromCharCode(digest[i]);
-  const b64 = typeof btoa === "function" ? btoa(binary) : Buffer.from(digest).toString("base64");
-  return `sha384-${b64}`;
-}
-async function loadWasmManifest(wasmPath) {
-  const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
-  const manifestPath = wasmPath.replace(/[^/\\]+$/, "wasm-manifest.json");
-  try {
-    if (isNode) {
-      const fs = await import("fs");
-      const { promisify } = await import("util");
-      const readFile = promisify(fs.readFile);
-      const text = await readFile(manifestPath, "utf8");
-      return JSON.parse(text);
-    } else {
-      const res = await fetch(manifestPath);
-      if (!res.ok) return null;
-      return await res.json();
-    }
-  } catch {
-    return null;
-  }
-}
-async function verifyWasmIntegrity(buffer, wasmPath, options = {}) {
-  const fileName = wasmPath.split(/[/\\]/).pop() || wasmPath;
-  const manifest = options.manifest !== void 0 ? options.manifest : await loadWasmManifest(wasmPath);
-  if (!manifest) {
-    if (options.required) {
-      throw new Error(
-        `WASM integrity check failed: no manifest at sibling of "${wasmPath}" (required)`
-      );
-    }
-    if (typeof console !== "undefined") {
-      console.warn(
-        `[wasm-integrity] no manifest found beside "${wasmPath}"; skipping SHA-384 verification (set options.required=true to fail closed)`
-      );
-    }
-    return;
-  }
-  const expected = manifest[fileName];
-  if (!expected) {
-    if (options.required) {
-      throw new Error(`WASM integrity check failed: manifest has no entry for "${fileName}"`);
-    }
-    return;
-  }
-  const actual = await sha384OfBuffer(buffer);
-  if (actual !== expected) {
-    throw new Error(
-      `WASM integrity check failed for "${fileName}": expected ${expected}, got ${actual}`
-    );
-  }
-}
-export {
-  loadWasmManifest,
-  sha384OfBuffer,
-  verifyWasmIntegrity
-};

package/dist/resolve-VJX3YTS6.js

@@ -1,43 +0,0 @@
-import "./chunk-4VMDO6W2.js";
-
-// src/backends/wasm/resolve.ts
-async function resolvePackagedWasm(metaUrl, wasmFile) {
-  const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
-  if (!isNode) return null;
-  const { fileURLToPath } = await import("url");
-  const { dirname, join } = await import("path");
-  const { existsSync } = await import("fs");
-  let dir = dirname(fileURLToPath(metaUrl));
-  for (let depth = 0; depth < 8; depth++) {
-    const candidates = [join(dir, "wasm", wasmFile), join(dir, "dist", "wasm", wasmFile)];
-    for (const candidate of candidates) {
-      if (existsSync(candidate)) return candidate;
-    }
-    const parent = dirname(dir);
-    if (parent === dir) break;
-    dir = parent;
-  }
-  return null;
-}
-async function defaultWasmLocation(metaUrl, wasmFile, opts) {
-  if (opts?.browser) {
-    return new URL(`./wasm/${wasmFile}`, metaUrl).href;
-  }
-  const { fileURLToPath } = await import("url");
-  const { dirname, join } = await import("path");
-  const { existsSync } = await import("fs");
-  let dir = dirname(fileURLToPath(metaUrl));
-  for (let depth = 0; depth < 8; depth++) {
-    if (existsSync(join(dir, "package.json"))) {
-      return join(dir, "dist", "wasm", wasmFile);
-    }
-    const parent = dirname(dir);
-    if (parent === dir) break;
-    dir = parent;
-  }
-  return join(dirname(fileURLToPath(metaUrl)), "wasm", wasmFile);
-}
-export {
-  defaultWasmLocation,
-  resolvePackagedWasm
-};