numbl 0.4.11 → 0.4.13

package/dist-cli/cli.js

@@ -312,7 +312,7 @@ import {
   delimiter,
   dirname as dirname3,
   join as join7,
-  relative,
+  relative as relative2,
   resolve as resolve2,
   sep
 } from "path";
@@ -21701,6 +21701,114 @@ var valuesAreEqual = (a, b) => {
   }
 };
 
+// src/numbl-core/runtime/datetime-format.ts
+var MONTH_ABBR = [
+  "Jan",
+  "Feb",
+  "Mar",
+  "Apr",
+  "May",
+  "Jun",
+  "Jul",
+  "Aug",
+  "Sep",
+  "Oct",
+  "Nov",
+  "Dec"
+];
+var MONTH_FULL = [
+  "January",
+  "February",
+  "March",
+  "April",
+  "May",
+  "June",
+  "July",
+  "August",
+  "September",
+  "October",
+  "November",
+  "December"
+];
+function pad(n, width) {
+  return String(n).padStart(width, "0");
+}
+function formatDatetimePattern(pattern, year, month, day, hour, minute, second) {
+  const wholeSec = Math.floor(second);
+  const frac = second - wholeSec;
+  const monthIdx = Math.max(0, Math.min(11, Math.floor(month) - 1));
+  const formatToken = (c, n) => {
+    switch (c) {
+      case "y":
+      case "u":
+        if (n === 2) return pad(Math.floor(year) % 100, 2);
+        return pad(Math.floor(year), n);
+      case "M":
+        if (n === 1) return String(Math.floor(month));
+        if (n === 2) return pad(Math.floor(month), 2);
+        if (n === 3) return MONTH_ABBR[monthIdx];
+        return MONTH_FULL[monthIdx];
+      case "d":
+        return pad(Math.floor(day), n);
+      case "H":
+        return pad(Math.floor(hour), n);
+      case "h":
+        return pad(Math.floor(hour) % 12 || 12, n);
+      case "m":
+        return pad(Math.floor(minute), n);
+      case "s":
+        return pad(wholeSec, n);
+      case "S":
+        return pad(Math.floor(frac * Math.pow(10, n)), n);
+      case "a":
+        return hour < 12 ? "AM" : "PM";
+      default:
+        throw new RuntimeError(
+          `datetime: unsupported format token '${c.repeat(n)}' in '${pattern}'`
+        );
+    }
+  };
+  let out = "";
+  let i = 0;
+  while (i < pattern.length) {
+    const c = pattern[i];
+    if (c === "'") {
+      i++;
+      let lit = "";
+      let closed = false;
+      while (i < pattern.length) {
+        if (pattern[i] === "'") {
+          if (pattern[i + 1] === "'") {
+            lit += "'";
+            i += 2;
+            continue;
+          }
+          i++;
+          closed = true;
+          break;
+        }
+        lit += pattern[i++];
+      }
+      if (!closed && lit !== "")
+        throw new RuntimeError(
+          `datetime: unterminated quoted literal in format '${pattern}'`
+        );
+      out += lit === "" ? "'" : lit;
+      continue;
+    }
+    if (!/[a-zA-Z]/.test(c)) {
+      out += c;
+      i++;
+      continue;
+    }
+    let j = i;
+    while (j < pattern.length && pattern[j] === c) j++;
+    out += formatToken(c, j - i);
+    i = j;
+  }
+  return out;
+}
+
 // src/numbl-core/runtime/display.ts
 function displayValue(v) {
   if (isRuntimeNumber(v)) {
@@ -21975,6 +22083,13 @@ function formatDatetimeOrDuration(v) {
   if (v.className === "duration") return formatDurationInstance(v);
   return null;
 }
+function textField(v, name) {
+  const fv = v.fields.get(name);
+  if (fv === void 0) return null;
+  if (isRuntimeString(fv)) return fv;
+  if (isRuntimeChar(fv)) return fv.value;
+  return null;
+}
 function formatDatetimeInstance(v) {
   const year = numField(v, "Year");
   const month = numField(v, "Month");
@@ -21984,6 +22099,10 @@ function formatDatetimeInstance(v) {
   const second = numField(v, "Second");
   if (year === null || month === null || day === null || hour === null || minute === null || second === null)
     return null;
+  const fmt = textField(v, "Format");
+  if (fmt !== null && fmt !== "") {
+    return formatDatetimePattern(fmt, year, month, day, hour, minute, second);
+  }
   const mIdx = Math.max(0, Math.min(11, Math.floor(month) - 1));
   const mon = DATETIME_MONTH_ABBR[mIdx];
   const dd = String(Math.floor(day)).padStart(2, "0");
@@ -23069,7 +23188,7 @@ function indexIntoTensorWithTensor(base, idx) {
     if (hasImag) imIndices.push(base.imag[k]);
   }
   const idxIs0x0 = idx.data.length === 0 && idx.shape.length === 2 && idx.shape[0] === 0 && idx.shape[1] === 0;
-  const baseIsVector = base.shape.length <= 2 && (base.shape[0] === 1 || base.shape[1] === 1 || base.shape.length === 1);
+  const baseIsVector = base.data.length !== 1 && base.shape.length <= 2 && (base.shape[0] === 1 || base.shape[1] === 1 || base.shape.length === 1);
   const idxIsVector = idx.shape.length <= 2 && (idx.shape[0] === 1 || idx.shape[1] === 1 || idx.shape.length === 1);
   const outShape = idxIs0x0 ? [0, 0] : baseIsVector && idxIsVector ? base.shape[0] === 1 ? [1, resultData.length] : [resultData.length, 1] : idx.shape;
   const imOut = hasImag && imIndices.some((x) => x !== 0) ? allocFloat64Array(imIndices) : void 0;
@@ -23131,6 +23250,17 @@ function indexIntoRTValue(base, indices) {
         return base.elements[i];
       }
       if (isRuntimeTensor(idx)) {
+        if (idx._isLogical) {
+          const newElements2 = [];
+          for (let k = 0; k < idx.data.length; k++) {
+            if (idx.data[k] !== 0) {
+              if (k >= base.elements.length)
+                throw new RuntimeError("Index exceeds struct array bounds");
+              newElements2.push(base.elements[k]);
+            }
+          }
+          return RTV.structArray(base.fieldNames, newElements2);
+        }
         const newElements = [];
         for (let k = 0; k < idx.data.length; k++) {
           const i = Math.round(idx.data[k]) - 1;
@@ -27482,8 +27612,8 @@ function sprintfFormat(fmt, args) {
               if (width > 0) {
                 const padChar = zeroPad ? "0" : " ";
                 const padLen = Math.max(0, width - sign.length - s.length);
-                const pad = padChar.repeat(padLen);
-                result += leftAlign ? sign + s + " ".repeat(padLen) : zeroPad ? sign + pad + s : pad + sign + s;
+                const pad2 = padChar.repeat(padLen);
+                result += leftAlign ? sign + s + " ".repeat(padLen) : zeroPad ? sign + pad2 + s : pad2 + sign + s;
               } else {
                 result += sign + s;
               }
@@ -37818,6 +37948,32 @@ function datetimeFromUtcMs(ms) {
 function datetimeFromDatenum(dn) {
   return datetimeFromUtcMs((dn - MATLAB_EPOCH_DAYS) * 864e5);
 }
+function withTimeZoneFormat(dt, timeZone, format) {
+  if (timeZone === void 0 && format === void 0) return dt;
+  const fields = new Map(dt.fields);
+  if (timeZone !== void 0) fields.set("TimeZone", RTV.char(timeZone));
+  if (format !== void 0) fields.set("Format", RTV.char(format));
+  return new RuntimeClassInstance("datetime", fields, false);
+}
+function todayParts(offsetDays, utc) {
+  const now = /* @__PURE__ */ new Date();
+  if (utc) {
+    const d2 = new Date(
+      Date.UTC(
+        now.getUTCFullYear(),
+        now.getUTCMonth(),
+        now.getUTCDate() + offsetDays
+      )
+    );
+    return [d2.getUTCFullYear(), d2.getUTCMonth() + 1, d2.getUTCDate()];
+  }
+  const d = new Date(
+    now.getFullYear(),
+    now.getMonth(),
+    now.getDate() + offsetDays
+  );
+  return [d.getFullYear(), d.getMonth() + 1, d.getDate()];
+}
 function makeDuration(totalSeconds) {
   const fields = /* @__PURE__ */ new Map([
     ["Seconds", RTV.num(totalSeconds)]
@@ -37867,9 +38023,10 @@ registerIBuiltin({
       "t = datetime('now' | 'today' | 'yesterday' | 'tomorrow')",
       "t = datetime(Y, M, D)",
       "t = datetime(Y, M, D, H, MI, S)",
-      "t = datetime(X, 'ConvertFrom', dateType)"
+      "t = datetime(X, 'ConvertFrom', dateType)",
+      "t = datetime(..., 'TimeZone', tz, 'Format', fmt)"
     ],
-    description: "Create a datetime value. Supported dateType values for ConvertFrom: 'datenum', 'posixtime', 'excel', 'excel1904'."
+    description: "Create a datetime value. Supported dateType values for ConvertFrom: 'datenum', 'posixtime', 'excel', 'excel1904'. Trailing 'TimeZone' (only '', 'local', 'UTC') and 'Format' name-value pairs are supported; the Format pattern is used by char() and display."
   },
   resolve: () => ({
     outputTypes: [
@@ -37881,75 +38038,83 @@ registerIBuiltin({
       }
     ],
     apply: (args) => {
-      if (args.length === 0) {
-        return datetimeFromDate(/* @__PURE__ */ new Date());
+      let timeZone;
+      let format;
+      let positional = args;
+      while (positional.length >= 2) {
+        const name = argToString(positional[positional.length - 2]);
+        if (name === null) break;
+        const lc = name.toLowerCase();
+        if (lc !== "timezone" && lc !== "format") break;
+        const value = argToString(positional[positional.length - 1]);
+        if (value === null)
+          throw new RuntimeError(`datetime: ${name} value must be text`);
+        if (lc === "timezone") timeZone = value;
+        else format = value;
+        positional = positional.slice(0, -2);
+      }
+      const tzLc = timeZone?.toLowerCase();
+      if (tzLc !== void 0 && tzLc !== "" && tzLc !== "local" && tzLc !== "utc") {
+        throw new RuntimeError(
+          `datetime: unsupported TimeZone '${timeZone}' (supported: '', 'local', 'UTC')`
+        );
       }
-      if (args.length === 1) {
-        const s = argToString(args[0]);
+      const utc = tzLc === "utc";
+      const finish = (dt) => withTimeZoneFormat(dt, timeZone, format);
+      if (positional.length === 0) {
+        return finish(
+          utc ? datetimeFromUtcMs(Date.now()) : datetimeFromDate(/* @__PURE__ */ new Date())
+        );
+      }
+      if (positional.length === 1) {
+        const s = argToString(positional[0]);
         if (s !== null) {
-          const now = /* @__PURE__ */ new Date();
           const lc = s.toLowerCase();
-          if (lc === "now") return datetimeFromDate(now);
-          if (lc === "today") {
-            return makeDatetime(
-              now.getFullYear(),
-              now.getMonth() + 1,
-              now.getDate(),
-              0,
-              0,
-              0
+          if (lc === "now") {
+            return finish(
+              utc ? datetimeFromUtcMs(Date.now()) : datetimeFromDate(/* @__PURE__ */ new Date())
             );
           }
-          if (lc === "yesterday" || lc === "tomorrow") {
-            const d = new Date(
-              now.getFullYear(),
-              now.getMonth(),
-              now.getDate() + (lc === "tomorrow" ? 1 : -1)
-            );
-            return makeDatetime(
-              d.getFullYear(),
-              d.getMonth() + 1,
-              d.getDate(),
-              0,
-              0,
-              0
-            );
+          if (lc === "today" || lc === "yesterday" || lc === "tomorrow") {
+            const offset = lc === "tomorrow" ? 1 : lc === "yesterday" ? -1 : 0;
+            const [Y, M, D] = todayParts(offset, utc);
+            return finish(makeDatetime(Y, M, D, 0, 0, 0));
           }
           throw new RuntimeError(
             `datetime: unrecognized relative-day string '${s}'`
           );
         }
       }
-      if (args.length === 3) {
-        const flag = argToString(args[1]);
-        const kind = argToString(args[2]);
+      if (positional.length === 3) {
+        const flag = argToString(positional[1]);
+        const kind = argToString(positional[2]);
         if (flag !== null && flag.toLowerCase() === "convertfrom") {
           if (kind === null)
             throw new RuntimeError("datetime: dateType must be text");
-          const x = toNumber(args[0]);
+          const x = toNumber(positional[0]);
           const k = kind.toLowerCase();
-          if (k === "datenum") return datetimeFromDatenum(x);
-          if (k === "posixtime") return datetimeFromUtcMs(x * 1e3);
+          if (k === "datenum") return finish(datetimeFromDatenum(x));
+          if (k === "posixtime") return finish(datetimeFromUtcMs(x * 1e3));
           if (k === "excel") {
-            return datetimeFromUtcMs((x - 25569) * 864e5);
+            return finish(datetimeFromUtcMs((x - 25569) * 864e5));
           }
           if (k === "excel1904") {
-            return datetimeFromUtcMs((x - 24107) * 864e5);
+            return finish(datetimeFromUtcMs((x - 24107) * 864e5));
           }
           throw new RuntimeError(
             `datetime: unsupported ConvertFrom type '${kind}'`
           );
         }
       }
-      if (args.length === 3 || args.length === 6 || args.length === 7) {
-        const Y = toNumber(args[0]);
-        const M = toNumber(args[1]);
-        const D = toNumber(args[2]);
-        const H2 = args.length >= 6 ? toNumber(args[3]) : 0;
-        const MI = args.length >= 6 ? toNumber(args[4]) : 0;
-        const S = args.length >= 6 ? toNumber(args[5]) : 0;
-        const MS = args.length >= 7 ? toNumber(args[6]) : 0;
-        return makeDatetime(Y, M, D, H2, MI, S + MS / 1e3);
+      if (positional.length === 3 || positional.length === 6 || positional.length === 7) {
+        const Y = toNumber(positional[0]);
+        const M = toNumber(positional[1]);
+        const D = toNumber(positional[2]);
+        const H2 = positional.length >= 6 ? toNumber(positional[3]) : 0;
+        const MI = positional.length >= 6 ? toNumber(positional[4]) : 0;
+        const S = positional.length >= 6 ? toNumber(positional[5]) : 0;
+        const MS = positional.length >= 7 ? toNumber(positional[6]) : 0;
+        return finish(makeDatetime(Y, M, D, H2, MI, S + MS / 1e3));
       }
       throw new RuntimeError(
         `datetime: unsupported argument pattern (nargin=${args.length})`
@@ -39906,6 +40071,7 @@ var SPECIAL_BUILTIN_NAMES = [
   "copyfile",
   "fileattrib",
   "unzip",
+  "zip",
   "dir",
   "warning",
   "input",
@@ -41018,6 +41184,26 @@ function registerSpecialBuiltins(rt) {
     }
     return void 0;
   });
+  registerSpecial("zip", (nargout, args) => {
+    const io = requireFileIO();
+    if (!io.zip)
+      throw new RuntimeError("zip is not available in this environment");
+    const margs = args.map((a) => ensureRuntimeValue(a));
+    if (margs.length < 2)
+      throw new RuntimeError("zip requires at least 2 arguments");
+    let zipfilename = toString(margs[0]);
+    const base = zipfilename.split(/[/\\]/).pop() ?? zipfilename;
+    if (!base.includes(".")) zipfilename += ".zip";
+    const fnames = margs[1];
+    const filenames = isRuntimeCell(fnames) ? fnames.data.map((v) => toString(ensureRuntimeValue(v))) : [toString(fnames)];
+    const rootdir = margs.length >= 3 ? toString(margs[2]) : ".";
+    const entryNames = io.zip(zipfilename, filenames, rootdir);
+    if (nargout >= 1) {
+      const cellData = entryNames.map((f) => RTV.char(f));
+      return RTV.cell(cellData, [1, cellData.length]);
+    }
+    return void 0;
+  });
   registerSpecial("dir", (nargout, args) => {
     const io = requireFileIO();
     if (!io.listDir)
@@ -43278,8 +43464,8 @@ function padExtent(e, tight) {
     return [min2, max2];
   }
   if (tight) return [min2, max2];
-  const pad = (max2 - min2) * 0.05;
-  return [min2 - pad, max2 + pad];
+  const pad2 = (max2 - min2) * 0.05;
+  return [min2 - pad2, max2 + pad2];
 }
 function applyExplicit(auto, explicit) {
   if (!explicit) return auto;
@@ -49484,13 +49670,19 @@ registerIBuiltin({
   resolve: (argTypes) => {
     if (argTypes.length !== 1) return null;
     const a = argTypes[0];
-    if (a.kind !== "char" && a.kind !== "string" && a.kind !== "number" && a.kind !== "tensor" && a.kind !== "cell" && a.kind !== "unknown")
+    if (a.kind !== "char" && a.kind !== "string" && a.kind !== "number" && a.kind !== "tensor" && a.kind !== "cell" && a.kind !== "unknown" && // datetime/duration only — other classes keep their own char methods
+    // (or the generic "does not support these argument types" error).
+    !(a.kind === "class_instance" && (a.className === "datetime" || a.className === "duration")))
       return null;
     return {
       outputTypes: [{ kind: "char" }],
       apply: (args) => {
         const v = args[0];
         if (isRuntimeChar(v)) return v;
+        if (isRuntimeClassInstance(v)) {
+          const s = formatDatetimeOrDuration(v);
+          if (s !== null) return RTV.char(s);
+        }
         if (isRuntimeString(v)) return RTV.char(v);
         if (isRuntimeNumber(v))
           return RTV.char(String.fromCharCode(Math.round(v)));
@@ -51878,7 +52070,7 @@ defineBuiltin({
           const result = invComplexJS(A.data, A.imag, n);
           return RTV.tensor(result.re, [n, n], result.im);
         }
-        const bridge = getEffectiveBridge("inv");
+        const bridge = getEffectiveBridge("inv", "inv");
         return RTV.tensor(bridge.inv(toF64(A.data), n), [n, n]);
       }
     }
@@ -54225,6 +54417,16 @@ defineBuiltin({
           parseReshapeDims(args, 1);
           return v;
         }
+        if (isRuntimeCell(v)) {
+          const total = v.data.length;
+          const s = parseReshapeDims(args, total);
+          while (s.length > 2 && s[s.length - 1] === 1) s.pop();
+          if (s.length > 2)
+            throw new RuntimeError("reshape: cell arrays must be 2-D");
+          const nr = s[0];
+          const nc = s.length >= 2 ? s[1] : 1;
+          return RTV.cell(v.data.slice(), [nr, nc]);
+        }
         if (isRuntimeClassInstanceArray(v)) {
           const total = v.elements.length;
           const s = parseReshapeDims(args, total);
@@ -64075,7 +64277,7 @@ function getSourceLine(getSource, file, line) {
 }
 
 // src/numbl-core/version.ts
-var NUMBL_VERSION = "0.4.11";
+var NUMBL_VERSION = "0.4.13";
 
 // src/cli-repl.ts
 import { createInterface } from "readline";
@@ -66709,6 +66911,12 @@ function evalMember(expr, nargout) {
     }
   }
   const base = this.evalExpr(expr.base);
+  const baseRv = ensureRuntimeValue(base);
+  if (isRuntimeStructArray(baseRv) && baseRv.elements.length !== 1) {
+    return baseRv.elements.map(
+      (el) => ensureRuntimeValue(this.rt.getMember(el, expr.name))
+    );
+  }
   return this.rt.getMember(base, expr.name);
 }
 function tryExtractDottedName(expr) {
@@ -66920,12 +67128,31 @@ function evalCellLiteral(expr) {
     }
     return RTV.cell(elements2, [1, elements2.length]);
   }
-  const numRows = expr.rows.length;
-  const numCols = expr.rows[0].length;
+  const rowValues = expr.rows.map((row) => {
+    const vals = [];
+    for (const e of row) {
+      const v = this.evalExpr(e);
+      if (Array.isArray(v)) {
+        for (const elem of v) vals.push(ensureRuntimeValue(elem));
+      } else {
+        vals.push(ensureRuntimeValue(v));
+      }
+    }
+    return vals;
+  });
+  const numRows = rowValues.length;
+  const numCols = rowValues[0].length;
+  for (const rv of rowValues) {
+    if (rv.length !== numCols) {
+      throw new RuntimeError(
+        "Cell array rows must have the same number of columns"
+      );
+    }
+  }
   const elements = [];
   for (let c = 0; c < numCols; c++) {
     for (let r = 0; r < numRows; r++) {
-      elements.push(ensureRuntimeValue(this.evalExpr(expr.rows[r][c])));
+      elements.push(rowValues[r][c]);
     }
   }
   return RTV.cell(elements, [numRows, numCols]);
@@ -71400,8 +71627,8 @@ function mtoc2_sprintf_format(fmt, args) {
               if (width > 0) {
                 const padChar = zeroPad ? "0" : " ";
                 const padLen = Math.max(0, width - sign.length - s.length);
-                const pad = padChar.repeat(padLen);
-                result += leftAlign ? sign + s + " ".repeat(padLen) : zeroPad ? sign + pad + s : pad + sign + s;
+                const pad2 = padChar.repeat(padLen);
+                result += leftAlign ? sign + s + " ".repeat(padLen) : zeroPad ? sign + pad2 + s : pad2 + sign + s;
               } else {
                 result += sign + s;
               }
@@ -71686,8 +71913,8 @@ function disp_complex_slice(re, im, offset, rows, cols) {
     for (let c = 0; c < cols; c++) {
       const idx = r + c * rows;
       const cell = cells[idx];
-      const pad = colWidths[c] - cell.length;
-      for (let i = 0; i < pad; i++) line += " ";
+      const pad2 = colWidths[c] - cell.length;
+      for (let i = 0; i < pad2; i++) line += " ";
       line += cell;
       if (c < cols - 1) line += "   ";
     }
@@ -73099,6 +73326,14 @@ function mtoc2_tensor_mtimes_real(a, b) {
       `mtoc2_tensor_mtimes_real: inner dimensions disagree (${k} vs ${k2})`
     );
   }
+  if (typeof $matmulAccel === "function") {
+    const accel = $matmulAccel(a.data, m, k, b.data, n);
+    if (accel) {
+      const acc = mtoc2_tensor_alloc(m, n);
+      acc.data.set(accel);
+      return acc;
+    }
+  }
   const r = mtoc2_tensor_alloc(m, n);
   for (let j = 0; j < n; j++) {
     for (let i = 0; i < m; i++) {
@@ -74310,7 +74545,7 @@ var JS_SNIPPETS = {
   "tensor_logspace.js": "// JS sibling of `tensor_logspace.h`. Build a 1\xD7n row tensor of n\n// logarithmically-spaced values from 10^a to 10^b. Byte-for-byte with\n// numbl's interpreter `logspace`, including the MATLAB special case\n// where an upper limit of exactly `pi` makes the last point `pi`.\n\n\nfunction mtoc2_tensor_logspace(a, b, n) {\n  if (n <= 0) return mtoc2_tensor_alloc(1, 0);\n  const isPi = b === Math.PI;\n  const endVal = isPi ? Math.PI : Math.pow(10, b);\n  const out = mtoc2_tensor_alloc(1, n);\n  if (n === 1) {\n    out.data[0] = endVal;\n    return out;\n  }\n  if (isPi) {\n    const logStart = Math.log10(Math.pow(10, a));\n    const logEnd = Math.log10(Math.PI);\n    for (let i = 0; i < n; i++) {\n      const t = logStart + ((logEnd - logStart) * i) / (n - 1);\n      out.data[i] = Math.pow(10, t);\n    }\n  } else {\n    for (let i = 0; i < n; i++) {\n      const t = a + ((b - a) * i) / (n - 1);\n      out.data[i] = Math.pow(10, t);\n    }\n  }\n  return out;\n}\n",
   "tensor_meshgrid.js": "// JS sibling of `tensor_meshgrid.h`. MATLAB-style coordinate grid.\n\n\nfunction numel(t) {\n  let n = 1;\n  for (const s of t.shape) n *= s;\n  return n;\n}\n\nfunction mtoc2_meshgrid_x(x, y) {\n  const M = numel(x);\n  const N = numel(y);\n  const X = mtoc2_tensor_alloc(N, M);\n  for (let j = 0; j < M; j++) {\n    const xj = x.data[j];\n    for (let i = 0; i < N; i++) X.data[i + j * N] = xj;\n  }\n  return X;\n}\n\nfunction mtoc2_meshgrid(x, y) {\n  const M = numel(x);\n  const N = numel(y);\n  const X = mtoc2_tensor_alloc(N, M);\n  const Y = mtoc2_tensor_alloc(N, M);\n  for (let j = 0; j < M; j++) {\n    const xj = x.data[j];\n    for (let i = 0; i < N; i++) {\n      X.data[i + j * N] = xj;\n      Y.data[i + j * N] = y.data[i];\n    }\n  }\n  return { X, Y };\n}\n\nfunction mtoc2_meshgrid_1arg(x) {\n  return mtoc2_meshgrid(x, x);\n}\n",
   "tensor_mtimes_complex.js": '// JS sibling of `tensor_mtimes_complex.h`. 2-D matrix multiplication\n// producing a complex result. Handles all three mixed cases\n// (complex \xD7 complex, complex \xD7 real, real \xD7 complex) by treating an\n// undefined imag lane as implicit zeros.\n\n\nfunction mtoc2_tensor_mtimes_complex_scalar(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_complex_scalar: inputs must be 2-D");\n  }\n  if (a.shape[0] !== 1 || b.shape[1] !== 1) {\n    throw new Error("mtoc2_tensor_mtimes_complex_scalar: requires 1\xD7k * k\xD71");\n  }\n  if (a.shape[1] !== b.shape[0]) {\n    throw new Error("mtoc2_tensor_mtimes_complex_scalar: inner-dim mismatch");\n  }\n  const k = a.shape[1];\n  const aim = a.imag;\n  const bim = b.imag;\n  let accR = 0;\n  let accI = 0;\n  for (let p = 0; p < k; p++) {\n    const ar = a.data[p];\n    const ai = aim !== undefined ? aim[p] : 0;\n    const br = b.data[p];\n    const bi = bim !== undefined ? bim[p] : 0;\n    accR += ar * br - ai * bi;\n    accI += ar * bi + ai * br;\n  }\n  return { re: accR, im: accI };\n}\n\nfunction mtoc2_tensor_mtimes_complex(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_complex: inputs must be 2-D");\n  }\n  const m = a.shape[0];\n  const k = a.shape[1];\n  const k2 = b.shape[0];\n  const n = b.shape[1];\n  if (k !== k2) {\n    throw new Error(\n      `mtoc2_tensor_mtimes_complex: inner-dim mismatch (${m}\xD7${k} * ${k2}\xD7${n})`\n    );\n  }\n  const out = mtoc2_tensor_alloc_complex(m, n);\n  const aim = a.imag;\n  const bim = b.imag;\n  for (let j = 0; j < n; j++) {\n    for (let i = 0; i < m; i++) {\n      let accR = 0;\n      let accI = 0;\n      for (let p = 0; p < k; p++) {\n        const ar = a.data[i + p * m];\n        const ai = aim !== undefined ? aim[i + p * m] : 0;\n        const br = b.data[p + j * k];\n        const bi = bim !== undefined ? bim[p + j * k] : 0;\n        accR += ar * br - ai * bi;\n        accI += ar * bi + ai * br;\n      }\n      out.data[i + j * m] = accR;\n      out.imag[i + j * m] = accI;\n    }\n  }\n  return out;\n}\n',
-  "tensor_mtimes_real.js": '// JS sibling of `tensor_mtimes_real.h`. Real 2-D matrix\n// multiplication A * B. Column-major in, column-major out.\n// `mtoc2_tensor_mtimes_real_scalar` mirrors the C 1\xD7k * k\xD71\n// inner-product path that returns a bare number instead of a\n// freshly-allocated 1\xD71 tensor.\n\n\nfunction mtoc2_tensor_mtimes_real(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_real: inputs must be 2-D");\n  }\n  const m = a.shape[0];\n  const k = a.shape[1];\n  const k2 = b.shape[0];\n  const n = b.shape[1];\n  if (k !== k2) {\n    throw new Error(\n      `mtoc2_tensor_mtimes_real: inner dimensions disagree (${k} vs ${k2})`\n    );\n  }\n  const r = mtoc2_tensor_alloc(m, n);\n  for (let j = 0; j < n; j++) {\n    for (let i = 0; i < m; i++) {\n      let s = 0;\n      for (let p = 0; p < k; p++) {\n        s += a.data[i + p * m] * b.data[p + j * k];\n      }\n      r.data[i + j * m] = s;\n    }\n  }\n  return r;\n}\n\nfunction mtoc2_tensor_mtimes_real_scalar(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_real_scalar: inputs must be 2-D");\n  }\n  // The translate-time type system has already proven a is 1\xD7k and b\n  // is k\xD71; this helper trusts that and just runs the inner product.\n  const k = a.shape[1];\n  let s = 0;\n  for (let p = 0; p < k; p++) {\n    s += a.data[p * a.shape[0]] * b.data[p];\n  }\n  return s;\n}\n',
+  "tensor_mtimes_real.js": '// JS sibling of `tensor_mtimes_real.h`. Real 2-D matrix\n// multiplication A * B. Column-major in, column-major out.\n// `mtoc2_tensor_mtimes_real_scalar` mirrors the C 1\xD7k * k\xD71\n// inner-product path that returns a bare number instead of a\n// freshly-allocated 1\xD71 tensor.\n\n\nfunction mtoc2_tensor_mtimes_real(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_real: inputs must be 2-D");\n  }\n  const m = a.shape[0];\n  const k = a.shape[1];\n  const k2 = b.shape[0];\n  const n = b.shape[1];\n  if (k !== k2) {\n    throw new Error(\n      `mtoc2_tensor_mtimes_real: inner dimensions disagree (${k} vs ${k2})`\n    );\n  }\n  // Optional WASM LAPACK accelerator. When a host installs `$matmulAccel`\n  // (browser worker, after loading a bridge), delegate the multiply; a\n  // falsy return means "declined" (too small / error) and we fall through\n  // to the JS loop below. `typeof` guard is safe even when unset.\n  if (typeof $matmulAccel === "function") {\n    const accel = $matmulAccel(a.data, m, k, b.data, n);\n    if (accel) {\n      const acc = mtoc2_tensor_alloc(m, n);\n      acc.data.set(accel);\n      return acc;\n    }\n  }\n  const r = mtoc2_tensor_alloc(m, n);\n  for (let j = 0; j < n; j++) {\n    for (let i = 0; i < m; i++) {\n      let s = 0;\n      for (let p = 0; p < k; p++) {\n        s += a.data[i + p * m] * b.data[p + j * k];\n      }\n      r.data[i + j * m] = s;\n    }\n  }\n  return r;\n}\n\nfunction mtoc2_tensor_mtimes_real_scalar(a, b) {\n  if (a.shape.length !== 2 || b.shape.length !== 2) {\n    throw new Error("mtoc2_tensor_mtimes_real_scalar: inputs must be 2-D");\n  }\n  // The translate-time type system has already proven a is 1\xD7k and b\n  // is k\xD71; this helper trusts that and just runs the inner product.\n  const k = a.shape[1];\n  let s = 0;\n  for (let p = 0; p < k; p++) {\n    s += a.data[p * a.shape[0]] * b.data[p];\n  }\n  return s;\n}\n',
   "tensor_norm.js": "// JS sibling of `tensor_norm.h`. Vector norms for real and complex\n// tensors. The complex variants tolerate a missing imag lane (real\n// tensor flowing through a complex-typed norm route) by treating it\n// as zero.\n\nfunction mtoc2_norm2_real(a) {\n  let acc = 0;\n  for (let i = 0; i < a.data.length; i++) {\n    const x = a.data[i];\n    acc += x * x;\n  }\n  return Math.sqrt(acc);\n}\n\nfunction mtoc2_norm2_complex(a) {\n  let acc = 0;\n  const im = a.imag;\n  for (let i = 0; i < a.data.length; i++) {\n    const re = a.data[i];\n    const imv = im !== undefined ? im[i] : 0;\n    acc += re * re + imv * imv;\n  }\n  return Math.sqrt(acc);\n}\n\nfunction mtoc2_norm_p_real(a, p) {\n  const n = a.data.length;\n  if (n === 0) return 0;\n  if (p === Infinity) {\n    let m = 0;\n    for (let i = 0; i < n; i++) {\n      const x = Math.abs(a.data[i]);\n      if (x > m) m = x;\n    }\n    return m;\n  }\n  if (p === -Infinity) {\n    let m = Math.abs(a.data[0]);\n    for (let i = 1; i < n; i++) {\n      const x = Math.abs(a.data[i]);\n      if (x < m) m = x;\n    }\n    return m;\n  }\n  if (p === 1) {\n    let acc = 0;\n    for (let i = 0; i < n; i++) acc += Math.abs(a.data[i]);\n    return acc;\n  }\n  if (p === 2) return mtoc2_norm2_real(a);\n  let acc = 0;\n  for (let i = 0; i < n; i++) acc += Math.pow(Math.abs(a.data[i]), p);\n  return Math.pow(acc, 1 / p);\n}\n\nfunction mtoc2_norm_p_complex(a, p) {\n  const n = a.data.length;\n  if (n === 0) return 0;\n  const im = a.imag;\n  const abs = i => {\n    const re = a.data[i];\n    const imv = im !== undefined ? im[i] : 0;\n    return Math.hypot(re, imv);\n  };\n  if (p === Infinity) {\n    let m = 0;\n    for (let i = 0; i < n; i++) {\n      const x = abs(i);\n      if (x > m) m = x;\n    }\n    return m;\n  }\n  if (p === -Infinity) {\n    let m = abs(0);\n    for (let i = 1; i < n; i++) {\n      const x = abs(i);\n      if (x < m) m = x;\n    }\n    return m;\n  }\n  if (p === 2) return mtoc2_norm2_complex(a);\n  if (p === 1) {\n    let acc = 0;\n    for (let i = 0; i < n; i++) acc += abs(i);\n    return acc;\n  }\n  let acc = 0;\n  for (let i = 0; i < n; i++) acc += Math.pow(abs(i), p);\n  return Math.pow(acc, 1 / p);\n}\n",
   "tensor_ones_nd.js": "// JS sibling of `tensor_ones_nd.h`. Fill the freshly-allocated\n// tensor with `1.0`.\n\n\nfunction mtoc2_tensor_ones_nd(ndim, dims) {\n  const t = mtoc2_tensor_alloc_nd(ndim, dims);\n  t.data.fill(1);\n  return t;\n}\n",
   "tensor_ones_square.js": "// JS sibling of `tensor_ones_square.h`. See `tensor_zeros_square.js`\n// for the rationale.\n\n\nfunction mtoc2_tensor_ones_square(n) {\n  return mtoc2_tensor_ones_nd(2, [n, n]);\n}\n",
@@ -82980,9 +83215,9 @@ function tileExact(src, padShape2, padReps) {
   }
   return out;
 }
-function padTo(arr, pad, len) {
+function padTo(arr, pad2, len) {
   const out = arr.slice();
-  while (out.length < len) out.push(pad);
+  while (out.length < len) out.push(pad2);
   return out;
 }
 function trimTrailingOnes(shape) {
@@ -96988,10 +97223,10 @@ import {
   chmodSync,
   realpathSync
 } from "fs";
-import { unzipSync } from "fflate";
+import { unzipSync, zipSync } from "fflate";
 import { execFileSync } from "child_process";
 import { homedir, tmpdir } from "os";
-import { join as join3, resolve, dirname, basename } from "path";
+import { join as join3, resolve, relative, dirname, basename } from "path";
 
 // src/cli-scan.ts
 import { readFileSync as readFileSync2, readdirSync, statSync } from "fs";
@@ -97474,6 +97709,69 @@ var NodeFileIOAdapter = class {
     }
     return extracted;
   }
+  zip(zipfilename, filenames, rootDir) {
+    const root = resolve(expandTilde(rootDir));
+    const dest = resolve(expandTilde(zipfilename));
+    const entries = {};
+    const names = [];
+    const toEntryName = (abs2) => {
+      const rel = relative(root, abs2);
+      if (rel.startsWith("..")) return basename(abs2);
+      return rel.split("\\").join("/");
+    };
+    const addFile = (abs2) => {
+      const name = toEntryName(abs2);
+      const buf = readFileSync3(abs2);
+      entries[name] = new Uint8Array(
+        buf.buffer,
+        buf.byteOffset,
+        buf.byteLength
+      );
+      names.push(name);
+    };
+    const addDir = (absDir) => {
+      for (const entry of readdirSync2(absDir)) {
+        const abs2 = join3(absDir, entry);
+        let st2;
+        try {
+          st2 = statSync2(abs2);
+        } catch {
+          continue;
+        }
+        if (st2.isDirectory()) addDir(abs2);
+        else addFile(abs2);
+      }
+    };
+    const addPath = (abs2) => {
+      const st2 = statSync2(abs2);
+      if (st2.isDirectory()) addDir(abs2);
+      else addFile(abs2);
+    };
+    for (const f of filenames) {
+      const p2 = expandTilde(f);
+      const abs2 = resolve(root, p2);
+      const last = basename(p2);
+      if (last.includes("*") || last.includes("?")) {
+        const dir = dirname(abs2);
+        const re = new RegExp(
+          "^" + last.replace(/[.+^${}()|[\]\\]/g, "\\$&").replace(/\*/g, ".*").replace(/\?/g, ".") + "$"
+        );
+        let dirEntries;
+        try {
+          dirEntries = readdirSync2(dir);
+        } catch {
+          continue;
+        }
+        for (const entry of dirEntries) {
+          if (re.test(entry)) addPath(join3(dir, entry));
+        }
+      } else {
+        addPath(abs2);
+      }
+    }
+    writeFileSync(dest, zipSync(entries));
+    return names;
+  }
   listDir(dirPath) {
     const p2 = expandTilde(dirPath);
     if (p2.includes("**")) {
@@ -97579,51 +97877,43 @@ var NodeFileIOAdapter = class {
     return results;
   }
   listDirRecursive(pattern) {
-    const idx = pattern.indexOf("**");
-    let baseDir = pattern.slice(0, idx);
-    if (baseDir.endsWith("/") || baseDir.endsWith("\\")) {
+    const norm2 = pattern.replace(/\\/g, "/");
+    const idx = norm2.indexOf("**");
+    let baseDir = norm2.slice(0, idx);
+    if (baseDir.endsWith("/")) {
       baseDir = baseDir.slice(0, -1);
     }
     if (!baseDir) baseDir = ".";
     const absBase = resolve(baseDir);
+    try {
+      if (!statSync2(absBase).isDirectory()) return [];
+    } catch {
+      return [];
+    }
+    const suffix = norm2.slice(idx);
+    const DIRS = "\0";
+    const ANY = "
";
+    const re = new RegExp(
+      "^" + suffix.replace(/[.+^${}()|[\]\\]/g, "\\$&").replace(/\*\*\//g, DIRS).replace(/\*\*/g, ANY).replace(/\*/g, "[^/]*").replace(/\?/g, "[^/]").replace(new RegExp(DIRS, "g"), "(?:.*/)?").replace(new RegExp(ANY, "g"), ".*") + "$"
+    );
     const results = [];
-    const walkDir = (dir) => {
+    const pushIfMatch = (rel, name, folder, bytes, isdir, mtimeMs) => {
+      if (re.test(rel)) results.push({ name, folder, bytes, isdir, mtimeMs });
+    };
+    const walkDir = (dir, rel) => {
+      const relPrefix = rel ? rel + "/" : "";
+      let dirMtime = 0;
       try {
-        const ds = statSync2(dir);
-        results.push({
-          name: ".",
-          folder: dir,
-          bytes: 0,
-          isdir: true,
-          mtimeMs: ds.mtimeMs
-        });
+        dirMtime = statSync2(dir).mtimeMs;
       } catch {
-        results.push({
-          name: ".",
-          folder: dir,
-          bytes: 0,
-          isdir: true,
-          mtimeMs: 0
-        });
       }
+      pushIfMatch(relPrefix + ".", ".", dir, 0, true, dirMtime);
+      let parentMtime = 0;
       try {
-        const ps = statSync2(resolve(dir, ".."));
-        results.push({
-          name: "..",
-          folder: dir,
-          bytes: 0,
-          isdir: true,
-          mtimeMs: ps.mtimeMs
-        });
+        parentMtime = statSync2(resolve(dir, "..")).mtimeMs;
       } catch {
-        results.push({
-          name: "..",
-          folder: dir,
-          bytes: 0,
-          isdir: true,
-          mtimeMs: 0
-        });
       }
+      pushIfMatch(relPrefix + "..", "..", dir, 0, true, parentMtime);
       let entries;
       try {
         entries = readdirSync2(dir);
@@ -97634,22 +97924,23 @@ var NodeFileIOAdapter = class {
       for (const entry of entries) {
         try {
           const es = statSync2(join3(dir, entry));
-          results.push({
-            name: entry,
-            folder: dir,
-            bytes: es.isDirectory() ? 0 : es.size,
-            isdir: es.isDirectory(),
-            mtimeMs: es.mtimeMs
-          });
+          pushIfMatch(
+            relPrefix + entry,
+            entry,
+            dir,
+            es.isDirectory() ? 0 : es.size,
+            es.isDirectory(),
+            es.mtimeMs
+          );
           if (es.isDirectory()) subdirs.push(entry);
         } catch {
         }
       }
       for (const sub2 of subdirs) {
-        walkDir(join3(dir, sub2));
+        walkDir(join3(dir, sub2), relPrefix + sub2);
       }
     };
-    walkDir(absBase);
+    walkDir(absBase, "");
     return results;
   }
   getEntry(fid) {
@@ -100069,7 +100360,7 @@ function unzipToFiles(zipData) {
 }
 
 // src/cli.ts
-import { zipSync } from "fflate";
+import { zipSync as zipSync2 } from "fflate";
 
 // src/numbl-core/native/qhull-node.ts
 var promise = null;
@@ -100094,8 +100385,8 @@ function variantAddonPath(fastMath) {
   const name = fastMath ? "numbl_addon.fastmath.node" : "numbl_addon.nofastmath.node";
   return join7(packageDir, "build", "Release", name);
 }
-function loadNativeAddon(fastMath) {
-  if (process.env.NUMBL_NO_NATIVE) return;
+function loadNativeAddon(fastMath, noNative = false) {
+  if (noNative || process.env.NUMBL_NO_NATIVE) return;
   const variantPath = variantAddonPath(fastMath);
   const flag = fastMath ? " --fast-math" : "";
   if (!existsSync3(variantPath)) {
@@ -100219,8 +100510,8 @@ function findTestFiles(dir) {
   walk2(dir);
   return results;
 }
-async function runTests(dir, optimization, fastMath) {
-  loadNativeAddon(fastMath ?? false);
+async function runTests(dir, optimization, fastMath, noNative) {
+  loadNativeAddon(fastMath ?? false, noNative ?? false);
   await loadQhullBackend();
   const absDir = resolve2(process.cwd(), dir);
   const testFiles = findTestFiles(absDir);
@@ -100232,7 +100523,7 @@ async function runTests(dir, optimization, fastMath) {
   let fail = 0;
   const failedScripts = [];
   for (const filepath of testFiles) {
-    const rel = relative(process.cwd(), filepath);
+    const rel = relative2(process.cwd(), filepath);
     const source = readFileSync6(filepath, "utf-8");
     const mainFileName = filepath;
     const scriptDir = dirname3(filepath);
@@ -100351,9 +100642,12 @@ Options (for run and eval):
                      across --opt levels). Off by default so all --opt
                      levels agree; --no-fast-math is accepted as the
                      (now default) opt-out.
+  --no-native        Do not load the native LAPACK addon; use the pure-JS
+                     fallbacks (same effect as NUMBL_NO_NATIVE=1)
 
 Environment variables:
-  NUMBL_PATH              Extra workspace directories (separated by ${delimiter})`);
+  NUMBL_PATH              Extra workspace directories (separated by ${delimiter})
+  NUMBL_NO_NATIVE         If set, do not load the native LAPACK addon`);
 }
 function parseOptions(args) {
   const opts = {
@@ -100372,7 +100666,8 @@ function parseOptions(args) {
     // and vectorizes transcendentals, so results drift from the JIT
     // kernels and diverge across --opt levels. Opt back in with
     // `--fast-math` (and a matching `build-addon --fast-math`).
-    fastMath: false
+    fastMath: false,
+    noNative: false
   };
   const envPath = process.env.NUMBL_PATH;
   if (envPath) {
@@ -100410,6 +100705,9 @@ function parseOptions(args) {
       case "--no-fast-math":
         opts.fastMath = false;
         break;
+      case "--no-native":
+        opts.noNative = true;
+        break;
       case "--dump-ast":
         opts.dumpAst = true;
         break;
@@ -100583,7 +100881,7 @@ async function cmdEval(args) {
   );
 }
 async function executeWithOptions(code, mainFileName, workspaceFiles, opts, searchPaths) {
-  loadNativeAddon(opts.fastMath);
+  loadNativeAddon(opts.fastMath, opts.noNative);
   await loadQhullBackend();
   const profiling = !!opts.profileOutput;
   const totalStart = performance.now();
@@ -100895,7 +101193,7 @@ function collectSiteFiles(projectDir, outDirAbs, ignore) {
     for (const e of entries) {
       const full = join7(dir, e.name);
       if (resolve2(full) === outDirAbs) continue;
-      const rel = relative(projectDir, full).split(sep).join("/");
+      const rel = relative2(projectDir, full).split(sep).join("/");
       if (matchesIgnore(rel)) continue;
       if (e.isDirectory()) {
         if (SITE_DEFAULT_IGNORE_DIRS.has(e.name)) continue;
@@ -101006,7 +101304,7 @@ async function cmdBuildSite(args) {
   zipEntries["numbl-project.json"] = new TextEncoder().encode(
     JSON.stringify(manifest, null, 2)
   );
-  const zip = zipSync(zipEntries, { level: 6 });
+  const zip = zipSync2(zipEntries, { level: 6 });
   writeFileSync4(join7(outDirAbs, "project.zip"), zip);
   const totalMb = zip.length / 1048576;
   console.log(
@@ -101071,7 +101369,7 @@ function cmdListBuiltins(flags) {
 }
 async function cmdRepl(args) {
   const opts = parseOptions(args);
-  loadNativeAddon(opts.fastMath);
+  loadNativeAddon(opts.fastMath, opts.noNative);
   await ensureMipCorePackage();
   const replSearchPaths = [...opts.extraPaths];
   const replFiles = opts.extraPaths.flatMap((d) => scanMFiles(d));
@@ -101104,7 +101402,7 @@ function cmdShowProfile(args) {
   }
   const data = JSON.parse(readFileSync6(filepath, "utf-8"));
   const fmt = (n, decimals = 1) => n.toFixed(decimals);
-  const pad = (s, w) => s.padStart(w);
+  const pad2 = (s, w) => s.padStart(w);
   const execTime = data.executionTimeMs || 1;
   const jitTime = data.jitCompileTimeMs ?? 0;
   const dispatches = data.dispatches ?? [];
@@ -101132,27 +101430,27 @@ function cmdShowProfile(args) {
   console.log("\nAll timings are self-time (exclusive) and mutually disjoint.");
   if (dispatches.length > 0) {
     const nameW = Math.max(20, ...dispatches.map((b) => b.name.length)) + 2;
-    const header = "Function".padEnd(nameW) + pad("Calls", 10) + pad("Self (ms)", 13) + pad("Avg (ms)", 11) + pad("% Exec", 9);
+    const header = "Function".padEnd(nameW) + pad2("Calls", 10) + pad2("Self (ms)", 13) + pad2("Avg (ms)", 11) + pad2("% Exec", 9);
     console.log("\nDispatched functions (user-defined / workspace):");
     console.log(header);
     console.log("\u2500".repeat(header.length));
     for (const b of dispatches) {
       const avg = b.callCount > 0 ? b.totalTimeMs / b.callCount : 0;
       const pct = b.totalTimeMs / execTime * 100;
-      const line = b.name.padEnd(nameW) + pad(String(b.callCount), 10) + pad(fmt(b.totalTimeMs, 2), 13) + pad(fmt(avg, 3), 11) + pad(fmt(pct, 1) + "%", 9);
+      const line = b.name.padEnd(nameW) + pad2(String(b.callCount), 10) + pad2(fmt(b.totalTimeMs, 2), 13) + pad2(fmt(avg, 3), 11) + pad2(fmt(pct, 1) + "%", 9);
       console.log(line);
     }
   }
   if (builtins.length > 0) {
     const nameW = Math.max(20, ...builtins.map((b) => b.name.length)) + 2;
-    const header = "Builtin".padEnd(nameW) + pad("Interp", 10) + pad("Fallback", 10) + pad("JIT", 10) + pad("Self (ms)", 13) + pad("% Exec", 9);
+    const header = "Builtin".padEnd(nameW) + pad2("Interp", 10) + pad2("Fallback", 10) + pad2("JIT", 10) + pad2("Self (ms)", 13) + pad2("% Exec", 9);
     console.log("\nBuiltin functions (call counts by source):");
     console.log(header);
     console.log("\u2500".repeat(header.length));
     for (const b of builtins) {
       const totalTime = b.fallback.totalTimeMs + b.interp.totalTimeMs + b.jit.totalTimeMs;
       const pct = totalTime / execTime * 100;
-      const line = b.name.padEnd(nameW) + pad(String(b.interp.callCount), 10) + pad(String(b.fallback.callCount), 10) + pad(String(b.jit.callCount), 10) + pad(fmt(totalTime, 2), 13) + pad(fmt(pct, 1) + "%", 9);
+      const line = b.name.padEnd(nameW) + pad2(String(b.interp.callCount), 10) + pad2(String(b.fallback.callCount), 10) + pad2(String(b.jit.callCount), 10) + pad2(fmt(totalTime, 2), 13) + pad2(fmt(pct, 1) + "%", 9);
       console.log(line);
     }
   } else {
@@ -101162,13 +101460,13 @@ function cmdShowProfile(args) {
   if (hotLoops.length > 0) {
     hotLoops.sort((a, b) => b.totalTimeMs - a.totalTimeMs);
     const locW = Math.max(20, ...hotLoops.map((l) => `${l.file}:${l.line}`.length)) + 2;
-    const header = "Location".padEnd(locW) + pad("Kind", 8) + pad("Calls", 8) + pad("Max Iters", 12) + pad("Total (ms)", 13);
+    const header = "Location".padEnd(locW) + pad2("Kind", 8) + pad2("Calls", 8) + pad2("Max Iters", 12) + pad2("Total (ms)", 13);
     console.log("\nInterpreted loops (>1000 iterations):");
     console.log(header);
     console.log("\u2500".repeat(header.length));
     for (const l of hotLoops) {
       const loc = `${l.file}:${l.line}`;
-      const line = loc.padEnd(locW) + pad(l.kind, 8) + pad(String(l.callCount ?? 1), 8) + pad(String(l.iterations), 12) + pad(fmt(l.totalTimeMs ?? 0, 2), 13);
+      const line = loc.padEnd(locW) + pad2(l.kind, 8) + pad2(String(l.callCount ?? 1), 8) + pad2(String(l.iterations), 12) + pad2(fmt(l.totalTimeMs ?? 0, 2), 13);
       console.log(line);
     }
   }
@@ -101230,7 +101528,12 @@ async function main() {
     case "run-tests": {
       const testOpts = parseOptions(rest);
       const dir = testOpts.positional.length > 0 ? testOpts.positional[0] : join7(packageDir, "numbl_test_scripts");
-      await runTests(dir, testOpts.optimization, testOpts.fastMath);
+      await runTests(
+        dir,
+        testOpts.optimization,
+        testOpts.fastMath,
+        testOpts.noNative
+      );
       break;
     }
     case "build-addon":