numbl 0.1.2 → 0.1.3

package/dist-cli/cli.js

@@ -195,6 +195,179 @@ var kstr = (value) => {
   return "unknown";
 };
 
+// src/numbl-core/native/lapack-bridge.ts
+var NATIVE_ADDON_EXPECTED_VERSION = 1;
+var _bridge = null;
+function setLapackBridge(bridge) {
+  _bridge = bridge;
+}
+function getLapackBridge() {
+  return _bridge;
+}
+
+// src/numbl-core/runtime/constructors.ts
+var RTV = {
+  num(value) {
+    return value;
+  },
+  tensor(data, shape, imag) {
+    const d = data instanceof FloatXArray ? data : new FloatXArray(data);
+    const im2 = imag ? imag instanceof FloatXArray ? imag : new FloatXArray(imag) : void 0;
+    const s = [...shape];
+    while (s.length > 2 && s[s.length - 1] === 1) s.pop();
+    return { kind: "tensor", data: d, imag: im2, shape: s, _rc: 1 };
+  },
+  /** Fast tensor constructor — data must be FloatXArray, shape already normalized (no trailing singletons). */
+  tensorRaw(data, shape) {
+    return { kind: "tensor", data, imag: void 0, shape, _rc: 1 };
+  },
+  /** Create a scalar tensor (1x1) */
+  scalar(value) {
+    return value;
+  },
+  /** Create a row vector [1 x n] */
+  row(data, imag) {
+    const im2 = imag ? new FloatXArray(imag) : void 0;
+    return {
+      kind: "tensor",
+      data: new FloatXArray(data),
+      imag: im2,
+      shape: [1, data.length],
+      _rc: 1
+    };
+  },
+  /** Create a column vector [n x 1] */
+  col(data, imag) {
+    const im2 = imag ? new FloatXArray(imag) : void 0;
+    return {
+      kind: "tensor",
+      data: new FloatXArray(data),
+      imag: im2,
+      shape: [data.length, 1],
+      _rc: 1
+    };
+  },
+  /** Create a matrix from row-major data */
+  matrix(rows, cols, data, imag) {
+    const d = data instanceof FloatXArray ? data : new FloatXArray(data);
+    const im2 = imag ? imag instanceof FloatXArray ? imag : new FloatXArray(imag) : void 0;
+    return { kind: "tensor", data: d, imag: im2, shape: [rows, cols], _rc: 1 };
+  },
+  string(value) {
+    return value;
+  },
+  char(value) {
+    return { kind: "char", value };
+  },
+  logical(value) {
+    return value;
+  },
+  cell(data, shape) {
+    return { kind: "cell", data, shape: [...shape], _rc: 1 };
+  },
+  struct(fields) {
+    const map = fields instanceof Map ? fields : new Map(Object.entries(fields));
+    return { kind: "struct", fields: map };
+  },
+  func(name, impl, captures = []) {
+    return { kind: "function", name, impl, captures };
+  },
+  classInstance(className, propertyNames, isHandleClass2 = false, defaults) {
+    const fields = /* @__PURE__ */ new Map();
+    for (const name of propertyNames) {
+      fields.set(
+        name,
+        defaults?.get(name) ?? RTV.tensor(new FloatXArray(0), [0, 0])
+      );
+    }
+    return {
+      kind: "class_instance",
+      className,
+      fields,
+      isHandleClass: isHandleClass2
+    };
+  },
+  complex(re2, im2) {
+    return { kind: "complex_number", re: re2, im: im2 };
+  },
+  dummyHandle() {
+    return { kind: "dummy_handle" };
+  },
+  structArray(fieldNames, elements) {
+    return { kind: "struct_array", fieldNames, elements };
+  },
+  sparseMatrix(m, n, ir, jc, pr, pi) {
+    return { kind: "sparse_matrix", m, n, ir, jc, pr, pi, _rc: 1 };
+  },
+  dictionary(entries, keyType, valueType) {
+    return {
+      kind: "dictionary",
+      entries: entries ?? /* @__PURE__ */ new Map(),
+      keyType,
+      valueType
+    };
+  }
+};
+var getItemTypeFromRuntimeValue = (value) => {
+  if (isRuntimeNumber(value)) {
+    return { kind: "Number" };
+  }
+  if (isRuntimeLogical(value)) {
+    return { kind: "Boolean" };
+  }
+  if (isRuntimeString(value)) {
+    return { kind: "String" };
+  }
+  switch (value.kind) {
+    case "complex_number":
+      return { kind: "ComplexNumber" };
+    case "char":
+      return { kind: "Char" };
+    case "tensor":
+      return {
+        kind: "Tensor",
+        isComplex: value.imag !== void 0,
+        isLogical: value._isLogical || void 0
+      };
+    case "cell":
+      return {
+        kind: "Cell",
+        elementType: value.data.length > 0 ? getItemTypeFromRuntimeValue(value.data[0]) : { kind: "Unknown" },
+        // length: value.shape.reduce((a, b) => a * b, 1),
+        length: "unknown"
+        // do not include length or we get too many jits
+      };
+    case "struct": {
+      const knownFields = {};
+      for (const [k, v] of value.fields) {
+        knownFields[k] = getItemTypeFromRuntimeValue(v);
+      }
+      return { kind: "Struct", knownFields };
+    }
+    case "function":
+      return {
+        kind: "Function",
+        params: [],
+        returns: { kind: "Unknown" }
+      };
+    case "class_instance":
+      return {
+        kind: "ClassInstance",
+        className: value.className
+      };
+    case "dummy_handle":
+      return { kind: "DummyHandle" };
+    case "struct_array":
+      return { kind: "Unknown" };
+    case "sparse_matrix":
+      return { kind: "SparseMatrix", isComplex: value.pi !== void 0 };
+    case "dictionary":
+      return { kind: "Dictionary" };
+    default:
+      return { kind: "Unknown" };
+  }
+};
+
 // src/numbl-core/interpreter/jit/jitTypes.ts
 function signFromNumber(v) {
   if (v > 0) return "positive";
@@ -637,169 +810,6 @@ function extractSnippet(source, offset, contextLines = 2) {
   return snippetLines.join("\n");
 }
 
-// src/numbl-core/runtime/constructors.ts
-var RTV = {
-  num(value) {
-    return value;
-  },
-  tensor(data, shape, imag) {
-    const d = data instanceof FloatXArray ? data : new FloatXArray(data);
-    const im2 = imag ? imag instanceof FloatXArray ? imag : new FloatXArray(imag) : void 0;
-    const s = [...shape];
-    while (s.length > 2 && s[s.length - 1] === 1) s.pop();
-    return { kind: "tensor", data: d, imag: im2, shape: s, _rc: 1 };
-  },
-  /** Fast tensor constructor — data must be FloatXArray, shape already normalized (no trailing singletons). */
-  tensorRaw(data, shape) {
-    return { kind: "tensor", data, imag: void 0, shape, _rc: 1 };
-  },
-  /** Create a scalar tensor (1x1) */
-  scalar(value) {
-    return value;
-  },
-  /** Create a row vector [1 x n] */
-  row(data, imag) {
-    const im2 = imag ? new FloatXArray(imag) : void 0;
-    return {
-      kind: "tensor",
-      data: new FloatXArray(data),
-      imag: im2,
-      shape: [1, data.length],
-      _rc: 1
-    };
-  },
-  /** Create a column vector [n x 1] */
-  col(data, imag) {
-    const im2 = imag ? new FloatXArray(imag) : void 0;
-    return {
-      kind: "tensor",
-      data: new FloatXArray(data),
-      imag: im2,
-      shape: [data.length, 1],
-      _rc: 1
-    };
-  },
-  /** Create a matrix from row-major data */
-  matrix(rows, cols, data, imag) {
-    const d = data instanceof FloatXArray ? data : new FloatXArray(data);
-    const im2 = imag ? imag instanceof FloatXArray ? imag : new FloatXArray(imag) : void 0;
-    return { kind: "tensor", data: d, imag: im2, shape: [rows, cols], _rc: 1 };
-  },
-  string(value) {
-    return value;
-  },
-  char(value) {
-    return { kind: "char", value };
-  },
-  logical(value) {
-    return value;
-  },
-  cell(data, shape) {
-    return { kind: "cell", data, shape: [...shape], _rc: 1 };
-  },
-  struct(fields) {
-    const map = fields instanceof Map ? fields : new Map(Object.entries(fields));
-    return { kind: "struct", fields: map };
-  },
-  func(name, impl, captures = []) {
-    return { kind: "function", name, impl, captures };
-  },
-  classInstance(className, propertyNames, isHandleClass2 = false, defaults) {
-    const fields = /* @__PURE__ */ new Map();
-    for (const name of propertyNames) {
-      fields.set(
-        name,
-        defaults?.get(name) ?? RTV.tensor(new FloatXArray(0), [0, 0])
-      );
-    }
-    return {
-      kind: "class_instance",
-      className,
-      fields,
-      isHandleClass: isHandleClass2
-    };
-  },
-  complex(re2, im2) {
-    return { kind: "complex_number", re: re2, im: im2 };
-  },
-  dummyHandle() {
-    return { kind: "dummy_handle" };
-  },
-  structArray(fieldNames, elements) {
-    return { kind: "struct_array", fieldNames, elements };
-  },
-  sparseMatrix(m, n, ir, jc, pr, pi) {
-    return { kind: "sparse_matrix", m, n, ir, jc, pr, pi, _rc: 1 };
-  },
-  dictionary(entries, keyType, valueType) {
-    return {
-      kind: "dictionary",
-      entries: entries ?? /* @__PURE__ */ new Map(),
-      keyType,
-      valueType
-    };
-  }
-};
-var getItemTypeFromRuntimeValue = (value) => {
-  if (isRuntimeNumber(value)) {
-    return { kind: "Number" };
-  }
-  if (isRuntimeLogical(value)) {
-    return { kind: "Boolean" };
-  }
-  if (isRuntimeString(value)) {
-    return { kind: "String" };
-  }
-  switch (value.kind) {
-    case "complex_number":
-      return { kind: "ComplexNumber" };
-    case "char":
-      return { kind: "Char" };
-    case "tensor":
-      return {
-        kind: "Tensor",
-        isComplex: value.imag !== void 0,
-        isLogical: value._isLogical || void 0
-      };
-    case "cell":
-      return {
-        kind: "Cell",
-        elementType: value.data.length > 0 ? getItemTypeFromRuntimeValue(value.data[0]) : { kind: "Unknown" },
-        // length: value.shape.reduce((a, b) => a * b, 1),
-        length: "unknown"
-        // do not include length or we get too many jits
-      };
-    case "struct": {
-      const knownFields = {};
-      for (const [k, v] of value.fields) {
-        knownFields[k] = getItemTypeFromRuntimeValue(v);
-      }
-      return { kind: "Struct", knownFields };
-    }
-    case "function":
-      return {
-        kind: "Function",
-        params: [],
-        returns: { kind: "Unknown" }
-      };
-    case "class_instance":
-      return {
-        kind: "ClassInstance",
-        className: value.className
-      };
-    case "dummy_handle":
-      return { kind: "DummyHandle" };
-    case "struct_array":
-      return { kind: "Unknown" };
-    case "sparse_matrix":
-      return { kind: "SparseMatrix", isComplex: value.pi !== void 0 };
-    case "dictionary":
-      return { kind: "Dictionary" };
-    default:
-      return { kind: "Unknown" };
-  }
-};
-
 // src/numbl-core/runtime/utils.ts
 function tensorSize2D(t) {
   const s = t.shape;
@@ -1718,6 +1728,28 @@ function binaryNumberOnly(argTypes) {
     return null;
   return [{ kind: "number" }];
 }
+var nativeUnaryOpCode = /* @__PURE__ */ new Map([
+  [Math.exp, 0],
+  [Math.log, 1],
+  [Math.log2, 2],
+  [Math.log10, 3],
+  [Math.sqrt, 4],
+  [Math.abs, 5],
+  [Math.floor, 6],
+  [Math.ceil, 7],
+  [Math.round, 8],
+  [Math.trunc, 9],
+  [Math.sin, 10],
+  [Math.cos, 11],
+  [Math.tan, 12],
+  [Math.asin, 13],
+  [Math.acos, 14],
+  [Math.atan, 15],
+  [Math.sinh, 16],
+  [Math.cosh, 17],
+  [Math.tanh, 18],
+  [Math.sign, 19]
+]);
 function applyUnaryElemwise(v, realFn, complexFn, name) {
   if (isRuntimeSparseMatrix(v))
     return applyUnaryElemwise(sparseToDense(v), realFn, complexFn, name);
@@ -1730,6 +1762,14 @@ function applyUnaryElemwise(v, realFn, complexFn, name) {
   if (isRuntimeTensor(v)) {
     const n = v.data.length;
     if (!v.imag) {
+      const opCode = nativeUnaryOpCode.get(realFn);
+      if (opCode !== void 0) {
+        const bridge = getLapackBridge();
+        if (bridge?.unaryElemwise && v.data instanceof Float64Array) {
+          const result = bridge.unaryElemwise(v.data, opCode);
+          return RTV.tensorRaw(result, v.shape.slice());
+        }
+      }
       const out = new FloatXArray(n);
       for (let i = 0; i < n; i++) out[i] = realFn(v.data[i]);
       return makeTensor(out, void 0, v.shape.slice());
@@ -3618,16 +3658,6 @@ function toString(v) {
   }
 }
 
-// src/numbl-core/native/lapack-bridge.ts
-var NATIVE_ADDON_EXPECTED_VERSION = 1;
-var _bridge = null;
-function setLapackBridge(bridge) {
-  _bridge = bridge;
-}
-function getLapackBridge() {
-  return _bridge;
-}
-
 // src/ts-lapack/src/utils/xerbla.ts
 function xerbla(srname, info) {
   throw new Error(
@@ -22461,6 +22491,18 @@ function tryNativeElemwiseReal(at, bt, opCode) {
   );
   return RTV.tensorRaw(result, at.shape);
 }
+function tryNativeElemwiseScalar(scalar, tensor, opCode, scalarOnLeft) {
+  if (tensor.imag) return null;
+  const bridge = getLapackBridge();
+  if (!bridge?.elemwiseScalar) return null;
+  const result = bridge.elemwiseScalar(
+    scalar,
+    tensor.data,
+    opCode,
+    scalarOnLeft
+  );
+  return RTV.tensorRaw(result, tensor.shape);
+}
 function tensorElemwiseComplex(at, bt, opCode, jsOp) {
   const bridge = getLapackBridge();
   if (bridge?.elemwiseComplex) {
@@ -22531,6 +22573,13 @@ function mAdd(a, b) {
       im: aIm + bIm
     }));
   }
+  if (isRuntimeNumber(a) && isRuntimeTensor(b)) {
+    const nr = tryNativeElemwiseScalar(a, b, ELEMWISE_ADD, true);
+    if (nr) return nr;
+  } else if (isRuntimeTensor(a) && isRuntimeNumber(b)) {
+    const nr = tryNativeElemwiseScalar(b, a, ELEMWISE_ADD, false);
+    if (nr) return nr;
+  }
   return binaryOp(a, b, (x, y) => x + y);
 }
 function mSub(a, b) {
@@ -22563,6 +22612,13 @@ function mSub(a, b) {
       im: aIm - bIm
     }));
   }
+  if (isRuntimeNumber(a) && isRuntimeTensor(b)) {
+    const nr = tryNativeElemwiseScalar(a, b, ELEMWISE_SUB, true);
+    if (nr) return nr;
+  } else if (isRuntimeTensor(a) && isRuntimeNumber(b)) {
+    const nr = tryNativeElemwiseScalar(b, a, ELEMWISE_SUB, false);
+    if (nr) return nr;
+  }
   return binaryOp(a, b, (x, y) => x - y);
 }
 function mMul(a, b) {
@@ -22579,6 +22635,13 @@ function mMul(a, b) {
       im: aRe * bIm + aIm * bRe
     }));
   }
