ftui-extras 0.3.1

//! Quake underwater warp distortion effect.
//!
//! Authentic port of Quake's `R_WarpScreen` / `turbsin` table: for each
//! output pixel, the source coordinates are displaced by sinusoidal waves
//! to create the classic underwater distortion.
//!
//! # Algorithm
//!
//! For each output pixel (x, y):
//! ```text
//! src_x = x + amplitude * sin(y * frequency + time)
//! src_y = y + amplitude * cos(x * frequency + time)
//! ```
//! Then sample from the inner effect's buffer at (src_x, src_y).
//!
//! # Determinism
//!
//! Uses pre-computed sin LUT; given same inputs, output is identical.
//!
//! # No Per-Frame Allocations
//!
//! Inner buffer and LUT are grow-only.

use crate::visual_fx::{BackdropFx, FxContext, FxQuality};
use ftui_render::cell::PackedRgba;

// ---------------------------------------------------------------------------
// Turbulence sin table (256 entries, matching Quake's turbsin)
// ---------------------------------------------------------------------------

/// Pre-computed sin table with 256 entries for fast warp displacement.
fn build_turbsin() -> [f64; 256] {
    let mut table = [0.0f64; 256];
    let mut i = 0;
    while i < 256 {
        table[i] = (i as f64 * std::f64::consts::TAU / 256.0).sin();
        i += 1;
    }
    table
}

// ---------------------------------------------------------------------------
// UnderwaterWarpFx
// ---------------------------------------------------------------------------

/// Quake-style underwater warp distortion.
///
/// Wraps an inner `BackdropFx` and applies sinusoidal displacement to
/// create the classic underwater distortion effect.
///
/// # Quality Degradation
///
/// - `Full`: Full-resolution warp with smooth sampling
/// - `Reduced`: Reduced amplitude (less distortion)
/// - `Minimal`: Very low amplitude, basic warp
/// - `Off`: No rendering
pub struct UnderwaterWarpFx {
    /// Pre-computed 256-entry sin lookup table.
    turbsin: [f64; 256],
    /// The inner effect being warped.
    inner: Box<dyn BackdropFx>,
    /// Scratch buffer for inner effect rendering.
    inner_buf: Vec<PackedRgba>,
    /// Per-column y displacement scratch (grow-only; no per-frame alloc).
    col_dy: Vec<f64>,
    /// Warp amplitude (in pixels). Default: 2.0.
    amplitude: f64,
    /// Warp frequency. Default: 0.3.
    frequency: f64,
    /// Cached width.
    last_width: u16,
    /// Cached height.
    last_height: u16,
}

impl UnderwaterWarpFx {
    /// Create a new underwater warp wrapping an inner effect.
    pub fn new(inner: Box<dyn BackdropFx>) -> Self {
        Self {
            turbsin: build_turbsin(),
            inner,
            inner_buf: Vec::new(),
            col_dy: Vec::new(),
            amplitude: 2.0,
            frequency: 0.3,
            last_width: 0,
            last_height: 0,
        }
    }

    /// Create with custom amplitude and frequency.
    pub fn with_params(inner: Box<dyn BackdropFx>, amplitude: f64, frequency: f64) -> Self {
        Self {
            turbsin: build_turbsin(),
            inner,
            inner_buf: Vec::new(),
            col_dy: Vec::new(),
            amplitude,
            frequency,
            last_width: 0,
            last_height: 0,
        }
    }

    /// Set the warp amplitude (in pixels).
    pub fn set_amplitude(&mut self, amplitude: f64) {
        self.amplitude = amplitude;
    }

    /// Set the warp frequency.
    pub fn set_frequency(&mut self, frequency: f64) {
        self.frequency = frequency;
    }

    /// Access the inner effect.
    pub fn inner(&self) -> &dyn BackdropFx {
        &*self.inner
    }

    /// Access the inner effect mutably.
    pub fn inner_mut(&mut self) -> &mut dyn BackdropFx {
        &mut *self.inner
    }

