ballistics-engine 0.13.3

package/LICENSE

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

package/README.md

@@ -0,0 +1,614 @@
+# Ballistics Engine
+
+A high-performance ballistics trajectory calculation engine with comprehensive physics modeling, automatic zeroing, and statistical analysis capabilities.
+
+**Project Website:** [https://ballistics.rs/](https://ballistics.rs/)
+
+## Features
+
+- **Full 3D Trajectory Integration** - Six-state ballistic modeling with adaptive RK45 and fixed-step RK4 integration methods
+- **Advanced Drag Models** - Support for G1, G7, and custom drag curves with automatic transonic corrections
+- **Automatic Zeroing** - Calculate sight adjustments and apply zero angles automatically
+- **Unit Conversion** - Seamless switching between Imperial (default) and Metric units
+- **BC Segmentation** - Velocity-dependent ballistic coefficient modeling with automatic estimation
+- **Atmospheric Modeling** - Temperature, pressure, humidity, and altitude effects with ICAO standard atmosphere
+- **Wind Effects** - 3D wind calculations with altitude-dependent wind shear modeling
+- **Monte Carlo Simulations** - Statistical analysis with parameter uncertainties
+- **BC Estimation** - Estimate ballistic coefficients from trajectory data
+- **Advanced Physics**:
+  - **Spin Effects**: Magnus effect, enhanced spin drift with decay modeling
+  - **Earth Effects**: Coriolis effect with latitude-dependent calculations
+  - **Angular Motion**: Gyroscopic precession and nutation physics
+  - **Transonic Analysis**: Pitch damping coefficients and stability warnings
+  - **Trajectory Sampling**: Regular interval data collection for analysis
+  - **Form Factor Corrections**: Bullet-specific drag adjustments
+- **Multiple Output Formats** - JSON, CSV, and formatted tables
+
+## Installation
+
+### From Source
+
+```bash
+git clone https://github.com/ajokela/ballistics-engine.git
+cd ballistics-engine
+cargo build --release
+```
+
+The binary will be at: `target/release/ballistics`
+
+## Quick Start
+
+### Basic Trajectory (Imperial Units - Default)
+
+```bash
+# .308 Winchester, 168gr bullet at 2700 fps
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 --max-range 1000
+
+# With automatic zeroing at 200 yards
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 --auto-zero 200 --max-range 500
+```
+
+### Metric Units
+
+```bash
+# Same bullet in metric units
+./ballistics trajectory --units metric -v 823 -b 0.475 -m 10.9 -d 7.82 --max-range 1000
+```
+
+## Unit Systems
+
+The engine supports two unit systems, selectable with the `--units` flag:
+
+### Imperial (Default)
+- **Velocity**: feet per second (fps)
+- **Mass**: grains
+- **Distance**: yards
+- **Diameter**: inches
+- **Temperature**: Fahrenheit
+- **Pressure**: inHg
+- **Wind**: mph
+
+### Metric
+- **Velocity**: meters per second (m/s)
+- **Mass**: grams
+- **Distance**: meters
+- **Diameter**: millimeters
+- **Temperature**: Celsius
+- **Pressure**: hPa (millibars)
+- **Wind**: m/s
+
+## Commands
+
+### Trajectory Calculation
+
+Calculate ballistic trajectory with environmental conditions:
+
+```bash
+# Imperial units (default)
+./ballistics trajectory \
+  -v 2700          # Velocity (fps)
+  -b 0.475         # Ballistic coefficient
+  -m 168           # Mass (grains)
+  -d 0.308         # Diameter (inches)
+  --drag-model g7  # G7 drag model
+  --angle 0        # Launch angle (degrees)
+  --max-range 1000 # Maximum range (yards)
+  --wind-speed 10  # Wind speed (mph)
+  --wind-direction 90  # Wind from right (degrees)
+  --temperature 59     # Temperature (Fahrenheit)
+  --pressure 29.92    # Pressure (inHg)
+  --humidity 50       # Relative humidity (%)
+  --altitude 0        # Altitude (feet)
+  --full             # Show all trajectory points
+```
+
+#### Auto-Zero Feature
+
+Automatically calculate and apply the zero angle for a specific distance:
+
+```bash
+# Zero at 200 yards and show trajectory to 500 yards
+./ballistics trajectory \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --auto-zero 200 \  # Automatically zero at 200 yards
+  --max-range 500 \
+  --full
+
+# Custom sight height for auto-zero
+./ballistics trajectory \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --auto-zero 100 \
+  --sight-height 0.055  # 2.2 inches in yards
+```
+
+#### Advanced BC Modeling
+
+Enable velocity-dependent BC modeling for more accurate long-range predictions:
+
+```bash
+# Enable BC segmentation (velocity-based BC changes)
+./ballistics trajectory \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --use-bc-segments \
+  --auto-zero 600 \
+  --max-range 1000
+```
+
+#### Advanced Physics - Magnus and Spin Drift
+
+Enable advanced gyroscopic and aerodynamic effects:
+
+```bash
+# Magnus effect and spin drift calculation
+./ballistics trajectory \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --twist-rate 10      # 1:10" barrel twist
+  --twist-right        # Right-hand twist
+  --enable-magnus      # Enable Magnus effect
+  --enable-spin-drift  # Enable enhanced spin drift
+  --wind-speed 10 \
+  --wind-direction 90 \
+  --max-range 1000
+
+# Coriolis effect for extreme long range
+./ballistics trajectory \
+  -v 3000 -b 0.750 -m 250 -d 0.338 \
+  --enable-coriolis \
+  --latitude 45        # Shooting latitude
+  --shooting-angle 90  # Azimuth (0=N, 90=E)
+  --max-range 2000
+```
+
+### Zero Calculation
+
+Calculate the sight adjustment needed to zero at a specific distance:
+
+```bash
+# Calculate zero for 200 yards
+./ballistics zero \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --target-distance 200
+
+# With custom sight height (default is 0.05 yards / 1.8 inches)
+./ballistics zero \
+  -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --target-distance 300 \
+  --sight-height 0.055  # 2.2 inches
+
+# Metric example
+./ballistics zero --units metric \
+  -v 823 -b 0.475 -m 10.9 -d 7.82 \
+  --target-distance 200  # 200 meters
+```
+
+Output includes:
+- Zero angle in degrees
+- Adjustment in MOA (Minutes of Angle)
+- Adjustment in mrad (milliradians)
+- Maximum ordinate (highest point of trajectory)
+
+### Monte Carlo Simulation
+
+Run statistical analysis with parameter variations:
+
+```bash
+./ballistics monte-carlo \
+  -v 2700         # Base velocity (fps)
+  -b 0.475        # Base BC
+  -m 168          # Mass (grains)
+  -d 0.308        # Diameter (inches)
+  -n 1000         # Number of simulations
+  --velocity-std 10    # Velocity std dev (fps)
+  --angle-std 0.5     # Angle std dev (degrees)
+  --bc-std 0.01       # BC std dev
+  --wind-std 2        # Wind speed std dev (mph)
+  --target-distance 300  # Target distance for hit probability
+```
+
+### BC Estimation
+
+Estimate ballistic coefficient from observed trajectory data:
+
+```bash
+./ballistics estimate-bc \
+  -v 2700 -m 168 -d 0.308 \
+  --distance1 100 --drop1 0.0   # First data point
+  --distance2 200 --drop2 0.023  # Second data point
+```
+
+## Advanced Features
+
+### Integration Methods
+
+The engine supports two numerical integration methods:
+
+- **RK45 (Dormand-Prince Adaptive)** - Default method, provides best accuracy with adaptive step sizing
+- **RK4 (Runge-Kutta 4th Order Fixed-Step)** - Available with `--use-rk4-fixed` flag for faster computation
+
+### Wind Shear Modeling
+
+Model altitude-dependent wind variations:
+
+```bash
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --wind-speed 10 --wind-direction 90 \
+  --enable-wind-shear \
+  --max-range 1000
+```
+
+### Transonic Stability Analysis
+
+Analyze projectile stability through the transonic regime:
+
+```bash
+./ballistics trajectory -v 3000 -b 0.475 -m 168 -d 0.308 \
+  --enable-pitch-damping \
+  --max-range 2000
+```
+
+Provides warnings about transonic instability and minimum pitch damping coefficients.
+
+### Trajectory Sampling
+
+Collect trajectory data at regular intervals for detailed analysis:
+
+```bash
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --sample-trajectory \
+  --sample-interval 25  # Sample every 25 yards
+  --max-range 1000 -o json
+```
+
+### Angular Motion Physics
+
+Model precession and nutation of spinning projectiles:
+
+```bash
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 \
+  --twist-rate 10 \
+  --enable-precession \
+  --max-range 1000
+```
+
+### Complete Advanced Physics Example
+
+```bash
+./ballistics trajectory \
+  -v 2850 -b 0.690 -m 230 -d 0.338 \
+  --drag-model g7 \
+  --twist-rate 8.5 --twist-right \
+  --enable-magnus \
+  --enable-coriolis \
+  --enable-spin-drift \
+  --enable-wind-shear \
+  --enable-pitch-damping \
+  --enable-precession \
+  --sample-trajectory \
+  --latitude 38.5 \
+  --shooting-angle 45 \
+  --wind-speed 15 --wind-direction 270 \
+  --altitude 6000 \
+  --max-range 2000
+```
+
+## Physics Modeling
+
+The ballistics engine implements comprehensive physics modeling for accurate trajectory prediction:
+
+### Aerodynamic Effects
+- **Drag Modeling** - Multiple drag functions (G1-G8, JBM, custom curves) with transonic flow corrections
+- **Form Factor** - Projectile efficiency corrections based on shape and design
+- **Reynolds Number Effects** - Viscosity and flow regime corrections
+
+### Gyroscopic Effects  
+- **Spin Drift** - Lateral deviation due to gyroscopic and Magnus effects
+- **Precession** - Gyroscopic precession of spinning projectile
+- **Nutation** - Oscillatory motion superimposed on precession
+- **Spin Decay** - Reduction in spin rate over time due to aerodynamic damping
+- **Pitch Damping** - Aerodynamic moments opposing angular motion
+
+### Environmental Effects
+- **Coriolis Effect** - Earth's rotation influence on long-range trajectories
+- **Magnus Effect** - Force from spinning projectile in crossflow
+- **Wind Shear** - Altitude-dependent wind variations
+- **Atmospheric Stratification** - Density and sound speed variations with altitude
+
+### Stability Modeling
+- **Dynamic Stability** - Gyroscopic and aerodynamic stability calculations
+- **Yaw of Repose** - Equilibrium yaw angle in crosswind
+- **Limit Cycle Yaw** - Bounded oscillatory motion analysis
+
+## FFI Layer
+
+The library includes a Foreign Function Interface (FFI) layer for integration with iOS, Android, and other platforms. The FFI provides C-compatible bindings for all major functionality.
+
+<img src="ios.png" alt="iOS Integration Example" width="35%">
+
+### FFI Features
+- **C-Compatible Structures** - All data structures use C-compatible layouts
+- **Safe Memory Management** - Proper handling of memory across language boundaries
+- **iOS/Swift Integration** - Ready for use with Swift through bridging headers
+- **Android/JNI Support** - Compatible with Java Native Interface
+- **Monte Carlo Simulation** - Statistical analysis with parameter variations
+- **Error Handling** - Graceful error propagation across FFI boundary
+
+### Example FFI Usage (C/Swift)
+```c
+// Create input parameters
+FFIBallisticInputs inputs = {
+    .muzzle_velocity = 823.0,        // m/s
+    .ballistic_coefficient = 0.475,
+    .mass = 0.0109,                  // kg
+    .diameter = 0.00782,             // meters
+    .drag_model = 0,                 // G1
+    .sight_height = 0.05,            // meters
+    .temperature = 15.0,             // Celsius
+    .altitude = 0.0
+};
+
+// Calculate trajectory
+FFITrajectoryResult* result = ballistics_calculate_trajectory(&inputs, NULL, NULL, 1000.0, 0.1);
+
+// Use results
+printf("Max range: %.2f meters\n", result->max_range);
+
+// Clean up
+ballistics_free_trajectory_result(result);
+```
+
+### Monte Carlo Simulation via FFI
+```c
+// Set up Monte Carlo parameters
+FFIMonteCarloParams params = {
+    .num_simulations = 1000,
+    .velocity_std_dev = 10.0,       // m/s variation
+    .angle_std_dev = 0.001,         // radian variation (elevation)
+    .bc_std_dev = 0.01,             // BC variation
+    .wind_speed_std_dev = 2.0,      // m/s wind variation
+    .target_distance = 600.0,       // Target at 600m
+    .azimuth_std_dev = 0.001        // radian variation (horizontal)
+};
+
+// Run simulation
+FFIMonteCarloResults* results = ballistics_monte_carlo(&inputs, NULL, &params);
+
+// Use statistical results
+printf("Mean range: %.2f m (σ=%.2f)\n", results->mean_range, results->std_dev_range);
+printf("Hit probability at 600m: %.1f%%\n", results->hit_probability * 100);
+
+// Access individual shots
+for (int i = 0; i < results->num_results; i++) {
+    printf("Shot %d: Range %.2f m, Impact velocity %.2f m/s\n", 
+           i, results->ranges[i], results->impact_velocities[i]);
+}
+
+// Clean up
+ballistics_free_monte_carlo_results(results);
+```
+
+## Output Formats
+
+All commands support three output formats via the `-o` flag:
+
+- **table** (default) - Formatted ASCII table for terminal display
+- **json** - Complete data in JSON format for programmatic use
+- **csv** - Comma-separated values for spreadsheet analysis
+
+## Practical Examples
+
+### Hunting Zero
+
+Zero a hunting rifle at 200 yards with environmental conditions:
+
+```bash
+# Calculate zero angle
+./ballistics zero \
+  -v 2650 -b 0.460 -m 180 -d 0.308 \
+  --target-distance 200
+
+# Verify trajectory with auto-zero
+./ballistics trajectory \
+  -v 2650 -b 0.460 -m 180 -d 0.308 \
+  --auto-zero 200 \
+  --max-range 400 \
+  --wind-speed 15 \
+  --wind-direction 270 \
+  --temperature 32 \
+  --humidity 30 \
+  --altitude 5000 \
+  --full
+```
+
+### Long Range Shooting
+
+Analyze trajectory for 1000-yard shot:
+
+```bash
+./ballistics trajectory \
+  -v 2850 -b 0.690 -m 230 -d 0.338 \
+  --drag-model g7 \
+  --auto-zero 100 \
+  --max-range 1100 \
+  --wind-speed 10 \
+  --wind-direction 45 \
+  --full \
+  -o json > trajectory.json
+```
+
+### Load Development
+
+Compare different loads using Monte Carlo:
+
+```bash
+# Load 1: Higher velocity, more variation
+./ballistics monte-carlo \
+  -v 2750 -b 0.475 -m 168 -d 0.308 \
+  -n 1000 \
+  --velocity-std 15 \
+  --target-distance 600
+
+# Load 2: Lower velocity, more consistent
+./ballistics monte-carlo \
+  -v 2680 -b 0.475 -m 168 -d 0.308 \
+  -n 1000 \
+  --velocity-std 8 \
+  --target-distance 600
+```
+
+## Advanced Features
+
+### BC Segmentation
+
+Velocity-dependent BC modeling accounts for how ballistic coefficient changes as the bullet slows down. Enable with `--use-bc-segments`:
+
+- Automatically estimates BC segments based on bullet characteristics
+- No external data required - uses caliber, weight, and BC
+- Identifies bullet type (Match, Hunting, VLD, etc.) from parameters
+- Applies physics-based BC degradation curves
+
+Example:
+```bash
+./ballistics trajectory -v 2700 -b 0.475 -m 168 -d 0.308 --use-bc-segments --max-range 1000
+```
+
+
+### Advanced Physics Modeling
+
+When enabled, the engine calculates:
+- **Magnus Effect** - Side force from spinning projectiles
+- **Spin Drift** - Lateral drift due to gyroscopic effects
+- **Coriolis Effect** - Earth rotation effects (with latitude input)
+- **Transonic Drag** - Enhanced drag modeling in transonic regime
+- **Reynolds Number Corrections** - Viscous effects at low velocities
+
+## Building from Source
+
+### Requirements
+
+- Rust 1.70 or later
+- Cargo build system
+
+### Build Commands
+
+```bash
+# Debug build
+cargo build
+
+# Release build (optimized)
+cargo build --release
+
+# Run tests
+cargo test
+
+# Build documentation
+cargo doc --open
+```
+
+## Library Usage
+
+Use as a Rust library in your own projects:
+
+```rust
+use ballistics_engine::{
+    BallisticInputs, TrajectorySolver, 
+    WindConditions, AtmosphericConditions
+};
+
+fn main() -> Result<(), Box<dyn std::error::Error>> {
+    let inputs = BallisticInputs {
+        muzzle_velocity: 823.0,  // m/s
+        launch_angle: 0.0,       // radians
+        ballistic_coefficient: 0.475,
+        mass: 0.0109,            // kg
+        diameter: 0.00782,       // meters
+        sight_height: 0.05,      // meters
+        ..Default::default()
+    };
+    
+    let wind = WindConditions {
+        speed: 5.0,              // m/s
+        direction: 1.5708,       // 90 degrees in radians
+        ..Default::default()
+    };
+    
+    let atmosphere = AtmosphericConditions {
+        temperature: 15.0,       // Celsius
+        pressure: 1013.25,       // hPa
+        humidity: 50.0,          // %
+        altitude: 0.0,           // meters
+        ..Default::default()
+    };
+    
+    let solver = TrajectorySolver::new(inputs, wind, atmosphere);
+    let result = solver.solve()?;
+    
+    println!("Max range: {:.2} m", result.max_range);
+    println!("Max height: {:.2} m", result.max_height);
+    println!("Time of flight: {:.3} s", result.time_of_flight);
+    
+    Ok(())
+}
+```
+
+## Performance
+
+Optimized Rust implementation provides:
+- Single trajectory (1000m): ~5ms
+- Monte Carlo (1000 runs): ~500ms  
+- BC estimation: ~50ms
+- Zero calculation: ~10ms
+
+## Common Ballistic Coefficients
+
+| Caliber | Weight | BC (G1) | BC (G7) | Description |
+|---------|--------|---------|---------|-------------|
+| .223    | 55gr   | 0.250   | -       | FMJ |
+| .223    | 77gr   | 0.362   | 0.182   | Match |
+| .308    | 168gr  | 0.475   | 0.224   | Match |
+| .308    | 175gr  | 0.505   | 0.253   | Match |
+| .308    | 180gr  | 0.480   | -       | Hunting |
+| .338    | 300gr  | 0.768   | 0.383   | Match |
+| 6.5mm   | 140gr  | 0.620   | 0.310   | Match |
+| .50     | 750gr  | 1.050   | 0.520   | Match |
+
+## Troubleshooting
+
+### Trajectory hits ground early
+- Check if you're using `--auto-zero` or setting `--angle` manually
+- Default angle is 0° (horizontal), which will hit ground quickly
+- Use `--auto-zero <distance>` to automatically calculate proper angle
+
+### Units confusion
+- Default is Imperial (fps, grains, yards)
+- Use `--units metric` for metric system
+- All inputs must match the selected unit system
+
+### Unexpected BC behavior
+- G1 and G7 models have different BC values for same bullet
+- G7 typically better for boat-tail bullets
+- BC segmentation automatically applied based on bullet type
+
+## Contributing
+
+Contributions are welcome! Please:
+1. Fork the repository
+2. Create a feature branch
+3. Add tests for new features
+4. Run `cargo test` and `cargo fmt`
+5. Submit a pull request
+
+## License
+
+This project is licensed under the MIT License - see LICENSE file for details.
+
+## Acknowledgments
+
+- Ballistics physics based on Robert McCoy's "Modern Exterior Ballistics"
+- Drag tables from military ballistics research
+- BC segmentation algorithms from Bryan Litz's research
+- Community contributions and testing
+
+## Support
+
+For issues, questions, or contributions:
+- GitHub Issues: [github.com/ajokela/ballistics-engine/issues](https://github.com/ajokela/ballistics-engine/issues)

package/ballistics_engine.d.ts

@@ -0,0 +1,49 @@
+/* tslint:disable */
+/* eslint-disable */
+/**
+ * Object-oriented calculator for programmatic use
+ * Provides a type-safe, fluent API alternative to the CLI interface
+ */
+export class Calculator {
+  free(): void;
+  [Symbol.dispose](): void;
+  setAltitude(altitude_ft: number): Calculator;
+  setDiameter(diameter_inches: number): Calculator;
+  setHumidity(humidity: number): Calculator;
+  setPressure(pressure_inhg: number): Calculator;
+  setVelocity(velocity_fps: number): Calculator;
+  setMaxRange(range_yards: number): Calculator;
+  setDragModel(model: string): Calculator;
+  setZeroRange(range_yards: number): Calculator;
+  setTemperature(temp_f: number): Calculator;
+  setSightHeight(height_inches: number): Calculator;
+  enableCoriolis(enabled: boolean, latitude?: number | null): Calculator;
+  /**
+   * Get full trajectory table as array of points
+   * Returns array of: [{ range_yards, drop_inches, windage_inches, velocity_fps, energy_ftlb, time_sec }, ...]
+   */
+  getFullTrajectory(): any;
+  /**
+   * Calculate trajectory and return result as JavaScript object
+   * Returns: { range_yards, drop_inches, windage_inches, velocity_fps, energy_ftlb, time_sec }
+   */
+  calculateTrajectory(range_yards: number): any;
+  enableSpinDrift(enabled: boolean, twist_rate?: number | null): Calculator;
+  /**
+   * Create a new calculator with default values
+   * Defaults: .308 Winchester 168gr at 2700 fps, standard atmosphere
+   */
+  constructor();
+  setBC(bc: number): Calculator;
+  setMass(mass_grains: number): Calculator;
+  setWind(speed_mph: number, direction_deg: number): Calculator;
+}
+export class WasmBallistics {
+  free(): void;
+  [Symbol.dispose](): void;
+  /**
+   * Run a command and return the output
+   */
+  runCommand(command: string): string;
+  constructor();
+}

package/ballistics_engine.js

@@ -0,0 +1,544 @@
+
+let imports = {};
+imports['__wbindgen_placeholder__'] = module.exports;
+
+let cachedUint8ArrayMemory0 = null;
+
+function getUint8ArrayMemory0() {
+    if (cachedUint8ArrayMemory0 === null || cachedUint8ArrayMemory0.byteLength === 0) {
+        cachedUint8ArrayMemory0 = new Uint8Array(wasm.memory.buffer);
+    }
+    return cachedUint8ArrayMemory0;
+}
+
+let cachedTextDecoder = new TextDecoder('utf-8', { ignoreBOM: true, fatal: true });
+
+cachedTextDecoder.decode();
+
+function decodeText(ptr, len) {
+    return cachedTextDecoder.decode(getUint8ArrayMemory0().subarray(ptr, ptr + len));
+}
+
+function getStringFromWasm0(ptr, len) {
+    ptr = ptr >>> 0;
+    return decodeText(ptr, len);
+}
+
+function addToExternrefTable0(obj) {
+    const idx = wasm.__externref_table_alloc();
+    wasm.__wbindgen_externrefs.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 getArrayU8FromWasm0(ptr, len) {
+    ptr = ptr >>> 0;
+    return getUint8ArrayMemory0().subarray(ptr / 1, ptr / 1 + len);
+}
+
+function isLikeNone(x) {
+    return x === undefined || x === null;
+}
+
+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;
+}
+
+function takeFromExternrefTable0(idx) {
+    const value = wasm.__wbindgen_externrefs.get(idx);
+    wasm.__externref_table_dealloc(idx);
+    return value;
+}
+
+const CalculatorFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_calculator_free(ptr >>> 0, 1));
+/**
+ * Object-oriented calculator for programmatic use
+ * Provides a type-safe, fluent API alternative to the CLI interface
+ */
+class Calculator {
+
+    static __wrap(ptr) {
+        ptr = ptr >>> 0;
+        const obj = Object.create(Calculator.prototype);
+        obj.__wbg_ptr = ptr;
+        CalculatorFinalization.register(obj, obj.__wbg_ptr, obj);
+        return obj;
+    }
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        CalculatorFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_calculator_free(ptr, 0);
+    }
+    /**
+     * @param {number} altitude_ft
+     * @returns {Calculator}
+     */
+    setAltitude(altitude_ft) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setAltitude(ptr, altitude_ft);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} diameter_inches
+     * @returns {Calculator}
+     */
+    setDiameter(diameter_inches) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setDiameter(ptr, diameter_inches);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} humidity
+     * @returns {Calculator}
+     */
+    setHumidity(humidity) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setHumidity(ptr, humidity);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} pressure_inhg
+     * @returns {Calculator}
+     */
+    setPressure(pressure_inhg) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setPressure(ptr, pressure_inhg);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} velocity_fps
+     * @returns {Calculator}
+     */
+    setVelocity(velocity_fps) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setVelocity(ptr, velocity_fps);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} range_yards
+     * @returns {Calculator}
+     */
+    setMaxRange(range_yards) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setMaxRange(ptr, range_yards);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {string} model
+     * @returns {Calculator}
+     */
+    setDragModel(model) {
+        const ptr = this.__destroy_into_raw();
+        const ptr0 = passStringToWasm0(model, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+        const len0 = WASM_VECTOR_LEN;
+        const ret = wasm.calculator_setDragModel(ptr, ptr0, len0);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} range_yards
+     * @returns {Calculator}
+     */
+    setZeroRange(range_yards) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setZeroRange(ptr, range_yards);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} temp_f
+     * @returns {Calculator}
+     */
+    setTemperature(temp_f) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setTemperature(ptr, temp_f);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} height_inches
+     * @returns {Calculator}
+     */
+    setSightHeight(height_inches) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setSightHeight(ptr, height_inches);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {boolean} enabled
+     * @param {number | null} [latitude]
+     * @returns {Calculator}
+     */
+    enableCoriolis(enabled, latitude) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_enableCoriolis(ptr, enabled, !isLikeNone(latitude), isLikeNone(latitude) ? 0 : latitude);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * Get full trajectory table as array of points
+     * Returns array of: [{ range_yards, drop_inches, windage_inches, velocity_fps, energy_ftlb, time_sec }, ...]
+     * @returns {any}
+     */
+    getFullTrajectory() {
+        const ret = wasm.calculator_getFullTrajectory(this.__wbg_ptr);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * Calculate trajectory and return result as JavaScript object
+     * Returns: { range_yards, drop_inches, windage_inches, velocity_fps, energy_ftlb, time_sec }
+     * @param {number} range_yards
+     * @returns {any}
+     */
+    calculateTrajectory(range_yards) {
+        const ret = wasm.calculator_calculateTrajectory(this.__wbg_ptr, range_yards);
+        if (ret[2]) {
+            throw takeFromExternrefTable0(ret[1]);
+        }
+        return takeFromExternrefTable0(ret[0]);
+    }
+    /**
+     * @param {boolean} enabled
+     * @param {number | null} [twist_rate]
+     * @returns {Calculator}
+     */
+    enableSpinDrift(enabled, twist_rate) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_enableSpinDrift(ptr, enabled, !isLikeNone(twist_rate), isLikeNone(twist_rate) ? 0 : twist_rate);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * Create a new calculator with default values
+     * Defaults: .308 Winchester 168gr at 2700 fps, standard atmosphere
+     */
+    constructor() {
+        const ret = wasm.calculator_new();
+        this.__wbg_ptr = ret >>> 0;
+        CalculatorFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+    /**
+     * @param {number} bc
+     * @returns {Calculator}
+     */
+    setBC(bc) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setBC(ptr, bc);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} mass_grains
+     * @returns {Calculator}
+     */
+    setMass(mass_grains) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setMass(ptr, mass_grains);
+        return Calculator.__wrap(ret);
+    }
+    /**
+     * @param {number} speed_mph
+     * @param {number} direction_deg
+     * @returns {Calculator}
+     */
+    setWind(speed_mph, direction_deg) {
+        const ptr = this.__destroy_into_raw();
+        const ret = wasm.calculator_setWind(ptr, speed_mph, direction_deg);
+        return Calculator.__wrap(ret);
+    }
+}
+if (Symbol.dispose) Calculator.prototype[Symbol.dispose] = Calculator.prototype.free;
+
+exports.Calculator = Calculator;
+
+const WasmBallisticsFinalization = (typeof FinalizationRegistry === 'undefined')
+    ? { register: () => {}, unregister: () => {} }
+    : new FinalizationRegistry(ptr => wasm.__wbg_wasmballistics_free(ptr >>> 0, 1));
+
+class WasmBallistics {
+
+    __destroy_into_raw() {
+        const ptr = this.__wbg_ptr;
+        this.__wbg_ptr = 0;
+        WasmBallisticsFinalization.unregister(this);
+        return ptr;
+    }
+
+    free() {
+        const ptr = this.__destroy_into_raw();
+        wasm.__wbg_wasmballistics_free(ptr, 0);
+    }
+    /**
+     * Run a command and return the output
+     * @param {string} command
+     * @returns {string}
+     */
+    runCommand(command) {
+        let deferred3_0;
+        let deferred3_1;
+        try {
+            const ptr0 = passStringToWasm0(command, wasm.__wbindgen_malloc, wasm.__wbindgen_realloc);
+            const len0 = WASM_VECTOR_LEN;
+            const ret = wasm.wasmballistics_runCommand(this.__wbg_ptr, ptr0, len0);
+            var ptr2 = ret[0];
+            var len2 = ret[1];
+            if (ret[3]) {
+                ptr2 = 0; len2 = 0;
+                throw takeFromExternrefTable0(ret[2]);
+            }
+            deferred3_0 = ptr2;
+            deferred3_1 = len2;
+            return getStringFromWasm0(ptr2, len2);
+        } finally {
+            wasm.__wbindgen_free(deferred3_0, deferred3_1, 1);
+        }
+    }
+    constructor() {
+        const ret = wasm.wasmballistics_new();
+        this.__wbg_ptr = ret >>> 0;
+        WasmBallisticsFinalization.register(this, this.__wbg_ptr, this);
+        return this;
+    }
+}
+if (Symbol.dispose) WasmBallistics.prototype[Symbol.dispose] = WasmBallistics.prototype.free;
+
+exports.WasmBallistics = WasmBallistics;
+
+exports.__wbg___wbindgen_is_function_ee8a6c5833c90377 = function(arg0) {
+    const ret = typeof(arg0) === 'function';
+    return ret;
+};
+
+exports.__wbg___wbindgen_is_object_c818261d21f283a4 = function(arg0) {
+    const val = arg0;
+    const ret = typeof(val) === 'object' && val !== null;
+    return ret;
+};
+
+exports.__wbg___wbindgen_is_string_fbb76cb2940daafd = function(arg0) {
+    const ret = typeof(arg0) === 'string';
+    return ret;
+};
+
+exports.__wbg___wbindgen_is_undefined_2d472862bd29a478 = function(arg0) {
+    const ret = arg0 === undefined;
+    return ret;
+};
+
+exports.__wbg___wbindgen_throw_b855445ff6a94295 = function(arg0, arg1) {
+    throw new Error(getStringFromWasm0(arg0, arg1));
+};
+
+exports.__wbg_call_525440f72fbfc0ea = function() { return handleError(function (arg0, arg1, arg2) {
+    const ret = arg0.call(arg1, arg2);
+    return ret;
+}, arguments) };
+
+exports.__wbg_call_e762c39fa8ea36bf = function() { return handleError(function (arg0, arg1) {
+    const ret = arg0.call(arg1);
+    return ret;
+}, arguments) };
+
+exports.__wbg_crypto_574e78ad8b13b65f = function(arg0) {
+    const ret = arg0.crypto;
+    return ret;
+};
+
+exports.__wbg_getRandomValues_b8f5dbd5f3995a9e = function() { return handleError(function (arg0, arg1) {
+    arg0.getRandomValues(arg1);
+}, arguments) };
+
+exports.__wbg_length_69bca3cb64fc8748 = function(arg0) {
+    const ret = arg0.length;
+    return ret;
+};
+
+exports.__wbg_msCrypto_a61aeb35a24c1329 = function(arg0) {
+    const ret = arg0.msCrypto;
+    return ret;
+};
+
+exports.__wbg_new_1acc0b6eea89d040 = function() {
+    const ret = new Object();
+    return ret;
+};
+
+exports.__wbg_new_e17d9f43105b08be = function() {
+    const ret = new Array();
+    return ret;
+};
+
+exports.__wbg_new_no_args_ee98eee5275000a4 = function(arg0, arg1) {
+    const ret = new Function(getStringFromWasm0(arg0, arg1));
+    return ret;
+};
+
+exports.__wbg_new_with_length_01aa0dc35aa13543 = function(arg0) {
+    const ret = new Uint8Array(arg0 >>> 0);
+    return ret;
+};
+
+exports.__wbg_node_905d3e251edff8a2 = function(arg0) {
+    const ret = arg0.node;
+    return ret;
+};
+
+exports.__wbg_process_dc0fbacc7c1c06f7 = function(arg0) {
+    const ret = arg0.process;
+    return ret;
+};
+
+exports.__wbg_prototypesetcall_2a6620b6922694b2 = function(arg0, arg1, arg2) {
+    Uint8Array.prototype.set.call(getArrayU8FromWasm0(arg0, arg1), arg2);
+};
+
+exports.__wbg_push_df81a39d04db858c = function(arg0, arg1) {
+    const ret = arg0.push(arg1);
+    return ret;
+};
+
+exports.__wbg_randomFillSync_ac0988aba3254290 = function() { return handleError(function (arg0, arg1) {
+    arg0.randomFillSync(arg1);
+}, arguments) };
+
+exports.__wbg_require_60cc747a6bc5215a = function() { return handleError(function () {
+    const ret = module.require;
+    return ret;
+}, arguments) };
+
+exports.__wbg_set_c2abbebe8b9ebee1 = function() { return handleError(function (arg0, arg1, arg2) {
+    const ret = Reflect.set(arg0, arg1, arg2);
+    return ret;
+}, arguments) };
+
+exports.__wbg_static_accessor_GLOBAL_89e1d9ac6a1b250e = function() {
+    const ret = typeof global === 'undefined' ? null : global;
+    return isLikeNone(ret) ? 0 : addToExternrefTable0(ret);
+};
+
+exports.__wbg_static_accessor_GLOBAL_THIS_8b530f326a9e48ac = function() {
+    const ret = typeof globalThis === 'undefined' ? null : globalThis;
+    return isLikeNone(ret) ? 0 : addToExternrefTable0(ret);
+};
+
+exports.__wbg_static_accessor_SELF_6fdf4b64710cc91b = function() {
+    const ret = typeof self === 'undefined' ? null : self;
+    return isLikeNone(ret) ? 0 : addToExternrefTable0(ret);
+};
+
+exports.__wbg_static_accessor_WINDOW_b45bfc5a37f6cfa2 = function() {
+    const ret = typeof window === 'undefined' ? null : window;
+    return isLikeNone(ret) ? 0 : addToExternrefTable0(ret);
+};
+
+exports.__wbg_subarray_480600f3d6a9f26c = function(arg0, arg1, arg2) {
+    const ret = arg0.subarray(arg1 >>> 0, arg2 >>> 0);
+    return ret;
+};
+
+exports.__wbg_versions_c01dfd4722a88165 = function(arg0) {
+    const ret = arg0.versions;
+    return ret;
+};
+
+exports.__wbindgen_cast_2241b6af4c4b2941 = function(arg0, arg1) {
+    // Cast intrinsic for `Ref(String) -> Externref`.
+    const ret = getStringFromWasm0(arg0, arg1);
+    return ret;
+};
+
+exports.__wbindgen_cast_cb9088102bce6b30 = function(arg0, arg1) {
+    // Cast intrinsic for `Ref(Slice(U8)) -> NamedExternref("Uint8Array")`.
+    const ret = getArrayU8FromWasm0(arg0, arg1);
+    return ret;
+};
+
+exports.__wbindgen_cast_d6cd19b81560fd6e = function(arg0) {
+    // Cast intrinsic for `F64 -> Externref`.
+    const ret = arg0;
+    return ret;
+};
+
+exports.__wbindgen_init_externref_table = function() {
+    const table = wasm.__wbindgen_externrefs;
+    const offset = table.grow(4);
+    table.set(0, undefined);
+    table.set(offset + 0, undefined);
+    table.set(offset + 1, null);
+    table.set(offset + 2, true);
+    table.set(offset + 3, false);
+    ;
+};
+
+const wasmPath = `${__dirname}/ballistics_engine_bg.wasm`;
+const wasmBytes = require('fs').readFileSync(wasmPath);
+const wasmModule = new WebAssembly.Module(wasmBytes);
+const wasm = exports.__wasm = new WebAssembly.Instance(wasmModule, imports).exports;
+
+wasm.__wbindgen_start();
+

package/package.json

@@ -0,0 +1,27 @@
+{
+  "name": "ballistics-engine",
+  "collaborators": [
+    "Alex Jokela <email@tinycomputers.io>"
+  ],
+  "description": "High-performance ballistics trajectory engine with professional physics",
+  "version": "0.13.3",
+  "license": "MIT OR Apache-2.0",
+  "repository": {
+    "type": "git",
+    "url": "https://github.com/ajokela/ballistics-engine"
+  },
+  "files": [
+    "ballistics_engine_bg.wasm",
+    "ballistics_engine.js",
+    "ballistics_engine.d.ts"
+  ],
+  "main": "ballistics_engine.js",
+  "homepage": "https://ballistics.rs/",
+  "types": "ballistics_engine.d.ts",
+  "keywords": [
+    "ballistics",
+    "trajectory",
+    "physics",
+    "simulation"
+  ]
+}