quantwave-core 0.9.0

//! Native O(1) streaming patterns for Batch D (three-bar patterns).

use crate::indicators::patterns::candle_settings::*;
use crate::traits::Next;

/// CDL3BLACKCROWS
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDL3BLACKCROWS {
    bars_seen: usize,
    window: CandleWindow,
    shadow_very_short_avg: RollingCandleAvg,
}

impl CDL3BLACKCROWS {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(4),
            shadow_very_short_avg: RollingCandleAvg::new(SHADOW_VERY_SHORT),
        }
    }
}

impl Default for CDL3BLACKCROWS {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDL3BLACKCROWS {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.shadow_very_short_avg.push(open, high, low, close);

        if self.bars_seen <= SHADOW_VERY_SHORT.avg_period + 3 {
            return 0.0;
        }

        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && candle_color(self.window.bar(1).open, self.window.bar(1).close) == -1
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
            && self.window.bar(1).close < self.window.bar(2).close
            && self.window.bar(0).close < self.window.bar(1).close
            && self.window.bar(2).open <= self.window.bar(3).open.max(self.window.bar(3).close)
            && self.window.bar(1).open <= self.window.bar(2).open
            && self.window.bar(1).open >= self.window.bar(2).close
            && self.window.bar(0).open <= self.window.bar(1).open
            && self.window.bar(0).open >= self.window.bar(1).close
            && lower_shadow(
                self.window.bar(2).open,
                self.window.bar(2).low,
                self.window.bar(2).close,
            ) < self.shadow_very_short_avg.val(2)
            && lower_shadow(
                self.window.bar(1).open,
                self.window.bar(1).low,
                self.window.bar(1).close,
            ) < self.shadow_very_short_avg.val(1)
            && lower_shadow(
                self.window.bar(0).open,
                self.window.bar(0).low,
                self.window.bar(0).close,
            ) < self.shadow_very_short_avg.val(0)
        {
            return -100.0;
        }
        0.0
    }
}

/// CDL3INSIDE
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDL3INSIDE {
    bars_seen: usize,
    window: CandleWindow,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDL3INSIDE {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDL3INSIDE {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDL3INSIDE {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen <= BODY_SHORT.avg_period.max(BODY_LONG.avg_period) + 2 {
            return 0.0;
        }

        if real_body(self.window.bar(2).open, self.window.bar(2).close) > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_short_avg.val(1)
            && self.window.bar(1).open.max(self.window.bar(1).close)
                < self.window.bar(2).open.max(self.window.bar(2).close)
            && self.window.bar(1).open.min(self.window.bar(1).close)
                > self.window.bar(2).open.min(self.window.bar(2).close)
        {
            if candle_color(self.window.bar(2).open, self.window.bar(2).close) == 1
                && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
                && self.window.bar(0).close < self.window.bar(2).open
            {
                return -100.0;
            } else if candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
                && candle_color(self.window.bar(0).open, self.window.bar(0).close) == 1
                && self.window.bar(0).close > self.window.bar(2).open
            {
                return 100.0;
            }
        }
        0.0
    }
}

/// CDL3OUTSIDE
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDL3OUTSIDE {
    bars_seen: usize,
    window: CandleWindow,
}

impl CDL3OUTSIDE {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
        }
    }
}

impl Default for CDL3OUTSIDE {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDL3OUTSIDE {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);

        if self.bars_seen <= 2 {
            return 0.0;
        }

        let bull = candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && candle_color(self.window.bar(1).open, self.window.bar(1).close) == 1
            && self.window.bar(1).close >= self.window.bar(2).open
            && self.window.bar(1).open <= self.window.bar(2).close
            && self.window.bar(0).close > self.window.bar(1).close;
        let bear = candle_color(self.window.bar(2).open, self.window.bar(2).close) == 1
            && candle_color(self.window.bar(1).open, self.window.bar(1).close) == -1
            && self.window.bar(1).open >= self.window.bar(2).close
            && self.window.bar(1).close <= self.window.bar(2).open
            && self.window.bar(0).close < self.window.bar(1).close;
        if bull {
            return 100.0;
        }
        if bear {
            return -100.0;
        }
        0.0
    }
}

/// CDL3STARSINSOUTH
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDL3STARSINSOUTH {
    bars_seen: usize,
    window: CandleWindow,
    body_long_avg: RollingCandleAvg,
    shadow_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
    shadow_vs_avg: RollingCandleAvg,
}

