@justinelliottcobb/amari-wasm 0.9.1 → 0.9.3

package/amari_wasm.js

@@ -1,5 +1,20 @@
 let wasm;
 
+function addToExternrefTable0(obj) {
+    const idx = wasm.__externref_table_alloc();
+    wasm.__wbindgen_export_2.set(idx, obj);
+    return idx;
+}
+
+function handleError(f, args) {
+    try {
+        return f.apply(this, args);
+    } catch (e) {
+        const idx = addToExternrefTable0(e);
+        wasm.__wbindgen_exn_store(idx);
+    }
+}
+
 let cachedUint8ArrayMemory0 = null;
 
 function getUint8ArrayMemory0() {
@@ -30,6 +45,195 @@ function getStringFromWasm0(ptr, len) {
     return decodeText(ptr, len);
 }
 
+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;
+}
+
+let WASM_VECTOR_LEN = 0;
+
+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
+        };
+    }
+}
+
+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;
+}
+
+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;
+}
+
+function isLikeNone(x) {
+    return x === undefined || x === null;
+}
+
+function takeFromExternrefTable0(idx) {
+    const value = wasm.__wbindgen_export_2.get(idx);
+    wasm.__externref_table_dealloc(idx);
+    return value;
+}
+
+function _assertClass(instance, klass) {
+    if (!(instance instanceof klass)) {
+        throw new Error(`expected instance of ${klass.name}`);
+    }
+}
+/**
+ * Convert velocity to Lorentz factor
+ * @param {number} velocity_magnitude
+ * @returns {number}
+ */
+export function velocity_to_gamma(velocity_magnitude) {
+    const ret = wasm.velocity_to_gamma(velocity_magnitude);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return ret[0];
+}
+
+/**
+ * Convert Lorentz factor to velocity
+ * @param {number} gamma
+ * @returns {number}
+ */
+export function gamma_to_velocity(gamma) {
+    const ret = wasm.gamma_to_velocity(gamma);
+    if (ret[2]) {
+        throw takeFromExternrefTable0(ret[1]);
+    }
+    return ret[0];
+}
+
+/**
+ * Calculate light deflection angle for photon grazing massive object
+ * @param {number} impact_parameter
+ * @param {number} mass
+ * @returns {number}
+ */
+export function light_deflection_angle(impact_parameter, mass) {
+    const ret = wasm.light_deflection_angle(impact_parameter, mass);
+    return ret;
+}
+
+/**
+ * Validate that this module loaded correctly
+ * @returns {boolean}
+ */
+export function validate_relativistic_module() {
+    const ret = wasm.validate_relativistic_module();
+    return ret !== 0;
+}
+
 let cachedFloat64ArrayMemory0 = null;
 
 function getFloat64ArrayMemory0() {
@@ -39,8 +243,6 @@ function getFloat64ArrayMemory0() {
     return cachedFloat64ArrayMemory0;
 }
 
-let WASM_VECTOR_LEN = 0;
-
 function passArrayF64ToWasm0(arg, malloc) {
     const ptr = malloc(arg.length * 8, 8) >>> 0;
     getFloat64ArrayMemory0().set(arg, ptr / 8);
@@ -48,22 +250,10 @@ function passArrayF64ToWasm0(arg, malloc) {
     return ptr;
 }
 
-function takeFromExternrefTable0(idx) {
-    const value = wasm.__wbindgen_export_0.get(idx);
-    wasm.__externref_table_dealloc(idx);
-    return value;
-}
-
 function getArrayF64FromWasm0(ptr, len) {
     ptr = ptr >>> 0;
     return getFloat64ArrayMemory0().subarray(ptr / 8, ptr / 8 + len);
 }
-
-function _assertClass(instance, klass) {
-    if (!(instance instanceof klass)) {
-        throw new Error(`expected instance of ${klass.name}`);
-    }
-}
 /**
  * Initialize the WASM module
  */
@@ -71,6 +261,116 @@ export function init() {
     wasm.init();
 }
 
+let cachedUint32ArrayMemory0 = null;
+
+function getUint32ArrayMemory0() {
+    if (cachedUint32ArrayMemory0 === null || cachedUint32ArrayMemory0.byteLength === 0) {
+        cachedUint32ArrayMemory0 = new Uint32Array(wasm.memory.buffer);
+    }
+    return cachedUint32ArrayMemory0;
+}
+
+function getArrayU32FromWasm0(ptr, len) {
+    ptr = ptr >>> 0;
+    return getUint32ArrayMemory0().subarray(ptr / 4, ptr / 4 + len);
+}
+
+function passArray32ToWasm0(arg, malloc) {
+    const ptr = malloc(arg.length * 4, 4) >>> 0;
+    getUint32ArrayMemory0().set(arg, ptr / 4);
+    WASM_VECTOR_LEN = arg.length;
+    return ptr;
+}
+
+const AutoDiffFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_autodiff_free(ptr >>> 0, 1));
+/**
+ * Automatic differentiation utilities
+ */
+export class AutoDiff {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        AutoDiffFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_autodiff_free(ptr, 0);
+    }
+    /**
+     * Linear layer forward pass (y = Wx + b) with automatic gradients
+     * @param {Float64Array} inputs
+     * @param {Float64Array} weights
+     * @param {Float64Array} bias
+     * @param {number} input_size
+     * @param {number} output_size
+     * @returns {Float64Array}
+     */
+    static linearLayer(inputs, weights, bias, input_size, output_size) {
+        const ptr0 = passArrayF64ToWasm0(inputs, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(weights, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ptr2 = passArrayF64ToWasm0(bias, wasm.__wbindgen_malloc);
+        const len2 = WASM_VECTOR_LEN;
+        const ret = wasm.autodiff_linearLayer(ptr0, len0, ptr1, len1, ptr2, len2, input_size, output_size);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v4 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v4;
+    }
+    /**
+     * Compute mean squared error with automatic gradients
+     * @param {Float64Array} predictions
+     * @param {Float64Array} targets
+     * @returns {WasmDualNumber}
+     */
+    static meanSquaredError(predictions, targets) {
+        const ptr0 = passArrayF64ToWasm0(predictions, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(targets, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.autodiff_meanSquaredError(ptr0, len0, ptr1, len1);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Create a dual number and evaluate a polynomial
+     * @param {number} x
+     * @param {Float64Array} coefficients
+     * @returns {WasmDualNumber}
+     */
+    static evaluatePolynomial(x, coefficients) {
+        const ptr0 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.autodiff_evaluatePolynomial(x, ptr0, len0);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Compute numerical derivative using finite differences (fallback implementation)
+     * @param {number} x
+     * @param {Function} f
+     * @param {number} h
+     * @returns {number}
+     */
+    static numericalDerivative(x, f, h) {
+        const ret = wasm.autodiff_numericalDerivative(x, f, h);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return ret[0];
+    }
+}
+if (Symbol.dispose) AutoDiff.prototype[Symbol.dispose] = AutoDiff.prototype.free;
+
 const BatchOperationsFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
     : new FinalizationRegistry(ptr => wasm.__wbg_batchoperations_free(ptr >>> 0, 1));
@@ -132,389 +432,2312 @@ export class BatchOperations {
 }
 if (Symbol.dispose) BatchOperations.prototype[Symbol.dispose] = BatchOperations.prototype.free;
 
-const PerformanceOperationsFinalization = (typeof FinalizationRegistry === 'undefined')
+const BatchOpsFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
-    : new FinalizationRegistry(ptr => wasm.__wbg_performanceoperations_free(ptr >>> 0, 1));
+    : new FinalizationRegistry(ptr => wasm.__wbg_batchops_free(ptr >>> 0, 1));
 /**
- * High-performance WASM operations with memory pooling
+ * Batch operations for efficient computation
  */
-export class PerformanceOperations {
+export class BatchOps {
 
     __destroy_into_raw() {
         const ptr = this.__wbg_ptr;
         this.__wbg_ptr = 0;
-        PerformanceOperationsFinalization.unregister(this);
+        BatchOpsFinalization.unregister(this);
         return ptr;
     }
 
     free() {
         const ptr = this.__destroy_into_raw();
-        wasm.__wbg_performanceoperations_free(ptr, 0);
+        wasm.__wbg_batchops_free(ptr, 0);
     }
     /**
-     * Batch normalize vectors for efficiency
-     * @param {Float64Array} vectors
-     * @param {number} vector_size
+     * Matrix multiplication with automatic gradients
+     * @param {Float64Array} a
+     * @param {Float64Array} b
+     * @param {number} a_rows
+     * @param {number} a_cols
+     * @param {number} b_rows
+     * @param {number} b_cols
      * @returns {Float64Array}
      */
-    static batchNormalize(vectors, vector_size) {
-        const ptr0 = passArrayF64ToWasm0(vectors, wasm.__wbindgen_malloc);
-        const len0 = WASM_VECTOR_LEN;
-        const ret = wasm.performanceoperations_batchNormalize(ptr0, len0, vector_size);
-        var v2 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
-        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
-        return v2;
-    }
-    /**
-     * Optimized vector dot product
-     * @param {Float64Array} v1
-     * @param {Float64Array} v2
-     * @returns {number}
-     */
-    static vectorDotProduct(v1, v2) {
-        const ptr0 = passArrayF64ToWasm0(v1, wasm.__wbindgen_malloc);
+    static matrixMultiply(a, b, a_rows, a_cols, b_rows, b_cols) {
+        const ptr0 = passArrayF64ToWasm0(a, wasm.__wbindgen_malloc);
         const len0 = WASM_VECTOR_LEN;
-        const ptr1 = passArrayF64ToWasm0(v2, wasm.__wbindgen_malloc);
+        const ptr1 = passArrayF64ToWasm0(b, wasm.__wbindgen_malloc);
         const len1 = WASM_VECTOR_LEN;
-        const ret = wasm.performanceoperations_vectorDotProduct(ptr0, len0, ptr1, len1);
-        return ret;
+        const ret = wasm.batchops_matrixMultiply(ptr0, len0, ptr1, len1, a_rows, a_cols, b_rows, b_cols);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v3;
     }
     /**
-     * Optimized vector operations for 3D space
-     * @param {Float64Array} v1
-     * @param {Float64Array} v2
+     * Batch evaluate polynomial with derivatives
+     * @param {Float64Array} x_values
+     * @param {Float64Array} coefficients
      * @returns {Float64Array}
      */
-    static vectorCrossProduct(v1, v2) {
-        const ptr0 = passArrayF64ToWasm0(v1, wasm.__wbindgen_malloc);
+    static batchPolynomial(x_values, coefficients) {
+        const ptr0 = passArrayF64ToWasm0(x_values, wasm.__wbindgen_malloc);
         const len0 = WASM_VECTOR_LEN;
-        const ptr1 = passArrayF64ToWasm0(v2, wasm.__wbindgen_malloc);
+        const ptr1 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
         const len1 = WASM_VECTOR_LEN;
-        const ret = wasm.performanceoperations_vectorCrossProduct(ptr0, len0, ptr1, len1);
+        const ret = wasm.batchops_batchPolynomial(ptr0, len0, ptr1, len1);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
         var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
         wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
         return v3;
     }
     /**
-     * Fast geometric product for hot paths with memory pooling
-     * @param {Float64Array} lhs
-     * @param {Float64Array} rhs
+     * Compute Jacobian matrix for vector function
+     * @param {Float64Array} input_values
+     * @param {string} function_name
      * @returns {Float64Array}
      */
-    static fastGeometricProduct(lhs, rhs) {
-        const ptr0 = passArrayF64ToWasm0(lhs, wasm.__wbindgen_malloc);
+    static computeJacobian(input_values, function_name) {
+        const ptr0 = passArrayF64ToWasm0(input_values, wasm.__wbindgen_malloc);
         const len0 = WASM_VECTOR_LEN;
-        const ptr1 = passArrayF64ToWasm0(rhs, wasm.__wbindgen_malloc);
+        const ptr1 = passStringToWasm0(function_name, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
         const len1 = WASM_VECTOR_LEN;
-        const ret = wasm.performanceoperations_fastGeometricProduct(ptr0, len0, ptr1, len1);
+        const ret = wasm.batchops_computeJacobian(ptr0, len0, ptr1, len1);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
         var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
         wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
         return v3;
     }
 }
-if (Symbol.dispose) PerformanceOperations.prototype[Symbol.dispose] = PerformanceOperations.prototype.free;
+if (Symbol.dispose) BatchOps.prototype[Symbol.dispose] = BatchOps.prototype.free;
 
-const WasmMultivectorFinalization = (typeof FinalizationRegistry === 'undefined')
+const MLOpsFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
-    : new FinalizationRegistry(ptr => wasm.__wbg_wasmmultivector_free(ptr >>> 0, 1));
+    : new FinalizationRegistry(ptr => wasm.__wbg_mlops_free(ptr >>> 0, 1));
 /**
- * WASM wrapper for Multivector with TypedArray support
+ * Advanced machine learning operations using dual numbers
  */
-export class WasmMultivector {
-
-    static __wrap(ptr) {
-        ptr = ptr >>> 0;
-        const obj = Object.create(WasmMultivector.prototype);
-        obj.__wbg_ptr = ptr;
-        WasmMultivectorFinalization.register(obj, obj.__wbg_ptr, obj);
-        return obj;
-    }
+export class MLOps {
 
     __destroy_into_raw() {
         const ptr = this.__wbg_ptr;
         this.__wbg_ptr = 0;
-        WasmMultivectorFinalization.unregister(this);
+        MLOpsFinalization.unregister(this);
         return ptr;
     }
 
     free() {
         const ptr = this.__destroy_into_raw();
-        wasm.__wbg_wasmmultivector_free(ptr, 0);
+        wasm.__wbg_mlops_free(ptr, 0);
     }
     /**
-     * Create a basis vector (0-indexed)
-     * @param {number} index
-     * @returns {WasmMultivector}
+     * Batch apply activation function with gradients
+     * @param {Float64Array} inputs
+     * @param {string} activation
+     * @returns {Float64Array}
      */
-    static basisVector(index) {
-        const ret = wasm.wasmmultivector_basisVector(index);
+    static batchActivation(inputs, activation) {
+        const ptr0 = passArrayF64ToWasm0(inputs, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passStringToWasm0(activation, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.mlops_batchActivation(ptr0, len0, ptr1, len1);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v3;
+    }
+    /**
+     * Compute cross-entropy loss with automatic gradients
+     * @param {Float64Array} predictions
+     * @param {Float64Array} targets
+     * @returns {WasmDualNumber}
+     */
+    static crossEntropyLoss(predictions, targets) {
+        const ptr0 = passArrayF64ToWasm0(predictions, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(targets, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.mlops_crossEntropyLoss(ptr0, len0, ptr1, len1);
         if (ret[2]) {
             throw takeFromExternrefTable0(ret[1]);
         }
-        return WasmMultivector.__wrap(ret[0]);
+        return WasmDualNumber.__wrap(ret[0]);
     }
     /**
-     * Inner product (dot product for vectors)
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Gradient descent step
+     * @param {Float64Array} parameters
+     * @param {Float64Array} gradients
+     * @param {number} learning_rate
+     * @returns {Float64Array}
      */
-    innerProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_innerProduct(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    static gradientDescentStep(parameters, gradients, learning_rate) {
+        const ptr0 = passArrayF64ToWasm0(parameters, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(gradients, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.mlops_gradientDescentStep(ptr0, len0, ptr1, len1, learning_rate);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v3;
     }
     /**
-     * Outer product (wedge product)
+     * Compute softmax with automatic gradients
+     * @param {Float64Array} inputs
+     * @returns {Float64Array}
+     */
+    static softmax(inputs) {
+        const ptr0 = passArrayF64ToWasm0(inputs, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.mlops_softmax(ptr0, len0);
+        var v2 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v2;
+    }
+}
+if (Symbol.dispose) MLOps.prototype[Symbol.dispose] = MLOps.prototype.free;
+
+const NetworkBatchOperationsFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_networkbatchoperations_free(ptr >>> 0, 1));
+/**
+ * Batch network operations for performance
+ */
+export class NetworkBatchOperations {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        NetworkBatchOperationsFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_networkbatchoperations_free(ptr, 0);
+    }
+    /**
+     * Fast geometric clustering using k-means in 3D space
+     * @param {Float64Array} positions
+     * @param {number} num_nodes
+     * @param {number} k
+     * @param {number} max_iterations
+     * @returns {Uint32Array}
+     */
+    static geometricKMeansClustering(positions, num_nodes, k, max_iterations) {
+        const ptr0 = passArrayF64ToWasm0(positions, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.networkbatchoperations_geometricKMeansClustering(ptr0, len0, num_nodes, k, max_iterations);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v2 = getArrayU32FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 4, 4);
+        return v2;
+    }
+    /**
+     * Compute all pairwise distances in a network
+     * @param {Float64Array} positions
+     * @param {number} num_nodes
+     * @returns {Float64Array}
+     */
+    static computeAllPairwiseDistances(positions, num_nodes) {
+        const ptr0 = passArrayF64ToWasm0(positions, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.networkbatchoperations_computeAllPairwiseDistances(ptr0, len0, num_nodes);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v2 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v2;
+    }
+}
+if (Symbol.dispose) NetworkBatchOperations.prototype[Symbol.dispose] = NetworkBatchOperations.prototype.free;
+
+const PerformanceOperationsFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_performanceoperations_free(ptr >>> 0, 1));
+/**
+ * High-performance WASM operations with memory pooling
+ */
+export class PerformanceOperations {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        PerformanceOperationsFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_performanceoperations_free(ptr, 0);
+    }
+    /**
+     * Batch normalize vectors for efficiency
+     * @param {Float64Array} vectors
+     * @param {number} vector_size
+     * @returns {Float64Array}
+     */
+    static batchNormalize(vectors, vector_size) {
+        const ptr0 = passArrayF64ToWasm0(vectors, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.performanceoperations_batchNormalize(ptr0, len0, vector_size);
+        var v2 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v2;
+    }
+    /**
+     * Optimized vector dot product
+     * @param {Float64Array} v1
+     * @param {Float64Array} v2
+     * @returns {number}
+     */
+    static vectorDotProduct(v1, v2) {
+        const ptr0 = passArrayF64ToWasm0(v1, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(v2, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.performanceoperations_vectorDotProduct(ptr0, len0, ptr1, len1);
+        return ret;
+    }
+    /**
+     * Optimized vector operations for 3D space
+     * @param {Float64Array} v1
+     * @param {Float64Array} v2
+     * @returns {Float64Array}
+     */
+    static vectorCrossProduct(v1, v2) {
+        const ptr0 = passArrayF64ToWasm0(v1, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(v2, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.performanceoperations_vectorCrossProduct(ptr0, len0, ptr1, len1);
+        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v3;
+    }
+    /**
+     * Fast geometric product for hot paths with memory pooling
+     * @param {Float64Array} lhs
+     * @param {Float64Array} rhs
+     * @returns {Float64Array}
+     */
+    static fastGeometricProduct(lhs, rhs) {
+        const ptr0 = passArrayF64ToWasm0(lhs, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(rhs, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.performanceoperations_fastGeometricProduct(ptr0, len0, ptr1, len1);
+        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v3;
+    }
+}
+if (Symbol.dispose) PerformanceOperations.prototype[Symbol.dispose] = PerformanceOperations.prototype.free;
+
+const TropicalBatchFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_tropicalbatch_free(ptr >>> 0, 1));
+/**
+ * Batch operations for tropical numbers
+ */
+export class TropicalBatch {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        TropicalBatchFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_tropicalbatch_free(ptr, 0);
+    }
+    /**
+     * Convert array of log probabilities to tropical numbers and find maximum
+     * @param {Float64Array} log_probs
+     * @returns {number}
+     */
+    static maxLogProb(log_probs) {
+        const ptr0 = passArrayF64ToWasm0(log_probs, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.tropicalbatch_maxLogProb(ptr0, len0);
+        return ret;
+    }
+    /**
+     * Viterbi algorithm helper: find best path through trellis
+     * @param {Float64Array} prev_scores
+     * @param {Float64Array} transition_scores
+     * @param {Float64Array} emission_scores
+     * @param {number} num_states
+     * @returns {Float64Array}
+     */
+    static viterbiStep(prev_scores, transition_scores, emission_scores, num_states) {
+        const ptr0 = passArrayF64ToWasm0(prev_scores, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(transition_scores, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ptr2 = passArrayF64ToWasm0(emission_scores, wasm.__wbindgen_malloc);
+        const len2 = WASM_VECTOR_LEN;
+        const ret = wasm.tropicalbatch_viterbiStep(ptr0, len0, ptr1, len1, ptr2, len2, num_states);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v4 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v4;
+    }
+    /**
+     * Batch tropical addition (max operation)
+     * @param {Float64Array} values
+     * @returns {number}
+     */
+    static batchTropicalAdd(values) {
+        const ptr0 = passArrayF64ToWasm0(values, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.tropicalbatch_batchTropicalAdd(ptr0, len0);
+        return ret;
+    }
+    /**
+     * Batch tropical multiplication (addition)
+     * @param {Float64Array} values
+     * @returns {number}
+     */
+    static batchTropicalMul(values) {
+        const ptr0 = passArrayF64ToWasm0(values, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.tropicalbatch_batchTropicalMul(ptr0, len0);
+        return ret;
+    }
+}
+if (Symbol.dispose) TropicalBatch.prototype[Symbol.dispose] = TropicalBatch.prototype.free;
+
+const TropicalMLOpsFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_tropicalmlops_free(ptr >>> 0, 1));
+/**
+ * Advanced tropical operations for machine learning and optimization
+ */
+export class TropicalMLOps {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        TropicalMLOpsFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_tropicalmlops_free(ptr, 0);
+    }
+    /**
+     * Compute shortest paths using tropical algebra (Floyd-Warshall)
+     * @param {any} distance_matrix
+     * @returns {Array<any>}
+     */
+    static shortestPaths(distance_matrix) {
+        const ret = wasm.tropicalmlops_shortestPaths(distance_matrix);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Compute tropical matrix multiplication for pathfinding
+     * @param {any} a
+     * @param {any} b
+     * @returns {Array<any>}
+     */
+    static matrixMultiply(a, b) {
+        const ret = wasm.tropicalmlops_matrixMultiply(a, b);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Compute tropical convex combination (used in optimization)
+     * @param {Float64Array} values
+     * @param {Float64Array} weights
+     * @returns {number}
+     */
+    static convexCombination(values, weights) {
+        const ptr0 = passArrayF64ToWasm0(values, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ptr1 = passArrayF64ToWasm0(weights, wasm.__wbindgen_malloc);
+        const len1 = WASM_VECTOR_LEN;
+        const ret = wasm.tropicalmlops_convexCombination(ptr0, len0, ptr1, len1);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return ret[0];
+    }
+}
+if (Symbol.dispose) TropicalMLOps.prototype[Symbol.dispose] = TropicalMLOps.prototype.free;
+
+const WasmDualNumberFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmdualnumber_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for single-variable dual numbers
+ */
+export class WasmDualNumber {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmDualNumber.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmDualNumberFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmDualNumberFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmdualnumber_free(ptr, 0);
+    }
+    /**
+     * Natural logarithm
+     * @returns {WasmDualNumber}
+     */
+    ln() {
+        const ret = wasm.wasmdualnumber_ln(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Addition
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    add(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_add(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Cosine function
+     * @returns {WasmDualNumber}
+     */
+    cos() {
+        const ret = wasm.wasmdualnumber_cos(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Division
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    div(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_div(this.__wbg_ptr, other.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Exponential function
+     * @returns {WasmDualNumber}
+     */
+    exp() {
+        const ret = wasm.wasmdualnumber_exp(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Maximum of two dual numbers
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    max(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_max(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Minimum of two dual numbers
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    min(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_min(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Multiplication
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    mul(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_mul(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Negation
+     * @returns {WasmDualNumber}
+     */
+    neg() {
+        const ret = wasm.wasmdualnumber_neg(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Create a new dual number with given real and dual parts
+     * @param {number} real
+     * @param {number} dual
+     */
+    constructor(real, dual) {
+        const ret = wasm.wasmdualnumber_new(real, dual);
+        this.__wbg_ptr = ret >>> 0;
+        WasmDualNumberFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Power function
+     * @param {number} exponent
+     * @returns {WasmDualNumber}
+     */
+    pow(exponent) {
+        const ret = wasm.wasmdualnumber_pow(this.__wbg_ptr, exponent);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Sine function
+     * @returns {WasmDualNumber}
+     */
+    sin() {
+        const ret = wasm.wasmdualnumber_sin(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Subtraction
+     * @param {WasmDualNumber} other
+     * @returns {WasmDualNumber}
+     */
+    sub(other) {
+        _assertClass(other, WasmDualNumber);
+        const ret = wasm.wasmdualnumber_sub(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Tangent function
+     * @returns {WasmDualNumber}
+     */
+    tan() {
+        const ret = wasm.wasmdualnumber_tan(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Hyperbolic cosine
+     * @returns {WasmDualNumber}
+     */
+    cosh() {
+        const ret = wasm.wasmdualnumber_cosh(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Integer power
+     * @param {number} n
+     * @returns {WasmDualNumber}
+     */
+    powi(n) {
+        const ret = wasm.wasmdualnumber_powi(this.__wbg_ptr, n);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * ReLU activation function
+     * @returns {WasmDualNumber}
+     */
+    relu() {
+        const ret = wasm.wasmdualnumber_relu(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Hyperbolic sine
+     * @returns {WasmDualNumber}
+     */
+    sinh() {
+        const ret = wasm.wasmdualnumber_sinh(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Square root
+     * @returns {WasmDualNumber}
+     */
+    sqrt() {
+        const ret = wasm.wasmdualnumber_sqrt(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Hyperbolic tangent
+     * @returns {WasmDualNumber}
+     */
+    tanh() {
+        const ret = wasm.wasmdualnumber_tanh(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Sigmoid activation function
+     * @returns {WasmDualNumber}
+     */
+    sigmoid() {
+        const ret = wasm.wasmdualnumber_sigmoid(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Create a constant (derivative = 0)
+     * @param {number} value
+     * @returns {WasmDualNumber}
+     */
+    static constant(value) {
+        const ret = wasm.wasmdualnumber_constant(value);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Get the dual part (derivative)
+     * @returns {number}
+     */
+    getDual() {
+        const ret = wasm.wasmdualnumber_getDual(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get the real part (function value)
+     * @returns {number}
+     */
+    getReal() {
+        const ret = wasm.wasmdualnumber_getReal(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Softplus activation function
+     * @returns {WasmDualNumber}
+     */
+    softplus() {
+        const ret = wasm.wasmdualnumber_softplus(this.__wbg_ptr);
+        return WasmDualNumber.__wrap(ret);
+    }
+    /**
+     * Create a variable (derivative = 1)
+     * @param {number} value
+     * @returns {WasmDualNumber}
+     */
+    static variable(value) {
+        const ret = wasm.wasmdualnumber_variable(value);
+        return WasmDualNumber.__wrap(ret);
+    }
+}
+if (Symbol.dispose) WasmDualNumber.prototype[Symbol.dispose] = WasmDualNumber.prototype.free;
+
+const WasmFourVelocityFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmfourvelocity_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for four-velocity
+ */
+export class WasmFourVelocity {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmFourVelocity.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmFourVelocityFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmFourVelocityFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmfourvelocity_free(ptr, 0);
+    }
+    /**
+     * Create four-velocity from 3-velocity components
+     * @param {number} vx
+     * @param {number} vy
+     * @param {number} vz
+     * @returns {WasmFourVelocity}
+     */
+    static from_velocity(vx, vy, vz) {
+        const ret = wasm.wasmfourvelocity_from_velocity(vx, vy, vz);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmFourVelocity.__wrap(ret[0]);
+    }
+    /**
+     * Check if normalized (u·u = c²)
+     * @returns {boolean}
+     */
+    is_normalized() {
+        const ret = wasm.wasmfourvelocity_is_normalized(this.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Get as spacetime vector
+     * @returns {WasmSpacetimeVector}
+     */
+    as_spacetime_vector() {
+        const ret = wasm.wasmfourvelocity_as_spacetime_vector(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
+    }
+    /**
+     * Get spatial velocity magnitude
+     * @returns {number}
+     */
+    spatial_velocity_magnitude() {
+        const ret = wasm.wasmfourvelocity_spatial_velocity_magnitude(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get Lorentz factor γ
+     * @returns {number}
+     */
+    gamma() {
+        const ret = wasm.wasmfourvelocity_gamma(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get rapidity
+     * @returns {number}
+     */
+    rapidity() {
+        const ret = wasm.wasmfourvelocity_rapidity(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get string representation
+     * @returns {string}
+     */
+    to_string() {
+        let deferred1_0;
+        let deferred1_1;
+        try {
+            const ret = wasm.wasmfourvelocity_to_string(this.__wbg_ptr);
+            deferred1_0 = ret[0];
+            deferred1_1 = ret[1];
+            return getStringFromWasm0(ret[0], ret[1]);
+        } finally {
+            wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+        }
+    }
+}
+if (Symbol.dispose) WasmFourVelocity.prototype[Symbol.dispose] = WasmFourVelocity.prototype.free;
+
+const WasmGeodesicIntegratorFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmgeodesicintegrator_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for geodesic integration
+ */
+export class WasmGeodesicIntegrator {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmGeodesicIntegrator.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmGeodesicIntegratorFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmGeodesicIntegratorFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmgeodesicintegrator_free(ptr, 0);
+    }
+    /**
+     * Propagate particle through spacetime
+     * @param {WasmRelativisticParticle} particle
+     * @param {number} integration_time
+     * @param {number} time_step
+     * @returns {Array<any>}
+     */
+    propagate_particle(particle, integration_time, time_step) {
+        _assertClass(particle, WasmRelativisticParticle);
+        const ret = wasm.wasmgeodesicintegrator_propagate_particle(this.__wbg_ptr, particle.__wbg_ptr, integration_time, time_step);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Create integrator with Schwarzschild metric
+     * @param {WasmSchwarzschildMetric} metric
+     * @returns {WasmGeodesicIntegrator}
+     */
+    static with_schwarzschild(metric) {
+        _assertClass(metric, WasmSchwarzschildMetric);
+        const ret = wasm.wasmgeodesicintegrator_with_schwarzschild(metric.__wbg_ptr);
+        return WasmGeodesicIntegrator.__wrap(ret);
+    }
+}
+if (Symbol.dispose) WasmGeodesicIntegrator.prototype[Symbol.dispose] = WasmGeodesicIntegrator.prototype.free;
+
+const WasmGeometricNetworkFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmgeometricnetwork_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for GeometricNetwork
+ */
+export class WasmGeometricNetwork {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmGeometricNetworkFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmgeometricnetwork_free(ptr, 0);
+    }
+    /**
+     * Get degree of a node
+     * @param {number} node
+     * @returns {number}
+     */
+    getDegree(node) {
+        const ret = wasm.wasmgeometricnetwork_getDegree(this.__wbg_ptr, node);
+        return ret >>> 0;
+    }
+    /**
+     * Get neighbors of a node
+     * @param {number} node
+     * @returns {Uint32Array}
+     */
+    getNeighbors(node) {
+        const ret = wasm.wasmgeometricnetwork_getNeighbors(this.__wbg_ptr, node);
+        var v1 = getArrayU32FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 4, 4);
+        return v1;
+    }
+    /**
+     * Find shortest path between two nodes
+     * Returns the path as an array of node indices, or null if no path exists
+     * @param {number} source
+     * @param {number} target
+     * @returns {Uint32Array | undefined}
+     */
+    shortestPath(source, target) {
+        const ret = wasm.wasmgeometricnetwork_shortestPath(this.__wbg_ptr, source, target);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        let v1;
+        if (ret[0] !== 0) {
+            v1 = getArrayU32FromWasm0(ret[0], ret[1]).slice();
+            wasm.__wbindgen_free(ret[0], ret[1] * 4, 4);
+        }
+        return v1;
+    }
+    /**
+     * Add a labeled node
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     * @param {string} label
+     * @returns {number}
+     */
+    addLabeledNode(x, y, z, label) {
+        const ptr0 = passStringToWasm0(label, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmgeometricnetwork_addLabeledNode(this.__wbg_ptr, x, y, z, ptr0, len0);
+        return ret >>> 0;
+    }
+    /**
+     * Find communities using geometric clustering
+     * @param {number} num_communities
+     * @returns {Uint32Array}
+     */
+    findCommunities(num_communities) {
+        const ret = wasm.wasmgeometricnetwork_findCommunities(this.__wbg_ptr, num_communities);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v1 = getArrayU32FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 4, 4);
+        return v1;
+    }
+    /**
+     * Get node position as [x, y, z]
+     * @param {number} node
+     * @returns {Float64Array | undefined}
+     */
+    getNodePosition(node) {
+        const ret = wasm.wasmgeometricnetwork_getNodePosition(this.__wbg_ptr, node);
+        let v1;
+        if (ret[0] !== 0) {
+            v1 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+            wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        }
+        return v1;
+    }
+    /**
+     * Compute geometric distance between two nodes
+     * @param {number} node1
+     * @param {number} node2
+     * @returns {number}
+     */
+    geometricDistance(node1, node2) {
+        const ret = wasm.wasmgeometricnetwork_geometricDistance(this.__wbg_ptr, node1, node2);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return ret[0];
+    }
+    /**
+     * Add an undirected edge
+     * @param {number} a
+     * @param {number} b
+     * @param {number} weight
+     */
+    addUndirectedEdge(a, b, weight) {
+        const ret = wasm.wasmgeometricnetwork_addUndirectedEdge(this.__wbg_ptr, a, b, weight);
+        if (ret[1]) {
+            throw takeFromExternrefTable0(ret[0]);
+        }
+    }
+    /**
+     * Find shortest path with distance between two nodes
+     * Returns an object with path and distance, or null if no path exists
+     * @param {number} source
+     * @param {number} target
+     * @returns {any}
+     */
+    shortestPathWithDistance(source, target) {
+        const ret = wasm.wasmgeometricnetwork_shortestPathWithDistance(this.__wbg_ptr, source, target);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Compute geometric centrality for all nodes
+     * @returns {Float64Array}
+     */
+    computeGeometricCentrality() {
+        const ret = wasm.wasmgeometricnetwork_computeGeometricCentrality(this.__wbg_ptr);
+        if (ret[3]) {
+            throw takeFromExternrefTable0(ret[2]);
+        }
+        var v1 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v1;
+    }
+    /**
+     * Create a new geometric network
+     */
+    constructor() {
+        const ret = wasm.wasmgeometricnetwork_new();
+        this.__wbg_ptr = ret >>> 0;
+        WasmGeometricNetworkFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Add an edge between two nodes
+     * @param {number} source
+     * @param {number} target
+     * @param {number} weight
+     */
+    addEdge(source, target, weight) {
+        const ret = wasm.wasmgeometricnetwork_addEdge(this.__wbg_ptr, source, target, weight);
+        if (ret[1]) {
+            throw takeFromExternrefTable0(ret[0]);
+        }
+    }
+    /**
+     * Add a node at the given position
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     * @returns {number}
+     */
+    addNode(x, y, z) {
+        const ret = wasm.wasmgeometricnetwork_addNode(this.__wbg_ptr, x, y, z);
+        return ret >>> 0;
+    }
+    /**
+     * Get number of edges
+     * @returns {number}
+     */
+    numEdges() {
+        const ret = wasm.wasmgeometricnetwork_numEdges(this.__wbg_ptr);
+        return ret >>> 0;
+    }
+    /**
+     * Get number of nodes
+     * @returns {number}
+     */
+    numNodes() {
+        const ret = wasm.wasmgeometricnetwork_numNodes(this.__wbg_ptr);
+        return ret >>> 0;
+    }
+}
+if (Symbol.dispose) WasmGeometricNetwork.prototype[Symbol.dispose] = WasmGeometricNetwork.prototype.free;
+
+const WasmMultiDualNumberFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmmultidualnumber_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for multi-variable dual numbers
+ */
+export class WasmMultiDualNumber {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmMultiDualNumber.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmMultiDualNumberFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmMultiDualNumberFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmmultidualnumber_free(ptr, 0);
+    }
+    /**
+     * Get a specific partial derivative
+     * @param {number} index
+     * @returns {number}
+     */
+    getPartial(index) {
+        const ret = wasm.wasmmultidualnumber_getPartial(this.__wbg_ptr, index);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return ret[0];
+    }
+    /**
+     * Get the gradient (all partial derivatives)
+     * @returns {Float64Array}
+     */
+    getGradient() {
+        const ret = wasm.wasmmultidualnumber_getGradient(this.__wbg_ptr);
+        var v1 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v1;
+    }
+    /**
+     * Get number of variables
+     * @returns {number}
+     */
+    getNumVars() {
+        const ret = wasm.wasmmultidualnumber_getNumVars(this.__wbg_ptr);
+        return ret >>> 0;
+    }
+    /**
+     * Addition
+     * @param {WasmMultiDualNumber} other
+     * @returns {WasmMultiDualNumber}
+     */
+    add(other) {
+        _assertClass(other, WasmMultiDualNumber);
+        const ret = wasm.wasmmultidualnumber_add(this.__wbg_ptr, other.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultiDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Multiplication
+     * @param {WasmMultiDualNumber} other
+     * @returns {WasmMultiDualNumber}
+     */
+    mul(other) {
+        _assertClass(other, WasmMultiDualNumber);
+        const ret = wasm.wasmmultidualnumber_mul(this.__wbg_ptr, other.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultiDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Create a new multi-dual number
+     * @param {number} real
+     * @param {Float64Array} duals
+     */
+    constructor(real, duals) {
+        const ptr0 = passArrayF64ToWasm0(duals, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmmultidualnumber_new(real, ptr0, len0);
+        this.__wbg_ptr = ret >>> 0;
+        WasmMultiDualNumberFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Square root
+     * @returns {WasmMultiDualNumber}
+     */
+    sqrt() {
+        const ret = wasm.wasmmultidualnumber_sqrt(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultiDualNumber.__wrap(ret[0]);
+    }
+    /**
+     * Create a constant (all derivatives are zero)
+     * @param {number} value
+     * @param {number} num_vars
+     * @returns {WasmMultiDualNumber}
+     */
+    static constant(value, num_vars) {
+        const ret = wasm.wasmmultidualnumber_constant(value, num_vars);
+        return WasmMultiDualNumber.__wrap(ret);
+    }
+    /**
+     * Get the real part (function value)
+     * @returns {number}
+     */
+    getReal() {
+        const ret = wasm.wasmdualnumber_getReal(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Create a variable with derivative 1 at the specified index
+     * @param {number} value
+     * @param {number} num_vars
+     * @param {number} var_index
+     * @returns {WasmMultiDualNumber}
+     */
+    static variable(value, num_vars, var_index) {
+        const ret = wasm.wasmmultidualnumber_variable(value, num_vars, var_index);
+        return WasmMultiDualNumber.__wrap(ret);
+    }
+}
+if (Symbol.dispose) WasmMultiDualNumber.prototype[Symbol.dispose] = WasmMultiDualNumber.prototype.free;
+
+const WasmMultivectorFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmmultivector_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for Multivector with TypedArray support
+ */
+export class WasmMultivector {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmMultivector.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmMultivectorFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmMultivectorFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmmultivector_free(ptr, 0);
+    }
+    /**
+     * Create a basis vector (0-indexed)
+     * @param {number} index
+     * @returns {WasmMultivector}
+     */
+    static basisVector(index) {
+        const ret = wasm.wasmmultivector_basisVector(index);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultivector.__wrap(ret[0]);
+    }
+    /**
+     * Inner product (dot product for vectors)
+     * @param {WasmMultivector} other
+     * @returns {WasmMultivector}
+     */
+    innerProduct(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_innerProduct(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Outer product (wedge product)
+     * @param {WasmMultivector} other
+     * @returns {WasmMultivector}
+     */
+    outerProduct(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_outerProduct(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Scalar product
+     * @param {WasmMultivector} other
+     * @returns {number}
+     */
+    scalarProduct(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_scalarProduct(this.__wbg_ptr, other.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get a specific coefficient
+     * @param {number} index
+     * @returns {number}
+     */
+    getCoefficient(index) {
+        const ret = wasm.wasmmultivector_getCoefficient(this.__wbg_ptr, index);
+        return ret;
+    }
+    /**
+     * Set a specific coefficient
+     * @param {number} index
+     * @param {number} value
+     */
+    setCoefficient(index, value) {
+        wasm.wasmmultivector_setCoefficient(this.__wbg_ptr, index, value);
+    }
+    /**
+     * Get coefficients as a Float64Array
+     * @returns {Float64Array}
+     */
+    getCoefficients() {
+        const ret = wasm.wasmmultivector_getCoefficients(this.__wbg_ptr);
+        var v1 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
+        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
+        return v1;
+    }
+    /**
+     * Grade projection
+     * @param {number} grade
+     * @returns {WasmMultivector}
+     */
+    gradeProjection(grade) {
+        const ret = wasm.wasmmultivector_gradeProjection(this.__wbg_ptr, grade);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Create from a Float64Array of coefficients
+     * @param {Float64Array} coefficients
+     * @returns {WasmMultivector}
+     */
+    static fromCoefficients(coefficients) {
+        const ptr0 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmmultivector_fromCoefficients(ptr0, len0);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultivector.__wrap(ret[0]);
+    }
+    /**
+     * Geometric product
+     * @param {WasmMultivector} other
+     * @returns {WasmMultivector}
+     */
+    geometricProduct(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_geometricProduct(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Add two multivectors
+     * @param {WasmMultivector} other
+     * @returns {WasmMultivector}
+     */
+    add(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_add(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Exponential (for bivectors to create rotors)
+     * @returns {WasmMultivector}
+     */
+    exp() {
+        const ret = wasm.wasmmultivector_exp(this.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Create a new zero multivector
+     */
+    constructor() {
+        const ret = wasm.wasmmultivector_new();
+        this.__wbg_ptr = ret >>> 0;
+        WasmMultivectorFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Subtract two multivectors
      * @param {WasmMultivector} other
      * @returns {WasmMultivector}
      */
-    outerProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_outerProduct(this.__wbg_ptr, other.__wbg_ptr);
+    sub(other) {
+        _assertClass(other, WasmMultivector);
+        const ret = wasm.wasmmultivector_sub(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Compute norm (alias for magnitude, maintained for compatibility)
+     * @returns {number}
+     */
+    norm() {
+        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Scale by a scalar
+     * @param {number} scalar
+     * @returns {WasmMultivector}
+     */
+    scale(scalar) {
+        const ret = wasm.wasmmultivector_scale(this.__wbg_ptr, scalar);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Create a scalar multivector
+     * @param {number} value
+     * @returns {WasmMultivector}
+     */
+    static scalar(value) {
+        const ret = wasm.wasmmultivector_scalar(value);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Compute inverse
+     * @returns {WasmMultivector}
+     */
+    inverse() {
+        const ret = wasm.wasmmultivector_inverse(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultivector.__wrap(ret[0]);
+    }
+    /**
+     * Reverse
+     * @returns {WasmMultivector}
+     */
+    reverse() {
+        const ret = wasm.wasmmultivector_reverse(this.__wbg_ptr);
+        return WasmMultivector.__wrap(ret);
+    }
+    /**
+     * Compute magnitude
+     * @returns {number}
+     */
+    magnitude() {
+        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Normalize
+     * @returns {WasmMultivector}
+     */
+    normalize() {
+        const ret = wasm.wasmmultivector_normalize(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmMultivector.__wrap(ret[0]);
+    }
+}
+if (Symbol.dispose) WasmMultivector.prototype[Symbol.dispose] = WasmMultivector.prototype.free;
+
+const WasmRelativisticConstantsFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmrelativisticconstants_free(ptr >>> 0, 1));
+/**
+ * Physical constants
+ */
+export class WasmRelativisticConstants {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmRelativisticConstantsFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmrelativisticconstants_free(ptr, 0);
+    }
+    /**
+     * Earth mass (kg)
+     * @returns {number}
+     */
+    static get earth_mass() {
+        const ret = wasm.wasmrelativisticconstants_earth_mass();
+        return ret;
+    }
+    /**
+     * Solar mass (kg)
+     * @returns {number}
+     */
+    static get solar_mass() {
+        const ret = wasm.wasmrelativisticconstants_solar_mass();
+        return ret;
+    }
+    /**
+     * Speed of light in vacuum (m/s)
+     * @returns {number}
+     */
+    static get speed_of_light() {
+        const ret = wasm.wasmrelativisticconstants_speed_of_light();
+        return ret;
+    }
+    /**
+     * Gravitational constant (m³/kg·s²)
+     * @returns {number}
+     */
+    static get gravitational_constant() {
+        const ret = wasm.wasmrelativisticconstants_gravitational_constant();
+        return ret;
+    }
+}
+if (Symbol.dispose) WasmRelativisticConstants.prototype[Symbol.dispose] = WasmRelativisticConstants.prototype.free;
+
+const WasmRelativisticParticleFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmrelativisticparticle_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for relativistic particles
+ */
+export class WasmRelativisticParticle {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmRelativisticParticle.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmRelativisticParticleFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmRelativisticParticleFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmrelativisticparticle_free(ptr, 0);
+    }
+    /**
+     * Get 3D position components
+     * @returns {Array<any>}
+     */
+    position_3d() {
+        const ret = wasm.wasmrelativisticparticle_position_3d(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get position as spacetime vector
+     * @returns {WasmSpacetimeVector}
+     */
+    position_4d() {
+        const ret = wasm.wasmrelativisticparticle_position_4d(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
+    }
+    /**
+     * Create particle with specified energy
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     * @param {number} direction_x
+     * @param {number} direction_y
+     * @param {number} direction_z
+     * @param {number} kinetic_energy
+     * @param {number} mass
+     * @param {number} charge
+     * @returns {WasmRelativisticParticle}
+     */
+    static with_energy(x, y, z, direction_x, direction_y, direction_z, kinetic_energy, mass, charge) {
+        const ret = wasm.wasmrelativisticparticle_with_energy(x, y, z, direction_x, direction_y, direction_z, kinetic_energy, mass, charge);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return WasmRelativisticParticle.__wrap(ret[0]);
+    }
+    /**
+     * Get total energy
+     * @returns {number}
+     */
+    total_energy() {
+        const ret = wasm.wasmrelativisticparticle_total_energy(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get four-velocity
+     * @returns {WasmFourVelocity}
+     */
+    four_velocity() {
+        const ret = wasm.wasmrelativisticparticle_four_velocity(this.__wbg_ptr);
+        return WasmFourVelocity.__wrap(ret);
+    }
+    /**
+     * Get kinetic energy
+     * @returns {number}
+     */
+    kinetic_energy() {
+        const ret = wasm.wasmrelativisticparticle_kinetic_energy(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get momentum magnitude
+     * @returns {number}
+     */
+    momentum_magnitude() {
+        const ret = wasm.wasmrelativisticparticle_momentum_magnitude(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Create a new relativistic particle
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     * @param {number} vx
+     * @param {number} vy
+     * @param {number} vz
+     * @param {number} spin
+     * @param {number} mass
+     * @param {number} charge
+     */
+    constructor(x, y, z, vx, vy, vz, spin, mass, charge) {
+        const ret = wasm.wasmrelativisticparticle_new(x, y, z, vx, vy, vz, spin, mass, charge);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        this.__wbg_ptr = ret[0] >>> 0;
+        WasmRelativisticParticleFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Get rest mass
+     * @returns {number}
+     */
+    mass() {
+        const ret = wasm.wasmrelativisticparticle_mass(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get electric charge
+     * @returns {number}
+     */
+    charge() {
+        const ret = wasm.wasmrelativisticparticle_charge(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get string representation
+     * @returns {string}
+     */
+    to_string() {
+        let deferred1_0;
+        let deferred1_1;
+        try {
+            const ret = wasm.wasmrelativisticparticle_to_string(this.__wbg_ptr);
+            deferred1_0 = ret[0];
+            deferred1_1 = ret[1];
+            return getStringFromWasm0(ret[0], ret[1]);
+        } finally {
+            wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+        }
+    }
+}
+if (Symbol.dispose) WasmRelativisticParticle.prototype[Symbol.dispose] = WasmRelativisticParticle.prototype.free;
+
+const WasmRotorFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmrotor_free(ptr >>> 0, 1));
+/**
+ * Rotor operations for WASM
+ */
+export class WasmRotor {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmRotor.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmRotorFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmRotorFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmrotor_free(ptr, 0);
+    }
+    /**
+     * Create a rotor from a bivector and angle
+     * @param {WasmMultivector} bivector
+     * @param {number} angle
+     * @returns {WasmRotor}
+     */
+    static fromBivector(bivector, angle) {
+        _assertClass(bivector, WasmMultivector);
+        const ret = wasm.wasmrotor_fromBivector(bivector.__wbg_ptr, angle);
+        return WasmRotor.__wrap(ret);
+    }
+    /**
+     * Apply rotor to a multivector
+     * @param {WasmMultivector} mv
+     * @returns {WasmMultivector}
+     */
+    apply(mv) {
+        _assertClass(mv, WasmMultivector);
+        const ret = wasm.wasmrotor_apply(this.__wbg_ptr, mv.__wbg_ptr);
         return WasmMultivector.__wrap(ret);
     }
     /**
-     * Scalar product
-     * @param {WasmMultivector} other
-     * @returns {number}
+     * Compose two rotors
+     * @param {WasmRotor} other
+     * @returns {WasmRotor}
+     */
+    compose(other) {
+        _assertClass(other, WasmRotor);
+        const ret = wasm.wasmrotor_compose(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmRotor.__wrap(ret);
+    }
+    /**
+     * Get inverse rotor
+     * @returns {WasmRotor}
+     */
+    inverse() {
+        const ret = wasm.wasmmultivector_reverse(this.__wbg_ptr);
+        return WasmRotor.__wrap(ret);
+    }
+}
+if (Symbol.dispose) WasmRotor.prototype[Symbol.dispose] = WasmRotor.prototype.free;
+
+const WasmSchwarzschildMetricFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmschwarzschildmetric_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for Schwarzschild metric
+ */
+export class WasmSchwarzschildMetric {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmSchwarzschildMetric.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmSchwarzschildMetricFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmSchwarzschildMetricFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmschwarzschildmetric_free(ptr, 0);
+    }
+    /**
+     * Check for singularity at given position
+     * @param {WasmSpacetimeVector} position
+     * @returns {boolean}
+     */
+    has_singularity(position) {
+        _assertClass(position, WasmSpacetimeVector);
+        const ret = wasm.wasmschwarzschildmetric_has_singularity(this.__wbg_ptr, position.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Compute effective potential for circular orbits
+     * @param {number} r
+     * @param {number} angular_momentum
+     * @returns {number}
+     */
+    effective_potential(r, angular_momentum) {
+        const ret = wasm.wasmschwarzschildmetric_effective_potential(this.__wbg_ptr, r, angular_momentum);
+        return ret;
+    }
+    /**
+     * Get Schwarzschild radius
+     * @returns {number}
+     */
+    schwarzschild_radius() {
+        const ret = wasm.wasmschwarzschildmetric_schwarzschild_radius(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Create Schwarzschild metric for the Sun
+     * @returns {WasmSchwarzschildMetric}
+     */
+    static sun() {
+        const ret = wasm.wasmschwarzschildmetric_sun();
+        return WasmSchwarzschildMetric.__wrap(ret);
+    }
+    /**
+     * Get central mass
+     * @returns {number}
+     */
+    mass() {
+        const ret = wasm.wasmdualnumber_getReal(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Create Schwarzschild metric for Earth
+     * @returns {WasmSchwarzschildMetric}
+     */
+    static earth() {
+        const ret = wasm.wasmschwarzschildmetric_earth();
+        return WasmSchwarzschildMetric.__wrap(ret);
+    }
+    /**
+     * Create Schwarzschild metric for custom mass
+     * @param {number} mass
+     * @returns {WasmSchwarzschildMetric}
+     */
+    static from_mass(mass) {
+        const ret = wasm.wasmschwarzschildmetric_from_mass(mass);
+        return WasmSchwarzschildMetric.__wrap(ret);
+    }
+}
+if (Symbol.dispose) WasmSchwarzschildMetric.prototype[Symbol.dispose] = WasmSchwarzschildMetric.prototype.free;
+
+const WasmSpacetimeVectorFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmspacetimevector_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for spacetime vectors
+ */
+export class WasmSpacetimeVector {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmSpacetimeVector.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmSpacetimeVectorFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmSpacetimeVectorFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmspacetimevector_free(ptr, 0);
+    }
+    /**
+     * Check if vector is timelike (massive particle)
+     * @returns {boolean}
+     */
+    is_timelike() {
+        const ret = wasm.wasmspacetimevector_is_timelike(this.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Check if vector is spacelike
+     * @returns {boolean}
+     */
+    is_spacelike() {
+        const ret = wasm.wasmspacetimevector_is_spacelike(this.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Compute Minkowski norm squared
+     * @returns {number}
+     */
+    norm_squared() {
+        const ret = wasm.wasmspacetimevector_norm_squared(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Compute Minkowski inner product with another spacetime vector
+     * @param {WasmSpacetimeVector} other
+     * @returns {number}
+     */
+    minkowski_dot(other) {
+        _assertClass(other, WasmSpacetimeVector);
+        const ret = wasm.wasmspacetimevector_minkowski_dot(this.__wbg_ptr, other.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get temporal component
+     * @returns {number}
+     */
+    get t() {
+        const ret = wasm.wasmfourvelocity_gamma(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get x component
+     * @returns {number}
+     */
+    get x() {
+        const ret = wasm.wasmdualnumber_getDual(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get y component
+     * @returns {number}
+     */
+    get y() {
+        const ret = wasm.wasmspacetimevector_y(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Get z component
+     * @returns {number}
+     */
+    get z() {
+        const ret = wasm.wasmspacetimevector_z(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Create a new spacetime vector with components (ct, x, y, z)
+     * @param {number} t
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     */
+    constructor(t, x, y, z) {
+        const ret = wasm.wasmspacetimevector_new(t, x, y, z);
+        this.__wbg_ptr = ret >>> 0;
+        WasmSpacetimeVectorFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Check if vector is null (lightlike)
+     * @returns {boolean}
+     */
+    is_null() {
+        const ret = wasm.wasmspacetimevector_is_null(this.__wbg_ptr);
+        return ret !== 0;
+    }
+    /**
+     * Create a timelike vector
+     * @param {number} t
+     * @returns {WasmSpacetimeVector}
+     */
+    static timelike(t) {
+        const ret = wasm.wasmspacetimevector_timelike(t);
+        return WasmSpacetimeVector.__wrap(ret);
+    }
+    /**
+     * Create a spacelike vector
+     * @param {number} x
+     * @param {number} y
+     * @param {number} z
+     * @returns {WasmSpacetimeVector}
+     */
+    static spacelike(x, y, z) {
+        const ret = wasm.wasmspacetimevector_spacelike(x, y, z);
+        return WasmSpacetimeVector.__wrap(ret);
+    }
+    /**
+     * Get string representation
+     * @returns {string}
+     */
+    to_string() {
+        let deferred1_0;
+        let deferred1_1;
+        try {
+            const ret = wasm.wasmspacetimevector_to_string(this.__wbg_ptr);
+            deferred1_0 = ret[0];
+            deferred1_1 = ret[1];
+            return getStringFromWasm0(ret[0], ret[1]);
+        } finally {
+            wasm.__wbindgen_free(deferred1_0, deferred1_1, 1);
+        }
+    }
+}
+if (Symbol.dispose) WasmSpacetimeVector.prototype[Symbol.dispose] = WasmSpacetimeVector.prototype.free;
+
+const WasmTrajectoryPointFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmtrajectorypoint_free(ptr >>> 0, 1));
+/**
+ * Trajectory point for particle propagation
+ */
+export class WasmTrajectoryPoint {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmTrajectoryPointFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmtrajectorypoint_free(ptr, 0);
+    }
+    /**
+     * Get position
+     * @returns {WasmSpacetimeVector}
      */
-    scalarProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_scalarProduct(this.__wbg_ptr, other.__wbg_ptr);
-        return ret;
+    get position() {
+        const ret = wasm.wasmfourvelocity_as_spacetime_vector(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
     }
     /**
-     * Get a specific coefficient
-     * @param {number} index
+     * Time coordinate
      * @returns {number}
      */
-    getCoefficient(index) {
-        const ret = wasm.wasmmultivector_getCoefficient(this.__wbg_ptr, index);
+    get time() {
+        const ret = wasm.__wbg_get_wasmtrajectorypoint_time(this.__wbg_ptr);
         return ret;
     }
     /**
-     * Set a specific coefficient
-     * @param {number} index
-     * @param {number} value
+     * Time coordinate
+     * @param {number} arg0
      */
-    setCoefficient(index, value) {
-        wasm.wasmmultivector_setCoefficient(this.__wbg_ptr, index, value);
+    set time(arg0) {
+        wasm.__wbg_set_wasmtrajectorypoint_time(this.__wbg_ptr, arg0);
+    }
+}
+if (Symbol.dispose) WasmTrajectoryPoint.prototype[Symbol.dispose] = WasmTrajectoryPoint.prototype.free;
+
+const WasmTropicalNumberFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmtropicalnumber_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for TropicalNumber
+ */
+export class WasmTropicalNumber {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(WasmTropicalNumber.prototype);
+        obj.__wbg_ptr = ptr;
+        WasmTropicalNumberFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmTropicalNumberFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmtropicalnumber_free(ptr, 0);
     }
     /**
-     * Get coefficients as a Float64Array
-     * @returns {Float64Array}
+     * Check if this is infinite
+     * @returns {boolean}
      */
-    getCoefficients() {
-        const ret = wasm.wasmmultivector_getCoefficients(this.__wbg_ptr);
-        var v1 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
-        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
-        return v1;
+    isInfinity() {
+        const ret = wasm.wasmtropicalnumber_isInfinity(this.__wbg_ptr);
+        return ret !== 0;
     }
     /**
-     * Grade projection
-     * @param {number} grade
-     * @returns {WasmMultivector}
+     * Tropical addition (max operation)
+     * @param {WasmTropicalNumber} other
+     * @returns {WasmTropicalNumber}
      */
-    gradeProjection(grade) {
-        const ret = wasm.wasmmultivector_gradeProjection(this.__wbg_ptr, grade);
-        return WasmMultivector.__wrap(ret);
+    tropicalAdd(other) {
+        _assertClass(other, WasmTropicalNumber);
+        const ret = wasm.wasmtropicalnumber_add(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Create from a Float64Array of coefficients
-     * @param {Float64Array} coefficients
-     * @returns {WasmMultivector}
+     * Tropical multiplication (addition)
+     * @param {WasmTropicalNumber} other
+     * @returns {WasmTropicalNumber}
      */
-    static fromCoefficients(coefficients) {
-        const ptr0 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
-        const len0 = WASM_VECTOR_LEN;
-        const ret = wasm.wasmmultivector_fromCoefficients(ptr0, len0);
-        if (ret[2]) {
-            throw takeFromExternrefTable0(ret[1]);
-        }
-        return WasmMultivector.__wrap(ret[0]);
+    tropicalMul(other) {
+        _assertClass(other, WasmTropicalNumber);
+        const ret = wasm.wasmtropicalnumber_mul(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Geometric product
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Tropical power (scalar multiplication)
+     * @param {number} n
+     * @returns {WasmTropicalNumber}
      */
-    geometricProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_geometricProduct(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    tropicalPow(n) {
+        const ret = wasm.wasmtropicalnumber_tropicalPow(this.__wbg_ptr, n);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Add two multivectors
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Create from log probability
+     * @param {number} log_p
+     * @returns {WasmTropicalNumber}
      */
-    add(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_add(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    static fromLogProb(log_p) {
+        const ret = wasm.wasmtropicalnumber_fromLogProb(log_p);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Exponential (for bivectors to create rotors)
-     * @returns {WasmMultivector}
+     * Standard addition (for convenience)
+     * @param {WasmTropicalNumber} other
+     * @returns {WasmTropicalNumber}
      */
-    exp() {
-        const ret = wasm.wasmmultivector_exp(this.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    add(other) {
+        _assertClass(other, WasmTropicalNumber);
+        const ret = wasm.wasmtropicalnumber_add(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Create a new zero multivector
+     * Standard multiplication (for convenience)
+     * @param {WasmTropicalNumber} other
+     * @returns {WasmTropicalNumber}
      */
-    constructor() {
-        const ret = wasm.wasmmultivector_new();
-        this.__wbg_ptr = ret >>> 0;
-        WasmMultivectorFinalization.register(this, this.__wbg_ptr, this);
-        return this;
+    mul(other) {
+        _assertClass(other, WasmTropicalNumber);
+        const ret = wasm.wasmtropicalnumber_mul(this.__wbg_ptr, other.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Subtract two multivectors
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Negation
+     * @returns {WasmTropicalNumber}
      */
-    sub(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_sub(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    neg() {
+        const ret = wasm.wasmtropicalnumber_neg(this.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Compute norm (alias for magnitude, maintained for compatibility)
-     * @returns {number}
+     * Create a new tropical number from a regular number
+     * @param {number} value
      */
-    norm() {
-        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
-        return ret;
+    constructor(value) {
+        const ret = wasm.wasmtropicalnumber_fromLogProb(value);
+        this.__wbg_ptr = ret >>> 0;
+        WasmTropicalNumberFinalization.register(this, this.__wbg_ptr, this);
+        return this;
     }
     /**
-     * Scale by a scalar
-     * @param {number} scalar
-     * @returns {WasmMultivector}
+     * Create tropical one (regular zero)
+     * @returns {WasmTropicalNumber}
      */
-    scale(scalar) {
-        const ret = wasm.wasmmultivector_scale(this.__wbg_ptr, scalar);
-        return WasmMultivector.__wrap(ret);
+    static one() {
+        const ret = wasm.wasmtropicalnumber_one();
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Create a scalar multivector
-     * @param {number} value
-     * @returns {WasmMultivector}
+     * Create tropical zero (negative infinity)
+     * @returns {WasmTropicalNumber}
      */
-    static scalar(value) {
-        const ret = wasm.wasmmultivector_scalar(value);
-        return WasmMultivector.__wrap(ret);
+    static zero() {
+        const ret = wasm.wasmtropicalnumber_zero();
+        return WasmTropicalNumber.__wrap(ret);
     }
     /**
-     * Compute inverse
-     * @returns {WasmMultivector}
+     * Check if this is tropical one (zero)
+     * @returns {boolean}
      */
-    inverse() {
-        const ret = wasm.wasmmultivector_inverse(this.__wbg_ptr);
-        if (ret[2]) {
-            throw takeFromExternrefTable0(ret[1]);
-        }
-        return WasmMultivector.__wrap(ret[0]);
+    isOne() {
+        const ret = wasm.wasmtropicalnumber_isOne(this.__wbg_ptr);
+        return ret !== 0;
     }
     /**
-     * Reverse
-     * @returns {WasmMultivector}
+     * Check if this is tropical zero (negative infinity)
+     * @returns {boolean}
      */
-    reverse() {
-        const ret = wasm.wasmmultivector_reverse(this.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    isZero() {
+        const ret = wasm.wasmtropicalnumber_isZero(this.__wbg_ptr);
+        return ret !== 0;
     }
     /**
-     * Compute magnitude
+     * Convert to probability (via exp)
      * @returns {number}
      */
-    magnitude() {
-        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
+    toProb() {
+        const ret = wasm.wasmtropicalnumber_toProb(this.__wbg_ptr);
         return ret;
     }
     /**
-     * Normalize
-     * @returns {WasmMultivector}
+     * Get the underlying value
+     * @returns {number}
      */
-    normalize() {
-        const ret = wasm.wasmmultivector_normalize(this.__wbg_ptr);
-        if (ret[2]) {
-            throw takeFromExternrefTable0(ret[1]);
-        }
-        return WasmMultivector.__wrap(ret[0]);
+    getValue() {
+        const ret = wasm.wasmdualnumber_getReal(this.__wbg_ptr);
+        return ret;
     }
 }
-if (Symbol.dispose) WasmMultivector.prototype[Symbol.dispose] = WasmMultivector.prototype.free;
+if (Symbol.dispose) WasmTropicalNumber.prototype[Symbol.dispose] = WasmTropicalNumber.prototype.free;
 
-const WasmRotorFinalization = (typeof FinalizationRegistry === 'undefined')
+const WasmTropicalPolynomialFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
-    : new FinalizationRegistry(ptr => wasm.__wbg_wasmrotor_free(ptr >>> 0, 1));
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmtropicalpolynomial_free(ptr >>> 0, 1));
 /**
- * Rotor operations for WASM
+ * WASM wrapper for TropicalPolynomial - polynomial operations in tropical algebra
  */
-export class WasmRotor {
-
-    static __wrap(ptr) {
-        ptr = ptr >>> 0;
-        const obj = Object.create(WasmRotor.prototype);
-        obj.__wbg_ptr = ptr;
-        WasmRotorFinalization.register(obj, obj.__wbg_ptr, obj);
-        return obj;
-    }
+export class WasmTropicalPolynomial {
 
     __destroy_into_raw() {
         const ptr = this.__wbg_ptr;
         this.__wbg_ptr = 0;
-        WasmRotorFinalization.unregister(this);
+        WasmTropicalPolynomialFinalization.unregister(this);
         return ptr;
     }
 
     free() {
         const ptr = this.__destroy_into_raw();
-        wasm.__wbg_wasmrotor_free(ptr, 0);
+        wasm.__wbg_wasmtropicalpolynomial_free(ptr, 0);
     }
     /**
-     * Create a rotor from a bivector and angle
-     * @param {WasmMultivector} bivector
-     * @param {number} angle
-     * @returns {WasmRotor}
+     * Find tropical roots of the polynomial
+     * @returns {Array<any>}
      */
-    static fromBivector(bivector, angle) {
-        _assertClass(bivector, WasmMultivector);
-        const ret = wasm.wasmrotor_fromBivector(bivector.__wbg_ptr, angle);
-        return WasmRotor.__wrap(ret);
+    tropical_roots() {
+        const ret = wasm.wasmtropicalpolynomial_tropical_roots(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Apply rotor to a multivector
-     * @param {WasmMultivector} mv
-     * @returns {WasmMultivector}
+     * Get the number of coefficients
+     * @returns {number}
      */
-    apply(mv) {
-        _assertClass(mv, WasmMultivector);
-        const ret = wasm.wasmrotor_apply(this.__wbg_ptr, mv.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    coefficients_count() {
+        const ret = wasm.wasmtropicalpolynomial_coefficients_count(this.__wbg_ptr);
+        return ret >>> 0;
     }
     /**
-     * Compose two rotors
-     * @param {WasmRotor} other
-     * @returns {WasmRotor}
+     * Create a new tropical polynomial from coefficients
+     * @param {Float64Array} coefficients
      */
-    compose(other) {
-        _assertClass(other, WasmRotor);
-        const ret = wasm.wasmrotor_compose(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmRotor.__wrap(ret);
+    constructor(coefficients) {
+        const ptr0 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmtropicalpolynomial_new(ptr0, len0);
+        this.__wbg_ptr = ret >>> 0;
+        WasmTropicalPolynomialFinalization.register(this, this.__wbg_ptr, this);
+        return this;
     }
     /**
-     * Get inverse rotor
-     * @returns {WasmRotor}
+     * Evaluate the polynomial at a given tropical number
+     * @param {WasmTropicalNumber} x
+     * @returns {WasmTropicalNumber}
      */
-    inverse() {
-        const ret = wasm.wasmmultivector_reverse(this.__wbg_ptr);
-        return WasmRotor.__wrap(ret);
+    evaluate(x) {
+        _assertClass(x, WasmTropicalNumber);
+        const ret = wasm.wasmtropicalpolynomial_evaluate(this.__wbg_ptr, x.__wbg_ptr);
+        return WasmTropicalNumber.__wrap(ret);
     }
 }
-if (Symbol.dispose) WasmRotor.prototype[Symbol.dispose] = WasmRotor.prototype.free;
+if (Symbol.dispose) WasmTropicalPolynomial.prototype[Symbol.dispose] = WasmTropicalPolynomial.prototype.free;
+
+const WasmTropicalViterbiFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmtropicalviterbi_free(ptr >>> 0, 1));
+/**
+ * WASM wrapper for TropicalViterbi - Hidden Markov Model decoding
+ */
+export class WasmTropicalViterbi {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmTropicalViterbiFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmtropicalviterbi_free(ptr, 0);
+    }
+    /**
+     * Compute forward probabilities for all states
+     * @param {Uint32Array} observations
+     * @returns {Array<any>}
+     */
+    forward_probabilities(observations) {
+        const ptr0 = passArray32ToWasm0(observations, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmtropicalviterbi_forward_probabilities(this.__wbg_ptr, ptr0, len0);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Create a new Viterbi decoder
+     *
+     * # Arguments
+     * * `transitions` - Transition probability matrix (2D array)
+     * * `emissions` - Emission probability matrix (2D array)
+     * @param {any} transitions
+     * @param {any} emissions
+     */
+    constructor(transitions, emissions) {
+        const ret = wasm.wasmtropicalviterbi_new(transitions, emissions);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        this.__wbg_ptr = ret[0] >>> 0;
+        WasmTropicalViterbiFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * Decode the most likely state sequence for given observations
+     *
+     * Returns an object with `states` (array of state indices) and `probability` (log probability)
+     * @param {Uint32Array} observations
+     * @returns {any}
+     */
+    decode(observations) {
+        const ptr0 = passArray32ToWasm0(observations, wasm.__wbindgen_malloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.wasmtropicalviterbi_decode(this.__wbg_ptr, ptr0, len0);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+}
+if (Symbol.dispose) WasmTropicalViterbi.prototype[Symbol.dispose] = WasmTropicalViterbi.prototype.free;
 
 const EXPECTED_RESPONSE_TYPES = new Set(['basic', 'cors', 'default']);
 
@@ -554,9 +2777,62 @@ async function __wbg_load(module, imports) {
 function __wbg_get_imports() {
     const imports = {};
     imports.wbg = {};
-    imports.wbg.__wbg_log_fb67a0d730fc146e = function(arg0, arg1) {
+    imports.wbg.__wbg_apply_dc50eb58583d2a57 = function() { return handleError(function (arg0, arg1, arg2) {
+        const ret = arg0.apply(arg1, arg2);
+        return ret;
+    }, arguments) };
+    imports.wbg.__wbg_from_88bc52ce20ba6318 = function(arg0) {
+        const ret = Array.from(arg0);
+        return ret;
+    };
+    imports.wbg.__wbg_get_0da715ceaecea5c8 = function(arg0, arg1) {
+        const ret = arg0[arg1 >>> 0];
+        return ret;
+    };
+    imports.wbg.__wbg_length_186546c51cd61acd = function(arg0) {
+        const ret = arg0.length;
+        return ret;
+    };
+    imports.wbg.__wbg_log_c123554c5b238740 = function(arg0, arg1) {
         console.log(getStringFromWasm0(arg0, arg1));
     };
+    imports.wbg.__wbg_new_19c25a3f2fa63a02 = function() {
+        const ret = new Object();
+        return ret;
+    };
+    imports.wbg.__wbg_new_1f3a344cf3123716 = function() {
+        const ret = new Array();
+        return ret;
+    };
+    imports.wbg.__wbg_push_330b2eb93e4e1212 = function(arg0, arg1) {
+        const ret = arg0.push(arg1);
+        return ret;
+    };
+    imports.wbg.__wbg_set_453345bcda80b89a = function() { return handleError(function (arg0, arg1, arg2) {
+        const ret = Reflect.set(arg0, arg1, arg2);
+        return ret;
+    }, arguments) };
+    imports.wbg.__wbg_wasmspacetimevector_new = function(arg0) {
+        const ret = WasmSpacetimeVector.__wrap(arg0);
+        return ret;
+    };
+    imports.wbg.__wbg_wasmtropicalnumber_new = function(arg0) {
+        const ret = WasmTropicalNumber.__wrap(arg0);
+        return ret;
+    };
+    imports.wbg.__wbg_wbindgendebugstring_99ef257a3ddda34d = function(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);
+    };
+    imports.wbg.__wbg_wbindgennumberget_f74b4c7525ac05cb = function(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);
+    };
     imports.wbg.__wbg_wbindgenthrow_451ec1a8469d7eb6 = function(arg0, arg1) {
         throw new Error(getStringFromWasm0(arg0, arg1));
     };
@@ -565,8 +2841,13 @@ function __wbg_get_imports() {
         const ret = getStringFromWasm0(arg0, arg1);
         return ret;
     };
+    imports.wbg.__wbindgen_cast_d6cd19b81560fd6e = function(arg0) {
+        // Cast intrinsic for `F64 -> Externref`.
+        const ret = arg0;
+        return ret;
+    };
     imports.wbg.__wbindgen_init_externref_table = function() {
-        const table = wasm.__wbindgen_export_0;
+        const table = wasm.__wbindgen_export_2;
         const offset = table.grow(4);
         table.set(0, undefined);
         table.set(offset + 0, undefined);
@@ -586,7 +2867,9 @@ function __wbg_init_memory(imports, memory) {
 function __wbg_finalize_init(instance, module) {
     wasm = instance.exports;
     __wbg_init.__wbindgen_wasm_module = module;
+    cachedDataViewMemory0 = null;
     cachedFloat64ArrayMemory0 = null;
+    cachedUint32ArrayMemory0 = null;
     cachedUint8ArrayMemory0 = null;