@chronon-sim/engine 0.1.0

package/dist/rng/distributions.d.ts

@@ -0,0 +1,98 @@
+import type { Rng } from "./rng.js";
+/** Anything that can produce a number on demand: constants and random distributions alike. */
+export interface SampleProvider {
+    nextSample(): number;
+}
+/** JSON-serialisable description of a distribution. Mirrors the JSON model format. */
+export type DistributionSpec = {
+    dist: "constant";
+    value: number;
+} | {
+    dist: "uniform";
+    min: number;
+    max: number;
+} | {
+    dist: "exponential";
+    mean: number;
+} | {
+    dist: "normal";
+    mean: number;
+    stdDev: number;
+} | {
+    dist: "triangular";
+    min: number;
+    mode: number;
+    max: number;
+} | {
+    dist: "lognormal";
+    mean: number;
+    stdDev: number;
+};
+/** A plain number (constant) or a distribution description. */
+export type SamplerSpec = number | DistributionSpec;
+/** Names of the supported distributions. */
+export declare const DISTRIBUTION_NAMES: readonly ["constant", "uniform", "exponential", "normal", "triangular", "lognormal"];
+/** Required numeric parameters of each distribution. */
+export declare const DISTRIBUTION_PARAMS: Record<(typeof DISTRIBUTION_NAMES)[number], readonly string[]>;
+/** Always returns the same value. */
+export declare class ConstantSampler implements SampleProvider {
+    readonly value: number;
+    constructor(value: number);
+    nextSample(): number;
+}
+/** Uniform on [min, max). */
+export declare class UniformSampler implements SampleProvider {
+    private readonly rng;
+    private readonly min;
+    private readonly max;
+    constructor(rng: Rng, min: number, max: number);
+    nextSample(): number;
+}
+/** Exponential with the given mean (rate = 1 / mean). */
+export declare class ExponentialSampler implements SampleProvider {
+    private readonly rng;
+    private readonly mean;
+    constructor(rng: Rng, mean: number);
+    nextSample(): number;
+}
+/** Normal (Gaussian) via the Box–Muller transform. Can return negative values. */
+export declare class NormalSampler implements SampleProvider {
+    private readonly rng;
+    private readonly mean;
+    private readonly stdDev;
+    private spare;
+    constructor(rng: Rng, mean: number, stdDev: number);
+    /** A standard-normal variate. */
+    nextStandard(): number;
+    nextSample(): number;
+}
+/** Triangular on [min, max] with the given mode. */
+export declare class TriangularSampler implements SampleProvider {
+    private readonly rng;
+    private readonly min;
+    private readonly mode;
+    private readonly max;
+    constructor(rng: Rng, min: number, mode: number, max: number);
+    nextSample(): number;
+}
+/**
+ * Lognormal parameterised by the mean and standard deviation of the lognormal variable itself
+ * (not of its logarithm), which is how service times are usually described.
+ */
+export declare class LognormalSampler implements SampleProvider {
+    private readonly normal;
+    private readonly mu;
+    private readonly sigma;
+    constructor(rng: Rng, mean: number, stdDev: number);
+    nextSample(): number;
+}
+/**
+ * Validate a raw (untrusted, e.g. parsed JSON) sampler spec.
+ * @returns human-readable problems; empty when the spec is valid
+ */
+export declare function validateSamplerSpec(spec: unknown): string[];
+/** Theoretical mean of a validated spec (used for documentation and sanity tests). */
+export declare function samplerMean(spec: SamplerSpec): number;
+/** Build a {@link SampleProvider} from a validated spec, drawing from `rng`. */
+export declare function createSampler(spec: SamplerSpec, rng: Rng): SampleProvider;
+//# sourceMappingURL=distributions.d.ts.map

