@justinelliottcobb/amari-wasm 0.9.0 → 0.9.2

package/amari_wasm.js

@@ -30,6 +30,206 @@ function getStringFromWasm0(ptr, len) {
     return decodeText(ptr, len);
 }
 
+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);
+    }
+}
+
+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 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 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];
+}
+
+/**
+ * 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];
+}
+
+/**
+ * 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 +239,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 +246,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 +257,20 @@ 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);
+}
+
 const BatchOperationsFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
     : new FinalizationRegistry(ptr => wasm.__wbg_batchoperations_free(ptr >>> 0, 1));
@@ -132,6 +332,64 @@ export class BatchOperations {
 }
 if (Symbol.dispose) BatchOperations.prototype[Symbol.dispose] = BatchOperations.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));
@@ -152,175 +410,464 @@ export class PerformanceOperations {
         wasm.__wbg_performanceoperations_free(ptr, 0);
     }
     /**
-     * Fast geometric product for hot paths with memory pooling
-     * @param {Float64Array} lhs
-     * @param {Float64Array} rhs
+     * Batch normalize vectors for efficiency
+     * @param {Float64Array} vectors
+     * @param {number} vector_size
      * @returns {Float64Array}
      */
-    static fastGeometricProduct(lhs, rhs) {
-        const ptr0 = passArrayF64ToWasm0(lhs, wasm.__wbindgen_malloc);
+    static batchNormalize(vectors, vector_size) {
+        const ptr0 = passArrayF64ToWasm0(vectors, 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();
+        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 v3;
+        return v2;
     }
     /**
-     * Optimized vector operations for 3D space
+     * Optimized vector dot product
      * @param {Float64Array} v1
      * @param {Float64Array} v2
-     * @returns {Float64Array}
+     * @returns {number}
      */
-    static vectorCrossProduct(v1, v2) {
+    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_vectorCrossProduct(ptr0, len0, ptr1, len1);
-        var v3 = getArrayF64FromWasm0(ret[0], ret[1]).slice();
-        wasm.__wbindgen_free(ret[0], ret[1] * 8, 8);
-        return v3;
+        const ret = wasm.performanceoperations_vectorDotProduct(ptr0, len0, ptr1, len1);
+        return ret;
     }
     /**
-     * Optimized vector dot product
+     * Optimized vector operations for 3D space
      * @param {Float64Array} v1
      * @param {Float64Array} v2
-     * @returns {number}
+     * @returns {Float64Array}
      */
-    static vectorDotProduct(v1, v2) {
+    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_vectorDotProduct(ptr0, len0, ptr1, len1);
-        return ret;
+        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;
     }
     /**
-     * Batch normalize vectors for efficiency
-     * @param {Float64Array} vectors
-     * @param {number} vector_size
+     * Fast geometric product for hot paths with memory pooling
+     * @param {Float64Array} lhs
+     * @param {Float64Array} rhs
      * @returns {Float64Array}
      */
-    static batchNormalize(vectors, vector_size) {
-        const ptr0 = passArrayF64ToWasm0(vectors, wasm.__wbindgen_malloc);
+    static fastGeometricProduct(lhs, rhs) {
+        const ptr0 = passArrayF64ToWasm0(lhs, 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();
+        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 v2;
+        return v3;
     }
 }
 if (Symbol.dispose) PerformanceOperations.prototype[Symbol.dispose] = PerformanceOperations.prototype.free;
 
-const WasmMultivectorFinalization = (typeof FinalizationRegistry === 'undefined')
+const WasmFourVelocityFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
-    : new FinalizationRegistry(ptr => wasm.__wbg_wasmmultivector_free(ptr >>> 0, 1));
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmfourvelocity_free(ptr >>> 0, 1));
 /**
- * WASM wrapper for Multivector with TypedArray support
+ * WASM wrapper for four-velocity
  */
-export class WasmMultivector {
+export class WasmFourVelocity {
 
     static __wrap(ptr) {
         ptr = ptr >>> 0;
-        const obj = Object.create(WasmMultivector.prototype);
+        const obj = Object.create(WasmFourVelocity.prototype);
         obj.__wbg_ptr = ptr;
-        WasmMultivectorFinalization.register(obj, obj.__wbg_ptr, obj);
+        WasmFourVelocityFinalization.register(obj, obj.__wbg_ptr, obj);
         return obj;
     }
 
     __destroy_into_raw() {
         const ptr = this.__wbg_ptr;
         this.__wbg_ptr = 0;
-        WasmMultivectorFinalization.unregister(this);
+        WasmFourVelocityFinalization.unregister(this);
         return ptr;
     }
 
     free() {
         const ptr = this.__destroy_into_raw();
-        wasm.__wbg_wasmmultivector_free(ptr, 0);
-    }
-    /**
-     * 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;
+        wasm.__wbg_wasmfourvelocity_free(ptr, 0);
     }
     /**
-     * Create from a Float64Array of coefficients
-     * @param {Float64Array} coefficients
-     * @returns {WasmMultivector}
+     * Create four-velocity from 3-velocity components
+     * @param {number} vx
+     * @param {number} vy
+     * @param {number} vz
+     * @returns {WasmFourVelocity}
      */
-    static fromCoefficients(coefficients) {
-        const ptr0 = passArrayF64ToWasm0(coefficients, wasm.__wbindgen_malloc);
-        const len0 = WASM_VECTOR_LEN;
-        const ret = wasm.wasmmultivector_fromCoefficients(ptr0, len0);
+    static from_velocity(vx, vy, vz) {
+        const ret = wasm.wasmfourvelocity_from_velocity(vx, vy, vz);
         if (ret[2]) {
             throw takeFromExternrefTable0(ret[1]);
         }
-        return WasmMultivector.__wrap(ret[0]);
+        return WasmFourVelocity.__wrap(ret[0]);
     }
     /**
-     * Create a scalar multivector
-     * @param {number} value
-     * @returns {WasmMultivector}
+     * Check if normalized (u·u = c²)
+     * @returns {boolean}
      */
-    static scalar(value) {
-        const ret = wasm.wasmmultivector_scalar(value);
-        return WasmMultivector.__wrap(ret);
+    is_normalized() {
+        const ret = wasm.wasmfourvelocity_is_normalized(this.__wbg_ptr);
+        return ret !== 0;
     }
     /**
-     * Create a basis vector (0-indexed)
-     * @param {number} index
-     * @returns {WasmMultivector}
+     * Get as spacetime vector
+     * @returns {WasmSpacetimeVector}
      */
-    static basisVector(index) {
-        const ret = wasm.wasmmultivector_basisVector(index);
-        if (ret[2]) {
-            throw takeFromExternrefTable0(ret[1]);
-        }
-        return WasmMultivector.__wrap(ret[0]);
+    as_spacetime_vector() {
+        const ret = wasm.wasmfourvelocity_as_spacetime_vector(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
     }
     /**
-     * Get coefficients as a Float64Array
-     * @returns {Float64Array}
+     * Get spatial velocity magnitude
+     * @returns {number}
      */
-    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;
+    spatial_velocity_magnitude() {
+        const ret = wasm.wasmfourvelocity_spatial_velocity_magnitude(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Get a specific coefficient
-     * @param {number} index
+     * Get Lorentz factor γ
      * @returns {number}
      */
-    getCoefficient(index) {
-        const ret = wasm.wasmmultivector_getCoefficient(this.__wbg_ptr, index);
+    gamma() {
+        const ret = wasm.wasmfourvelocity_gamma(this.__wbg_ptr);
         return ret;
     }
     /**
-     * Set a specific coefficient
-     * @param {number} index
-     * @param {number} value
+     * Get rapidity
+     * @returns {number}
      */
-    setCoefficient(index, value) {
-        wasm.wasmmultivector_setCoefficient(this.__wbg_ptr, index, value);
+    rapidity() {
+        const ret = wasm.wasmfourvelocity_rapidity(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Geometric product
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Get string representation
+     * @returns {string}
      */
-    geometricProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_geometricProduct(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    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 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)
@@ -332,120 +879,394 @@ export class 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);
+    /**
+     * 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}
+     */
+    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);
     }
     /**
-     * Scalar product
-     * @param {WasmMultivector} other
-     * @returns {number}
+     * Get 3D position components
+     * @returns {Array<any>}
      */
-    scalarProduct(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_scalarProduct(this.__wbg_ptr, other.__wbg_ptr);
+    position_3d() {
+        const ret = wasm.wasmrelativisticparticle_position_3d(this.__wbg_ptr);
         return ret;
     }
     /**
-     * Reverse
-     * @returns {WasmMultivector}
+     * Get position as spacetime vector
+     * @returns {WasmSpacetimeVector}
      */
-    reverse() {
-        const ret = wasm.wasmmultivector_reverse(this.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    position_4d() {
+        const ret = wasm.wasmrelativisticparticle_position_4d(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
     }
     /**
-     * Grade projection
-     * @param {number} grade
-     * @returns {WasmMultivector}
+     * 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}
      */
-    gradeProjection(grade) {
-        const ret = wasm.wasmmultivector_gradeProjection(this.__wbg_ptr, grade);
-        return WasmMultivector.__wrap(ret);
+    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]);
     }
     /**
-     * Exponential (for bivectors to create rotors)
-     * @returns {WasmMultivector}
+     * Get total energy
+     * @returns {number}
      */
-    exp() {
-        const ret = wasm.wasmmultivector_exp(this.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    total_energy() {
+        const ret = wasm.wasmrelativisticparticle_total_energy(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Compute magnitude
-     * @returns {number}
+     * Get four-velocity
+     * @returns {WasmFourVelocity}
      */
-    magnitude() {
-        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
-        return ret;
+    four_velocity() {
+        const ret = wasm.wasmrelativisticparticle_four_velocity(this.__wbg_ptr);
+        return WasmFourVelocity.__wrap(ret);
     }
     /**
-     * Compute norm (alias for magnitude, maintained for compatibility)
+     * Get kinetic energy
      * @returns {number}
      */
-    norm() {
-        const ret = wasm.wasmmultivector_magnitude(this.__wbg_ptr);
+    kinetic_energy() {
+        const ret = wasm.wasmrelativisticparticle_kinetic_energy(this.__wbg_ptr);
         return ret;
     }
     /**
-     * Normalize
-     * @returns {WasmMultivector}
+     * Get momentum magnitude
+     * @returns {number}
      */
-    normalize() {
-        const ret = wasm.wasmmultivector_normalize(this.__wbg_ptr);
-        if (ret[2]) {
-            throw takeFromExternrefTable0(ret[1]);
-        }
-        return WasmMultivector.__wrap(ret[0]);
+    momentum_magnitude() {
+        const ret = wasm.wasmrelativisticparticle_momentum_magnitude(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Compute inverse
-     * @returns {WasmMultivector}
+     * 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
      */
-    inverse() {
-        const ret = wasm.wasmmultivector_inverse(this.__wbg_ptr);
+    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]);
         }
-        return WasmMultivector.__wrap(ret[0]);
+        this.__wbg_ptr = ret[0] >>> 0;
+        WasmRelativisticParticleFinalization.register(this, this.__wbg_ptr, this);
+        return this;
     }
     /**
-     * Add two multivectors
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Get rest mass
+     * @returns {number}
      */
-    add(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_add(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    mass() {
+        const ret = wasm.wasmrelativisticparticle_mass(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Subtract two multivectors
-     * @param {WasmMultivector} other
-     * @returns {WasmMultivector}
+     * Get electric charge
+     * @returns {number}
      */
-    sub(other) {
-        _assertClass(other, WasmMultivector);
-        const ret = wasm.wasmmultivector_sub(this.__wbg_ptr, other.__wbg_ptr);
-        return WasmMultivector.__wrap(ret);
+    charge() {
+        const ret = wasm.wasmrelativisticparticle_charge(this.__wbg_ptr);
+        return ret;
     }
     /**
-     * Scale by a scalar
-     * @param {number} scalar
-     * @returns {WasmMultivector}
+     * Get string representation
+     * @returns {string}
      */
-    scale(scalar) {
-        const ret = wasm.wasmmultivector_scale(this.__wbg_ptr, scalar);
-        return WasmMultivector.__wrap(ret);
+    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) WasmMultivector.prototype[Symbol.dispose] = WasmMultivector.prototype.free;
+if (Symbol.dispose) WasmRelativisticParticle.prototype[Symbol.dispose] = WasmRelativisticParticle.prototype.free;
 
 const WasmRotorFinalization = (typeof FinalizationRegistry === 'undefined')
     ? { register: () => {}, unregister: () => {} }
@@ -516,6 +1337,295 @@ export class WasmRotor {
 }
 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.wasmschwarzschildmetric_mass(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.wasmspacetimevector_x(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}
+     */
+    get position() {
+        const ret = wasm.wasmfourvelocity_as_spacetime_vector(this.__wbg_ptr);
+        return WasmSpacetimeVector.__wrap(ret);
+    }
+    /**
+     * Time coordinate
+     * @returns {number}
+     */
+    get time() {
+        const ret = wasm.__wbg_get_wasmtrajectorypoint_time(this.__wbg_ptr);
+        return ret;
+    }
+    /**
+     * Time coordinate
+     * @param {number} arg0
+     */
+    set time(arg0) {
+        wasm.__wbg_set_wasmtrajectorypoint_time(this.__wbg_ptr, arg0);
+    }
+}
+if (Symbol.dispose) WasmTrajectoryPoint.prototype[Symbol.dispose] = WasmTrajectoryPoint.prototype.free;
+
 const EXPECTED_RESPONSE_TYPES = new Set(['basic', 'cors', 'default']);
 
 async function __wbg_load(module, imports) {
@@ -554,9 +1664,36 @@ async function __wbg_load(module, imports) {
 function __wbg_get_imports() {
     const imports = {};
     imports.wbg = {};
-    imports.wbg.__wbg_log_6acc76af826c364a = function(arg0, arg1) {
+    imports.wbg.__wbg_log_92e4f1a348990392 = 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_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_wbindgenthrow_451ec1a8469d7eb6 = function(arg0, arg1) {
         throw new Error(getStringFromWasm0(arg0, arg1));
     };
@@ -565,8 +1702,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 +1728,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;