impl CDL3STARSINSOUTH {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            shadow_long_avg: RollingCandleAvg::new(SHADOW_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
            shadow_vs_avg: RollingCandleAvg::new(SHADOW_VERY_SHORT),
        }
    }
}

impl Default for CDL3STARSINSOUTH {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDL3STARSINSOUTH {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.shadow_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);
        self.shadow_vs_avg.push(open, high, low, close);

        // BODY_LONG, BODY_SHORT and SHADOW_VERY_SHORT all use avg_period 10;
        // SHADOW_LONG.avg_period is 0 and cannot raise the max.
        if self.bars_seen <= BODY_LONG.avg_period + 2 {
            return 0.0;
        }

        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && candle_color(self.window.bar(1).open, self.window.bar(1).close) == -1
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
            && real_body(self.window.bar(2).open, self.window.bar(2).close)
                > self.body_long_avg.val(2)
            && lower_shadow(
                self.window.bar(2).open,
                self.window.bar(2).low,
                self.window.bar(2).close,
            ) > self.shadow_long_avg.val(2)
            && self.window.bar(1).open.min(self.window.bar(1).close)
                > self.window.bar(2).open.min(self.window.bar(2).close)
            && self.window.bar(1).open.max(self.window.bar(1).close)
                < self.window.bar(2).open.max(self.window.bar(2).close)
            && self.window.bar(1).low < self.window.bar(2).low
            // 3rd: short body, both shadows short, range inside the 2nd's range
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                < self.body_short_avg.val(0)
            && upper_shadow(
                self.window.bar(0).open,
                self.window.bar(0).high,
                self.window.bar(0).close,
            ) < self.shadow_vs_avg.val(0)
            && lower_shadow(
                self.window.bar(0).open,
                self.window.bar(0).low,
                self.window.bar(0).close,
            ) < self.shadow_vs_avg.val(0)
            && self.window.bar(0).low > self.window.bar(1).low
            && self.window.bar(0).high < self.window.bar(1).high
        {
            return 100.0;
        }
        0.0
    }
}

/// CDL3WHITESOLDIERS
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDL3WHITESOLDIERS {
    bars_seen: usize,
    window: CandleWindow,
    near_avg: RollingCandleAvg,
    far_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
    shadow_very_short_avg: RollingCandleAvg,
}

impl CDL3WHITESOLDIERS {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(4),
            near_avg: RollingCandleAvg::new(NEAR),
            far_avg: RollingCandleAvg::new(FAR),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
            shadow_very_short_avg: RollingCandleAvg::new(SHADOW_VERY_SHORT),
        }
    }
}

impl Default for CDL3WHITESOLDIERS {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDL3WHITESOLDIERS {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.near_avg.push(open, high, low, close);
        self.far_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);
        self.shadow_very_short_avg.push(open, high, low, close);

        // Oracle lookback: max(SVS, BODY_SHORT, FAR, NEAR) + 2
        let lookback = SHADOW_VERY_SHORT
            .avg_period
            .max(BODY_SHORT.avg_period)
            .max(FAR.avg_period)
            .max(NEAR.avg_period)
            + 2;
        if self.bars_seen <= lookback {
            return 0.0;
        }

        let b2 = self.window.bar(2);
        let b1 = self.window.bar(1);
        let b0 = self.window.bar(0);

        if candle_color(b2.open, b2.close) == 1
            && candle_color(b1.open, b1.close) == 1
            && candle_color(b0.open, b0.close) == 1
            && b1.close > b2.close
            && b0.close > b1.close
            // Short upper shadows
            && upper_shadow(b2.open, b2.high, b2.close) < self.shadow_very_short_avg.val(2)
            && upper_shadow(b1.open, b1.high, b1.close) < self.shadow_very_short_avg.val(1)
            && upper_shadow(b0.open, b0.high, b0.close) < self.shadow_very_short_avg.val(0)
            // Opens within or near the previous body
            && b1.open > b2.open
            && b1.open <= b2.close + self.near_avg.val(1)
            && b0.open > b1.open
            && b0.open <= b1.close + self.near_avg.val(0)
            // Bodies not far shorter than the prior one
            && real_body(b1.open, b1.close) > real_body(b2.open, b2.close) - self.far_avg.val(1)
            && real_body(b0.open, b0.close) > real_body(b1.open, b1.close) - self.far_avg.val(0)
            // Last body not short
            && real_body(b0.open, b0.close) > self.body_short_avg.val(0)
        {
            return 100.0;
        }
        0.0
    }
}

