@iyulab/u-doe 0.2.0

package/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 iyulab
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

package/README.md

@@ -0,0 +1,127 @@
+# u-doe
+
+Design of Experiments (DOE) library for Rust — classical design generation, effect estimation, RSM, and multi-response optimization.
+
+## Features
+
+- **Design Generation**: Full factorial (2^k), fractional factorial (2^(k-p)), Plackett-Burman screening, Central Composite Design (CCD), Box-Behnken, Taguchi orthogonal arrays, Definitive Screening Design (DSD)
+- **Analysis**: Effect estimation (main effects + 2FI), half-normal plot data, DOE ANOVA with F-statistics and p-values, Taguchi S/N ratio analysis
+- **Response Surface Methodology**: Second-order OLS model fitting, steepest ascent path
+- **Power & Sample Size**: Two-level factorial power, required replicates, power curves
+- **Multi-response Optimization**: Derringer-Suich desirability functions
+- **Coding**: Actual ↔ coded value transformation (Encoder)
+
+## Quick Start
+
+```toml
+[dependencies]
+u-doe = "0.2"
+```
+
+## Examples
+
+### 2^4 Full Factorial Design
+
+```rust
+use u_doe::design::factorial::full_factorial;
+use u_doe::analysis::effects::estimate_effects;
+
+let design = full_factorial(4).unwrap(); // 16 runs, Yates order
+let responses = vec![
+    45.0, 71.0, 48.0, 65.0, 68.0, 60.0, 80.0, 65.0,
+    43.0, 100.0, 45.0, 104.0, 75.0, 86.0, 70.0, 96.0,
+];
+
+let effects = estimate_effects(&design, &responses, 2).unwrap();
+// Prints A = 21.625, AC = -18.125, ...
+for e in &effects {
+    println!("{}: {:.3}", e.name, e.estimate);
+}
+```
+
+### Response Surface Methodology
+
+```rust
+use u_doe::design::ccd::{ccd, AlphaType};
+use u_doe::analysis::rsm::{fit_rsm, steepest_ascent};
+
+let design = ccd(2, AlphaType::Rotatable, 3).unwrap();
+let responses = vec![/* measured values */];
+// let model = fit_rsm(&design, &responses).unwrap();
+// let path = steepest_ascent(&model, 5, 0.5);
+```
+
+### Multi-response Desirability
+
+```rust
+use u_doe::optimization::desirability::{ResponseSpec, overall_desirability};
+
+let specs = vec![
+    ResponseSpec::maximize(50.0, 100.0, 100.0, 1.0),
+    ResponseSpec::minimize(0.0, 0.0, 30.0, 1.0),
+];
+let d = overall_desirability(&specs, &[80.0, 10.0]);
+println!("Overall desirability: {d:.3}");
+```
+
+### Taguchi S/N Ratio
+
+```rust
+use u_doe::analysis::taguchi_sn::{signal_to_noise, sn_factor_effects, SnGoal};
+use u_doe::design::factorial::full_factorial;
+
+// Replicated responses per run (e.g., 3 replicates each)
+let responses = vec![
+    vec![10.0, 11.0, 10.5],  // run 1
+    vec![20.0, 21.0, 19.5],  // run 2
+    vec![15.0, 14.5, 15.5],  // run 3
+    vec![25.0, 26.0, 24.5],  // run 4
+];
+let sn_values = signal_to_noise(&responses, SnGoal::LargerIsBetter).unwrap();
+
+let design = full_factorial(2).unwrap();
+let effects = sn_factor_effects(&design, &sn_values).unwrap();
+for e in &effects {
+    println!("Factor {}: optimal level {}, Δ = {:.2} dB",
+             e.factor_index, e.optimal_level, e.delta);
+}
+```
+
+### Power & Sample Size
+
+```rust
+use u_doe::power::{two_level_factorial_power, required_replicates};
+
+// Power for a 2^3 full factorial, detecting effect of 2σ with α=0.05
+let power = two_level_factorial_power(3, 0, 2, 2.0, 1.0, 0.05);
+println!("Power = {power:.3}");  // ~0.95 at n=2 replicates
+
+// Minimum replicates to achieve 80% power
+let n = required_replicates(3, 0, 2.0, 1.0, 0.05, 0.80, 10);
+println!("Required replicates: {n}");
+```
+
+## Design Types
+
+| Design | Function | Use Case |
+|--------|----------|----------|
+| Full Factorial | `full_factorial(k)` | k = 2..7, all factor combinations |
+| Fractional Factorial | `fractional_factorial(k, p)` | Screening, k=3..7, p=1..3 |
+| Plackett-Burman | `plackett_burman(k)` | Screening, k ≤ 19, N = 8/12/16/20 |
+| CCD | `ccd(k, alpha, n_center)` | RSM, k = 2..6 |
+| Box-Behnken | `box_behnken(k, n_center)` | RSM, k = 3/4/5 |
+| Taguchi | `taguchi_array(name, k)` | Robust design, L4–L27 |
+| DSD | `definitive_screening(k)` | Screening + quadratic, k = 2..12, 2k+1 runs |
+
+## References
+
+- Montgomery, D.C. (2019). *Introduction to Statistical Quality Control*, 8th ed. Wiley.
+- Derringer, G. & Suich, R. (1980). Simultaneous Optimization of Several Response Variables. *Journal of Quality Technology* 12(4), 214–219.
+- Plackett, R.L. & Burman, J.P. (1946). The Design of Optimum Multifactorial Experiments. *Biometrika* 33(4), 305–325.
+- Jones, B. & Nachtsheim, C.J. (2011). A Class of Three-Level Designs for Definitive Screening in the Presence of Second-Order Effects. *Journal of Quality Technology* 43(1), 1–15.
+- Taguchi, G. (1986). *Introduction to Quality Engineering*. Asian Productivity Organization.
+
+## Related
+
+- [`u-analytics`](https://crates.io/crates/u-analytics) — SPC, process capability, statistical analysis
+- [`u-numflow`](https://crates.io/crates/u-numflow) — Math primitives (used internally)

package/package.json

@@ -0,0 +1,33 @@
+{
+  "name": "@iyulab/u-doe",
+  "type": "module",
+  "collaborators": [
+    "iyulab"
+  ],
+  "description": "Design of Experiments (DOE) framework: factorial, Plackett-Burman, CCD, Box-Behnken, Taguchi, effects analysis, RSM, and desirability optimization.",
+  "version": "0.2.0",
+  "license": "MIT",
+  "repository": {
+    "type": "git",
+    "url": "https://github.com/iyulab/u-doe"
+  },
+  "files": [
+    "u_doe_bg.wasm",
+    "u_doe.js",
+    "u_doe_bg.js",
+    "u_doe.d.ts"
+  ],
+  "main": "u_doe.js",
+  "types": "u_doe.d.ts",
+  "sideEffects": [
+    "./u_doe.js",
+    "./snippets/*"
+  ],
+  "keywords": [
+    "doe",
+    "design-of-experiments",
+    "rsm",
+    "factorial",
+    "statistics"
+  ]
+}

package/u_doe.d.ts

@@ -0,0 +1,110 @@
+/* tslint:disable */
+/* eslint-disable */
+
+/**
+ * Generate a Box-Behnken Design (BBD).
+ *
+ * Supported: k = 3, 4, or 5.
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is not 3, 4, or 5, or `n_center == 0`.
+ */
+export function box_behnken(k: number, n_center: number): any;
+
+/**
+ * Generate a Central Composite Design (CCD).
+ *
+ * `design_type`: `"FaceCentered"` | `"Rotatable"` | `"Inscribed"`
+ *
+ * `n_center`: number of center point replicates (≥ 1, default 3 if 0 is passed is rejected).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of range (2..=6), `n_center == 0`,
+ * or `design_type` is unrecognised.
+ */
+export function ccd(k: number, design_type: string, n_center: number): any;
+
+/**
+ * Generate a Definitive Screening Design (DSD).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of the supported range.
+ */
+export function definitive_screening(k: number): any;
+
+/**
+ * Perform DOE ANOVA.
+ *
+ * `design_json`: JSON array-of-arrays `[[f64]]` — the coded design matrix (rows = runs).
+ * `responses`: flat array of response values, one per run.
+ * `factor_names_json`: JSON array of factor name strings used as column labels.
+ * `effect_names_json`: JSON array of effect names to include (e.g. `["A","B","AB"]`).
+ *
+ * Returns an ANOVA result object with `effects`, `residual_ss`, `residual_df`,
+ * `total_ss`, `r_squared`, `r_squared_adj`.
+ *
+ * # Errors
+ * Returns an error string if dimensions do not match or JSON is malformed.
+ */
+export function doe_anova(design_json: any, responses: Float64Array, factor_names_json: any, effect_names_json: any): any;
+
+/**
+ * Generate a 2^k full factorial design.
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of range (1..=7).
+ */
+export function full_factorial(k: number): any;
+
+/**
+ * Compute Taguchi Signal-to-Noise ratios.
+ *
+ * `responses_json`: JSON array-of-arrays `[[f64]]` — one inner array per run,
+ *   containing the replicate measurements for that run.
+ * `goal`: `"LargerIsBetter"` | `"SmallerIsBetter"` | `"NominalIsBest"`
+ *
+ * Returns a flat `[f64]` of SN values in dB, one per run.
+ *
+ * # Errors
+ * Returns an error string if the goal string is unrecognised, or if the
+ * response data violates the requirements for the chosen goal.
+ */
+export function signal_to_noise(responses_json: any, goal: string): any;
+
+/**
+ * Get a Taguchi orthogonal array.
+ *
+ * `name`: `"L4"` | `"L8"` | `"L9"` | `"L12"` | `"L16"` | `"L18"` | `"L27"`
+ *
+ * `k`: number of factors to use (must be ≤ max columns for the array).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if the array name is unknown or `k` exceeds capacity.
+ */
+export function taguchi_array(name: string, k: number): any;
+
+/**
+ * Compute the statistical power of a 2^(k-p) factorial design.
+ *
+ * Uses the normal approximation. Returns power in [0, 1], or 0.0 for invalid
+ * inputs (k=0, p≥k, n_replicates=0, effect_size≤0, sigma≤0).
+ *
+ * # Arguments
+ * * `k` — total number of factors
+ * * `p` — number of generators (p=0 = full factorial)
+ * * `n_replicates` — number of replicates
+ * * `effect_size` — detectable effect δ (in response units)
+ * * `sigma` — process standard deviation σ
+ * * `alpha` — type-I error rate (e.g. 0.05)
+ */
+export function two_level_factorial_power(k: number, p: number, n_replicates: number, effect_size: number, sigma: number, alpha: number): number;

package/u_doe.js

@@ -0,0 +1,9 @@
+/* @ts-self-types="./u_doe.d.ts" */
+
+import * as wasm from "./u_doe_bg.wasm";
+import { __wbg_set_wasm } from "./u_doe_bg.js";
+__wbg_set_wasm(wasm);
+wasm.__wbindgen_start();
+export {
+    box_behnken, ccd, definitive_screening, doe_anova, full_factorial, signal_to_noise, taguchi_array, two_level_factorial_power
+} from "./u_doe_bg.js";

package/u_doe_bg.js

@@ -0,0 +1,553 @@
+/**
+ * Generate a Box-Behnken Design (BBD).
+ *
+ * Supported: k = 3, 4, or 5.
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is not 3, 4, or 5, or `n_center == 0`.
+ * @param {number} k
+ * @param {number} n_center
+ * @returns {any}
+ */
+export function box_behnken(k, n_center) {
+    const ret = wasm.box_behnken(k, n_center);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Generate a Central Composite Design (CCD).
+ *
+ * `design_type`: `"FaceCentered"` | `"Rotatable"` | `"Inscribed"`
+ *
+ * `n_center`: number of center point replicates (≥ 1, default 3 if 0 is passed is rejected).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of range (2..=6), `n_center == 0`,
+ * or `design_type` is unrecognised.
+ * @param {number} k
+ * @param {string} design_type
+ * @param {number} n_center
+ * @returns {any}
+ */
+export function ccd(k, design_type, n_center) {
+    const ptr0 = passStringToWasm0(design_type, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.ccd(k, ptr0, len0, n_center);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Generate a Definitive Screening Design (DSD).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of the supported range.
+ * @param {number} k
+ * @returns {any}
+ */
+export function definitive_screening(k) {
+    const ret = wasm.definitive_screening(k);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Perform DOE ANOVA.
+ *
+ * `design_json`: JSON array-of-arrays `[[f64]]` — the coded design matrix (rows = runs).
+ * `responses`: flat array of response values, one per run.
+ * `factor_names_json`: JSON array of factor name strings used as column labels.
+ * `effect_names_json`: JSON array of effect names to include (e.g. `["A","B","AB"]`).
+ *
+ * Returns an ANOVA result object with `effects`, `residual_ss`, `residual_df`,
+ * `total_ss`, `r_squared`, `r_squared_adj`.
+ *
+ * # Errors
+ * Returns an error string if dimensions do not match or JSON is malformed.
+ * @param {any} design_json
+ * @param {Float64Array} responses
+ * @param {any} factor_names_json
+ * @param {any} effect_names_json
+ * @returns {any}
+ */
+export function doe_anova(design_json, responses, factor_names_json, effect_names_json) {
+    const ptr0 = passArrayF64ToWasm0(responses, wasm.__wbindgen_malloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.doe_anova(design_json, ptr0, len0, factor_names_json, effect_names_json);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Generate a 2^k full factorial design.
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if `k` is out of range (1..=7).
+ * @param {number} k
+ * @returns {any}
+ */
+export function full_factorial(k) {
+    const ret = wasm.full_factorial(k);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Compute Taguchi Signal-to-Noise ratios.
+ *
+ * `responses_json`: JSON array-of-arrays `[[f64]]` — one inner array per run,
+ *   containing the replicate measurements for that run.
+ * `goal`: `"LargerIsBetter"` | `"SmallerIsBetter"` | `"NominalIsBest"`
+ *
+ * Returns a flat `[f64]` of SN values in dB, one per run.
+ *
+ * # Errors
+ * Returns an error string if the goal string is unrecognised, or if the
+ * response data violates the requirements for the chosen goal.
+ * @param {any} responses_json
+ * @param {string} goal
+ * @returns {any}
+ */
+export function signal_to_noise(responses_json, goal) {
+    const ptr0 = passStringToWasm0(goal, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.signal_to_noise(responses_json, ptr0, len0);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Get a Taguchi orthogonal array.
+ *
+ * `name`: `"L4"` | `"L8"` | `"L9"` | `"L12"` | `"L16"` | `"L18"` | `"L27"`
+ *
+ * `k`: number of factors to use (must be ≤ max columns for the array).
+ *
+ * Returns `{ data: [[f64]], factor_names: [str], run_count: usize, factor_count: usize }`.
+ *
+ * # Errors
+ * Returns an error string if the array name is unknown or `k` exceeds capacity.
+ * @param {string} name
+ * @param {number} k
+ * @returns {any}
+ */
+export function taguchi_array(name, k) {
+    const ptr0 = passStringToWasm0(name, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len0 = WASM_VECTOR_LEN;
+    const ret = wasm.taguchi_array(ptr0, len0, k);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return takeFromExternrefTable0(ret[0]);
+}
+
+/**
+ * Compute the statistical power of a 2^(k-p) factorial design.
+ *
+ * Uses the normal approximation. Returns power in [0, 1], or 0.0 for invalid
+ * inputs (k=0, p≥k, n_replicates=0, effect_size≤0, sigma≤0).
+ *
+ * # Arguments
+ * * `k` — total number of factors
+ * * `p` — number of generators (p=0 = full factorial)
+ * * `n_replicates` — number of replicates
+ * * `effect_size` — detectable effect δ (in response units)
+ * * `sigma` — process standard deviation σ
+ * * `alpha` — type-I error rate (e.g. 0.05)
+ * @param {number} k
+ * @param {number} p
+ * @param {number} n_replicates
+ * @param {number} effect_size
+ * @param {number} sigma
+ * @param {number} alpha
+ * @returns {number}
+ */
+export function two_level_factorial_power(k, p, n_replicates, effect_size, sigma, alpha) {
+    const ret = wasm.two_level_factorial_power(k, p, n_replicates, effect_size, sigma, alpha);
+    return ret;
+}
+export function __wbg_Error_83742b46f01ce22d(arg0, arg1) {
+    const ret = Error(getStringFromWasm0(arg0, arg1));
+    return ret;
+}
+export function __wbg_String_8564e559799eccda(arg0, arg1) {
+    const ret = String(arg1);
+    const ptr1 = passStringToWasm0(ret, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len1 = WASM_VECTOR_LEN;
+    getDataViewMemory0().setInt32(arg0 + 4 * 1, len1, true);
+    getDataViewMemory0().setInt32(arg0 + 4 * 0, ptr1, true);
+}
+export function __wbg___wbindgen_boolean_get_c0f3f60bac5a78d1(arg0) {
+    const v = arg0;
+    const ret = typeof(v) === 'boolean' ? v : undefined;
+    return isLikeNone(ret) ? 0xFFFFFF : ret ? 1 : 0;
+}
+export function __wbg___wbindgen_debug_string_5398f5bb970e0daa(arg0, arg1) {
+    const ret = debugString(arg1);
+    const ptr1 = passStringToWasm0(ret, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    const len1 = WASM_VECTOR_LEN;
+    getDataViewMemory0().setInt32(arg0 + 4 * 1, len1, true);
+    getDataViewMemory0().setInt32(arg0 + 4 * 0, ptr1, true);
+}
+export function __wbg___wbindgen_is_function_3c846841762788c1(arg0) {
+    const ret = typeof(arg0) === 'function';
+    return ret;
+}
+export function __wbg___wbindgen_is_object_781bc9f159099513(arg0) {
+    const val = arg0;
+    const ret = typeof(val) === 'object' && val !== null;
+    return ret;
+}
+export function __wbg___wbindgen_jsval_loose_eq_5bcc3bed3c69e72b(arg0, arg1) {
+    const ret = arg0 == arg1;
+    return ret;
+}
+export function __wbg___wbindgen_number_get_34bb9d9dcfa21373(arg0, arg1) {
+    const obj = arg1;
+    const ret = typeof(obj) === 'number' ? obj : undefined;
+    getDataViewMemory0().setFloat64(arg0 + 8 * 1, isLikeNone(ret) ? 0 : ret, true);
+    getDataViewMemory0().setInt32(arg0 + 4 * 0, !isLikeNone(ret), true);
+}
+export function __wbg___wbindgen_string_get_395e606bd0ee4427(arg0, arg1) {
+    const obj = arg1;
+    const ret = typeof(obj) === 'string' ? obj : undefined;
+    var ptr1 = isLikeNone(ret) ? 0 : passStringToWasm0(ret, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+    var len1 = WASM_VECTOR_LEN;
+    getDataViewMemory0().setInt32(arg0 + 4 * 1, len1, true);
+    getDataViewMemory0().setInt32(arg0 + 4 * 0, ptr1, true);
+}
+export function __wbg___wbindgen_throw_6ddd609b62940d55(arg0, arg1) {
+    throw new Error(getStringFromWasm0(arg0, arg1));
+}
+export function __wbg_call_e133b57c9155d22c() { return handleError(function (arg0, arg1) {
+    const ret = arg0.call(arg1);
+    return ret;
+}, arguments); }
+export function __wbg_done_08ce71ee07e3bd17(arg0) {
+    const ret = arg0.done;
+    return ret;
+}
+export function __wbg_get_326e41e095fb2575() { return handleError(function (arg0, arg1) {
+    const ret = Reflect.get(arg0, arg1);
+    return ret;
+}, arguments); }
+export function __wbg_get_unchecked_329cfe50afab7352(arg0, arg1) {
+    const ret = arg0[arg1 >>> 0];
+    return ret;
+}
+export function __wbg_instanceof_ArrayBuffer_101e2bf31071a9f6(arg0) {
+    let result;
+    try {
+        result = arg0 instanceof ArrayBuffer;
+    } catch (_) {
+        result = false;
+    }
+    const ret = result;
+    return ret;
+}
+export function __wbg_instanceof_Uint8Array_740438561a5b956d(arg0) {
+    let result;
+    try {
+        result = arg0 instanceof Uint8Array;
+    } catch (_) {
+        result = false;
+    }
+    const ret = result;
+    return ret;
+}
+export function __wbg_isArray_33b91feb269ff46e(arg0) {
+    const ret = Array.isArray(arg0);
+    return ret;
+}
+export function __wbg_iterator_d8f549ec8fb061b1() {
+    const ret = Symbol.iterator;
+    return ret;
+}
+export function __wbg_length_b3416cf66a5452c8(arg0) {
+    const ret = arg0.length;
+    return ret;
+}
+export function __wbg_length_ea16607d7b61445b(arg0) {
+    const ret = arg0.length;
+    return ret;
+}
+export function __wbg_new_5f486cdf45a04d78(arg0) {
+    const ret = new Uint8Array(arg0);
+    return ret;
+}
+export function __wbg_new_a70fbab9066b301f() {
+    const ret = new Array();
+    return ret;
+}
+export function __wbg_new_ab79df5bd7c26067() {
+    const ret = new Object();
+    return ret;
+}
+export function __wbg_next_11b99ee6237339e3() { return handleError(function (arg0) {
+    const ret = arg0.next();
+    return ret;
+}, arguments); }
+export function __wbg_next_e01a967809d1aa68(arg0) {
+    const ret = arg0.next;
+    return ret;
+}
+export function __wbg_prototypesetcall_d62e5099504357e6(arg0, arg1, arg2) {
+    Uint8Array.prototype.set.call(getArrayU8FromWasm0(arg0, arg1), arg2);
+}
+export function __wbg_set_282384002438957f(arg0, arg1, arg2) {
+    arg0[arg1 >>> 0] = arg2;
+}
+export function __wbg_set_6be42768c690e380(arg0, arg1, arg2) {
+    arg0[arg1] = arg2;
+}
+export function __wbg_value_21fc78aab0322612(arg0) {
+    const ret = arg0.value;
+    return ret;
+}
+export function __wbindgen_cast_0000000000000001(arg0) {
+    // Cast intrinsic for `F64 -> Externref`.
+    const ret = arg0;
+    return ret;
+}
+export function __wbindgen_cast_0000000000000002(arg0, arg1) {
+    // Cast intrinsic for `Ref(String) -> Externref`.
+    const ret = getStringFromWasm0(arg0, arg1);
+    return ret;
+}
+export function __wbindgen_cast_0000000000000003(arg0) {
+    // Cast intrinsic for `U64 -> Externref`.
+    const ret = BigInt.asUintN(64, arg0);
+    return ret;
+}
+export function __wbindgen_init_externref_table() {
+    const table = wasm.__wbindgen_externrefs;
+    const offset = table.grow(4);
+    table.set(0, undefined);
+    table.set(offset + 0, undefined);
+    table.set(offset + 1, null);
+    table.set(offset + 2, true);
+    table.set(offset + 3, false);
+}
+function addToExternrefTable0(obj) {
+    const idx = wasm.__externref_table_alloc();
+    wasm.__wbindgen_externrefs.set(idx, obj);
+    return idx;
+}
+
+function debugString(val) {
+    // primitive types
+    const type = typeof val;
+    if (type == 'number' || type == 'boolean' || val == null) {
+        return  `${val}`;
+    }
+    if (type == 'string') {
+        return `"${val}"`;
+    }
+    if (type == 'symbol') {
+        const description = val.description;
+        if (description == null) {
+            return 'Symbol';
+        } else {
+            return `Symbol(${description})`;
+        }
+    }
+    if (type == 'function') {
+        const name = val.name;
+        if (typeof name == 'string' && name.length > 0) {
+            return `Function(${name})`;
+        } else {
+            return 'Function';
+        }
+    }
+    // objects
+    if (Array.isArray(val)) {
+        const length = val.length;
+        let debug = '[';
+        if (length > 0) {
+            debug += debugString(val[0]);
+        }
+        for(let i = 1; i < length; i++) {
+            debug += ', ' + debugString(val[i]);
+        }
+        debug += ']';
+        return debug;
+    }
+    // Test for built-in
+    const builtInMatches = /\[object ([^\]]+)\]/.exec(toString.call(val));
+    let className;
+    if (builtInMatches && builtInMatches.length > 1) {
+        className = builtInMatches[1];
+    } else {
+        // Failed to match the standard '[object ClassName]'
+        return toString.call(val);
+    }
+    if (className == 'Object') {
+        // we're a user defined class or Object
+        // JSON.stringify avoids problems with cycles, and is generally much
+        // easier than looping through ownProperties of `val`.
+        try {
+            return 'Object(' + JSON.stringify(val) + ')';
+        } catch (_) {
+            return 'Object';
+        }
+    }
+    // errors
+    if (val instanceof Error) {
+        return `${val.name}: ${val.message}\n${val.stack}`;
+    }
+    // TODO we could test for more things here, like `Set`s and `Map`s.
+    return className;
+}
+
+function getArrayU8FromWasm0(ptr, len) {
+    ptr = ptr >>> 0;
+    return getUint8ArrayMemory0().subarray(ptr / 1, ptr / 1 + len);
+}
+
+let cachedDataViewMemory0 = null;
+function getDataViewMemory0() {
+    if (cachedDataViewMemory0 === null || cachedDataViewMemory0.buffer.detached === true || (cachedDataViewMemory0.buffer.detached === undefined && cachedDataViewMemory0.buffer !== wasm.memory.buffer)) {
+        cachedDataViewMemory0 = new DataView(wasm.memory.buffer);
+    }
+    return cachedDataViewMemory0;
+}
+
+let cachedFloat64ArrayMemory0 = null;
+function getFloat64ArrayMemory0() {
+    if (cachedFloat64ArrayMemory0 === null || cachedFloat64ArrayMemory0.byteLength === 0) {
+        cachedFloat64ArrayMemory0 = new Float64Array(wasm.memory.buffer);
+    }
+    return cachedFloat64ArrayMemory0;
+}
+
+function getStringFromWasm0(ptr, len) {
+    ptr = ptr >>> 0;
+    return decodeText(ptr, len);
+}
+
+let cachedUint8ArrayMemory0 = null;
+function getUint8ArrayMemory0() {
+    if (cachedUint8ArrayMemory0 === null || cachedUint8ArrayMemory0.byteLength === 0) {
+        cachedUint8ArrayMemory0 = new Uint8Array(wasm.memory.buffer);
+    }
+    return cachedUint8ArrayMemory0;
+}
+
+function handleError(f, args) {
+    try {
+        return f.apply(this, args);
+    } catch (e) {
+        const idx = addToExternrefTable0(e);
+        wasm.__wbindgen_exn_store(idx);
+    }
+}
+
+function isLikeNone(x) {
+    return x === undefined || x === null;
+}
+
+function passArrayF64ToWasm0(arg, malloc) {
+    const ptr = malloc(arg.length * 8, 8) >>> 0;
+    getFloat64ArrayMemory0().set(arg, ptr / 8);
+    WASM_VECTOR_LEN = arg.length;
+    return ptr;
+}
+
+function passStringToWasm0(arg, malloc, realloc) {
+    if (realloc === undefined) {
+        const buf = cachedTextEncoder.encode(arg);
+        const ptr = malloc(buf.length, 1) >>> 0;
+        getUint8ArrayMemory0().subarray(ptr, ptr + buf.length).set(buf);
+        WASM_VECTOR_LEN = buf.length;
+        return ptr;
+    }
+
+    let len = arg.length;
+    let ptr = malloc(len, 1) >>> 0;
+
+    const mem = getUint8ArrayMemory0();
+
+    let offset = 0;
+
+    for (; offset < len; offset++) {
+        const code = arg.charCodeAt(offset);
+        if (code > 0x7F) break;
+        mem[ptr + offset] = code;
+    }
+    if (offset !== len) {
+        if (offset !== 0) {
+            arg = arg.slice(offset);
+        }
+        ptr = realloc(ptr, len, len = offset + arg.length * 3, 1) >>> 0;
+        const view = getUint8ArrayMemory0().subarray(ptr + offset, ptr + len);
+        const ret = cachedTextEncoder.encodeInto(arg, view);
+
+        offset += ret.written;
+        ptr = realloc(ptr, len, offset, 1) >>> 0;
+    }
+
+    WASM_VECTOR_LEN = offset;
+    return ptr;
+}
+
+function takeFromExternrefTable0(idx) {
+    const value = wasm.__wbindgen_externrefs.get(idx);
+    wasm.__externref_table_dealloc(idx);
+    return value;
+}
+
+let cachedTextDecoder = new TextDecoder('utf-8', { ignoreBOM: true, fatal: true });
+cachedTextDecoder.decode();
+const MAX_SAFARI_DECODE_BYTES = 2146435072;
+let numBytesDecoded = 0;
+function decodeText(ptr, len) {
+    numBytesDecoded += len;
+    if (numBytesDecoded >= MAX_SAFARI_DECODE_BYTES) {
+        cachedTextDecoder = new TextDecoder('utf-8', { ignoreBOM: true, fatal: true });
+        cachedTextDecoder.decode();
+        numBytesDecoded = len;
+    }
+    return cachedTextDecoder.decode(getUint8ArrayMemory0().subarray(ptr, ptr + len));
+}
+
+const cachedTextEncoder = new TextEncoder();
+
+if (!('encodeInto' in cachedTextEncoder)) {
+    cachedTextEncoder.encodeInto = function (arg, view) {
+        const buf = cachedTextEncoder.encode(arg);
+        view.set(buf);
+        return {
+            read: arg.length,
+            written: buf.length
+        };
+    };
+}
+
+let WASM_VECTOR_LEN = 0;
+
+
+let wasm;
+export function __wbg_set_wasm(val) {
+    wasm = val;
+}