numbl 0.0.11 → 0.0.12

package/dist-cli/cli.js

@@ -20504,6 +20504,82 @@ var erfcxScalar = (x) => {
 };
 
 // src/numbl-core/builtins/math.ts
+function cRecip(re, im) {
+  const d = re * re + im * im;
+  return { re: re / d, im: -im / d };
+}
+function composeRecip(complexFn) {
+  return (re, im) => {
+    const r = cRecip(re, im);
+    return complexFn(r.re, r.im);
+  };
+}
+function cAcos(re, im) {
+  const w2re = 1 - re * re + im * im, w2im = -2 * re * im;
+  const w2r = Math.sqrt(w2re * w2re + w2im * w2im);
+  const w3re = Math.sqrt((w2r + w2re) / 2);
+  const w3im = w2im >= 0 ? Math.sqrt((w2r - w2re) / 2) : -Math.sqrt((w2r - w2re) / 2);
+  const w4re = re - w3im, w4im = im + w3re;
+  const w5re = Math.log(Math.sqrt(w4re * w4re + w4im * w4im));
+  const w5im = Math.atan2(w4im, w4re);
+  return { re: w5im, im: -w5re };
+}
+function cAsin(re, im) {
+  const w1re = -im, w1im = re;
+  const w2re = 1 - re * re + im * im, w2im = -2 * re * im;
+  const w2r = Math.sqrt(w2re * w2re + w2im * w2im);
+  const w3re = Math.sqrt((w2r + w2re) / 2);
+  const w3im = w2im >= 0 ? Math.sqrt((w2r - w2re) / 2) : -Math.sqrt((w2r - w2re) / 2);
+  const w4re = w1re + w3re, w4im = w1im + w3im;
+  const w5re = Math.log(Math.sqrt(w4re * w4re + w4im * w4im));
+  const w5im = Math.atan2(w4im, w4re);
+  return { re: w5im, im: -w5re };
+}
+function cAcosh(re, im) {
+  const a1re = re - 1, a1im = im;
+  const a1r = Math.sqrt(a1re * a1re + a1im * a1im);
+  const s1re = Math.sqrt((a1r + a1re) / 2);
+  const s1im = a1im >= 0 ? Math.sqrt((a1r - a1re) / 2) : -Math.sqrt((a1r - a1re) / 2);
+  const a2re = re + 1, a2im = im;
+  const a2r = Math.sqrt(a2re * a2re + a2im * a2im);
+  const s2re = Math.sqrt((a2r + a2re) / 2);
+  const s2im = a2im >= 0 ? Math.sqrt((a2r - a2re) / 2) : -Math.sqrt((a2r - a2re) / 2);
+  const w2re = s1re * s2re - s1im * s2im;
+  const w2im = s1re * s2im + s1im * s2re;
+  const w3re = re + w2re, w3im = im + w2im;
+  return {
+    re: Math.log(Math.sqrt(w3re * w3re + w3im * w3im)),
+    im: Math.atan2(w3im, w3re)
+  };
+}
+function cAsinh(re, im) {
+  const z2re = re * re - im * im, z2im = 2 * re * im;
+  const w1re = z2re + 1, w1im = z2im;
+  const w1r = Math.sqrt(w1re * w1re + w1im * w1im);
+  const w2re = Math.sqrt((w1r + w1re) / 2);
+  const w2im = w1im >= 0 ? Math.sqrt((w1r - w1re) / 2) : -Math.sqrt((w1r - w1re) / 2);
+  const w3re = re + w2re, w3im = im + w2im;
+  return {
+    re: Math.log(Math.sqrt(w3re * w3re + w3im * w3im)),
+    im: Math.atan2(w3im, w3re)
+  };
+}
+function cAtanh(re, im) {
+  if (im === 0 && (re > 1 || re < -1)) {
+    return {
+      re: 0.5 * Math.log(Math.abs((1 + re) / (1 - re))),
+      im: re > 0 ? Math.PI / 2 : -Math.PI / 2
+    };
+  }
+  const w1re = 1 + re, w1im = im;
+  const w2re = 1 - re, w2im = -im;
+  const denom = w2re * w2re + w2im * w2im;
+  const w3re = (w1re * w2re + w1im * w2im) / denom;
+  const w3im = (w1im * w2re - w1re * w2im) / denom;
+  const w4re = Math.log(Math.sqrt(w3re * w3re + w3im * w3im));
+  const w4im = Math.atan2(w3im, w3re);
+  return { re: w4re / 2, im: w4im / 2 };
+}
 function applyToComplex(v, fn) {
   if (isRuntimeComplexNumber(v)) {
     const result = fn(v.re, v.im);
@@ -20783,41 +20859,16 @@ function registerMathFunctions() {
   );
   register(
     "asin",
-    complexElemwise(
-      Math.asin,
-      (re, im) => {
-        const w1re = -im, w1im = re;
-        const w2re = 1 - re * re + im * im, w2im = -2 * re * im;
-        const w2r = Math.sqrt(w2re * w2re + w2im * w2im);
-        const w3re = Math.sqrt((w2r + w2re) / 2);
-        const w3im = w2im >= 0 ? Math.sqrt((w2r - w2re) / 2) : -Math.sqrt((w2r - w2re) / 2);
-        const w4re = w1re + w3re, w4im = w1im + w3im;
-        const w5re = Math.log(Math.sqrt(w4re * w4re + w4im * w4im));
-        const w5im = Math.atan2(w4im, w4re);
-        return { re: w5im, im: -w5re };
-      },
-      "ComplexNumber",
-      { nativeJs: "Math.asin" }
-    ),
+    complexElemwise(Math.asin, cAsin, "ComplexNumber", {
+      nativeJs: "Math.asin"
+    }),
     1
   );
   register(
     "acos",
-    complexElemwise(
-      Math.acos,
-      (re, im) => {
-        const w2re = 1 - re * re + im * im, w2im = -2 * re * im;
-        const w2r = Math.sqrt(w2re * w2re + w2im * w2im);
-        const w3re = Math.sqrt((w2r + w2re) / 2);
-        const w3im = w2im >= 0 ? Math.sqrt((w2r - w2re) / 2) : -Math.sqrt((w2r - w2re) / 2);
-        const w4re = re - w3im, w4im = im + w3re;
-        const w5re = Math.log(Math.sqrt(w4re * w4re + w4im * w4im));
-        const w5im = Math.atan2(w4im, w4re);
-        return { re: w5im, im: -w5re };
-      },
-      "ComplexNumber",
-      { nativeJs: "Math.acos" }
-    ),
+    complexElemwise(Math.acos, cAcos, "ComplexNumber", {
+      nativeJs: "Math.acos"
+    }),
     1
   );
   register(
@@ -20883,63 +20934,23 @@ function registerMathFunctions() {
   );
   register(
     "asinh",
-    complexElemwise(
-      Math.asinh,
-      (re, im) => {
-        const z2re = re * re - im * im, z2im = 2 * re * im;
-        const w1re = z2re + 1, w1im = z2im;
-        const w1r = Math.sqrt(w1re * w1re + w1im * w1im);
-        const w2re = Math.sqrt((w1r + w1re) / 2);
-        const w2im = w1im >= 0 ? Math.sqrt((w1r - w1re) / 2) : -Math.sqrt((w1r - w1re) / 2);
-        const w3re = re + w2re, w3im = im + w2im;
-        return {
-          re: Math.log(Math.sqrt(w3re * w3re + w3im * w3im)),
-          im: Math.atan2(w3im, w3re)
-        };
-      },
-      "ComplexNumber",
-      { nativeJs: "Math.asinh" }
-    ),
+    complexElemwise(Math.asinh, cAsinh, "ComplexNumber", {
+      nativeJs: "Math.asinh"
+    }),
     1
   );
   register(
     "acosh",
-    complexElemwise(
-      Math.acosh,
-      (re, im) => {
-        const z2re = re * re - im * im, z2im = 2 * re * im;
-        const w1re = z2re - 1, w1im = z2im;
-        const w1r = Math.sqrt(w1re * w1re + w1im * w1im);
-        const w2re = Math.sqrt((w1r + w1re) / 2);
-        const w2im = w1im >= 0 ? Math.sqrt((w1r - w1re) / 2) : -Math.sqrt((w1r - w1re) / 2);
-        const w3re = re + w2re, w3im = im + w2im;
-        return {
-          re: Math.log(Math.sqrt(w3re * w3re + w3im * w3im)),
-          im: Math.atan2(w3im, w3re)
-        };
-      },
-      "ComplexNumber",
-      { nativeJs: "Math.acosh" }
-    ),
+    complexElemwise(Math.acosh, cAcosh, "ComplexNumber", {
+      nativeJs: "Math.acosh"
+    }),
     1
   );
   register(
     "atanh",
-    complexElemwise(
-      Math.atanh,
-      (re, im) => {
-        const w1re = 1 + re, w1im = im;
-        const w2re = 1 - re, w2im = -im;
-        const denom = w2re * w2re + w2im * w2im;
-        const w3re = (w1re * w2re + w1im * w2im) / denom;
-        const w3im = (w1im * w2re - w1re * w2im) / denom;
-        const w4re = Math.log(Math.sqrt(w3re * w3re + w3im * w3im));
-        const w4im = Math.atan2(w3im, w3re);
-        return { re: w4re / 2, im: w4im / 2 };
-      },
-      "ComplexNumber",
-      { nativeJs: "Math.atanh" }
-    ),
+    complexElemwise(Math.atanh, cAtanh, "ComplexNumber", {
+      nativeJs: "Math.atanh"
+    }),
     1
   );
   const deg2rad = Math.PI / 180;
@@ -21474,17 +21485,29 @@ function registerMathFunctions() {
   );
   register(
     "asech",
-    elemwise((x) => Math.acosh(1 / x)),
+    complexElemwise(
+      (x) => Math.acosh(1 / x),
+      composeRecip(cAcosh),
+      "ComplexNumber"
+    ),
     1
   );
   register(
     "acsch",
-    elemwise((x) => Math.asinh(1 / x)),
+    complexElemwise(
+      (x) => Math.asinh(1 / x),
+      composeRecip(cAsinh),
+      "ComplexNumber"
+    ),
     1
   );
   register(
     "acoth",
-    elemwise((x) => Math.atanh(1 / x)),
+    complexElemwise(
+      (x) => Math.atanh(1 / x),
+      composeRecip(cAtanh),
+      "ComplexNumber"
+    ),
     1
   );
   register(
@@ -21504,12 +21527,20 @@ function registerMathFunctions() {
   );
   register(
     "asec",
-    elemwise((x) => Math.acos(1 / x)),
+    complexElemwise(
+      (x) => Math.acos(1 / x),
+      composeRecip(cAcos),
+      "ComplexNumber"
+    ),
     1
   );
   register(
     "acsc",
-    elemwise((x) => Math.asin(1 / x)),
+    complexElemwise(
+      (x) => Math.asin(1 / x),
+      composeRecip(cAsin),
+      "ComplexNumber"
+    ),
     1
   );
   register(
@@ -21517,14 +21548,26 @@ function registerMathFunctions() {
     elemwise((x) => Math.atan(1 / x)),
     1
   );
+  const toDeg = (fn) => (re, im) => {
+    const r = fn(re, im);
+    return { re: r.re * rad2deg, im: r.im * rad2deg };
+  };
   register(
     "asecd",
-    elemwise((x) => Math.acos(1 / x) * rad2deg),
+    complexElemwise(
+      (x) => Math.acos(1 / x) * rad2deg,
+      toDeg(composeRecip(cAcos)),
+      "ComplexNumber"
+    ),
     1
   );
   register(
     "acscd",
-    elemwise((x) => Math.asin(1 / x) * rad2deg),
+    complexElemwise(
+      (x) => Math.asin(1 / x) * rad2deg,
+      toDeg(composeRecip(cAsin)),
+      "ComplexNumber"
+    ),
     1
   );
   register(
@@ -28669,6 +28712,9 @@ var Runtime = class _Runtime {
   RTV = RTV;
   // Persistent variable storage: funcJsId → varName → value
   persistentStore = /* @__PURE__ */ new Map();
+  // Global variable storage: varName → value (shared across all functions)
+  // Public object so generated code can read/write as $rt.$g["name"]
+  $g = {};
   // Accessor guard: prevents recursive getter/setter/subsref calls
   activeAccessors = /* @__PURE__ */ new Set();
   // Callbacks into the compiler for class information (set by executeCode).
@@ -38967,8 +39013,15 @@ function genStmt(cg, stmt) {
       }
       break;
     }
-    case "Global":
+    case "Global": {
+      for (const v of stmt.vars) {
+        cg.globalVarRefs.set(
+          v.variable.id.id,
+          `$rt.$g[${JSON.stringify(v.name)}]`
+        );
+      }
       break;
+    }
     case "Persistent": {
       if (cg.currentFunctionJsId) {
         const funcId = JSON.stringify(cg.currentFunctionJsId);
@@ -39472,6 +39525,9 @@ var Codegen = class {
   generationInProgress = /* @__PURE__ */ new Map();
   // VarId → JS name cache
   varNameCache = /* @__PURE__ */ new Map();
+  // Global variable overrides: VarId → $rt.$g["name"] expression.
+  // Set by Global statement codegen so all reads/writes go through the shared store.
+  globalVarRefs = /* @__PURE__ */ new Map();
   // Stack of variable reference overrides (for closure capture in anonymous functions).
   // Each entry maps varId → replacement JS expression.
   varRefOverrides = [];
@@ -39552,6 +39608,8 @@ var Codegen = class {
       const override = this.varRefOverrides[i].get(id);
       if (override !== void 0) return override;
     }
+    const globalRef = this.globalVarRefs.get(id);
+    if (globalRef !== void 0) return globalRef;
     const cached = this.varNameCache.get(id);
     if (cached) return cached;
     const lastUnderscore = id.lastIndexOf("_");
@@ -40242,8 +40300,27 @@ var stdlibFiles = [
     source: "function parse(obj, varargin)\n    obj.parse(varargin{:});\nend\n"
   }
 ];
+var shimFiles = [
+  {
+    name: "+matlab/+internal/+decomposition/DenseLU.m",
+    source: "classdef DenseLU\n    properties\n        L\n        U\n        P\n    end\n    methods\n        function obj = DenseLU(A)\n            [obj.L, obj.U, obj.P] = lu(A);\n        end\n        function x = solve(obj, b, ~)\n            x = obj.U \\ (obj.L \\ (obj.P * b));\n        end\n    end\nend"
+  },
+  {
+    name: "+matlab/+internal/+math/blkdiag.m",
+    source: "function B = blkdiag(varargin)\n    B = builtin('blkdiag', varargin{:});\nend\n"
+  },
+  {
+    name: "decomposition.m",
+    source: "classdef decomposition\n    properties\n        A\n        checkCondition\n        isIllCond\n    end\n    methods\n        function obj = decomposition(A, varargin)\n            obj.A = A;\n            obj.checkCondition = true;\n            for i = 1:2:length(varargin)\n                if strcmpi(varargin{i}, 'CheckCondition')\n                    obj.checkCondition = varargin{i+1};\n                end\n            end\n            if obj.checkCondition\n                obj.isIllCond = (rcond(A) < eps);\n            else\n                obj.isIllCond = false;\n            end\n        end\n        function x = mldivide(obj, b)\n            x = obj.A \\ b;\n        end\n    end\nend\n"
+  },
+  {
+    name: "isIllConditioned.m",
+    source: "function result = isIllConditioned(dA)\n    result = dA.isIllCond;\nend\n"
+  }
+];
 
 // src/numbl-core/codegen/generateMainScriptCode.ts
+var SHIM_SEARCH_PATH = "__numbl_shims__";
 function generateMainScriptCode(source, mainFileName = "script.m", workspaceFiles, initialVariableNames, searchPaths) {
   let _t0 = performance.now();
   const ast = parseMFile(source, mainFileName);
@@ -40269,7 +40346,6 @@ function generateMainScriptCode(source, mainFileName = "script.m", workspaceFile
   }
   const _parseMainMs = performance.now() - _t0;
   const ctx = new LoweringContext(source, mainFileName);
-  ctx.registry.searchPaths = searchPaths ?? [];
   const mWorkspaceFiles = [];
   const jsWorkspaceFiles = [];
   if (workspaceFiles) {
@@ -40284,6 +40360,13 @@ function generateMainScriptCode(source, mainFileName = "script.m", workspaceFile
   for (const f of stdlibFiles) {
     mWorkspaceFiles.push(f);
   }
+  for (const f of shimFiles) {
+    mWorkspaceFiles.push({
+      name: `${SHIM_SEARCH_PATH}/${f.name}`,
+      source: f.source
+    });
+  }
+  ctx.registry.searchPaths = [...searchPaths ?? [], SHIM_SEARCH_PATH];
   _t0 = performance.now();
   const jsUserFunctions = loadJsUserFunctions(jsWorkspaceFiles);
   const _loadJsUserFunctionsMs = performance.now() - _t0;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "numbl",
-  "version": "0.0.11",
+  "version": "0.0.12",
   "description": "Run .m source files in the browser and on the command line by compiling to JavaScript",
   "license": "Apache-2.0",
   "type": "module",
@@ -31,7 +31,8 @@
     "dev": "vite",
     "build": "npm run build:web && npm run build:cli && npm run build:plot-viewer",
     "build:plot-viewer": "vite build --config vite.plot-viewer.config.ts",
-    "build:web": "tsc -b && vite build",
+    "bundle:stdlib": "tsx scripts/bundle-stdlib.ts",
+    "build:web": "npm run bundle:stdlib && tsc -b && vite build",
     "build:cli": "esbuild src/cli.ts --bundle --platform=node --format=esm --outfile=dist-cli/cli.js --packages=external --loader:.wasm=file",
     "build:server": "esbuild server/execution-service.ts --bundle --platform=node --format=esm --outfile=server/dist/execution-service.js --packages=external",
     "server:dev": "node --env-file=server/.env node_modules/.bin/tsx server/execution-service.ts",