/// CDLABANDONEDBABY
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLABANDONEDBABY {
    bars_seen: usize,
    window: CandleWindow,
    body_doji_avg: RollingCandleAvg,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDLABANDONEDBABY {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_doji_avg: RollingCandleAvg::new(BODY_DOJI),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDLABANDONEDBABY {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLABANDONEDBABY {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_doji_avg.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen
            <= BODY_DOJI
                .avg_period
                .max(BODY_LONG.avg_period)
                .max(BODY_SHORT.avg_period)
                + 2
        {
            return 0.0;
        }

        let penetration = 0.3;
        let base = real_body(self.window.bar(2).open, self.window.bar(2).close)
            > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_doji_avg.val(1)
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                > self.body_short_avg.val(0);
        let bull = base
            && candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == 1
            && self.window.bar(1).high < self.window.bar(2).low
            && self.window.bar(0).low > self.window.bar(1).high
            && self.window.bar(0).close
                > self.window.bar(2).close
                    + real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration;
        let bear = base
            && candle_color(self.window.bar(2).open, self.window.bar(2).close) == 1
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
            && self.window.bar(1).low > self.window.bar(2).high
            && self.window.bar(0).high < self.window.bar(1).low
            && self.window.bar(0).close
                < self.window.bar(2).close
                    - real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration;
        if bull {
            return 100.0;
        }
        if bear {
            return -100.0;
        }
        0.0
    }
}

/// CDLADVANCEBLOCK
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLADVANCEBLOCK {
    bars_seen: usize,
    window: CandleWindow,
    far_avg: RollingCandleAvg,
    near_avg: RollingCandleAvg,
    shadow_long_avg: RollingCandleAvg,
    shadow_short_avg: RollingCandleAvg,
    body_long_avg: RollingCandleAvg,
}

impl CDLADVANCEBLOCK {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            far_avg: RollingCandleAvg::new(FAR),
            near_avg: RollingCandleAvg::new(NEAR),
            shadow_long_avg: RollingCandleAvg::new(SHADOW_LONG),
            shadow_short_avg: RollingCandleAvg::new(SHADOW_SHORT),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
        }
    }
}

impl Default for CDLADVANCEBLOCK {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLADVANCEBLOCK {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.far_avg.push(open, high, low, close);
        self.near_avg.push(open, high, low, close);
        self.shadow_long_avg.push(open, high, low, close);
        self.shadow_short_avg.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);

        // Oracle lookback: max(SHADOW_LONG, SHADOW_SHORT, FAR, NEAR, BODY_LONG) + 2
        // SHADOW_LONG.avg_period is 0, so it cannot raise the max.
        let lookback = SHADOW_SHORT
            .avg_period
            .max(FAR.avg_period)
            .max(NEAR.avg_period)
            .max(BODY_LONG.avg_period)
            + 2;
        if self.bars_seen <= lookback {
            return 0.0;
        }

        let b2 = self.window.bar(2);
        let b1 = self.window.bar(1);
        let b0 = self.window.bar(0);

        let rb2 = real_body(b2.open, b2.close);
        let rb1 = real_body(b1.open, b1.close);
        let rb0 = real_body(b0.open, b0.close);

        // Three white candles with rising closes, opens within/near the prior body,
        // first candle long-bodied with a short upper shadow.
        let base = candle_color(b2.open, b2.close) == 1
            && candle_color(b1.open, b1.close) == 1
            && candle_color(b0.open, b0.close) == 1
            && b1.close > b2.close
            && b0.close > b1.close
            && b1.open > b2.open
            && b1.open <= b2.close + self.near_avg.val(1)
            && b0.open > b1.open
            && b0.open <= b1.close + self.near_avg.val(0)
            && rb2 > self.body_long_avg.val(2)
            && upper_shadow(b2.open, b2.high, b2.close) < self.shadow_short_avg.val(2);

        // Weakness: bodies shrinking and/or upper shadows lengthening.
        let weakness = base
            && ((rb1 < rb2 - self.far_avg.val(1) && rb0 < rb1 + self.near_avg.val(0))
                || (rb0 < rb1
                    && rb1 < rb2
                    && (upper_shadow(b0.open, b0.high, b0.close) > self.shadow_long_avg.val(0)
                        || upper_shadow(b1.open, b1.high, b1.close)
                            > self.shadow_long_avg.val(1)))
                || (rb0 < rb1 - self.far_avg.val(0)));

        if weakness { -100.0 } else { 0.0 }
    }
}