    /// Look up sin value from pre-computed table with interpolation.
    #[inline]
    fn turbsin_lookup(&self, phase: f64) -> f64 {
        // Map phase to 0..256 range
        let idx_f = phase.rem_euclid(256.0);
        let idx0 = idx_f as usize & 255;
        let idx1 = (idx0 + 1) & 255;
        let frac = idx_f - idx_f.floor();
        self.turbsin[idx0] * (1.0 - frac) + self.turbsin[idx1] * frac
    }
}

impl BackdropFx for UnderwaterWarpFx {
    fn name(&self) -> &'static str {
        "Underwater Warp"
    }

    fn resize(&mut self, width: u16, height: u16) {
        self.inner.resize(width, height);
        self.last_width = width;
        self.last_height = height;
    }

    fn render(&mut self, ctx: FxContext<'_>, out: &mut [PackedRgba]) {
        let w = ctx.width as usize;
        let h = ctx.height as usize;
        if w == 0 || h == 0 {
            return;
        }

        if !ctx.quality.is_enabled() {
            return;
        }

        let len = w * h;

        // Ensure inner buffer
        if self.inner_buf.len() < len {
            self.inner_buf.resize(len, PackedRgba::rgb(0, 0, 0));
        }

        // Render inner effect to scratch buffer
        self.inner.render(ctx, &mut self.inner_buf[..len]);

        self.last_width = ctx.width;
        self.last_height = ctx.height;

        // Quality-adjusted amplitude
        let amp = match ctx.quality {
            FxQuality::Full => self.amplitude,
            FxQuality::Reduced => self.amplitude * 0.6,
            FxQuality::Minimal => self.amplitude * 0.3,
            FxQuality::Off => return,
        };

        let time = ctx.time_seconds * 40.0; // Speed factor matching Quake's time scale
        let freq = self.frequency;
        let _w_f = w as f64;
        let _h_f = h as f64;
        let freq10 = freq * 10.0;
        let time_y = time * 0.9;
        let w_max = (w - 1) as f64;
        let h_max = (h - 1) as f64;

        // Precompute per-column y-displacements (phase_y depends only on x).
        // Reuse grow-only scratch to avoid per-frame allocation.
        if self.col_dy.len() != w {
            self.col_dy.resize(w, 0.0);
        }
        for x in 0..w {
            let phase_y = (x as f64) * freq10 + time_y;
            self.col_dy[x] = amp * self.turbsin_lookup(phase_y);
        }
        let col_dy = &self.col_dy;

        // Apply warp displacement
        for y in 0..h {
            // phase_x depends only on y — hoist out of inner x-loop
            let phase_x = (y as f64) * freq10 + time;
            let row_dx = amp * self.turbsin_lookup(phase_x);
            let y_f = y as f64;
            let dst_base = y * w;

            for x in 0..w {
                let dy = col_dy[x];

                // Source coordinates with clamping
                let src_x = ((x as f64 + row_dx).round()).clamp(0.0, w_max) as usize;
                let src_y = ((y_f + dy).round()).clamp(0.0, h_max) as usize;

                out[dst_base + x] = self.inner_buf[src_y * w + src_x];
            }
        }
    }
}