package/dist/rng/distributions.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"distributions.d.ts","sourceRoot":"","sources":["../../src/rng/distributions.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,GAAG,EAAE,MAAM,UAAU,CAAC;AAEpC,8FAA8F;AAC9F,MAAM,WAAW,cAAc;IAC7B,UAAU,IAAI,MAAM,CAAC;CACtB;AAED,sFAAsF;AACtF,MAAM,MAAM,gBAAgB,GACxB;IAAE,IAAI,EAAE,UAAU,CAAC;IAAC,KAAK,EAAE,MAAM,CAAA;CAAE,GACnC;IAAE,IAAI,EAAE,SAAS,CAAC;IAAC,GAAG,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,MAAM,CAAA;CAAE,GAC7C;IAAE,IAAI,EAAE,aAAa,CAAC;IAAC,IAAI,EAAE,MAAM,CAAA;CAAE,GACrC;IAAE,IAAI,EAAE,QAAQ,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,GAChD;IAAE,IAAI,EAAE,YAAY,CAAC;IAAC,GAAG,EAAE,MAAM,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,GAAG,EAAE,MAAM,CAAA;CAAE,GAC9D;IAAE,IAAI,EAAE,WAAW,CAAC;IAAC,IAAI,EAAE,MAAM,CAAC;IAAC,MAAM,EAAE,MAAM,CAAA;CAAE,CAAC;AAExD,+DAA+D;AAC/D,MAAM,MAAM,WAAW,GAAG,MAAM,GAAG,gBAAgB,CAAC;AAEpD,4CAA4C;AAC5C,eAAO,MAAM,kBAAkB,sFAOrB,CAAC;AAEX,wDAAwD;AACxD,eAAO,MAAM,mBAAmB,EAAE,MAAM,CAAC,CAAC,OAAO,kBAAkB,CAAC,CAAC,MAAM,CAAC,EAAE,SAAS,MAAM,EAAE,CAO9F,CAAC;AAEF,qCAAqC;AACrC,qBAAa,eAAgB,YAAW,cAAc;IACxC,QAAQ,CAAC,KAAK,EAAE,MAAM;gBAAb,KAAK,EAAE,MAAM;IAClC,UAAU,IAAI,MAAM;CAGrB;AAED,6BAA6B;AAC7B,qBAAa,cAAe,YAAW,cAAc;IAEjD,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,GAAG;gBAFH,GAAG,EAAE,GAAG,EACR,GAAG,EAAE,MAAM,EACX,GAAG,EAAE,MAAM;IAE9B,UAAU,IAAI,MAAM;CAGrB;AAED,yDAAyD;AACzD,qBAAa,kBAAmB,YAAW,cAAc;IAErD,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,IAAI;gBADJ,GAAG,EAAE,GAAG,EACR,IAAI,EAAE,MAAM;IAE/B,UAAU,IAAI,MAAM;CAGrB;AAED,kFAAkF;AAClF,qBAAa,aAAc,YAAW,cAAc;IAGhD,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,IAAI;IACrB,OAAO,CAAC,QAAQ,CAAC,MAAM;IAJzB,OAAO,CAAC,KAAK,CAAuB;gBAEjB,GAAG,EAAE,GAAG,EACR,IAAI,EAAE,MAAM,EACZ,MAAM,EAAE,MAAM;IAEjC,iCAAiC;IACjC,YAAY,IAAI,MAAM;IAWtB,UAAU,IAAI,MAAM;CAGrB;AAED,oDAAoD;AACpD,qBAAa,iBAAkB,YAAW,cAAc;IAEpD,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,GAAG;IACpB,OAAO,CAAC,QAAQ,CAAC,IAAI;IACrB,OAAO,CAAC,QAAQ,CAAC,GAAG;gBAHH,GAAG,EAAE,GAAG,EACR,GAAG,EAAE,MAAM,EACX,IAAI,EAAE,MAAM,EACZ,GAAG,EAAE,MAAM;IAE9B,UAAU,IAAI,MAAM;CASrB;AAED;;;GAGG;AACH,qBAAa,gBAAiB,YAAW,cAAc;IACrD,OAAO,CAAC,QAAQ,CAAC,MAAM,CAAgB;IACvC,OAAO,CAAC,QAAQ,CAAC,EAAE,CAAS;IAC5B,OAAO,CAAC,QAAQ,CAAC,KAAK,CAAS;gBACnB,GAAG,EAAE,GAAG,EAAE,IAAI,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM;IAUlD,UAAU,IAAI,MAAM;CAGrB;AAUD;;;GAGG;AACH,wBAAgB,mBAAmB,CAAC,IAAI,EAAE,OAAO,GAAG,MAAM,EAAE,CAyC3D;AAED,sFAAsF;AACtF,wBAAgB,WAAW,CAAC,IAAI,EAAE,WAAW,GAAG,MAAM,CAcrD;AAED,gFAAgF;AAChF,wBAAgB,aAAa,CAAC,IAAI,EAAE,WAAW,EAAE,GAAG,EAAE,GAAG,GAAG,cAAc,CAgBzE"}

package/dist/rng/distributions.js