/// CDLEVENINGDOJISTAR
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLEVENINGDOJISTAR {
    bars_seen: usize,
    window: CandleWindow,
    body_doji_avg: RollingCandleAvg,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDLEVENINGDOJISTAR {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_doji_avg: RollingCandleAvg::new(BODY_DOJI),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDLEVENINGDOJISTAR {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLEVENINGDOJISTAR {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_doji_avg.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen
            <= BODY_DOJI
                .avg_period
                .max(BODY_LONG.avg_period)
                .max(BODY_SHORT.avg_period)
                + 2
        {
            return 0.0;
        }

        let penetration = 0.3;
        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == 1
            && real_body(self.window.bar(2).open, self.window.bar(2).close)
                > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_doji_avg.val(1)
            && (self.window.bar(1).open.min(self.window.bar(1).close)
                > self.window.bar(2).open.max(self.window.bar(2).close))
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                > self.body_short_avg.val(0)
            && self.window.bar(0).close
                < self.window.bar(2).close
                    - real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration
        {
            return -100.0;
        }
        0.0
    }
}

/// CDLEVENINGSTAR
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLEVENINGSTAR {
    bars_seen: usize,
    window: CandleWindow,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDLEVENINGSTAR {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDLEVENINGSTAR {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLEVENINGSTAR {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen <= BODY_LONG.avg_period.max(BODY_SHORT.avg_period) + 2 {
            return 0.0;
        }

        let penetration = 0.3;
        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == 1
            && real_body(self.window.bar(2).open, self.window.bar(2).close)
                > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_short_avg.val(1)
            && (self.window.bar(1).open.min(self.window.bar(1).close)
                > self.window.bar(2).open.max(self.window.bar(2).close))
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == -1
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                > self.body_short_avg.val(0)
            && self.window.bar(0).close
                < self.window.bar(2).close
                    - real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration
        {
            return -100.0;
        }
        0.0
    }
}

/// CDLMORNINGDOJISTAR
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLMORNINGDOJISTAR {
    bars_seen: usize,
    window: CandleWindow,
    body_doji_avg: RollingCandleAvg,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDLMORNINGDOJISTAR {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_doji_avg: RollingCandleAvg::new(BODY_DOJI),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDLMORNINGDOJISTAR {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLMORNINGDOJISTAR {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_doji_avg.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen
            <= BODY_DOJI
                .avg_period
                .max(BODY_LONG.avg_period)
                .max(BODY_SHORT.avg_period)
                + 2
        {
            return 0.0;
        }

        let penetration = 0.3;
        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && real_body(self.window.bar(2).open, self.window.bar(2).close)
                > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_doji_avg.val(1)
            && (self.window.bar(1).open.max(self.window.bar(1).close)
                < self.window.bar(2).open.min(self.window.bar(2).close))
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == 1
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                > self.body_short_avg.val(0)
            && self.window.bar(0).close
                > self.window.bar(2).close
                    + real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration
        {
            return 100.0;
        }
        0.0
    }
}

/// CDLMORNINGSTAR
#[derive(Debug, Clone)]
#[allow(non_camel_case_types)]
pub struct CDLMORNINGSTAR {
    bars_seen: usize,
    window: CandleWindow,
    body_long_avg: RollingCandleAvg,
    body_short_avg: RollingCandleAvg,
}

impl CDLMORNINGSTAR {
    pub fn new() -> Self {
        Self {
            bars_seen: 0,
            window: CandleWindow::new(3),
            body_long_avg: RollingCandleAvg::new(BODY_LONG),
            body_short_avg: RollingCandleAvg::new(BODY_SHORT),
        }
    }
}

impl Default for CDLMORNINGSTAR {
    fn default() -> Self {
        Self::new()
    }
}

impl Next<(f64, f64, f64, f64)> for CDLMORNINGSTAR {
    type Output = f64;

    fn next(&mut self, (open, high, low, close): (f64, f64, f64, f64)) -> Self::Output {
        self.bars_seen += 1;
        self.window.push(open, high, low, close);
        self.body_long_avg.push(open, high, low, close);
        self.body_short_avg.push(open, high, low, close);

        if self.bars_seen <= BODY_LONG.avg_period.max(BODY_SHORT.avg_period) + 2 {
            return 0.0;
        }

        let penetration = 0.3;
        if candle_color(self.window.bar(2).open, self.window.bar(2).close) == -1
            && real_body(self.window.bar(2).open, self.window.bar(2).close)
                > self.body_long_avg.val(2)
            && real_body(self.window.bar(1).open, self.window.bar(1).close)
                <= self.body_short_avg.val(1)
            && (self.window.bar(1).open.max(self.window.bar(1).close)
                < self.window.bar(2).open.min(self.window.bar(2).close))
            && candle_color(self.window.bar(0).open, self.window.bar(0).close) == 1
            && real_body(self.window.bar(0).open, self.window.bar(0).close)
                > self.body_short_avg.val(0)
            && self.window.bar(0).close
                > self.window.bar(2).close
                    + real_body(self.window.bar(2).open, self.window.bar(2).close) * penetration
        {
            return 100.0;
        }
        0.0
    }
}

#[cfg(test)]
mod tests {
    use super::*;
    use crate::indicators::patterns::fixtures as fx;
    // talib-rs's cdl_3blackcrows is defective; the native implementation is correct
    // and is deliberately NOT made bug-compatible. Evidence:
    //
    //  * talib-rs seeds `shadow_sum[k]` at `bar_offset = start - 3 + k` but consumes
    //    that slot as the average for bar `i - 2 + k` — off by one — and then
    //    adds/subtracts against the skewed set, so the window becomes non-contiguous
    //    and corresponds to no principled rolling average.
    //  * Over 200,000 bars: a from-scratch model of C TA-Lib semantics gives 8
    //    mismatches against the oracle; a one-bar-shifted model gives 11; the native
    //    implementation gives exactly the same 8 as the C-correct model. So
    //    native == C TA-Lib and the oracle is the outlier.
    //
    // Measured divergence: 8 bars / 200,000 on the deterministic walk, and never
    // more than 1 per run over 40,000 random walks of length <= 200. Bound of 12
    // leaves headroom while staying well below the oracle's 33 non-zero signals.
    crate::test_pattern_parity!(
        test_cdl3blackcrows_parity,
        CDL3BLACKCROWS,
        talib_rs::pattern::cdl_3blackcrows,
        fx::three_black_crows,
        oracle_exempt = "talib-rs's cdl_3blackcrows seeds shadow_sum[k] at bar `start-3+k` but \
                         uses it for bar `i-2+k`, then rolls that skewed set, producing a \
                         non-contiguous window matching no principled rolling average. A \
                         from-scratch C-TA-Lib-semantics model and the native implementation \
                         produce the identical 8 mismatches vs the oracle over 200,000 bars; a \
                         shifted model produces 11. Native == C TA-Lib; the oracle is wrong.",
        max_mismatches = 12,
        reference_mismatches = 8
    );
    crate::test_pattern_parity!(
        test_cdl3inside_parity,
        CDL3INSIDE,
        talib_rs::pattern::cdl_3inside,
        fx::three_inside
    );
    crate::test_pattern_parity!(
        test_cdl3outside_parity,
        CDL3OUTSIDE,
        talib_rs::pattern::cdl_3outside
    );
    crate::test_pattern_parity!(
        test_cdl3starsinsouth_parity,
        CDL3STARSINSOUTH,
        talib_rs::pattern::cdl_3starsinsouth,
        fx::three_stars_in_south
    );
    crate::test_pattern_parity!(
        test_cdl3whitesoldiers_parity,
        CDL3WHITESOLDIERS,
        talib_rs::pattern::cdl_3whitesoldiers
    );
    crate::test_pattern_parity!(
        test_cdlabandonedbaby_parity,
        CDLABANDONEDBABY,
        talib_rs::pattern::cdl_abandonedbaby,
        fx::abandonedbaby
    );
    crate::test_pattern_parity!(
        test_cdladvanceblock_parity,
        CDLADVANCEBLOCK,
        talib_rs::pattern::cdl_advanceblock
    );
    crate::test_pattern_parity!(
        test_cdleveningdojistar_parity,
        CDLEVENINGDOJISTAR,
        talib_rs::pattern::cdl_eveningdojistar,
        fx::eveningdojistar
    );
    crate::test_pattern_parity!(
        test_cdleveningstar_parity,
        CDLEVENINGSTAR,
        talib_rs::pattern::cdl_eveningstar,
        fx::eveningstar
    );
    crate::test_pattern_parity!(
        test_cdlmorningdojistar_parity,
        CDLMORNINGDOJISTAR,
        talib_rs::pattern::cdl_morningdojistar,
        fx::morningdojistar
    );
    crate::test_pattern_parity!(
        test_cdlmorningstar_parity,
        CDLMORNINGSTAR,
        talib_rs::pattern::cdl_morningstar,
        fx::morningstar
    );
}