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
119
120

pub type CIResult<T> = Result<T, CIError>;

pub use crate::interval::IntervalError;

use num_traits::Float;

#[derive(thiserror::Error, Debug)]
pub enum CIError {
    #[error("Too few samples to compute: {0}")]
    TooFewSamples(usize),

    #[error("Too few successes: {0} (population: {1}; n*p={2}))")]
    TooFewSuccesses(usize, usize, f64),

    #[error("Too few failures: {0} (population: {1}; n*q={2}))")]
    TooFewFailures(usize, usize, f64),

    #[error("Invalid confidence level (must be )): {0}")]
    InvalidConfidenceLevel(f64),

    #[error("Invalid quantile (must be in (0, 1)): {0}")]
    InvalidQuantile(f64),

    #[error("Invalid number of successes: {0} (population: {1})")]
    InvalidSuccesses(usize, usize),

    #[error("Geometric/harmonic mean require strictly positive values: found {0}")]
    NonPositiveValue(f64),

    #[error("Invalid input data found")]
    InvalidInputData,

    #[error("Float type conversion error: {0}")]
    FloatConversionError(String),

    #[error("Index error: {0} should be in [0, {1})")]
    IndexError(f64, usize),

    // wrapper errors
    #[error("String error: {0}")]
    Error(String),

    #[error("Interval error: {0}")]
    IntervalError(#[from] IntervalError),
}

impl From<&str> for CIError {
    fn from(s: &str) -> Self {
        CIError::Error(s.to_string())
    }
}

#[derive(thiserror::Error, Debug)]
pub enum ConversionError {
    #[error("Empty interval has no concrete bounds")]
    NoConcreteBoundsError,

    #[error("Degenerate interval has single bound")]
    SingleBoundError,
}

///
/// Decorator trait used to convert from a generic [`Float`] type to a [`CIResult<f64>`]
///
pub(crate) trait FloatConversion<F: Float> {
    fn try_f64(&self, var_name: &str) -> CIResult<f64>;
}

impl<F: Float> FloatConversion<F> for F {
    #[inline]
    fn try_f64(&self, var_name: &str) -> CIResult<f64> {
        self.to_f64().ok_or_else(|| {
            CIError::FloatConversionError(format!(
                "Error converting {} ({}) to f64",
                var_name,
                std::any::type_name::<F>()
            ))
        })
    }
}

///
/// Decorator trait used to convert from an [`Option<F>`] to a [`CIResult<F>`]
///
pub(crate) trait FloatReverseConversion<F: num_traits::Float> {
    fn convert(&self, var_name: &str) -> CIResult<F>;
}

impl<F: Float> FloatReverseConversion<F> for Option<F> {
    #[inline]
    fn convert(&self, var_name: &str) -> CIResult<F> {
        self.ok_or_else(|| {
            CIError::FloatConversionError(format!(
                "Error converting {} to {}",
                var_name,
                std::any::type_name::<F>()
            ))
        })
    }
}

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

    fn string_to_error() -> CIResult<()> {
        Err("This is a string error")?
    }

    #[test]
    fn test_string_to_error() {
        let err = string_to_error();
        assert!(err.is_err());
        match err {
            Err(CIError::Error(s)) => assert_eq!(s, "This is a string error"),
            Err(e) => panic!("Unexpected error type: {:?}", e),
            Ok(_) => panic!("Unexpected success"),
        }
    }
}