mantis-ta 0.6.5

Composable technical analysis and strategy engine for Rust
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118

use crate::indicators::Indicator;
use crate::types::Candle;
use crate::utils::ringbuf::RingBuf;

/// Exponential Moving Average over closing prices.
///
/// # Examples
/// ```rust
/// use mantis_ta::indicators::{Indicator, EMA};
/// use mantis_ta::types::Candle;
///
/// let prices = [1.0, 2.0, 3.0, 4.0];
/// let candles: Vec<Candle> = prices
///     .iter()
///     .enumerate()
///     .map(|(i, p)| Candle {
///         timestamp: i as i64,
///         open: *p,
///         high: *p,
///         low: *p,
///         close: *p,
///         volume: 0.0,
///     })
///     .collect();
///
/// let out = EMA::new(3).calculate(&candles);
/// // Warmup: first 2 bars None, third is seeded SMA
/// assert!(out[0].is_none());
/// assert!(out[1].is_none());
/// assert!(out[2].is_some());
/// ```
#[derive(Debug, Clone)]
pub struct EMA {
    period: usize,
    multiplier: f64,
    warmup: RingBuf<f64>,
    ema: Option<f64>,
}

impl EMA {
    pub fn new(period: usize) -> Self {
        assert!(period > 0, "period must be > 0");
        Self {
            period,
            multiplier: 2.0 / (period as f64 + 1.0),
            warmup: RingBuf::new(period, 0.0),
            ema: None,
        }
    }

    #[inline]
    fn update(&mut self, value: f64) -> Option<f64> {
        if let Some(prev) = self.ema {
            let next = (value - prev) * self.multiplier + prev;
            self.ema = Some(next);
            return self.ema;
        }

        // Warmup using simple average.
        self.warmup.push(value);
        if self.warmup.len() < self.period {
            return None;
        }
        let sum: f64 = self.warmup.iter().copied().sum();
        let sma = sum / self.period as f64;
        self.ema = Some(sma);
        self.ema
    }
}

impl Indicator for EMA {
    type Output = f64;

    fn next(&mut self, candle: &Candle) -> Option<Self::Output> {
        self.update(candle.close)
    }

    fn reset(&mut self) {
        self.ema = None;
        self.warmup = RingBuf::new(self.period, 0.0);
    }

    fn warmup_period(&self) -> usize {
        self.period
    }

    fn clone_boxed(&self) -> Box<dyn Indicator<Output = Self::Output>> {
        Box::new(self.clone())
    }
}

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

    #[test]
    fn computes_ema_after_warmup() {
        let mut ema = EMA::new(3);
        let prices = [1.0, 2.0, 3.0, 4.0];
        let candles: Vec<Candle> = prices
            .iter()
            .map(|p| Candle {
                timestamp: 0,
                open: *p,
                high: *p,
                low: *p,
                close: *p,
                volume: 0.0,
            })
            .collect();

        let outputs: Vec<_> = candles.iter().map(|c| ema.next(c)).collect();
        assert_eq!(outputs[0], None);
        assert_eq!(outputs[1], None);
        assert!(outputs[2].is_some());
        assert!(outputs[3].is_some());
    }
}