@@ -0,0 +1,224 @@
+/** Names of the supported distributions. */
+export const DISTRIBUTION_NAMES = [
+    "constant",
+    "uniform",
+    "exponential",
+    "normal",
+    "triangular",
+    "lognormal",
+];
+/** Required numeric parameters of each distribution. */
+export const DISTRIBUTION_PARAMS = {
+    constant: ["value"],
+    uniform: ["min", "max"],
+    exponential: ["mean"],
+    normal: ["mean", "stdDev"],
+    triangular: ["min", "mode", "max"],
+    lognormal: ["mean", "stdDev"],
+};
+/** Always returns the same value. */
+export class ConstantSampler {
+    value;
+    constructor(value) {
+        this.value = value;
+    }
+    nextSample() {
+        return this.value;
+    }
+}
+/** Uniform on [min, max). */
+export class UniformSampler {
+    rng;
+    min;
+    max;
+    constructor(rng, min, max) {
+        this.rng = rng;
+        this.min = min;
+        this.max = max;
+    }
+    nextSample() {
+        return this.min + (this.max - this.min) * this.rng.nextFloat();
+    }
+}
+/** Exponential with the given mean (rate = 1 / mean). */
+export class ExponentialSampler {
+    rng;
+    mean;
+    constructor(rng, mean) {
+        this.rng = rng;
+        this.mean = mean;
+    }
+    nextSample() {
+        return -this.mean * Math.log(this.rng.nextFloatOpen());
+    }
+}
+/** Normal (Gaussian) via the Box–Muller transform. Can return negative values. */
+export class NormalSampler {
+    rng;
+    mean;
+    stdDev;
+    spare = null;
+    constructor(rng, mean, stdDev) {
+        this.rng = rng;
+        this.mean = mean;
+        this.stdDev = stdDev;
+    }
+    /** A standard-normal variate. */
+    nextStandard() {
+        if (this.spare !== null) {
+            const z = this.spare;
+            this.spare = null;
+            return z;
+        }
+        const r = Math.sqrt(-2 * Math.log(this.rng.nextFloatOpen()));
+        const theta = 2 * Math.PI * this.rng.nextFloat();
+        this.spare = r * Math.sin(theta);
+        return r * Math.cos(theta);
+    }
+    nextSample() {
+        return this.mean + this.stdDev * this.nextStandard();
+    }
+}
+/** Triangular on [min, max] with the given mode. */
+export class TriangularSampler {
+    rng;
+    min;
+    mode;
+    max;
+    constructor(rng, min, mode, max) {
+        this.rng = rng;
+        this.min = min;
+        this.mode = mode;
+        this.max = max;
+    }
+    nextSample() {
+        const { min, mode, max } = this;
+        if (max === min)
+            return min;
+        const u = this.rng.nextFloat();
+        const fc = (mode - min) / (max - min);
+        return u < fc
+            ? min + Math.sqrt(u * (max - min) * (mode - min))
+            : max - Math.sqrt((1 - u) * (max - min) * (max - mode));
+    }
+}
+/**
+ * Lognormal parameterised by the mean and standard deviation of the lognormal variable itself
+ * (not of its logarithm), which is how service times are usually described.
+ */
+export class LognormalSampler {
+    normal;
+    mu;
+    sigma;
+    constructor(rng, mean, stdDev) {
+        // sigma^2 = ln(1 + (stdDev/mean)^2). log1p keeps precision when the ratio is tiny. When the squared ratio
+        // overflows (a subnormal mean, or an enormous stdDev) 1 + ratio^2 is just ratio^2, so use 2 ln(stdDev/mean)
+        // computed from the two logs; the direct formula would give ln(Infinity) - Infinity = NaN.
+        const ratioSquared = (stdDev / mean) ** 2;
+        const sigmaSquared = Number.isFinite(ratioSquared) ? Math.log1p(ratioSquared) : 2 * (Math.log(stdDev) - Math.log(mean));
+        this.sigma = Math.sqrt(sigmaSquared);
+        this.mu = Math.log(mean) - 0.5 * sigmaSquared;
+        this.normal = new NormalSampler(rng, 0, 1);
+    }
+    nextSample() {
+        return Math.exp(this.mu + this.sigma * this.normal.nextStandard());
+    }
+}
+function isRecord(v) {
+    return typeof v === "object" && v !== null && !Array.isArray(v);
+}
+function finite(v) {
+    return typeof v === "number" && Number.isFinite(v);
+}
+/**
+ * Validate a raw (untrusted, e.g. parsed JSON) sampler spec.
+ * @returns human-readable problems; empty when the spec is valid
+ */
+export function validateSamplerSpec(spec) {
+    if (finite(spec))
+        return [];
+    if (!isRecord(spec)) {
+        return ["expected a number or a distribution object such as {\"dist\": \"exponential\", \"mean\": 1}"];
+    }
+    const name = spec["dist"];
+    if (typeof name !== "string" || !DISTRIBUTION_NAMES.includes(name)) {
+        return [`"dist" must be one of ${DISTRIBUTION_NAMES.join(", ")}`];
+    }
+    const params = DISTRIBUTION_PARAMS[name];
+    const problems = [];
+    for (const key of params) {
+        if (!finite(spec[key]))
+            problems.push(`${name}: "${key}" must be a finite number`);
+    }
+    for (const key of Object.keys(spec)) {
+        if (key !== "dist" && !params.includes(key))
+            problems.push(`${name}: unknown parameter "${key}"`);
+    }
+    if (problems.length > 0)
+        return problems;
+    const p = spec;
+    const get = (k) => p[k];
+    switch (name) {
+        case "uniform":
+            if (get("min") > get("max"))
+                problems.push("uniform: min must be <= max");
+            break;
+        case "exponential":
+            if (!(get("mean") > 0))
+                problems.push("exponential: mean must be > 0");
+            break;
+        case "normal":
+            if (get("stdDev") < 0)
+                problems.push("normal: stdDev must be >= 0");
+            break;
+        case "triangular":
+            if (!(get("min") <= get("mode") && get("mode") <= get("max"))) {
+                problems.push("triangular: requires min <= mode <= max");
+            }
+            break;
+        case "lognormal":
+            if (!(get("mean") > 0))
+                problems.push("lognormal: mean must be > 0");
+            if (get("stdDev") < 0)
+                problems.push("lognormal: stdDev must be >= 0");
+            break;
+    }
+    return problems;
+}
+/** Theoretical mean of a validated spec (used for documentation and sanity tests). */
+export function samplerMean(spec) {
+    if (typeof spec === "number")
+        return spec;
+    switch (spec.dist) {
+        case "constant":
+            return spec.value;
+        case "uniform":
+            return (spec.min + spec.max) / 2;
+        case "exponential":
+        case "normal":
+        case "lognormal":
+            return spec.mean;
+        case "triangular":
+            return (spec.min + spec.mode + spec.max) / 3;
+    }
+}
+/** Build a {@link SampleProvider} from a validated spec, drawing from `rng`. */
+export function createSampler(spec, rng) {
+    if (typeof spec === "number")
+        return new ConstantSampler(spec);
+    switch (spec.dist) {
+        case "constant":
+            return new ConstantSampler(spec.value);
+        case "uniform":
+            return new UniformSampler(rng, spec.min, spec.max);
+        case "exponential":
+            return new ExponentialSampler(rng, spec.mean);
+        case "normal":
+            return new NormalSampler(rng, spec.mean, spec.stdDev);
+        case "triangular":
+            return new TriangularSampler(rng, spec.min, spec.mode, spec.max);
+        case "lognormal":
+            return new LognormalSampler(rng, spec.mean, spec.stdDev);
+    }
+}
+//# sourceMappingURL=distributions.js.map