impl std::fmt::Debug for UnderwaterWarpFx {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        f.debug_struct("UnderwaterWarpFx")
            .field("amplitude", &self.amplitude)
            .field("frequency", &self.frequency)
            .field("last_width", &self.last_width)
            .field("last_height", &self.last_height)
            .finish()
    }
}

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

    /// Gradient effect for testing warp displacement.
    struct GradientFx;
    impl BackdropFx for GradientFx {
        fn name(&self) -> &'static str {
            "Gradient"
        }
        fn render(&mut self, ctx: FxContext<'_>, out: &mut [PackedRgba]) {
            let w = ctx.width as usize;
            let h = ctx.height as usize;
            for y in 0..h {
                for x in 0..w {
                    let r = ((x * 255) / w.max(1)) as u8;
                    let g = ((y * 255) / h.max(1)) as u8;
                    out[y * w + x] = PackedRgba::rgb(r, g, 128);
                }
            }
        }
    }

    fn default_theme() -> ThemeInputs {
        ThemeInputs::default_dark()
    }

    fn make_ctx(width: u16, height: u16, frame: u64) -> FxContext<'static> {
        let theme = Box::leak(Box::new(default_theme()));
        FxContext {
            width,
            height,
            frame,
            time_seconds: frame as f64 / 60.0,
            quality: FxQuality::Full,
            theme,
        }
    }

    #[test]
    fn warp_produces_distortion() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(20, 15, 5);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp.render(ctx, &mut buf);

        // Output should differ from a plain gradient
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx, &mut plain_buf);

        // At least some pixels should be displaced
        let different = buf
            .iter()
            .zip(plain_buf.iter())
            .filter(|(a, b)| a != b)
            .count();
        assert!(different > 0, "Warp should displace at least some pixels");
    }

    #[test]
    fn warp_zero_dimensions() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(0, 0, 0);
        let mut buf = vec![];
        warp.render(ctx, &mut buf);
        // Should not panic
    }

    #[test]
    fn warp_deterministic() {
        let inner1 = Box::new(GradientFx);
        let inner2 = Box::new(GradientFx);
        let mut warp1 = UnderwaterWarpFx::new(inner1);
        let mut warp2 = UnderwaterWarpFx::new(inner2);

        let ctx = make_ctx(20, 15, 10);
        let mut buf1 = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf2 = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp1.render(ctx, &mut buf1);
        warp2.render(ctx, &mut buf2);
        assert_eq!(buf1, buf2);
    }

    #[test]
    fn turbsin_table_range() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        for i in 0..256 {
            let v = warp.turbsin[i];
            assert!((-1.0..=1.0).contains(&v), "turbsin[{i}] = {v} out of range");
        }
    }

    #[test]
    fn custom_params() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::with_params(inner, 5.0, 0.5);
        assert!((warp.amplitude - 5.0).abs() < f64::EPSILON);
        assert!((warp.frequency - 0.5).abs() < f64::EPSILON);

        warp.set_amplitude(3.0);
        warp.set_frequency(0.2);
        assert!((warp.amplitude - 3.0).abs() < f64::EPSILON);
        assert!((warp.frequency - 0.2).abs() < f64::EPSILON);
    }

    // ── Turbsin table ──────────────────────────────────────────────

    #[test]
    fn turbsin_has_256_entries() {
        let table = build_turbsin();
        assert_eq!(table.len(), 256);
    }

    #[test]
    fn turbsin_values_in_sine_range() {
        let table = build_turbsin();
        for (i, &v) in table.iter().enumerate() {
            assert!(
                (-1.0..=1.0).contains(&v),
                "turbsin[{i}] = {v} out of [-1, 1]"
            );
        }
    }

    #[test]
    fn turbsin_starts_at_zero() {
        let table = build_turbsin();
        assert!(
            table[0].abs() < 1e-10,
            "sin(0) should be ~0, got {}",
            table[0]
        );
    }

    #[test]
    fn turbsin_quarter_is_one() {
        let table = build_turbsin();
        // index 64 = 256/4 = 90 degrees -> sin(pi/2) = 1
        assert!(
            (table[64] - 1.0).abs() < 1e-10,
            "sin(pi/2) should be ~1, got {}",
            table[64]
        );
    }

    #[test]
    fn turbsin_half_is_zero() {
        let table = build_turbsin();
        // index 128 = 256/2 = 180 degrees -> sin(pi) ≈ 0
        assert!(
            table[128].abs() < 1e-10,
            "sin(pi) should be ~0, got {}",
            table[128]
        );
    }

    #[test]
    fn turbsin_three_quarters_is_neg_one() {
        let table = build_turbsin();
        // index 192 = 256*3/4 = 270 degrees -> sin(3*pi/2) = -1
        assert!(
            (table[192] + 1.0).abs() < 1e-10,
            "sin(3pi/2) should be ~-1, got {}",
            table[192]
        );
    }

    // ── Turbsin lookup interpolation ───────────────────────────────

    #[test]
    fn turbsin_lookup_at_integer_indices() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        for i in 0..256 {
            let looked_up = warp.turbsin_lookup(i as f64);
            let direct = warp.turbsin[i];
            assert!(
                (looked_up - direct).abs() < 1e-10,
                "lookup({i}) = {looked_up}, direct = {direct}"
            );
        }
    }

    #[test]
    fn turbsin_lookup_interpolates_between_entries() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        // Midpoint between entries 0 and 1
        let mid = warp.turbsin_lookup(0.5);
        let expected = (warp.turbsin[0] + warp.turbsin[1]) / 2.0;
        assert!(
            (mid - expected).abs() < 1e-10,
            "midpoint: got {mid}, expected {expected}"
        );
    }

    #[test]
    fn turbsin_lookup_wraps_around() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        let at_0 = warp.turbsin_lookup(0.0);
        let at_256 = warp.turbsin_lookup(256.0);
        assert!(
            (at_0 - at_256).abs() < 1e-10,
            "should wrap at 256: {at_0} vs {at_256}"
        );
    }

    #[test]
    fn turbsin_lookup_negative_phase_wraps() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        let pos = warp.turbsin_lookup(64.0);
        let neg = warp.turbsin_lookup(-192.0); // -192 mod 256 = 64
        assert!(
            (pos - neg).abs() < 1e-10,
            "negative phase should wrap: {pos} vs {neg}"
        );
    }

    // ── Quality degradation ────────────────────────────────────────

    #[test]
    fn quality_off_skips_rendering() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let theme = Box::leak(Box::new(default_theme()));
        let ctx = FxContext {
            width: 10,
            height: 10,
            frame: 5,
            time_seconds: 0.1,
            quality: FxQuality::Off,
            theme,
        };
        let sentinel = PackedRgba::rgb(42, 42, 42);
        let mut buf = vec![sentinel; 100];
        warp.render(ctx, &mut buf);
        // Buffer should remain untouched
        assert_eq!(buf[0], sentinel);
        assert_eq!(buf[50], sentinel);
    }

    #[test]
    fn reduced_quality_less_distortion() {
        let theme = Box::leak(Box::new(default_theme()));

        let inner1 = Box::new(GradientFx);
        let mut warp_full = UnderwaterWarpFx::with_params(inner1, 4.0, 0.3);
        let ctx_full = FxContext {
            width: 20,
            height: 15,
            frame: 5,
            time_seconds: 5.0 / 60.0,
            quality: FxQuality::Full,
            theme,
        };

        let inner2 = Box::new(GradientFx);
        let mut warp_reduced = UnderwaterWarpFx::with_params(inner2, 4.0, 0.3);
        let ctx_reduced = FxContext {
            width: 20,
            height: 15,
            frame: 5,
            time_seconds: 5.0 / 60.0,
            quality: FxQuality::Reduced,
            theme,
        };

        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx_full, &mut plain_buf);

        let mut buf_full = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf_reduced = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp_full.render(ctx_full, &mut buf_full);
        warp_reduced.render(ctx_reduced, &mut buf_reduced);

        // Count displaced pixels for each
        let displaced_full = buf_full
            .iter()
            .zip(plain_buf.iter())
            .filter(|(a, b)| a != b)
            .count();
        let displaced_reduced = buf_reduced
            .iter()
            .zip(plain_buf.iter())
            .filter(|(a, b)| a != b)
            .count();

        assert!(
            displaced_full >= displaced_reduced,
            "full quality should have >= displaced pixels: full={displaced_full}, reduced={displaced_reduced}"
        );
    }

    // ── Warp properties ────────────────────────────────────────────

    #[test]
    fn zero_amplitude_produces_no_distortion() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::with_params(inner, 0.0, 0.3);
        let ctx = make_ctx(20, 15, 5);

        let mut warp_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp.render(ctx, &mut warp_buf);

        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx, &mut plain_buf);

        assert_eq!(
            warp_buf, plain_buf,
            "zero amplitude should produce no distortion"
        );
    }

    #[test]
    fn warp_varies_over_time() {
        let inner1 = Box::new(GradientFx);
        let inner2 = Box::new(GradientFx);
        let mut warp1 = UnderwaterWarpFx::new(inner1);
        let mut warp2 = UnderwaterWarpFx::new(inner2);

        let mut buf1 = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf2 = vec![PackedRgba::rgb(0, 0, 0); 300];

        let ctx_t0 = make_ctx(20, 15, 0);
        let ctx_t5 = make_ctx(20, 15, 30); // different time
        warp1.render(ctx_t0, &mut buf1);
        warp2.render(ctx_t5, &mut buf2);

        assert_ne!(buf1, buf2, "warp should vary over time");
    }

    #[test]
    fn name_returns_expected() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        assert_eq!(warp.name(), "Underwater Warp");
    }

    #[test]
    fn debug_format_includes_fields() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        let debug = format!("{warp:?}");
        assert!(debug.contains("UnderwaterWarpFx"));
        assert!(debug.contains("amplitude"));
        assert!(debug.contains("frequency"));
    }

    #[test]
    fn inner_accessor() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        assert_eq!(warp.inner().name(), "Gradient");
    }

    #[test]
    fn inner_mut_accessor() {
        let mut warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        assert_eq!(warp.inner_mut().name(), "Gradient");
    }

    #[test]
    fn resize_updates_cached_dimensions() {
        let mut warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        warp.resize(80, 40);
        assert_eq!(warp.last_width, 80);
        assert_eq!(warp.last_height, 40);
    }

    #[test]
    fn default_amplitude_and_frequency() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        assert!((warp.amplitude - 2.0).abs() < f64::EPSILON);
        assert!((warp.frequency - 0.3).abs() < f64::EPSILON);
    }

    #[test]
    fn all_output_pixels_from_inner_source() {
        // With warp, every output pixel should be a pixel from the inner effect
        // (just displaced, not interpolated or invented)
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(20, 15, 5);

        let mut warp_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp.render(ctx, &mut warp_buf);

        // Build set of all possible inner pixels
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx, &mut plain_buf);

        let inner_set: std::collections::HashSet<PackedRgba> = plain_buf.iter().copied().collect();

        for (i, px) in warp_buf.iter().enumerate() {
            assert!(
                inner_set.contains(px),
                "output pixel {i} ({px:?}) is not from inner effect"
            );
        }
    }

    // =========================================================================
    // Quality Tier Ordering (bd-3g6pi)
    // =========================================================================

    #[test]
    fn minimal_quality_less_distortion_than_reduced() {
        let theme = Box::leak(Box::new(default_theme()));

        let inner_r = Box::new(GradientFx);
        let mut warp_reduced = UnderwaterWarpFx::with_params(inner_r, 4.0, 0.3);
        let ctx_reduced = FxContext {
            width: 20,
            height: 15,
            frame: 5,
            time_seconds: 5.0 / 60.0,
            quality: FxQuality::Reduced,
            theme,
        };

        let inner_m = Box::new(GradientFx);
        let mut warp_minimal = UnderwaterWarpFx::with_params(inner_m, 4.0, 0.3);
        let ctx_minimal = FxContext {
            quality: FxQuality::Minimal,
            ..ctx_reduced
        };

        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx_reduced, &mut plain_buf);

        let mut buf_reduced = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf_minimal = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp_reduced.render(ctx_reduced, &mut buf_reduced);
        warp_minimal.render(ctx_minimal, &mut buf_minimal);

        let displaced_reduced = buf_reduced
            .iter()
            .zip(plain_buf.iter())
            .filter(|(a, b)| a != b)
            .count();
        let displaced_minimal = buf_minimal
            .iter()
            .zip(plain_buf.iter())
            .filter(|(a, b)| a != b)
            .count();

        assert!(
            displaced_reduced >= displaced_minimal,
            "Reduced should displace >= Minimal: reduced={displaced_reduced}, minimal={displaced_minimal}"
        );
    }

    #[test]
    fn quality_tier_distortion_ordering() {
        // Full >= Reduced >= Minimal > Off (0)
        let theme = Box::leak(Box::new(default_theme()));
        let base_ctx = FxContext {
            width: 20,
            height: 15,
            frame: 5,
            time_seconds: 5.0 / 60.0,
            quality: FxQuality::Full,
            theme,
        };

        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(base_ctx, &mut plain_buf);

        let mut counts = Vec::new();
        for quality in [
            FxQuality::Full,
            FxQuality::Reduced,
            FxQuality::Minimal,
            FxQuality::Off,
        ] {
            let inner = Box::new(GradientFx);
            let mut warp = UnderwaterWarpFx::with_params(inner, 4.0, 0.3);
            let ctx = FxContext {
                quality,
                ..base_ctx
            };
            let sentinel = PackedRgba::rgb(42, 42, 42);
            let mut buf = vec![sentinel; 300];
            warp.render(ctx, &mut buf);

            let displaced = if quality == FxQuality::Off {
                // Off doesn't render, all sentinels remain
                0
            } else {
                buf.iter()
                    .zip(plain_buf.iter())
                    .filter(|(a, b)| a != b)
                    .count()
            };
            counts.push((quality, displaced));
        }

        // Full >= Reduced >= Minimal >= Off
        for i in 0..counts.len() - 1 {
            assert!(
                counts[i].1 >= counts[i + 1].1,
                "{:?} ({}) should displace >= {:?} ({})",
                counts[i].0,
                counts[i].1,
                counts[i + 1].0,
                counts[i + 1].1,
            );
        }
    }

    #[test]
    fn minimal_quality_produces_output() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let theme = Box::leak(Box::new(default_theme()));
        let ctx = FxContext {
            width: 10,
            height: 10,
            frame: 5,
            time_seconds: 0.1,
            quality: FxQuality::Minimal,
            theme,
        };
        let sentinel = PackedRgba::rgb(99, 99, 99);
        let mut buf = vec![sentinel; 100];
        warp.render(ctx, &mut buf);

        // Output should not be all sentinels (i.e., rendering happened)
        let changed = buf.iter().filter(|c| **c != sentinel).count();
        assert!(
            changed > 0,
            "Minimal quality should produce output, not leave buffer untouched"
        );
    }

    // =========================================================================
    // High Amplitude Clamping (bd-3g6pi)
    // =========================================================================

    #[test]
    fn high_amplitude_does_not_panic_or_oob() {
        // Very large amplitude should still produce valid output (clamped coords)
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::with_params(inner, 100.0, 0.5);
        let ctx = make_ctx(20, 15, 5);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        warp.render(ctx, &mut buf);
        // Should not panic, and all pixels should be from inner source
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 300];
        plain.render(ctx, &mut plain_buf);

        let inner_set: std::collections::HashSet<PackedRgba> = plain_buf.iter().copied().collect();
        for (i, px) in buf.iter().enumerate() {
            assert!(
                inner_set.contains(px),
                "high-amp pixel {i} ({px:?}) not from inner"
            );
        }
    }

    #[test]
    fn very_high_amplitude_clamps_to_edges() {
        // With extremely high amplitude, most source coords will clamp to edges.
        // This means many output pixels will be edge pixels.
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::with_params(inner, 1000.0, 0.5);
        let ctx = make_ctx(10, 10, 3);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 100];
        warp.render(ctx, &mut buf);

        // The gradient FX produces (0,0,128) at top-left and various values at edges.
        // With extreme amplitude, many pixels should be clamped to boundary values.
        // Just verify no panic and all pixels are from the inner set.
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 100];
        plain.render(ctx, &mut plain_buf);

        let inner_set: std::collections::HashSet<PackedRgba> = plain_buf.iter().copied().collect();
        for px in &buf {
            assert!(inner_set.contains(px));
        }
    }

    // =========================================================================
    // Edge Dimensions (bd-3g6pi)
    // =========================================================================

    #[test]
    fn warp_1x1_dimension() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(1, 1, 5);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 1];
        warp.render(ctx, &mut buf);
        // With 1x1, there's only one pixel and it must come from inner (0,0)
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 1];
        plain.render(ctx, &mut plain_buf);
        assert_eq!(
            buf[0], plain_buf[0],
            "1x1 should produce the only inner pixel"
        );
    }

    #[test]
    fn warp_zero_width_nonzero_height() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(0, 10, 5);
        let mut buf = vec![];
        warp.render(ctx, &mut buf);
        // Should not panic
    }

    #[test]
    fn warp_nonzero_width_zero_height() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(10, 0, 5);
        let mut buf = vec![];
        warp.render(ctx, &mut buf);
        // Should not panic
    }

    #[test]
    fn warp_2x2_minimal_grid() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(2, 2, 5);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 4];
        warp.render(ctx, &mut buf);
        // All pixels must come from inner
        let mut plain = GradientFx;
        let mut plain_buf = vec![PackedRgba::rgb(0, 0, 0); 4];
        plain.render(ctx, &mut plain_buf);
        let inner_set: std::collections::HashSet<PackedRgba> = plain_buf.iter().copied().collect();
        for px in &buf {
            assert!(inner_set.contains(px), "2x2 pixel not from inner");
        }
    }

    // =========================================================================
    // Grow-Only Buffer Verification (bd-3g6pi)
    // =========================================================================

    #[test]
    fn inner_buf_grow_only_on_shrink() {
        // Render at large size, then small size. inner_buf should not shrink.
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);

        let ctx_large = make_ctx(30, 20, 1);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 600];
        warp.render(ctx_large, &mut buf);
        let cap_after_large = warp.inner_buf.capacity();

        let ctx_small = make_ctx(10, 5, 2);
        let mut buf_small = vec![PackedRgba::rgb(0, 0, 0); 50];
        warp.render(ctx_small, &mut buf_small);
        let cap_after_small = warp.inner_buf.capacity();

        assert!(
            cap_after_small >= cap_after_large,
            "inner_buf should be grow-only: {} < {}",
            cap_after_small,
            cap_after_large
        );
    }

    #[test]
    fn no_per_frame_allocation_on_repeated_renders() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        let ctx = make_ctx(20, 15, 0);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 300];

        // First render allocates
        warp.render(ctx, &mut buf);
        let cap_inner = warp.inner_buf.capacity();
        let cap_col = warp.col_dy.capacity();

        // Subsequent renders at same size should not grow
        for frame in 1..10 {
            let ctx_n = make_ctx(20, 15, frame);
            warp.render(ctx_n, &mut buf);
        }

        assert_eq!(
            warp.inner_buf.capacity(),
            cap_inner,
            "inner_buf grew during repeated renders"
        );
        assert_eq!(
            warp.col_dy.capacity(),
            cap_col,
            "col_dy grew during repeated renders"
        );
    }

    // =========================================================================
    // Turbsin Symmetry (bd-3g6pi)
    // =========================================================================

    #[test]
    fn turbsin_antisymmetry() {
        // sin(x) = -sin(-x) => turbsin[i] ≈ -turbsin[256-i] for i in 1..256
        let table = build_turbsin();
        for i in 1..256 {
            let pos = table[i];
            let neg = table[256 - i];
            assert!(
                (pos + neg).abs() < 1e-10,
                "turbsin[{i}] ({pos}) + turbsin[{}] ({neg}) should ≈ 0",
                256 - i
            );
        }
    }

    #[test]
    fn turbsin_matches_actual_sin() {
        // Each entry should match sin(i * 2*pi / 256) to high precision
        let table = build_turbsin();
        for (i, &actual) in table.iter().enumerate() {
            let expected = (i as f64 * std::f64::consts::TAU / 256.0).sin();
            assert!(
                (actual - expected).abs() < 1e-14,
                "turbsin[{i}] = {actual}, expected {expected}"
            );
        }
    }

    // =========================================================================
    // Render Updates Dimensions (bd-3g6pi)
    // =========================================================================

    #[test]
    fn render_updates_last_width_height() {
        let inner = Box::new(GradientFx);
        let mut warp = UnderwaterWarpFx::new(inner);
        assert_eq!(warp.last_width, 0);
        assert_eq!(warp.last_height, 0);

        let ctx = make_ctx(25, 18, 1);
        let mut buf = vec![PackedRgba::rgb(0, 0, 0); 450];
        warp.render(ctx, &mut buf);

        assert_eq!(warp.last_width, 25);
        assert_eq!(warp.last_height, 18);
    }

    // =========================================================================
    // Determinism with Different Parameters (bd-3g6pi)
    // =========================================================================

    #[test]
    fn deterministic_with_custom_params() {
        let inner1 = Box::new(GradientFx);
        let inner2 = Box::new(GradientFx);
        let mut w1 = UnderwaterWarpFx::with_params(inner1, 3.5, 0.7);
        let mut w2 = UnderwaterWarpFx::with_params(inner2, 3.5, 0.7);

        let ctx = make_ctx(15, 10, 7);
        let mut buf1 = vec![PackedRgba::rgb(0, 0, 0); 150];
        let mut buf2 = vec![PackedRgba::rgb(0, 0, 0); 150];
        w1.render(ctx, &mut buf1);
        w2.render(ctx, &mut buf2);
        assert_eq!(buf1, buf2, "Custom params should be deterministic");
    }

    #[test]
    fn different_amplitude_produces_different_output() {
        let inner1 = Box::new(GradientFx);
        let inner2 = Box::new(GradientFx);
        let mut w_low = UnderwaterWarpFx::with_params(inner1, 1.0, 0.3);
        let mut w_high = UnderwaterWarpFx::with_params(inner2, 5.0, 0.3);

        let ctx = make_ctx(20, 15, 5);
        let mut buf_low = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf_high = vec![PackedRgba::rgb(0, 0, 0); 300];
        w_low.render(ctx, &mut buf_low);
        w_high.render(ctx, &mut buf_high);

        assert_ne!(
            buf_low, buf_high,
            "Different amplitudes should produce different output"
        );
    }

    #[test]
    fn different_frequency_produces_different_output() {
        let inner1 = Box::new(GradientFx);
        let inner2 = Box::new(GradientFx);
        let mut w_low = UnderwaterWarpFx::with_params(inner1, 2.0, 0.1);
        let mut w_high = UnderwaterWarpFx::with_params(inner2, 2.0, 0.9);

        let ctx = make_ctx(20, 15, 5);
        let mut buf_low = vec![PackedRgba::rgb(0, 0, 0); 300];
        let mut buf_high = vec![PackedRgba::rgb(0, 0, 0); 300];
        w_low.render(ctx, &mut buf_low);
        w_high.render(ctx, &mut buf_high);

        assert_ne!(
            buf_low, buf_high,
            "Different frequencies should produce different output"
        );
    }

    // =========================================================================
    // Turbsin Lookup Edge Cases (bd-3g6pi)
    // =========================================================================

    #[test]
    fn turbsin_lookup_quarter_fractions() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        // At 0.25 between entries 0 and 1
        let at_quarter = warp.turbsin_lookup(0.25);
        let expected = warp.turbsin[0] * 0.75 + warp.turbsin[1] * 0.25;
        assert!(
            (at_quarter - expected).abs() < 1e-10,
            "quarter interpolation: got {at_quarter}, expected {expected}"
        );

        // At 0.75 between entries 0 and 1
        let at_three_quarter = warp.turbsin_lookup(0.75);
        let expected_tq = warp.turbsin[0] * 0.25 + warp.turbsin[1] * 0.75;
        assert!(
            (at_three_quarter - expected_tq).abs() < 1e-10,
            "3/4 interpolation: got {at_three_quarter}, expected {expected_tq}"
        );
    }

    #[test]
    fn turbsin_lookup_large_positive_phase() {
        let warp = UnderwaterWarpFx::new(Box::new(GradientFx));
        // 512.0 mod 256 = 0.0, should equal lookup(0)
        let at_0 = warp.turbsin_lookup(0.0);
        let at_512 = warp.turbsin_lookup(512.0);
        assert!(
            (at_0 - at_512).abs() < 1e-10,
            "512 should wrap to 0: {at_0} vs {at_512}"
        );
    }
}