+  if (isRuntimeNumber(a) && isRuntimeTensor(b)) {
+    const nr = tryNativeElemwiseScalar(a, b, ELEMWISE_MUL, true);
+    if (nr) return nr;
+  } else if (isRuntimeTensor(a) && isRuntimeNumber(b)) {
+    const nr = tryNativeElemwiseScalar(b, a, ELEMWISE_MUL, false);
+    if (nr) return nr;
+  }
   return binaryOp(a, b, (x, y) => x * y);
 }
 function mElemMul(a, b) {
@@ -22611,6 +22674,13 @@ function mElemMul(a, b) {
       im: aRe * bIm + aIm * bRe
     }));
   }
+  if (isRuntimeNumber(a) && isRuntimeTensor(b)) {
+    const nr = tryNativeElemwiseScalar(a, b, ELEMWISE_MUL, true);
+    if (nr) return nr;
+  } else if (isRuntimeTensor(a) && isRuntimeNumber(b)) {
+    const nr = tryNativeElemwiseScalar(b, a, ELEMWISE_MUL, false);
+    if (nr) return nr;
+  }
   return binaryOp(a, b, (x, y) => x * y);
 }
 function mDiv(a, b) {
@@ -22649,6 +22719,13 @@ function mElemDiv(a, b) {
   if (isComplexOrMixed(a, b)) {
     return complexBinaryOp(a, b, complexDivide);
   }
+  if (isRuntimeNumber(a) && isRuntimeTensor(b)) {
+    const nr = tryNativeElemwiseScalar(a, b, ELEMWISE_DIV, true);
+    if (nr) return nr;
+  } else if (isRuntimeTensor(a) && isRuntimeNumber(b)) {
+    const nr = tryNativeElemwiseScalar(b, a, ELEMWISE_DIV, false);
+    if (nr) return nr;
+  }
   return binaryOp(a, b, (x, y) => x / y);
 }
 function mLeftDiv(a, b) {
@@ -28565,6 +28642,44 @@ registerIBuiltin({
     }
   })
 });