package/dist/rng/distributions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"distributions.js","sourceRoot":"","sources":["../../src/rng/distributions.ts"],"names":[],"mappings":"AAmBA,4CAA4C;AAC5C,MAAM,CAAC,MAAM,kBAAkB,GAAG;IAChC,UAAU;IACV,SAAS;IACT,aAAa;IACb,QAAQ;IACR,YAAY;IACZ,WAAW;CACH,CAAC;AAEX,wDAAwD;AACxD,MAAM,CAAC,MAAM,mBAAmB,GAAmE;IACjG,QAAQ,EAAE,CAAC,OAAO,CAAC;IACnB,OAAO,EAAE,CAAC,KAAK,EAAE,KAAK,CAAC;IACvB,WAAW,EAAE,CAAC,MAAM,CAAC;IACrB,MAAM,EAAE,CAAC,MAAM,EAAE,QAAQ,CAAC;IAC1B,UAAU,EAAE,CAAC,KAAK,EAAE,MAAM,EAAE,KAAK,CAAC;IAClC,SAAS,EAAE,CAAC,MAAM,EAAE,QAAQ,CAAC;CAC9B,CAAC;AAEF,qCAAqC;AACrC,MAAM,OAAO,eAAe;IACL;IAArB,YAAqB,KAAa;QAAb,UAAK,GAAL,KAAK,CAAQ;IAAG,CAAC;IACtC,UAAU;QACR,OAAO,IAAI,CAAC,KAAK,CAAC;IACpB,CAAC;CACF;AAED,6BAA6B;AAC7B,MAAM,OAAO,cAAc;IAEN;IACA;IACA;IAHnB,YACmB,GAAQ,EACR,GAAW,EACX,GAAW;QAFX,QAAG,GAAH,GAAG,CAAK;QACR,QAAG,GAAH,GAAG,CAAQ;QACX,QAAG,GAAH,GAAG,CAAQ;IAC3B,CAAC;IACJ,UAAU;QACR,OAAO,IAAI,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC;IACjE,CAAC;CACF;AAED,yDAAyD;AACzD,MAAM,OAAO,kBAAkB;IAEV;IACA;IAFnB,YACmB,GAAQ,EACR,IAAY;QADZ,QAAG,GAAH,GAAG,CAAK;QACR,SAAI,GAAJ,IAAI,CAAQ;IAC5B,CAAC;IACJ,UAAU;QACR,OAAO,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,CAAC,CAAC;IACzD,CAAC;CACF;AAED,kFAAkF;AAClF,MAAM,OAAO,aAAa;IAGL;IACA;IACA;IAJX,KAAK,GAAkB,IAAI,CAAC;IACpC,YACmB,GAAQ,EACR,IAAY,EACZ,MAAc;QAFd,QAAG,GAAH,GAAG,CAAK;QACR,SAAI,GAAJ,IAAI,CAAQ;QACZ,WAAM,GAAN,MAAM,CAAQ;IAC9B,CAAC;IACJ,iCAAiC;IACjC,YAAY;QACV,IAAI,IAAI,CAAC,KAAK,KAAK,IAAI,EAAE,CAAC;YACxB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC;YACrB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC;YAClB,OAAO,CAAC,CAAC;QACX,CAAC;QACD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,aAAa,EAAE,CAAC,CAAC,CAAC;QAC7D,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC;QACjD,IAAI,CAAC,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QACjC,OAAO,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;IAC7B,CAAC;IACD,UAAU;QACR,OAAO,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,MAAM,GAAG,IAAI,CAAC,YAAY,EAAE,CAAC;IACvD,CAAC;CACF;AAED,oDAAoD;AACpD,MAAM,OAAO,iBAAiB;IAET;IACA;IACA;IACA;IAJnB,YACmB,GAAQ,EACR,GAAW,EACX,IAAY,EACZ,GAAW;QAHX,QAAG,GAAH,GAAG,CAAK;QACR,QAAG,GAAH,GAAG,CAAQ;QACX,SAAI,GAAJ,IAAI,CAAQ;QACZ,QAAG,GAAH,GAAG,CAAQ;IAC3B,CAAC;IACJ,UAAU;QACR,MAAM,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC;QAChC,IAAI,GAAG,KAAK,GAAG;YAAE,OAAO,GAAG,CAAC;QAC5B,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC;QAC/B,MAAM,EAAE,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;QACtC,OAAO,CAAC,GAAG,EAAE;YACX,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC;YACjD,CAAC,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC;IAC5D,CAAC;CACF;AAED;;;GAGG;AACH,MAAM,OAAO,gBAAgB;IACV,MAAM,CAAgB;IACtB,EAAE,CAAS;IACX,KAAK,CAAS;IAC/B,YAAY,GAAQ,EAAE,IAAY,EAAE,MAAc;QAChD,0GAA0G;QAC1G,4GAA4G;QAC5G,2FAA2F;QAC3F,MAAM,YAAY,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC;QAC1C,MAAM,YAAY,GAAG,MAAM,CAAC,QAAQ,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC;QACxH,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,YAAY,CAAC,CAAC;QACrC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,GAAG,GAAG,YAAY,CAAC;QAC9C,IAAI,CAAC,MAAM,GAAG,IAAI,aAAa,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,CAAC;IACD,UAAU;QACR,OAAO,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,YAAY,EAAE,CAAC,CAAC;IACrE,CAAC;CACF;AAED,SAAS,QAAQ,CAAC,CAAU;IAC1B,OAAO,OAAO,CAAC,KAAK,QAAQ,IAAI,CAAC,KAAK,IAAI,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;AAClE,CAAC;AAED,SAAS,MAAM,CAAC,CAAU;IACxB,OAAO,OAAO,CAAC,KAAK,QAAQ,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;AACrD,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,mBAAmB,CAAC,IAAa;IAC/C,IAAI,MAAM,CAAC,IAAI,CAAC;QAAE,OAAO,EAAE,CAAC;IAC5B,IAAI,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QACpB,OAAO,CAAC,6FAA6F,CAAC,CAAC;IACzG,CAAC;IACD,MAAM,IAAI,GAAG,IAAI,CAAC,MAAM,CAAC,CAAC;IAC1B,IAAI,OAAO,IAAI,KAAK,QAAQ,IAAI,CAAE,kBAAwC,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE,CAAC;QAC1F,OAAO,CAAC,yBAAyB,kBAAkB,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACpE,CAAC;IACD,MAAM,MAAM,GAAG,mBAAmB,CAAC,IAAwC,CAAC,CAAC;IAC7E,MAAM,QAAQ,GAAa,EAAE,CAAC;IAC9B,KAAK,MAAM,GAAG,IAAI,MAAM,EAAE,CAAC;QACzB,IAAI,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;YAAE,QAAQ,CAAC,IAAI,CAAC,GAAG,IAAI,MAAM,GAAG,2BAA2B,CAAC,CAAC;IACrF,CAAC;IACD,KAAK,MAAM,GAAG,IAAI,MAAM,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC;QACpC,IAAI,GAAG,KAAK,MAAM,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,GAAG,CAAC;YAAE,QAAQ,CAAC,IAAI,CAAC,GAAG,IAAI,wBAAwB,GAAG,GAAG,CAAC,CAAC;IACpG,CAAC;IACD,IAAI,QAAQ,CAAC,MAAM,GAAG,CAAC;QAAE,OAAO,QAAQ,CAAC;IACzC,MAAM,CAAC,GAAG,IAA8B,CAAC;IACzC,MAAM,GAAG,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,CAAW,CAAC;IAC1C,QAAQ,IAAI,EAAE,CAAC;QACb,KAAK,SAAS;YACZ,IAAI,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,KAAK,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,6BAA6B,CAAC,CAAC;YAC1E,MAAM;QACR,KAAK,aAAa;YAChB,IAAI,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,+BAA+B,CAAC,CAAC;YACvE,MAAM;QACR,KAAK,QAAQ;YACX,IAAI,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,6BAA6B,CAAC,CAAC;YACpE,MAAM;QACR,KAAK,YAAY;YACf,IAAI,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,KAAK,CAAC,CAAC,EAAE,CAAC;gBAC9D,QAAQ,CAAC,IAAI,CAAC,yCAAyC,CAAC,CAAC;YAC3D,CAAC;YACD,MAAM;QACR,KAAK,WAAW;YACd,IAAI,CAAC,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,6BAA6B,CAAC,CAAC;YACrE,IAAI,GAAG,CAAC,QAAQ,CAAC,GAAG,CAAC;gBAAE,QAAQ,CAAC,IAAI,CAAC,gCAAgC,CAAC,CAAC;YACvE,MAAM;IACV,CAAC;IACD,OAAO,QAAQ,CAAC;AAClB,CAAC;AAED,sFAAsF;AACtF,MAAM,UAAU,WAAW,CAAC,IAAiB;IAC3C,IAAI,OAAO,IAAI,KAAK,QAAQ;QAAE,OAAO,IAAI,CAAC;IAC1C,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,UAAU;YACb,OAAO,IAAI,CAAC,KAAK,CAAC;QACpB,KAAK,SAAS;YACZ,OAAO,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QACnC,KAAK,aAAa,CAAC;QACnB,KAAK,QAAQ,CAAC;QACd,KAAK,WAAW;YACd,OAAO,IAAI,CAAC,IAAI,CAAC;QACnB,KAAK,YAAY;YACf,OAAO,CAAC,IAAI,CAAC,GAAG,GAAG,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IACjD,CAAC;AACH,CAAC;AAED,gFAAgF;AAChF,MAAM,UAAU,aAAa,CAAC,IAAiB,EAAE,GAAQ;IACvD,IAAI,OAAO,IAAI,KAAK,QAAQ;QAAE,OAAO,IAAI,eAAe,CAAC,IAAI,CAAC,CAAC;IAC/D,QAAQ,IAAI,CAAC,IAAI,EAAE,CAAC;QAClB,KAAK,UAAU;YACb,OAAO,IAAI,eAAe,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;QACzC,KAAK,SAAS;YACZ,OAAO,IAAI,cAAc,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACrD,KAAK,aAAa;YAChB,OAAO,IAAI,kBAAkB,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC;QAChD,KAAK,QAAQ;YACX,OAAO,IAAI,aAAa,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;QACxD,KAAK,YAAY;YACf,OAAO,IAAI,iBAAiB,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC;QACnE,KAAK,WAAW;YACd,OAAO,IAAI,gBAAgB,CAAC,GAAG,EAAE,IAAI,CAAC,IAAI,EAAE,IAAI,CAAC,MAAM,CAAC,CAAC;IAC7D,CAAC;AACH,CAAC"}

