accurate 0.4.1

(more or less) accurate floating point algorithms
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

//! Naive floating point dot product

use std::ops::Add;

#[cfg(doc)]
use embed_doc_image::embed_doc_image;

use num_traits::Float;

use super::traits::DotAccumulator;

/// Naive floating point dot product
///
/// ![][NaiveDot]
///
/// # Examples
///
/// ```
/// use accurate::traits::*;
/// use accurate::dot::NaiveDot;
///
/// let d = NaiveDot::zero() + (1.0, 1.0) + (2.0, 2.0) + (3.0, 3.0);
/// assert_eq!(14.0f64, d.dot());
/// ```
#[cfg_attr(doc, embed_doc_image("NaiveDot", "images/NaiveDot.svg"))]
#[derive(Copy, Clone, Debug)]
pub struct NaiveDot<F>(F);

impl<F> DotAccumulator<F> for NaiveDot<F>
where
    F: Float,
{
    #[inline]
    fn dot(self) -> F {
        self.0
    }
}

impl<F> Add<(F, F)> for NaiveDot<F>
where
    F: Float,
{
    type Output = Self;

    #[inline]
    fn add(self, rhs: (F, F)) -> Self::Output {
        NaiveDot(self.0 + rhs.0 * rhs.1)
    }
}

impl<F> From<F> for NaiveDot<F>
where
    F: Float,
{
    fn from(x: F) -> Self {
        NaiveDot(x)
    }
}

impl<F> Add for NaiveDot<F>
where
    F: Float,
{
    type Output = Self;

    #[inline]
    fn add(self, rhs: Self) -> Self::Output {
        NaiveDot(self.0 + rhs.0)
    }
}

unsafe impl<F> Send for NaiveDot<F> where F: Send {}