@danielsimonjr/mathts-matrix 0.1.14 → 0.2.1

package/dist/index.d.ts

@@ -1125,7 +1125,7 @@ declare function getCachedFeatures(): WasmFeatures | null;
  * WASM Matrix Backend (AssemblyScript)
  *
  * Implements MatrixBackend using the AssemblyScript-compiled WebAssembly
- * module (`lib/wasm/mathts-as.wasm`, ~42 KB). Falls back to JSBackend when
+ * module (`dist/wasm/mathts-as.wasm`, ~65 KB). Falls back to JSBackend when
  * the WASM module is unavailable or for small matrices.
  *
  * ABI:
@@ -1178,7 +1178,7 @@ declare function getCachedFeatures(): WasmFeatures | null;
 interface WASMBackendConfig {
     /** Minimum elements to use WASM (default: 100) */
     minElements?: number;
-    /** Path to WASM file (defaults to the AS artifact, `lib/wasm/mathts-as.wasm`) */
+    /** Path to WASM file (defaults to the packaged AS artifact, `dist/wasm/mathts-as.wasm`) */
     wasmPath?: string;
     /** Enable SIMD code paths if the AS module exposes them (currently a no-op
      * — AS export surface does not include `simd*` ops; kept for API stability). */
@@ -1203,11 +1203,7 @@ declare class WASMBackend implements MatrixBackend {
     initialize(): Promise<void>;
     private doInitialize;
     /**
-     * Resolve the path to the AssemblyScript artifact (`lib/wasm/mathts-as.wasm`).
-     *
-     * On Windows, `URL.pathname` returns "/C:/foo/bar.wasm" which fs.readFile
-     * interprets as drive-relative ("C:\C:\foo\..."). `fileURLToPath` does the
-     * platform-correct conversion.
+     * Resolve the path to the AssemblyScript artifact (`dist/wasm/mathts-as.wasm`).
      */
     private resolveAsWasmPath;
     /**
@@ -2478,18 +2474,17 @@ interface SVDOptions {
     maxIterations?: number;
     /** Convergence tolerance */
     tolerance?: number;
-    /** Whether to compute full U and V matrices */
+    /**
+     * Complete the thin factor to a square orthonormal basis (U → m×m for tall
+     * inputs, V → n×n for wide ones), numpy `full_matrices=True` style, so the
+     * rectangular-Σ reconstruction `A = U·Σ·Vᵀ` holds. Default `false` = thin
+     * factors (the library's long-standing behavior; the option was previously
+     * accepted but ignored — B-4).
+     */
     fullMatrices?: boolean;
     /** Threshold for rank determination */
     rankTolerance?: number;
 }
-/**
- * Compute SVD of a matrix
- *
- * @param matrix - Input matrix (m x n)
- * @param options - Computation options
- * @returns SVD decomposition
- */
 declare function svd(matrix: number[][] | Float64Array, options?: SVDOptions): SVDResult;
 /**
  * Compute only singular values (faster than full SVD)

package/dist/index.js

@@ -2481,22 +2481,17 @@ var WASMBackend = class {
     }
   }
   /**
-   * Resolve the path to the AssemblyScript artifact (`lib/wasm/mathts-as.wasm`).
-   *
-   * On Windows, `URL.pathname` returns "/C:/foo/bar.wasm" which fs.readFile
-   * interprets as drive-relative ("C:\C:\foo\..."). `fileURLToPath` does the
-   * platform-correct conversion.
+   * Resolve the path to the AssemblyScript artifact (`dist/wasm/mathts-as.wasm`).
    */
   async resolveAsWasmPath() {
     const isNode = typeof process !== "undefined" && process.versions?.node !== void 0;
+    const { resolvePackagedWasm, defaultWasmLocation } = await import("./resolve-VJX3YTS6.js");
     if (isNode) {
-      const { resolvePackagedWasm } = await import("./resolve-JRHDDNVQ.js");
       const found = await resolvePackagedWasm(import.meta.url, "mathts-as.wasm");
       if (found) return found;
-      const { fileURLToPath } = await import("url");
-      return fileURLToPath(new URL(`../../../lib/wasm/mathts-as.wasm`, import.meta.url));
+      return defaultWasmLocation(import.meta.url, "mathts-as.wasm");
     }
-    return new URL(`../../../lib/wasm/mathts-as.wasm`, import.meta.url).href;
+    return defaultWasmLocation(import.meta.url, "mathts-as.wasm", { browser: true });
   }
   /**
    * Compile + instantiate the AssemblyScript WASM artifact. Each
@@ -5483,27 +5478,22 @@ var WasmLoader = class _WasmLoader {
    * (`mathts-as.wasm`); the legacy second toolchain was removed in the
    * WASM-backend migration (complete 2026-06-26).
    *
-   * The path is resolved relative to this source file's location so it is
-   * CWD-independent.  This file lives at:
-   *   <repo-root>/matrix/src/backends/WasmLoader.ts
-   * so the repo root is three directories up, and the artifact is at
-   *   <repo-root>/lib/wasm/mathts-as.wasm
-   *
-   * Both branches use `new URL(relative, import.meta.url)` to resolve a
-   * file: URL. In Node we convert via `fileURLToPath` (which correctly
-   * strips the leading slash on Windows so the drive letter is not
-   * doubled); in the browser we keep the full `.href` for fetch().
+   * Resolution is package-relative and CWD-independent: the packaged artifact
+   * (`dist/wasm/`) is preferred via `resolvePackagedWasm`; when absent, the
+   * canonical expected location is returned so the missing-binary warning is
+   * actionable (see `defaultWasmLocation` in wasm/resolve.ts — B-3).
    */
   async getDefaultWasmPath() {
     const wasmFile = "mathts-as.wasm";
     if (this.isNode) {
-      const { resolvePackagedWasm } = await import("./resolve-JRHDDNVQ.js");
+      const { resolvePackagedWasm } = await import("./resolve-VJX3YTS6.js");
       const found = await resolvePackagedWasm(import.meta.url, wasmFile);
       if (found) return found;
-      const { fileURLToPath } = await import("url");
-      return fileURLToPath(new URL(`../../../lib/wasm/${wasmFile}`, import.meta.url));
+      const { defaultWasmLocation: defaultWasmLocation2 } = await import("./resolve-VJX3YTS6.js");
+      return defaultWasmLocation2(import.meta.url, wasmFile);
     }
-    return new URL(`../../../lib/wasm/${wasmFile}`, import.meta.url).href;
+    const { defaultWasmLocation } = await import("./resolve-VJX3YTS6.js");
+    return defaultWasmLocation(import.meta.url, wasmFile, { browser: true });
   }
   async loadNodeWasm(path, totalStart) {
     const fs = await import("fs");
@@ -6622,11 +6612,32 @@ function jacobiSVD(A) {
   }
   return { d, U, V };
 }
+function completeOrthonormalBasis(cols, m) {
+  const k = cols[0]?.length ?? 0;
+  const basis = Array.from({ length: k }, (_, j) => cols.map((row2) => row2[j]));
+  for (let cand = 0; cand < m && basis.length < m; cand++) {
+    const v = Array.from({ length: m }, (_, i) => i === cand ? 1 : 0);
+    for (let pass = 0; pass < 2; pass++) {
+      for (const b of basis) {
+        let dot = 0;
+        for (let i = 0; i < m; i++) dot += v[i] * b[i];
+        for (let i = 0; i < m; i++) v[i] -= dot * b[i];
+      }
+    }
+    let norm4 = 0;
+    for (let i = 0; i < m; i++) norm4 += v[i] * v[i];
+    norm4 = Math.sqrt(norm4);
+    if (norm4 > 0.1) {
+      basis.push(v.map((x) => x / norm4));
+    }
+  }
+  return Array.from({ length: m }, (_, i) => basis.map((b) => b[i]));
+}
 function svd(matrix2, options = {}) {
   const {
     maxIterations = DEFAULT_MAX_ITERATIONS,
     tolerance = DEFAULT_TOLERANCE,
-    fullMatrices: _fullMatrices = true,
+    fullMatrices = false,
     rankTolerance = 1e-10
   } = options;
   let A;
@@ -6727,12 +6738,12 @@ function svd(matrix2, options = {}) {
     return {
       U: sortedV,
       S: sortedS,
-      V: sortedU,
+      V: fullMatrices && sortedU.length > (sortedU[0]?.length ?? 0) ? completeOrthonormalBasis(sortedU, sortedU.length) : sortedU,
       rank
     };
   }
   return {
-    U: sortedU,
+    U: fullMatrices && sortedU.length > (sortedU[0]?.length ?? 0) ? completeOrthonormalBasis(sortedU, sortedU.length) : sortedU,
     S: sortedS,
     V: sortedV,
     rank

package/dist/{resolve-JRHDDNVQ.js → resolve-VJX3YTS6.js}

@@ -19,6 +19,25 @@ async function resolvePackagedWasm(metaUrl, wasmFile) {
   }
   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
 };

package/dist/wasm/wasm-manifest.json

@@ -1,3 +1,3 @@
 {
-  "mathts-as.wasm": "sha384-eSck1dw9Rtu8oUAQ/0GMpC1L4bhAPNVS72581QGAyd1xDlRvonQT/XA2RZjpxghB"
+  "mathts-as.wasm": "sha384-q4RGhtVhsLmeCfG5lYpL70f51laygEpfcKGuxKWFnoC9iSXZPXMUEYE9V8W8gYZu"
 }

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@danielsimonjr/mathts-matrix",
-  "version": "0.1.14",
+  "version": "0.2.1",
   "description": "Matrix operations for MathTS with WASM/WebGPU backend support",
   "author": "Daniel Simon Jr.",
   "license": "MIT",
@@ -32,8 +32,8 @@
     "build:prod": "tsup src/index.ts --format esm --dts --clean --minify --treeshake && node scripts/copy-wasm.mjs"
   },
   "dependencies": {
-    "@danielsimonjr/mathts-core": "^0.4.0",
-    "@danielsimonjr/mathts-parallel": "^0.3.0"
+    "@danielsimonjr/mathts-core": "^0.5.0",
+    "@danielsimonjr/mathts-parallel": "^0.3.3"
   },
   "devDependencies": {
     "@types/node": "^25.5.2",