package/dist/rng/index.d.ts

@@ -0,0 +1,4 @@
+export { Rng, deriveSeed } from "./rng.js";
+export { ConstantSampler, UniformSampler, ExponentialSampler, NormalSampler, TriangularSampler, LognormalSampler, createSampler, validateSamplerSpec, samplerMean, DISTRIBUTION_NAMES, DISTRIBUTION_PARAMS, } from "./distributions.js";
+export type { SampleProvider, DistributionSpec, SamplerSpec } from "./distributions.js";
+//# sourceMappingURL=index.d.ts.map

package/dist/rng/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/rng/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,GAAG,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC;AAC3C,OAAO,EACL,eAAe,EACf,cAAc,EACd,kBAAkB,EAClB,aAAa,EACb,iBAAiB,EACjB,gBAAgB,EAChB,aAAa,EACb,mBAAmB,EACnB,WAAW,EACX,kBAAkB,EAClB,mBAAmB,GACpB,MAAM,oBAAoB,CAAC;AAC5B,YAAY,EAAE,cAAc,EAAE,gBAAgB,EAAE,WAAW,EAAE,MAAM,oBAAoB,CAAC"}

package/dist/rng/index.js

@@ -0,0 +1,3 @@
+export { Rng, deriveSeed } from "./rng.js";
+export { ConstantSampler, UniformSampler, ExponentialSampler, NormalSampler, TriangularSampler, LognormalSampler, createSampler, validateSamplerSpec, samplerMean, DISTRIBUTION_NAMES, DISTRIBUTION_PARAMS, } from "./distributions.js";
+//# sourceMappingURL=index.js.map

package/dist/rng/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/rng/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,GAAG,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC;AAC3C,OAAO,EACL,eAAe,EACf,cAAc,EACd,kBAAkB,EAClB,aAAa,EACb,iBAAiB,EACjB,gBAAgB,EAChB,aAAa,EACb,mBAAmB,EACnB,WAAW,EACX,kBAAkB,EAClB,mBAAmB,GACpB,MAAM,oBAAoB,CAAC"}

package/dist/rng/rng.d.ts

@@ -0,0 +1,26 @@
+/**
+ * Derive a child seed from a base seed and an index (e.g. replication number).
+ * The result is a non-negative integer below 2^53.
+ */
+export declare function deriveSeed(baseSeed: number, index: number): number;
+/**
+ * Seeded pseudo-random number generator (xoshiro128**, 128-bit state).
+ *
+ * A generator is identified by `(seed, streamId)`: the same pair always yields the same sequence,
+ * and different stream ids yield statistically independent sequences. Give every component its
+ * own stream so that changing one component never shifts another component's random numbers.
+ */
+export declare class Rng {
+    private s0;
+    private s1;
+    private s2;
+    private s3;
+    constructor(seed: number, streamId?: string | number);
+    /** Next uniformly distributed unsigned 32-bit integer. */
+    nextUint32(): number;
+    /** Next uniform float in [0, 1) with 53 bits of resolution. */
+    nextFloat(): number;
+    /** Next uniform float in (0, 1] — safe as an argument to Math.log. */
+    nextFloatOpen(): number;
+}
+//# sourceMappingURL=rng.d.ts.map

package/dist/rng/rng.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"rng.d.ts","sourceRoot":"","sources":["../../src/rng/rng.ts"],"names":[],"mappings":"AA6BA;;;GAGG;AACH,wBAAgB,UAAU,CAAC,QAAQ,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,MAAM,CAOlE;AAeD;;;;;;GAMG;AACH,qBAAa,GAAG;IACd,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,EAAE,CAAS;IACnB,OAAO,CAAC,EAAE,CAAS;gBAEP,IAAI,EAAE,MAAM,EAAE,QAAQ,GAAE,MAAM,GAAG,MAAU;IAYvD,0DAA0D;IAC1D,UAAU,IAAI,MAAM;IAcpB,+DAA+D;IAC/D,SAAS,IAAI,MAAM;IAMnB,sEAAsE;IACtE,aAAa,IAAI,MAAM;CAGxB"}

