Function calculate_damped_yaw_of_repose 

pub fn calculate_damped_yaw_of_repose(
    stability_factor: f64,
    velocity_mps: f64,
    spin_rate_rad_s: f64,
    _wind_velocity_mps: f64,
    pitch_rate_rad_s: f64,
    air_density_kg_m3: f64,
    caliber_inches: f64,
    length_inches: f64,
    mass_grains: f64,
    mach: f64,
    bullet_type: &str,
) -> (f64, f64)

Calculate yaw of repose with pitch damping effects.

Returns the equilibrium yaw and a signed convergence rate in s^-1: positive values converge toward equilibrium, while negative values identify a divergent pitch mode.