+var WEBOPTIONS_DEFAULTS = {
+  CharacterEncoding: RTV.char("auto"),
+  UserAgent: RTV.char("numbl"),
+  Timeout: 5,
+  Username: RTV.char(""),
+  Password: RTV.char(""),
+  KeyName: RTV.char(""),
+  KeyValue: RTV.char(""),
+  ContentType: RTV.char("auto"),
+  MediaType: RTV.char("auto"),
+  RequestMethod: RTV.char("auto"),
+  ArrayFormat: RTV.char("csv"),
+  CertificateFilename: RTV.char("default")
+};
+registerIBuiltin({
+  name: "weboptions",
+  resolve: () => ({
+    outputTypes: [{ kind: "struct" }],
+    apply: (args) => {
+      const fields = /* @__PURE__ */ new Map();
+      for (const [k, v] of Object.entries(WEBOPTIONS_DEFAULTS)) {
+        fields.set(k, v);
+      }
+      if (args.length % 2 !== 0)
+        throw new RuntimeError("weboptions: expected name-value pairs");
+      for (let i = 0; i < args.length; i += 2) {
+        const key = args[i];
+        if (!isRuntimeChar(key) && !isRuntimeString(key))
+          throw new RuntimeError("weboptions: option name must be a string");
+        const name = isRuntimeChar(key) ? key.value : key;
+        if (!(name in WEBOPTIONS_DEFAULTS))
+          throw new RuntimeError(`weboptions: unknown option '${name}'`);
+        fields.set(name, args[i + 1]);
+      }
+      return RTV.struct(Object.fromEntries(fields));
+    }
+  })
+});
 registerIBuiltin({
   name: "odeget",
   resolve: () => ({
@@ -29119,6 +29234,130 @@ registerIBuiltin({
   }
 });
 
+// src/numbl-core/interpreter/builtins/time-system.ts
+var ticTime = 0;
+function getTicTime() {
+  return ticTime;
+}
+defineBuiltin({
+  name: "tic",
+  cases: [
+    {
+      match: (argTypes) => argTypes.length === 0 ? [{ kind: "number" }] : null,
+      apply: () => {
+        ticTime = performance.now();
+        return RTV.num(ticTime / 1e3);
+      }
+    }
+  ]
+});
+defineBuiltin({
+  name: "clock",
+  cases: [
+    {
+      match: (argTypes) => argTypes.length === 0 ? [{ kind: "tensor", isComplex: false, shape: [1, 6] }] : null,
+      apply: () => {
+        const now = /* @__PURE__ */ new Date();
+        return RTV.tensor(
+          new FloatXArray([
+            now.getFullYear(),
+            now.getMonth() + 1,
+            now.getDate(),
+            now.getHours(),
+            now.getMinutes(),
+            now.getSeconds() + now.getMilliseconds() / 1e3
+          ]),
+          [1, 6]
+        );
+      }
+    }
+  ]
+});
+defineBuiltin({
+  name: "etime",
+  cases: [
+    {
+      match: (argTypes) => argTypes.length === 2 ? [{ kind: "number" }] : null,
+      apply: (args) => {
+        const t1 = args[0];
+        const t0 = args[1];
+        if (!isRuntimeTensor(t1) || !isRuntimeTensor(t0))
+          throw new RuntimeError("etime: arguments must be clock vectors");
+        const toMs = (t) => {
+          const d = new Date(
+            t.data[0],
+            t.data[1] - 1,
+            t.data[2],
+            t.data[3],
+            t.data[4],
+            Math.floor(t.data[5]),
+            t.data[5] % 1 * 1e3
+          );
+          return d.getTime();
+        };
+        return RTV.num((toMs(t1) - toMs(t0)) / 1e3);
+      }
+    }
+  ]
+});
+defineBuiltin({
+  name: "version",
+  cases: [
+    {
+      match: (argTypes) => argTypes.length === 0 ? [{ kind: "char" }] : null,
+      apply: () => RTV.char("9.14.0")
+    }
+  ]
+});
+function getComputerStrings() {
+  if (typeof process === "undefined") {
+    return { str: "BROWSER", arch: "browser" };
+  }
+  const platform = process.platform;
+  const cpuArch = process.arch;
+  if (platform === "win32") return { str: "PCWIN64", arch: "win64" };
+  if (platform === "darwin") {
+    if (cpuArch === "arm64") return { str: "MACA64", arch: "maca64" };
+    return { str: "MACI64", arch: "maci64" };
+  }
+  return { str: "GLNXA64", arch: "glnxa64" };
+}
+registerIBuiltin({
+  name: "computer",
+  resolve: () => ({
+    outputTypes: [{ kind: "unknown" }],
+    apply: (args, nargout) => {
+      const info = getComputerStrings();
+      if (args.length >= 1 && isRuntimeChar(args[0])) {
+        const arg = toString(args[0]);
+        if (arg === "arch") return RTV.char(info.arch);
+        throw new RuntimeError(`computer: unknown argument '${arg}'`);
+      }
+      const maxsize = 2 ** 48 - 1;
+      const endian = RTV.char("L");
+      if (nargout <= 1) return RTV.char(info.str);
+      if (nargout === 2) return [RTV.char(info.str), RTV.num(maxsize)];
+      return [RTV.char(info.str), RTV.num(maxsize), endian];
+    }
+  })
+});
+var _platform = typeof process !== "undefined" ? process.platform : "linux";
+for (const [name, val] of [
+  ["ismac", _platform === "darwin"],
+  ["ispc", _platform === "win32"],
+  ["isunix", _platform !== "win32"]
+]) {
+  defineBuiltin({
+    name,
+    cases: [
+      {
+        match: (argTypes) => argTypes.length === 0 ? [{ kind: "boolean" }] : null,
+        apply: () => RTV.logical(val)
+      }
+    ]
+  });
+}
+
 // src/numbl-core/runtime/plotUtils.ts
 var COLOR_SHORT = {
   r: [1, 0, 0],
@@ -31451,7 +31690,8 @@ var SPECIAL_BUILTIN_NAMES = [
   "cd",
   "ode45",
   "ode23",
-  "deval"
+  "deval",
+  "toc"
 ];
 
 // src/numbl-core/runtime/specialBuiltins.ts
@@ -31502,6 +31742,14 @@ function registerSpecialBuiltins(rt) {
     }
     return 0;
   });
+  registerSpecial("toc", (nargout) => {
+    const elapsed = (performance.now() - getTicTime()) / 1e3;
+    if (nargout === 0) {
+      rt.output(`Elapsed time is ${elapsed.toFixed(6)} seconds.
+`);
+    }
+    return RTV.num(elapsed);
+  });
   registerSpecial("warning", (_nargout, args) => {
     if (args.length === 0) return RTV.num(0);
     const margs = args.map((a) => ensureRuntimeValue(a));
@@ -32140,6 +32388,40 @@ function registerSpecialBuiltins(rt) {
       RTV.char("")
     ];
   });
+  function extractWebOptions(arg) {
+    if (!isRuntimeStruct(arg)) return void 0;
+    const s = arg;
+    if (!s.fields.has("Timeout")) return void 0;
+    const opts = {};
+    const t = s.fields.get("Timeout");
+    if (t !== void 0) opts.timeout = toNumber(t);
+    const rm = s.fields.get("RequestMethod");
+    if (rm !== void 0) {
+      const v = toString(rm);
+      if (v !== "auto") opts.requestMethod = v;
+    }
+    const un = s.fields.get("Username");
+    if (un !== void 0) {
+      const v = toString(un);
+      if (v) opts.username = v;
+    }
+    const pw = s.fields.get("Password");
+    if (pw !== void 0) {
+      const v = toString(pw);
+      if (v) opts.password = v;
+    }
+    const kn = s.fields.get("KeyName");
+    if (kn !== void 0) {
+      const v = toString(kn);
+      if (v) opts.keyName = v;
+    }
+    const kv = s.fields.get("KeyValue");
+    if (kv !== void 0) {
+      const v = toString(kv);
+      if (v) opts.keyValue = v;
+    }
+    return opts;
+  }
   registerSpecial("websave", (_nargout, args) => {
     const io = requireFileIO();
     const margs = args.map((a) => ensureRuntimeValue(a));
@@ -32149,8 +32431,14 @@ function registerSpecialBuiltins(rt) {
       throw new RuntimeError("websave is not available in this environment");
     const filename = toString(margs[0]);
     let url = toString(margs[1]);
+    let webOpts;
+    let queryEnd = margs.length;
+    if (margs.length > 2) {
+      webOpts = extractWebOptions(margs[margs.length - 1]);
+      if (webOpts) queryEnd = margs.length - 1;
+    }
     const queryParts = [];
-    for (let i = 2; i + 1 < margs.length; i += 2) {
+    for (let i = 2; i + 1 < queryEnd; i += 2) {
       const name = encodeURIComponent(toString(margs[i]));
       const value = encodeURIComponent(toString(margs[i + 1]));
       queryParts.push(`${name}=${value}`);
@@ -32159,7 +32447,7 @@ function registerSpecialBuiltins(rt) {
       const sep = url.includes("?") ? "&" : "?";
       url += sep + queryParts.join("&");
     }
-    io.websave(url, filename);
+    io.websave(url, filename, webOpts);
     return RTV.char(filename);
   });
   registerSpecial("webread", (_nargout, args) => {
@@ -32170,8 +32458,14 @@ function registerSpecialBuiltins(rt) {
     if (!io.webread)
       throw new RuntimeError("webread is not available in this environment");
     let url = toString(margs[0]);
+    let webOpts;
+    let queryEnd = margs.length;
+    if (margs.length > 1) {
+      webOpts = extractWebOptions(margs[margs.length - 1]);
+      if (webOpts) queryEnd = margs.length - 1;
+    }
     const queryParts = [];
-    for (let i = 1; i + 1 < margs.length; i += 2) {
+    for (let i = 1; i + 1 < queryEnd; i += 2) {
       const name = encodeURIComponent(toString(margs[i]));
       const value = encodeURIComponent(toString(margs[i + 1]));
       queryParts.push(`${name}=${value}`);
@@ -32180,7 +32474,7 @@ function registerSpecialBuiltins(rt) {
       const sep = url.includes("?") ? "&" : "?";
       url += sep + queryParts.join("&");
     }
-    const text = io.webread(url);
+    const text = io.webread(url, webOpts);
     try {
       const parsed = JSON.parse(text);
       return convertJsonValue(parsed);
@@ -38960,7 +39254,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length < 1 || argTypes.length > 2)
+        if (nargout > 1 || argTypes.length < 1 || argTypes.length > 2)
           return null;
         if (!isNumericJitType(argTypes[0])) return null;
         return [NUM];
@@ -39077,7 +39371,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length < 1 || argTypes.length > 3)
+        if (nargout > 1 || argTypes.length < 1 || argTypes.length > 3)
           return null;
         if (!isNumericJitType(argTypes[0])) return null;
         const a = argTypes[0];
@@ -39156,7 +39450,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length !== 2) return null;
+        if (nargout > 1 || argTypes.length !== 2) return null;
         if (!isNumericJitType(argTypes[0]) || !isNumericJitType(argTypes[1]))
           return null;
         const hasTensor = argTypes.some((t) => t.kind === "tensor");
@@ -39240,7 +39534,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (argTypes.length !== 1 || nargout !== 1) return null;
+        if (argTypes.length !== 1 || nargout > 1) return null;
         if (!isNumericJitType(argTypes[0])) return null;
         const hasComplex = argTypes[0].kind === "complex_or_number" || argTypes[0].kind === "tensor" && argTypes[0].isComplex;
         return [hasComplex ? COMPLEX_OR_NUM : NUM];
@@ -39352,7 +39646,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (argTypes.length !== 1 || nargout !== 1) return null;
+        if (argTypes.length !== 1 || nargout > 1) return null;
         if (!isNumericJitType(argTypes[0])) return null;
         return [NUM];
       },
@@ -39377,7 +39671,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1) return null;
+        if (nargout > 1) return null;
         if (argTypes.length !== 2 && argTypes.length !== 3) return null;
         if (!isNumericJitType(argTypes[0]) || !isNumericJitType(argTypes[1]))
           return null;
@@ -39443,7 +39737,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (argTypes.length !== 1 || nargout !== 1) return null;
+        if (argTypes.length !== 1 || nargout > 1) return null;
         if (!isNumericJitType(argTypes[0])) return null;
         const a = argTypes[0];
         if (a.kind === "number" || a.kind === "boolean") return [NUM];
@@ -39556,11 +39850,11 @@ function invComplexJS(dataRe, dataIm, n) {
 registerIBuiltin({
   name: "svd",
   resolve: (argTypes, nargout) => {
-    if (nargout < 1 || nargout > 3 || argTypes.length < 1 || argTypes.length > 2)
+    if (nargout < 0 || nargout > 3 || argTypes.length < 1 || argTypes.length > 2)
       return null;
     if (!isNumericJitType(argTypes[0])) return null;
     const c = tensorType();
-    if (nargout === 1)
+    if (nargout <= 1)
       return { outputTypes: [c], apply: (args, n) => svdApply(args, n) };
     return { outputTypes: [c, c, c], apply: (args, n) => svdApply(args, n) };
   }
@@ -39571,7 +39865,7 @@ function svdApply(args, nargout) {
   const A = args[0];
   if (isRuntimeNumber(A)) {
     const val = Math.abs(A);
-    if (nargout === 1) return RTV.tensor(new FloatXArray([val]), [1, 1]);
+    if (nargout <= 1) return RTV.tensor(new FloatXArray([val]), [1, 1]);
     return [
       RTV.tensor(new FloatXArray([A >= 0 ? 1 : -1]), [1, 1]),
       RTV.tensor(new FloatXArray([val]), [1, 1]),
@@ -39598,7 +39892,7 @@ function svdApply(args, nargout) {
       nargout === 3
     );
     if (!result) throw new RuntimeError("svd: complex SVD failed");
-    if (nargout === 1) return RTV.tensor(new FloatXArray(result.S), [k, 1]);
+    if (nargout <= 1) return RTV.tensor(new FloatXArray(result.S), [k, 1]);
     const uCols = econ ? k : m;
     const vCols = econ ? k : n;
     return [
@@ -39619,7 +39913,7 @@ function svdApply(args, nargout) {
   if (bridge?.svd) {
     const result = bridge.svd(toF64(A.data), m, n, econ, nargout === 3);
     if (result) {
-      if (nargout === 1) return RTV.tensor(new FloatXArray(result.S), [k, 1]);
+      if (nargout <= 1) return RTV.tensor(new FloatXArray(result.S), [k, 1]);
       const uCols = econ ? k : m;
       const vCols = econ ? k : n;
       return [
@@ -39689,12 +39983,12 @@ function powerIterationEigenvalues(A_data, n, numEigenvalues) {
 registerIBuiltin({
   name: "qr",
   resolve: (argTypes, nargout) => {
-    if (nargout < 1 || nargout > 3 || argTypes.length < 1 || argTypes.length > 2)
+    if (nargout < 0 || nargout > 3 || argTypes.length < 1 || argTypes.length > 2)
       return null;
     if (!isNumericJitType(argTypes[0])) return null;
     const isComplex2 = argTypes[0].kind === "tensor" && argTypes[0].isComplex;
     const t = tensorType(isComplex2 || void 0);
-    if (nargout === 1)
+    if (nargout <= 1)
       return { outputTypes: [t], apply: (args, n) => qrApply(args, n) };
     if (nargout === 3)
       return {
@@ -39711,7 +40005,7 @@ function qrApply(args, nargout) {
   if (isRuntimeNumber(A)) {
     const val = A;
     const s = val >= 0 ? 1 : -1;
-    if (nargout === 1) return RTV.tensor(new FloatXArray([s * val]), [1, 1]);
+    if (nargout <= 1) return RTV.tensor(new FloatXArray([s * val]), [1, 1]);
     if (nargout === 3)
       return [
         RTV.tensor(new FloatXArray([s]), [1, 1]),
@@ -39744,7 +40038,7 @@ function qrApply(args, nargout) {
       nargout === 2
     );
     if (!result) throw new RuntimeError("qr: complex QR failed");
-    if (nargout === 1) {
+    if (nargout <= 1) {
       const rRows = econ ? k : m;
       return RTV.tensor(
         new FloatXArray(result.RRe),
@@ -39770,7 +40064,7 @@ function qrApply(args, nargout) {
   if (bridge?.qr) {
     const result = bridge.qr(toF64(A.data), m, n, econ, nargout === 2);
     if (result) {
-      if (nargout === 1)
+      if (nargout <= 1)
         return RTV.tensor(new FloatXArray(result.R), [econ ? k : m, n]);
       const qCols = econ ? k : m;
       return [
@@ -39813,7 +40107,7 @@ function qrApply(args, nargout) {
         R_data[colMajorIndex(j + i, c, m)] -= scale * v[i];
     }
   }
-  if (nargout === 1) {
+  if (nargout <= 1) {
     if (econ) {
       const R_econ = new FloatXArray(k * n);
       for (let r = 0; r < k; r++)
@@ -40029,11 +40323,11 @@ function qrPivotApply(A, m, n, k, econ) {
 registerIBuiltin({
   name: "lu",
   resolve: (argTypes, nargout) => {
-    if (nargout < 1 || nargout > 3) return null;
+    if (nargout < 0 || nargout > 3) return null;
     if (argTypes.length < 1 || argTypes.length > 2) return null;
     if (!isNumericJitType(argTypes[0])) return null;
     const t = tensorType();
-    if (nargout === 1)
+    if (nargout <= 1)
       return { outputTypes: [t], apply: (args, n) => luApply(args, n) };
     if (nargout === 2)
       return { outputTypes: [t, t], apply: (args, n) => luApply(args, n) };
@@ -40052,7 +40346,7 @@ function luApply(args, nargout) {
     if (nargout <= 2) {
       const L2 = RTV.tensor(new FloatXArray([1]), [1, 1]);
       const U2 = RTV.tensor(new FloatXArray([val]), [1, 1]);
-      if (nargout === 1) return L2;
+      if (nargout <= 1) return L2;
       return [L2, U2];
     }
     const L = RTV.tensor(new FloatXArray([1]), [1, 1]);
@@ -40090,7 +40384,7 @@ function luApply(args, nargout) {
       if (U_im && LU_im) U_im[i + j * k] = LU_im[i + j * m];
     }
   }
-  if (nargout === 1)
+  if (nargout <= 1)
     return RTV.tensor(
       new FloatXArray(LU_re),
       [m, n],
@@ -40166,11 +40460,11 @@ function ipivToPermVector(ipiv, m) {
 registerIBuiltin({
   name: "eig",
   resolve: (argTypes, nargout) => {
-    if (nargout < 1 || nargout > 3 || argTypes.length < 1 || argTypes.length > 3)
+    if (nargout < 0 || nargout > 3 || argTypes.length < 1 || argTypes.length > 3)
       return null;
     if (!isNumericJitType(argTypes[0])) return null;
     const c = tensorType(true);
-    if (nargout === 1)
+    if (nargout <= 1)
       return { outputTypes: [c], apply: (args, n) => eigApply(args, n) };
     if (nargout === 2)
       return { outputTypes: [c, c], apply: (args, n) => eigApply(args, n) };
@@ -40184,7 +40478,7 @@ function eigApply(args, nargout) {
   const { balance, outputForm } = parseEigOptionsRuntime(args);
   if (isRuntimeNumber(A)) {
     const val = A;
-    if (nargout === 1) return RTV.num(val);
+    if (nargout <= 1) return RTV.num(val);
     const V = RTV.tensor(new FloatXArray([1]), [1, 1]);
     if (nargout === 2) {
       if (outputForm === "vector") return [V, RTV.num(val)];
@@ -40213,7 +40507,7 @@ function eigApply(args, nargout) {
     );
     if (!result2) throw new RuntimeError("eig: complex eig failed");
     const { wRe, wIm, VLRe, VLIm, VRRe, VRIm } = result2;
-    if (nargout === 1) return maybeComplexTensor(wRe, [n, 1], wIm);
+    if (nargout <= 1) return maybeComplexTensor(wRe, [n, 1], wIm);
     const Vout2 = computeVR && VRRe && VRIm ? maybeComplexTensor(VRRe, [n, n], VRIm) : RTV.tensor(new FloatXArray(n * n), [n, n]);
     const Dout2 = outputForm === "vector" ? maybeComplexTensor(wRe, [n, 1], wIm) : buildDiagMatrix(wRe, wIm, n);
     if (nargout === 2) return [Vout2, Dout2];
@@ -40226,7 +40520,7 @@ function eigApply(args, nargout) {
   if (!result) throw new RuntimeError("eig: LAPACK eig failed");
   const { wr, wi, VL, VR } = result;
   const hasComplex = wi.some((v) => v !== 0);
-  if (nargout === 1) return maybeComplexTensor(wr, [n, 1], wi);
+  if (nargout <= 1) return maybeComplexTensor(wr, [n, 1], wi);
   const Vout = computeVR && VR ? buildEigenvectorMatrix(VR, wi, n, hasComplex) : RTV.tensor(new FloatXArray(n * n), [n, n]);
   const Dout = outputForm === "vector" ? maybeComplexTensor(wr, [n, 1], wi) : buildDiagMatrix(wr, wi, n);
   if (nargout === 2) return [Vout, Dout];
@@ -40250,11 +40544,11 @@ function parseEigOptionsRuntime(args) {
 registerIBuiltin({
   name: "chol",
   resolve: (argTypes, nargout) => {
-    if (nargout < 1 || nargout > 3 || argTypes.length < 1 || argTypes.length > 3)
+    if (nargout < 0 || nargout > 3 || argTypes.length < 1 || argTypes.length > 3)
       return null;
     if (!isNumericJitType(argTypes[0])) return null;
     const t = tensorType();
-    if (nargout === 1)
+    if (nargout <= 1)
       return { outputTypes: [t], apply: (args, n) => cholApply(args, n) };
     if (nargout === 2)
       return { outputTypes: [t, NUM], apply: (args, n) => cholApply(args, n) };
@@ -40415,7 +40709,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (argTypes.length !== 2 || nargout !== 1) return null;
+        if (argTypes.length !== 2 || nargout > 1) return null;
         if (!isNumericJitType(argTypes[0]) || !isNumericJitType(argTypes[1]))
           return null;
         return [tensorType()];
@@ -40458,7 +40752,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length < 1 || argTypes.length > 2)
+        if (nargout > 1 || argTypes.length < 1 || argTypes.length > 2)
           return null;
         if (!isNumericJitType(argTypes[0])) return null;
         return [NUM];
@@ -40519,7 +40813,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length < 1 || argTypes.length > 2)
+        if (nargout > 1 || argTypes.length < 1 || argTypes.length > 2)
           return null;
         if (!isNumericJitType(argTypes[0])) return null;
         return [NUM];
@@ -40567,7 +40861,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length < 1 || argTypes.length > 2)
+        if (nargout > 1 || argTypes.length < 1 || argTypes.length > 2)
           return null;
         if (!isNumericJitType(argTypes[0])) return null;
         const a = argTypes[0];
@@ -40664,7 +40958,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length !== 2) return null;
+        if (nargout > 1 || argTypes.length !== 2) return null;
         if (!isNumericJitType(argTypes[0]) || !isNumericJitType(argTypes[1]))
           return null;
         return [tensorType()];
@@ -40700,7 +40994,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length === 0) return null;
+        if (nargout > 1 || argTypes.length === 0) return null;
         if (!argTypes.every(isNumericJitType)) return null;
         return [tensorType()];
       },
@@ -40743,7 +41037,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1) return null;
+        if (nargout > 1) return null;
         if (argTypes.length !== 2 && argTypes.length !== 4) return null;
         const xIdx = 0, yIdx = argTypes.length === 4 ? 2 : 1;
         if (!isNumericJitType(argTypes[xIdx]) || !isNumericJitType(argTypes[yIdx]))
@@ -40872,7 +41166,7 @@ defineBuiltin({
   cases: [
     {
       match: (argTypes, nargout) => {
-        if (nargout !== 1 || argTypes.length !== 1) return null;
+        if (nargout > 1 || argTypes.length !== 1) return null;
         if (argTypes[0].kind !== "tensor") return null;
         return [tensorType(argTypes[0].isComplex)];
       },
@@ -44734,9 +45028,11 @@ registerIBuiltin({
 // src/numbl-core/helpers/prng.ts
 var _rngState = null;
 var _rngSeed = 0;
+var _bmSpare = null;
 function setRngShuffle() {
   _rngState = null;
   _rngSeed = 0;
+  _bmSpare = null;
 }
 function setRngSeed(seed) {
   _rngSeed = seed;
@@ -44752,6 +45048,7 @@ function splitmix32(seed) {
   };
 }
 function seedRng(seed) {
+  _bmSpare = null;
   const sm = splitmix32(seed);
   _rngState = new Uint32Array([sm(), sm(), sm(), sm()]);
   if (_rngState[0] === 0 && _rngState[1] === 0 && _rngState[2] === 0 && _rngState[3] === 0) {
@@ -44778,10 +45075,60 @@ function rngRandom() {
   return (xoshiro128ss() >>> 0) / 4294967296;
 }
 function boxMullerRandom() {
-  let u = 0, v = 0;
-  while (u === 0) u = rngRandom();
-  while (v === 0) v = rngRandom();
-  return Math.sqrt(-2 * Math.log(u)) * Math.cos(2 * Math.PI * v);
+  if (_bmSpare !== null) {
+    const s2 = _bmSpare;
+    _bmSpare = null;
+    return s2;
+  }
+  let u, v, s;
+  do {
+    u = 2 * rngRandom() - 1;
+    v = 2 * rngRandom() - 1;
+    s = u * u + v * v;
+  } while (s >= 1 || s === 0);
+  const mul = Math.sqrt(-2 * Math.log(s) / s);
+  _bmSpare = v * mul;
+  return u * mul;
+}
+function fillRandn(data) {
+  if (_rngState !== null) {
+    const bridge = getLapackBridge();
+    if (bridge?.fillRandn) {
+      const r = bridge.fillRandn(
+        _rngState,
+        data.length,
+        _bmSpare ?? 0,
+        _bmSpare !== null
+      );
+      _bmSpare = r.hasSpare ? r.spare : null;
+      if (data instanceof Float64Array) {
+        data.set(r.data);
+      } else {
+        for (let i2 = 0; i2 < data.length; i2++) data[i2] = r.data[i2];
+      }
+      return;
+    }
+  }
+  const n = data.length;
+  let i = 0;
+  if (_bmSpare !== null && i < n) {
+    data[i++] = _bmSpare;
+    _bmSpare = null;
+  }
+  for (; i + 1 < n; i += 2) {
+    let u, v, s;
+    do {
+      u = 2 * rngRandom() - 1;
+      v = 2 * rngRandom() - 1;
+      s = u * u + v * v;
+    } while (s >= 1 || s === 0);
+    const mul = Math.sqrt(-2 * Math.log(s) / s);
+    data[i] = u * mul;
+    data[i + 1] = v * mul;
+  }
+  if (i < n) {
+    data[i] = boxMullerRandom();
+  }
 }
 function getRngStateStruct() {
   const stateArray = _rngState ? RTV.tensor(new FloatXArray(Array.from(_rngState).map((v) => v)), [4, 1]) : RTV.tensor(new FloatXArray(0), [0, 1]);
@@ -44849,7 +45196,7 @@ registerIBuiltin({
     };
   }
 });
-function registerRandBuiltin(name, gen) {
+function registerRandBuiltin(name, gen, bulkFill) {
   defineBuiltin({
     name,
     cases: [
@@ -44888,7 +45235,11 @@ function registerRandBuiltin(name, gen) {
           if (shape.length === 1) shape.push(shape[0]);
           const n = numel(shape);
           const data = new FloatXArray(n);
-          for (let i = 0; i < n; i++) data[i] = gen();
+          if (bulkFill) {
+            bulkFill(data);
+          } else {
+            for (let i = 0; i < n; i++) data[i] = gen();
+          }
           return RTV.tensor(data, shape);
         }
       }
@@ -44896,7 +45247,7 @@ function registerRandBuiltin(name, gen) {
   });
 }
 registerRandBuiltin("rand", rngRandom);
-registerRandBuiltin("randn", boxMullerRandom);
+registerRandBuiltin("randn", boxMullerRandom, fillRandn);
 defineBuiltin({
   name: "randi",
   cases: [
@@ -45341,136 +45692,6 @@ registerIBuiltin({
   }
 });
 
-// src/numbl-core/interpreter/builtins/time-system.ts
-var ticTime = 0;
-defineBuiltin({
-  name: "tic",
-  cases: [
-    {
-      match: (argTypes) => argTypes.length === 0 ? [{ kind: "number" }] : null,
-      apply: () => {
-        ticTime = performance.now();
-        return RTV.num(ticTime / 1e3);
-      }
-    }
-  ]
-});
-defineBuiltin({
-  name: "toc",
-  cases: [
-    {
-      match: (argTypes) => argTypes.length === 0 ? [{ kind: "number" }] : null,
-      apply: () => RTV.num((performance.now() - ticTime) / 1e3)
-    }
-  ]
-});
-defineBuiltin({
-  name: "clock",
-  cases: [
-    {
-      match: (argTypes) => argTypes.length === 0 ? [{ kind: "tensor", isComplex: false, shape: [1, 6] }] : null,
-      apply: () => {
-        const now = /* @__PURE__ */ new Date();
-        return RTV.tensor(
-          new FloatXArray([
-            now.getFullYear(),
-            now.getMonth() + 1,
-            now.getDate(),
-            now.getHours(),
-            now.getMinutes(),
-            now.getSeconds() + now.getMilliseconds() / 1e3
-          ]),
-          [1, 6]
-        );
-      }
-    }
-  ]
-});
-defineBuiltin({
-  name: "etime",
-  cases: [
-    {
-      match: (argTypes) => argTypes.length === 2 ? [{ kind: "number" }] : null,
-      apply: (args) => {
-        const t1 = args[0];
-        const t0 = args[1];
-        if (!isRuntimeTensor(t1) || !isRuntimeTensor(t0))
-          throw new RuntimeError("etime: arguments must be clock vectors");
-        const toMs = (t) => {
-          const d = new Date(
-            t.data[0],
-            t.data[1] - 1,
-            t.data[2],
-            t.data[3],
-            t.data[4],
-            Math.floor(t.data[5]),
-            t.data[5] % 1 * 1e3
-          );
-          return d.getTime();
-        };
-        return RTV.num((toMs(t1) - toMs(t0)) / 1e3);
-      }
-    }
-  ]
-});
-defineBuiltin({
-  name: "version",
-  cases: [
-    {
-      match: (argTypes) => argTypes.length === 0 ? [{ kind: "char" }] : null,
-      apply: () => RTV.char("9.14.0")
-    }
-  ]
-});
-function getComputerStrings() {
-  if (typeof process === "undefined") {
-    return { str: "BROWSER", arch: "browser" };
-  }
-  const platform = process.platform;
-  const cpuArch = process.arch;
-  if (platform === "win32") return { str: "PCWIN64", arch: "win64" };
-  if (platform === "darwin") {
-    if (cpuArch === "arm64") return { str: "MACA64", arch: "maca64" };
-    return { str: "MACI64", arch: "maci64" };
-  }
-  return { str: "GLNXA64", arch: "glnxa64" };
-}
-registerIBuiltin({
-  name: "computer",
-  resolve: () => ({
-    outputTypes: [{ kind: "unknown" }],
-    apply: (args, nargout) => {
-      const info = getComputerStrings();
-      if (args.length >= 1 && isRuntimeChar(args[0])) {
-        const arg = toString(args[0]);
-        if (arg === "arch") return RTV.char(info.arch);
-        throw new RuntimeError(`computer: unknown argument '${arg}'`);
-      }
-      const maxsize = 2 ** 48 - 1;
-      const endian = RTV.char("L");
-      if (nargout <= 1) return RTV.char(info.str);
-      if (nargout === 2) return [RTV.char(info.str), RTV.num(maxsize)];
-      return [RTV.char(info.str), RTV.num(maxsize), endian];
-    }
-  })
-});
-var _platform = typeof process !== "undefined" ? process.platform : "linux";
-for (const [name, val] of [
-  ["ismac", _platform === "darwin"],
-  ["ispc", _platform === "win32"],
-  ["isunix", _platform !== "win32"]
-]) {
-  defineBuiltin({
-    name,
-    cases: [
-      {
-        match: (argTypes) => argTypes.length === 0 ? [{ kind: "boolean" }] : null,
-        apply: () => RTV.logical(val)
-      }
-    ]
-  });
-}
-
 // src/numbl-core/interpreter/builtins/sparse.ts
 function toNumericArray(v, name) {
   if (isRuntimeNumber(v)) return [v];
@@ -48330,7 +48551,7 @@ function getSourceLine(getSource, file, line) {
 }
 
 // src/numbl-core/version.ts
-var NUMBL_VERSION = "0.1.2";
+var NUMBL_VERSION = "0.1.3";
 
 // src/cli-repl.ts
 import { createInterface } from "readline";
@@ -49592,9 +49813,7 @@ function lowerStmt(ctx, stmt) {
       result = lowerAssign(ctx, stmt);
       break;
     case "ExprStmt":
-      if (!stmt.suppressed) return null;
-      result = lowerExprStmt(ctx, stmt);
-      break;
+      return null;
     case "If":
       result = lowerIf(ctx, stmt);
       break;
@@ -49669,11 +49888,6 @@ function lowerMultiAssign(ctx, stmt) {
     }
   ];
 }
-function lowerExprStmt(ctx, stmt) {
-  const expr = lowerExpr(ctx, stmt.expr);
-  if (!expr) return null;
-  return [{ tag: "ExprStmt", expr }];
-}
 function lowerIf(ctx, stmt) {
   const cond = lowerExpr(ctx, stmt.cond);
   if (!cond) return null;
@@ -50742,10 +50956,11 @@ function execStmt(stmt) {
   }
   switch (stmt.type) {
     case "ExprStmt": {
-      const val = this.evalExpr(stmt.expr);
+      const val = this.evalExprNargout(stmt.expr, 0);
       const singleVal = Array.isArray(val) ? val[0] : val;
       const rv = ensureRuntimeValue(singleVal);
       this.ans = rv;
+      this.env.set("ans", rv);
       if (!stmt.suppressed && !this.isOutputExpr(stmt.expr)) {
         this.rt.displayResult(rv);
       }
@@ -50795,10 +51010,12 @@ function execStmt(stmt) {
         const rv = this.rt.share(i < values.length ? values[i] : void 0);
         this.assignLValue(lv, rv);
       }
-      if (!stmt.suppressed && values.length > 0) {
-        const firstLv = stmt.lvalues[0];
-        if (firstLv.type === "Var") {
-          this.rt.displayAssign(firstLv.name, ensureRuntimeValue(values[0]));
+      if (!stmt.suppressed) {
+        for (let i = 0; i < stmt.lvalues.length; i++) {
+          const lv = stmt.lvalues[i];
+          if (lv.type === "Var" && i < values.length) {
+            this.rt.displayAssign(lv.name, ensureRuntimeValue(values[i]));
+          }
         }
       }
       return null;
@@ -51155,6 +51372,18 @@ function evalArgs(argExprs) {
   }
   return args;
 }
+var binopProfileName = {
+  ["Add" /* Add */]: "plus",
+  ["Sub" /* Sub */]: "minus",
+  ["Mul" /* Mul */]: "mtimes",
+  ["ElemMul" /* ElemMul */]: "times",
+  ["Div" /* Div */]: "mrdivide",
+  ["ElemDiv" /* ElemDiv */]: "rdivide",
+  ["LeftDiv" /* LeftDiv */]: "mldivide",
+  ["ElemLeftDiv" /* ElemLeftDiv */]: "ldivide",
+  ["Pow" /* Pow */]: "mpower",
+  ["ElemPow" /* ElemPow */]: "power"
+};
 function evalBinary(expr) {
   if (expr.op === "AndAnd" /* AndAnd */) {
     const left2 = this.evalExpr(expr.left);
@@ -51180,6 +51409,13 @@ function evalBinary(expr) {
     if (!isNaN(r)) return r;
     return mPow(ensureRuntimeValue(left), ensureRuntimeValue(right));
   }
+  if (this.rt.profilingEnabled && (typeof left !== "number" || typeof right !== "number")) {
+    const opName = binopProfileName[expr.op] ?? expr.op;
+    this.rt.profileEnter("builtin:interp:" + opName);
+    const result = binop(expr.op, left, right);
+    this.rt.profileLeave();
+    return result;
+  }
   return binop(expr.op, left, right);
 }
 function evalUnary(expr) {
@@ -51578,7 +51814,6 @@ function switchMatch2(switchVal, caseVal) {
   return switchMatch(switchVal, caseVal);
 }
 function isOutputExpr(expr) {
-  if (expr.type !== "FuncCall") return false;
   const outputFunctions = [
     "disp",
     "display",
@@ -51586,9 +51821,12 @@ function isOutputExpr(expr) {
     "warning",
     "assert",
     "tic",
+    "toc",
     "help"
   ];
-  return outputFunctions.includes(expr.name);
+  if (expr.type === "FuncCall") return outputFunctions.includes(expr.name);
+  if (expr.type === "Ident") return outputFunctions.includes(expr.name);
+  return false;
 }
 
 // src/numbl-core/interpreter/interpreterFunctions.ts
@@ -52256,8 +52494,9 @@ function callUserFunction(fn, args, nargout, narginOverride) {
     }
     const hasVarargout = fn.outputs.length > 0 && fn.outputs[fn.outputs.length - 1] === "varargout";
     const regularOutputs = hasVarargout ? fn.outputs.slice(0, -1) : fn.outputs;
+    const collectCount = nargout === 0 && regularOutputs.length > 0 ? 1 : Math.min(regularOutputs.length, nargout);
     const outputs = [];
-    for (let i = 0; i < Math.min(regularOutputs.length, nargout); i++) {
+    for (let i = 0; i < collectCount; i++) {
       const val = this.env.get(regularOutputs[i]);
       if (val === void 0 && nargout >= i + 1) {
         throw new RuntimeError(
@@ -53827,10 +54066,17 @@ function executeCode(source, options = {}, workspaceFiles, mainFileName = "scrip
   if (options.log) {
     options.log("AST parsed, starting interpretation");
   }
+  const isRepl = mainFileName === "repl";
   const localFunctions = [];
   const localClasses = [];
   for (const stmt of ast.body) {
     if (stmt.type === "Function") {
+      if (isRepl) {
+        throw new RuntimeError(
+          "Function definitions are not supported in the REPL. Save the function to a .m file instead.",
+          stmt.span
+        );
+      }
       localFunctions.push(stmt);
     } else if (stmt.type === "ClassDef") {
       localClasses.push(stmt);
@@ -54056,7 +54302,7 @@ ${jsCode}`
       output: rt.outputLines,
       generatedJS: jitSections.length > 0 ? `// Interpreter mode \u2014 JIT compiled sections:
 
-${jitSections.join("\n\n")}` : "// interpreted mode \u2014 no JS generated",
+${jitSections.join("\n\n")}` : "// No JS generated",
       plotInstructions: rt.plotInstructions,
       returnValue: interpreter.ans ?? RTV.num(0),
       variableValues: interpreter.getVariableValues(),
@@ -54082,7 +54328,7 @@ ${jitSections.join("\n\n")}` : "// interpreted mode \u2014 no JS generated",
   } catch (e) {
     const generatedJS = jitSections.length > 0 ? `// Interpreter mode \u2014 JIT compiled sections:
 
-${jitSections.join("\n\n")}` : "// interpreted mode \u2014 no JS generated";
+${jitSections.join("\n\n")}` : "// No JS generated";
     if (e instanceof RuntimeError) {
       if (e.line === null && rt.$file && rt.$line > 0) {
         e.file = rt.$file;
@@ -54132,7 +54378,7 @@ import {
   rmSync,
   rmdirSync
 } from "fs";
-import { inflateRawSync } from "zlib";
+import { unzipSync } from "fflate";
 import { execFileSync } from "child_process";
 import { homedir, tmpdir } from "os";
 import { join as join3, resolve, dirname, basename } from "path";
@@ -54208,6 +54454,21 @@ function permissionToFlags(permission) {
   }
 }
 var READ_CHUNK_SIZE = 8192;
+function buildFetchOptions(options) {
+  const init = {};
+  if (options?.requestMethod) {
+    init.method = options.requestMethod.toUpperCase();
+  }
+  const headers = {};
+  if (options?.username && options?.password) {
+    headers["Authorization"] = "Basic " + Buffer.from(`${options.username}:${options.password}`).toString("base64");
+  }
+  if (options?.keyName && options?.keyValue) {
+    headers[options.keyName] = options.keyValue;
+  }
+  if (Object.keys(headers).length > 0) init.headers = headers;
+  return init;
+}
 var NodeFileIOAdapter = class {
   nextFid = 3;
   // 0=stdin, 1=stdout, 2=stderr reserved
@@ -54381,35 +54642,43 @@ var NodeFileIOAdapter = class {
   tempdir() {
     return tmpdir();
   }
-  websave(url, filename) {
+  websave(url, filename, options) {
     const dest = expandTilde(filename);
+    const timeoutMs = options?.timeout ? Math.round(options.timeout * 1e3) : 3e4;
+    const fetchOpts = buildFetchOptions(options);
     const script = `
-      fetch(${JSON.stringify(url)})
-        .then(r => { if (!r.ok) throw new Error('HTTP ' + r.status); return r.arrayBuffer(); })
+      const ac = new AbortController();
+      const t = setTimeout(() => ac.abort(), ${timeoutMs});
+      fetch(${JSON.stringify(url)}, { ...${JSON.stringify(fetchOpts)}, signal: ac.signal })
+        .then(r => { clearTimeout(t); if (!r.ok) throw new Error('HTTP ' + r.status); return r.arrayBuffer(); })
         .then(buf => require('fs').writeFileSync(${JSON.stringify(dest)}, Buffer.from(buf)))
         .catch(e => { process.stderr.write(e.message + '\\n'); process.exit(1); });
     `;
     try {
       execFileSync(process.execPath, ["-e", script], {
         stdio: ["ignore", "ignore", "pipe"],
-        timeout: 3e4
+        timeout: timeoutMs + 5e3
       });
     } catch (e) {
       const msg = e instanceof Error && "stderr" in e ? e.stderr.toString().trim() : String(e);
       throw new Error(`websave: failed to download ${url}: ${msg}`);
     }
   }
-  webread(url) {
+  webread(url, options) {
+    const timeoutMs = options?.timeout ? Math.round(options.timeout * 1e3) : 3e4;
+    const fetchOpts = buildFetchOptions(options);
     const script = `
-      fetch(${JSON.stringify(url)})
-        .then(r => { if (!r.ok) throw new Error('HTTP ' + r.status); return r.text(); })
+      const ac = new AbortController();
+      const t = setTimeout(() => ac.abort(), ${timeoutMs});
+      fetch(${JSON.stringify(url)}, { ...${JSON.stringify(fetchOpts)}, signal: ac.signal })
+        .then(r => { clearTimeout(t); if (!r.ok) throw new Error('HTTP ' + r.status); return r.text(); })
         .then(t => process.stdout.write(t))
         .catch(e => { process.stderr.write(e.message + '\\n'); process.exit(1); });
     `;
     try {
       const result = execFileSync(process.execPath, ["-e", script], {
         stdio: ["ignore", "pipe", "pipe"],
-        timeout: 3e4,
+        timeout: timeoutMs + 5e3,
         maxBuffer: 50 * 1024 * 1024
       });
       return result.toString();
@@ -54466,56 +54735,18 @@ var NodeFileIOAdapter = class {
     const dest = resolve(expandTilde(outputfolder));
     mkdirSync(dest, { recursive: true });
     const buf = readFileSync3(src);
-    let eocdOffset = -1;
-    for (let i = buf.length - 22; i >= 0; i--) {
-      if (buf[i] === 80 && buf[i + 1] === 75 && buf[i + 2] === 5 && buf[i + 3] === 6) {
-        eocdOffset = i;
-        break;
-      }
-    }
-    if (eocdOffset === -1) throw new Error("unzip: invalid ZIP file");
-    const cdOffset = buf.readUInt32LE(eocdOffset + 16);
-    const cdEntries = buf.readUInt16LE(eocdOffset + 10);
+    const files = unzipSync(
+      new Uint8Array(buf.buffer, buf.byteOffset, buf.byteLength)
+    );
     const extracted = [];
-    let pos = cdOffset;
-    for (let i = 0; i < cdEntries; i++) {
-      if (buf.readUInt32LE(pos) !== 33639248)
-        throw new Error("unzip: corrupt central directory");
-      const method = buf.readUInt16LE(pos + 10);
-      const nameLen = buf.readUInt16LE(pos + 28);
-      const extraLen = buf.readUInt16LE(pos + 30);
-      const commentLen = buf.readUInt16LE(pos + 32);
-      const localHeaderOffset = buf.readUInt32LE(pos + 42);
-      const entryName = buf.toString("utf-8", pos + 46, pos + 46 + nameLen);
-      pos += 46 + nameLen + extraLen + commentLen;
+    for (const [entryName, data] of Object.entries(files)) {
       if (entryName.endsWith("/")) {
         mkdirSync(join3(dest, entryName), { recursive: true });
         continue;
       }
-      const localPos = localHeaderOffset;
-      if (buf.readUInt32LE(localPos) !== 67324752)
-        throw new Error("unzip: corrupt local header");
-      const localNameLen = buf.readUInt16LE(localPos + 26);
-      const localExtraLen = buf.readUInt16LE(localPos + 28);
-      const compressedSize = buf.readUInt32LE(localPos + 18);
-      const dataStart = localPos + 30 + localNameLen + localExtraLen;
-      const compressedData = buf.subarray(
-        dataStart,
-        dataStart + compressedSize
-      );
-      let fileData;
-      if (method === 0) {
-        fileData = Buffer.from(compressedData);
-      } else if (method === 8) {
-        fileData = inflateRawSync(compressedData);
-      } else {
-        throw new Error(
-          `unzip: unsupported compression method ${method} for ${entryName}`
-        );
-      }
       const outPath = join3(dest, entryName);
       mkdirSync(dirname(outPath), { recursive: true });
-      writeFileSync(outPath, fileData);
+      writeFileSync(outPath, data);
       extracted.push(outPath);
     }
     return extracted;
@@ -55241,65 +55472,15 @@ async function runRepl(initialWorkspaceFiles, onDrawnow, initialSearchPaths, nat
 }
 
 // src/vfs/unzipToFiles.ts
-import { inflateSync } from "fflate";
-var TEXT_DECODER = new TextDecoder("utf-8");
+import { unzipSync as unzipSync2 } from "fflate";
 function unzipToFiles(zipData) {
-  let eocdOffset = -1;
-  for (let i = zipData.length - 22; i >= 0; i--) {
-    if (zipData[i] === 80 && zipData[i + 1] === 75 && zipData[i + 2] === 5 && zipData[i + 3] === 6) {
-      eocdOffset = i;
-      break;
-    }
-  }
-  if (eocdOffset === -1) throw new Error("Invalid ZIP file");
-  const view = new DataView(
-    zipData.buffer,
-    zipData.byteOffset,
-    zipData.byteLength
-  );
-  const cdOffset = view.getUint32(eocdOffset + 16, true);
-  const cdEntries = view.getUint16(eocdOffset + 10, true);
-  const files = [];
-  let pos = cdOffset;
-  for (let i = 0; i < cdEntries; i++) {
-    if (view.getUint32(pos, true) !== 33639248) {
-      throw new Error("Corrupt ZIP central directory");
-    }
-    const method = view.getUint16(pos + 10, true);
-    const nameLen = view.getUint16(pos + 28, true);
-    const extraLen = view.getUint16(pos + 30, true);
-    const commentLen = view.getUint16(pos + 32, true);
-    const localHeaderOffset = view.getUint32(pos + 42, true);
-    const entryName = TEXT_DECODER.decode(
-      zipData.subarray(pos + 46, pos + 46 + nameLen)
-    );
-    pos += 46 + nameLen + extraLen + commentLen;
-    if (entryName.endsWith("/")) continue;
-    const localPos = localHeaderOffset;
-    if (view.getUint32(localPos, true) !== 67324752) {
-      throw new Error("Corrupt ZIP local header");
-    }
-    const localNameLen = view.getUint16(localPos + 26, true);
-    const localExtraLen = view.getUint16(localPos + 28, true);
-    const compressedSize = view.getUint32(localPos + 18, true);
-    const dataStart = localPos + 30 + localNameLen + localExtraLen;
-    const compressedData = zipData.subarray(
-      dataStart,
-      dataStart + compressedSize
-    );
-    let content;
-    if (method === 0) {
-      content = new Uint8Array(compressedData);
-    } else if (method === 8) {
-      content = inflateSync(compressedData);
-    } else {
-      throw new Error(
-        `Unsupported ZIP compression method ${method} for ${entryName}`
-      );
-    }
-    files.push({ path: entryName, content });
+  const files = unzipSync2(zipData);
+  const out = [];
+  for (const [path, content] of Object.entries(files)) {
+    if (path.endsWith("/")) continue;
+    out.push({ path, content });
   }
-  return files;
+  return out;
 }
 
 // src/cli.ts