package/dist/rng/rng.js

@@ -0,0 +1,101 @@
+/** murmur3 32-bit finalizer: a strong integer mixer. */
+function fmix32(h) {
+    h ^= h >>> 16;
+    h = Math.imul(h, 0x85ebca6b);
+    h ^= h >>> 13;
+    h = Math.imul(h, 0xc2b2ae35);
+    h ^= h >>> 16;
+    return h >>> 0;
+}
+/** Two independent 32-bit string hashes (FNV-1a and a multiplicative variant) = a 64-bit stream hash. */
+function hashString(s) {
+    let a = 0x811c9dc5;
+    let b = 0x9747b28c;
+    for (let i = 0; i < s.length; i++) {
+        const c = s.charCodeAt(i);
+        a = Math.imul(a ^ c, 0x01000193);
+        b = Math.imul(b + c + 1, 0x5bd1e995);
+        b ^= b >>> 15;
+    }
+    return [fmix32(a), fmix32(b ^ s.length)];
+}
+function splitSeed(seed) {
+    if (!Number.isFinite(seed))
+        throw new RangeError(`seed must be a finite number, got ${seed}`);
+    const s = Math.abs(Math.trunc(seed));
+    return [s >>> 0, Math.floor(s / 4294967296) >>> 0];
+}
+/**
+ * Derive a child seed from a base seed and an index (e.g. replication number).
+ * The result is a non-negative integer below 2^53.
+ */
+export function deriveSeed(baseSeed, index) {
+    const [lo, hi] = splitSeed(baseSeed);
+    const [ilo, ihi] = splitSeed(index);
+    const a = fmix32(lo ^ fmix32(ilo + 0x9e3779b9));
+    const b = fmix32(hi ^ fmix32(ihi + 0x7f4a7c15) ^ a);
+    const c = fmix32(a + b + 0x165667b1);
+    return (c & 0x1fffff) * 4294967296 + fmix32(a ^ (b + 0x27d4eb2f));
+}
+function splitmix32(seed) {
+    let a = seed | 0;
+    return () => {
+        a = (a + 0x9e3779b9) | 0;
+        let t = a ^ (a >>> 16);
+        t = Math.imul(t, 0x21f0aaad);
+        t ^= t >>> 15;
+        t = Math.imul(t, 0x735a2d97);
+        t ^= t >>> 15;
+        return t >>> 0;
+    };
+}
+/**
+ * Seeded pseudo-random number generator (xoshiro128**, 128-bit state).
+ *
+ * A generator is identified by `(seed, streamId)`: the same pair always yields the same sequence,
+ * and different stream ids yield statistically independent sequences. Give every component its
+ * own stream so that changing one component never shifts another component's random numbers.
+ */
+export class Rng {
+    s0;
+    s1;
+    s2;
+    s3;
+    constructor(seed, streamId = 0) {
+        const [lo, hi] = splitSeed(seed);
+        const [h1, h2] = hashString(String(streamId));
+        const g1 = splitmix32(fmix32(lo ^ fmix32(h1)));
+        const g2 = splitmix32(fmix32(hi + 0x632be5ab) ^ fmix32(h2 + 0x1b873593) ^ g1());
+        this.s0 = g1();
+        this.s1 = g1();
+        this.s2 = g2();
+        this.s3 = g2();
+        if ((this.s0 | this.s1 | this.s2 | this.s3) === 0)
+            this.s0 = 1;
+    }
+    /** Next uniformly distributed unsigned 32-bit integer. */
+    nextUint32() {
+        const { s1 } = this;
+        const x = Math.imul(s1, 5);
+        const result = Math.imul((x << 7) | (x >>> 25), 9) >>> 0;
+        const t = s1 << 9;
+        this.s2 ^= this.s0;
+        this.s3 ^= s1;
+        this.s1 ^= this.s2;
+        this.s0 ^= this.s3;
+        this.s2 ^= t;
+        this.s3 = (this.s3 << 11) | (this.s3 >>> 21);
+        return result;
+    }
+    /** Next uniform float in [0, 1) with 53 bits of resolution. */
+    nextFloat() {
+        const a = this.nextUint32() >>> 5; // 27 bits
+        const b = this.nextUint32() >>> 6; // 26 bits
+        return (a * 67108864 + b) / 9007199254740992;
+    }
+    /** Next uniform float in (0, 1] — safe as an argument to Math.log. */
+    nextFloatOpen() {
+        return 1 - this.nextFloat();
+    }
+}
+//# sourceMappingURL=rng.js.map

