alice-sdf 6.0.0

//! Moon (crescent) SDF (Deep Fried Edition)
//!
//! 2D crescent/moon shape in XZ plane, extruded along Y-axis.
//!
//! Based on Inigo Quilez's sdMoon formula.
//!
//! Author: Moroya Sakamoto

use crate::compiled::real::LawMinMax;
use glam::{Vec2, Vec3};

/// SDF for a crescent moon shape, extruded along Y-axis
///
/// - `d`: distance between the two circle centers
/// - `ra`: outer circle radius
/// - `rb`: inner (subtracted) circle radius
/// - `half_height`: half the extrusion height along Y
#[inline(always)]
pub fn sdf_moon(p: Vec3, d: f32, ra: f32, rb: f32, half_height: f32) -> f32 {
    // 2D moon in XZ plane (IQ's sdMoon)
    let qx = p.x.abs();
    let qz = p.z;
    let q = Vec2::new(qx, qz);

    let d_outer = q.length() - ra;
    let d_inner = (q - Vec2::new(d, 0.0)).length() - rb;
    let d_2d = d_outer.law_max(-d_inner);

    // Extrude along Y
    let d_y = p.y.abs() - half_height;
    let w = Vec2::new(d_2d.law_max(0.0), d_y.law_max(0.0));
    d_2d.law_max(d_y).law_min(0.0) + w.length()
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_moon_inside() {
        // Point on the crescent body (opposite side from inner circle)
        let d = sdf_moon(Vec3::new(0.0, 0.0, -0.8), 0.5, 1.0, 0.8, 0.5);
        assert!(d < 0.0, "Should be inside the crescent, got {}", d);
    }

    #[test]
    fn test_moon_far_outside() {
        let d = sdf_moon(Vec3::new(5.0, 0.0, 0.0), 0.5, 1.0, 0.8, 0.5);
        assert!(d > 0.0, "Far point should be outside, got {}", d);
    }

    #[test]
    fn test_moon_symmetry_x() {
        let d1 = sdf_moon(Vec3::new(0.3, 0.1, -0.5), 0.5, 1.0, 0.8, 0.5);
        let d2 = sdf_moon(Vec3::new(-0.3, 0.1, -0.5), 0.5, 1.0, 0.8, 0.5);
        assert!((d1 - d2).abs() < 0.001, "Should be symmetric in X");
    }

    #[test]
    fn test_moon_symmetry_y() {
        let d1 = sdf_moon(Vec3::new(0.1, 0.2, -0.5), 0.5, 1.0, 0.8, 0.5);
        let d2 = sdf_moon(Vec3::new(0.1, -0.2, -0.5), 0.5, 1.0, 0.8, 0.5);
        assert!((d1 - d2).abs() < 0.001, "Should be symmetric in Y");
    }
}