package/dist/rng/rng.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"rng.js","sourceRoot":"","sources":["../../src/rng/rng.ts"],"names":[],"mappings":"AAAA,wDAAwD;AACxD,SAAS,MAAM,CAAC,CAAS;IACvB,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;IACd,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IAC7B,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;IACd,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;IAC7B,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;IACd,OAAO,CAAC,KAAK,CAAC,CAAC;AACjB,CAAC;AAED,yGAAyG;AACzG,SAAS,UAAU,CAAC,CAAS;IAC3B,IAAI,CAAC,GAAG,UAAU,CAAC;IACnB,IAAI,CAAC,GAAG,UAAU,CAAC;IACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QAClC,MAAM,CAAC,GAAG,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,CAAC;QAC1B,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,UAAU,CAAC,CAAC;QACjC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,UAAU,CAAC,CAAC;QACrC,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;IAChB,CAAC;IACD,OAAO,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC;AAC3C,CAAC;AAED,SAAS,SAAS,CAAC,IAAY;IAC7B,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,IAAI,CAAC;QAAE,MAAM,IAAI,UAAU,CAAC,qCAAqC,IAAI,EAAE,CAAC,CAAC;IAC9F,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC;IACrC,OAAO,CAAC,CAAC,KAAK,CAAC,EAAE,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC,CAAC;AACrD,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,UAAU,CAAC,QAAgB,EAAE,KAAa;IACxD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC,QAAQ,CAAC,CAAC;IACrC,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;IACpC,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,GAAG,MAAM,CAAC,GAAG,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC,CAAC;IACpD,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC;IACrC,OAAO,CAAC,CAAC,GAAG,QAAQ,CAAC,GAAG,UAAU,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;AACpE,CAAC;AAED,SAAS,UAAU,CAAC,IAAY;IAC9B,IAAI,CAAC,GAAG,IAAI,GAAG,CAAC,CAAC;IACjB,OAAO,GAAG,EAAE;QACV,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;QACzB,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC;QACvB,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QAC7B,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;QACd,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC;QAC7B,CAAC,IAAI,CAAC,KAAK,EAAE,CAAC;QACd,OAAO,CAAC,KAAK,CAAC,CAAC;IACjB,CAAC,CAAC;AACJ,CAAC;AAED;;;;;;GAMG;AACH,MAAM,OAAO,GAAG;IACN,EAAE,CAAS;IACX,EAAE,CAAS;IACX,EAAE,CAAS;IACX,EAAE,CAAS;IAEnB,YAAY,IAAY,EAAE,WAA4B,CAAC;QACrD,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC;QACjC,MAAM,CAAC,EAAE,EAAE,EAAE,CAAC,GAAG,UAAU,CAAC,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC;QAC9C,MAAM,EAAE,GAAG,UAAU,CAAC,MAAM,CAAC,EAAE,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;QAC/C,MAAM,EAAE,GAAG,UAAU,CAAC,MAAM,CAAC,EAAE,GAAG,UAAU,CAAC,GAAG,MAAM,CAAC,EAAE,GAAG,UAAU,CAAC,GAAG,EAAE,EAAE,CAAC,CAAC;QAChF,IAAI,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACf,IAAI,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACf,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC;YAAE,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IACjE,CAAC;IAED,0DAA0D;IAC1D,UAAU;QACR,MAAM,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;QACpB,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC;QAC3B,MAAM,MAAM,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC;QACzD,MAAM,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC;QAClB,IAAI,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,IAAI,EAAE,CAAC;QACd,IAAI,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;QACnB,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;QACb,IAAI,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,IAAI,EAAE,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,KAAK,EAAE,CAAC,CAAC;QAC7C,OAAO,MAAM,CAAC;IAChB,CAAC;IAED,+DAA+D;IAC/D,SAAS;QACP,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC,CAAC,UAAU;QAC7C,MAAM,CAAC,GAAG,IAAI,CAAC,UAAU,EAAE,KAAK,CAAC,CAAC,CAAC,UAAU;QAC7C,OAAO,CAAC,CAAC,GAAG,QAAQ,GAAG,CAAC,CAAC,GAAG,gBAAgB,CAAC;IAC/C,CAAC;IAED,sEAAsE;IACtE,aAAa;QACX,OAAO,CAAC,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;IAC9B,CAAC;CACF"}

package/dist/run/assertions.d.ts

@@ -0,0 +1,9 @@
+import type { AssertionDefinition } from "../format/index.js";
+import type { AssertionResult, RunResults } from "./results.js";
+/**
+ * Check assertions against run results. An assertion on an undefined value (an unknown output, no
+ * observations, or a confidence interval with a single replication) fails, because a gate that
+ * silently passes when its data is missing is worse than no gate.
+ */
+export declare function evaluateAssertions(results: RunResults, assertions: readonly AssertionDefinition[]): AssertionResult[];
+//# sourceMappingURL=assertions.d.ts.map

package/dist/run/assertions.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"assertions.d.ts","sourceRoot":"","sources":["../../src/run/assertions.ts"],"names":[],"mappings":"AAAA,OAAO,KAAK,EAAE,mBAAmB,EAAE,MAAM,oBAAoB,CAAC;AAC9D,OAAO,KAAK,EAAE,eAAe,EAAE,UAAU,EAAE,MAAM,cAAc,CAAC;AAyChE;;;;GAIG;AACH,wBAAgB,kBAAkB,CAAC,OAAO,EAAE,UAAU,EAAE,UAAU,EAAE,SAAS,mBAAmB,EAAE,GAAG,eAAe,EAAE,CAYrH"}

package/dist/run/assertions.js

@@ -0,0 +1,57 @@
+function compare(actual, op, value) {
+    switch (op) {
+        case "<":
+            return actual < value;
+        case "<=":
+            return actual <= value;
+        case ">":
+            return actual > value;
+        case ">=":
+            return actual >= value;
+        case "==":
+            return actual === value;
+    }
+}
+function fmt(v) {
+    return Number.isInteger(v) ? String(v) : String(Number(v.toPrecision(5)));
+}
+/** The number an assertion compares, or null when it is undefined for these results. */
+function pick(results, a) {
+    const summary = results.outputs.find((o) => o.id === a.output);
+    if (summary === undefined)
+        return null;
+    switch (a.statistic ?? "mean") {
+        case "mean":
+            return summary.mean;
+        case "ci95Low":
+            return summary.ci95?.low ?? null;
+        case "ci95High":
+            return summary.ci95?.high ?? null;
+        case "min":
+        case "max": {
+            const values = results.replications.map((r) => r.outputs[a.output]).filter((v) => typeof v === "number");
+            if (values.length === 0)
+                return null;
+            return a.statistic === "min" ? Math.min(...values) : Math.max(...values);
+        }
+    }
+}
+/**
+ * Check assertions against run results. An assertion on an undefined value (an unknown output, no
+ * observations, or a confidence interval with a single replication) fails, because a gate that
+ * silently passes when its data is missing is worse than no gate.
+ */
+export function evaluateAssertions(results, assertions) {
+    return assertions.map((assertion) => {
+        const stat = assertion.statistic ?? "mean";
+        const label = assertion.name ? `${assertion.name}: ` : "";
+        const actual = pick(results, assertion);
+        const subject = `${assertion.output} (${stat})`;
+        if (actual === null) {
+            return { assertion, actual: null, passed: false, message: `${label}${subject} is undefined, cannot check ${assertion.op} ${fmt(assertion.value)}` };
+        }
+        const passed = compare(actual, assertion.op, assertion.value);
+        return { assertion, actual, passed, message: `${label}${subject} = ${fmt(actual)} ${passed ? assertion.op : `violates ${assertion.op}`} ${fmt(assertion.value)}` };
+    });
+}
+//# sourceMappingURL=assertions.js.map

package/dist/run/assertions.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"assertions.js","sourceRoot":"","sources":["../../src/run/assertions.ts"],"names":[],"mappings":"AAGA,SAAS,OAAO,CAAC,MAAc,EAAE,EAA6B,EAAE,KAAa;IAC3E,QAAQ,EAAE,EAAE,CAAC;QACX,KAAK,GAAG;YACN,OAAO,MAAM,GAAG,KAAK,CAAC;QACxB,KAAK,IAAI;YACP,OAAO,MAAM,IAAI,KAAK,CAAC;QACzB,KAAK,GAAG;YACN,OAAO,MAAM,GAAG,KAAK,CAAC;QACxB,KAAK,IAAI;YACP,OAAO,MAAM,IAAI,KAAK,CAAC;QACzB,KAAK,IAAI;YACP,OAAO,MAAM,KAAK,KAAK,CAAC;IAC5B,CAAC;AACH,CAAC;AAED,SAAS,GAAG,CAAC,CAAS;IACpB,OAAO,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC5E,CAAC;AAED,wFAAwF;AACxF,SAAS,IAAI,CAAC,OAAmB,EAAE,CAAsB;IACvD,MAAM,OAAO,GAAG,OAAO,CAAC,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,MAAM,CAAC,CAAC;IAC/D,IAAI,OAAO,KAAK,SAAS;QAAE,OAAO,IAAI,CAAC;IACvC,QAAQ,CAAC,CAAC,SAAS,IAAI,MAAM,EAAE,CAAC;QAC9B,KAAK,MAAM;YACT,OAAO,OAAO,CAAC,IAAI,CAAC;QACtB,KAAK,SAAS;YACZ,OAAO,OAAO,CAAC,IAAI,EAAE,GAAG,IAAI,IAAI,CAAC;QACnC,KAAK,UAAU;YACb,OAAO,OAAO,CAAC,IAAI,EAAE,IAAI,IAAI,IAAI,CAAC;QACpC,KAAK,KAAK,CAAC;QACX,KAAK,KAAK,CAAC,CAAC,CAAC;YACX,MAAM,MAAM,GAAG,OAAO,CAAC,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,OAAO,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAe,EAAE,CAAC,OAAO,CAAC,KAAK,QAAQ,CAAC,CAAC;YACtH,IAAI,MAAM,CAAC,MAAM,KAAK,CAAC;gBAAE,OAAO,IAAI,CAAC;YACrC,OAAO,CAAC,CAAC,SAAS,KAAK,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,CAAC,CAAC;QAC3E,CAAC;IACH,CAAC;AACH,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,kBAAkB,CAAC,OAAmB,EAAE,UAA0C;IAChG,OAAO,UAAU,CAAC,GAAG,CAAC,CAAC,SAAS,EAAE,EAAE;QAClC,MAAM,IAAI,GAAG,SAAS,CAAC,SAAS,IAAI,MAAM,CAAC;QAC3C,MAAM,KAAK,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC;QAC1D,MAAM,MAAM,GAAG,IAAI,CAAC,OAAO,EAAE,SAAS,CAAC,CAAC;QACxC,MAAM,OAAO,GAAG,GAAG,SAAS,CAAC,MAAM,KAAK,IAAI,GAAG,CAAC;QAChD,IAAI,MAAM,KAAK,IAAI,EAAE,CAAC;YACpB,OAAO,EAAE,SAAS,EAAE,MAAM,EAAE,IAAI,EAAE,MAAM,EAAE,KAAK,EAAE,OAAO,EAAE,GAAG,KAAK,GAAG,OAAO,+BAA+B,SAAS,CAAC,EAAE,IAAI,GAAG,CAAC,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC;QACtJ,CAAC;QACD,MAAM,MAAM,GAAG,OAAO,CAAC,MAAM,EAAE,SAAS,CAAC,EAAE,EAAE,SAAS,CAAC,KAAK,CAAC,CAAC;QAC9D,OAAO,EAAE,SAAS,EAAE,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,GAAG,KAAK,GAAG,OAAO,MAAM,GAAG,CAAC,MAAM,CAAC,IAAI,MAAM,CAAC,CAAC,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,CAAC,YAAY,SAAS,CAAC,EAAE,EAAE,IAAI,GAAG,CAAC,SAAS,CAAC,KAAK,CAAC,EAAE,EAAE,CAAC;IACrK,CAAC,CAAC,CAAC;AACL,CAAC"}

package/dist/run/index.d.ts

@@ -0,0 +1,6 @@
+export { runModel, defaultTimeSeries, replicationSeed } from "./run.js";
+export type { RunOptions } from "./run.js";
+export { evaluateAssertions } from "./assertions.js";
+export { Simulation } from "./simulation.js";
+export type { RunResults, RunSettingsUsed, ReplicationResult, OutputSummary, TimeSeriesResult, TimeSeriesReplication, AssertionResult, } from "./results.js";
+//# sourceMappingURL=index.d.ts.map

package/dist/run/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../src/run/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,QAAQ,EAAE,iBAAiB,EAAE,eAAe,EAAE,MAAM,UAAU,CAAC;AACxE,YAAY,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC;AAC3C,OAAO,EAAE,kBAAkB,EAAE,MAAM,iBAAiB,CAAC;AACrD,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,YAAY,EACV,UAAU,EACV,eAAe,EACf,iBAAiB,EACjB,aAAa,EACb,gBAAgB,EAChB,qBAAqB,EACrB,eAAe,GAChB,MAAM,cAAc,CAAC"}

package/dist/run/index.js

@@ -0,0 +1,4 @@
+export { runModel, defaultTimeSeries, replicationSeed } from "./run.js";
+export { evaluateAssertions } from "./assertions.js";
+export { Simulation } from "./simulation.js";
+//# sourceMappingURL=index.js.map

package/dist/run/index.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.js","sourceRoot":"","sources":["../../src/run/index.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,QAAQ,EAAE,iBAAiB,EAAE,eAAe,EAAE,MAAM,UAAU,CAAC;AAExE,OAAO,EAAE,kBAAkB,EAAE,MAAM,iBAAiB,CAAC;AACrD,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC"}