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
73
74
75
76
77
78
79
80
81
82
83
84
85
86

//! `OnlineExactSum` for dot product

use std::ops::Add;

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

use super::traits::DotAccumulator;
use crate::sum::traits::SumAccumulator;
use crate::sum::OnlineExactSum;
use crate::util::traits::TwoProduct;
use crate::util::two_product;

/// Calculates the dot product using product transformation and `OnlineExactSum`
///
/// ![][OnlineExactDot]
///
/// # Examples
///
/// ```
/// use accurate::traits::*;
/// use accurate::dot::OnlineExactDot;
///
/// let d = OnlineExactDot::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("OnlineExactDot", "images/OnlineExactDot.svg"))]
#[derive(Clone, Debug)]
pub struct OnlineExactDot<F> {
    s: OnlineExactSum<F>,
}

impl<F> DotAccumulator<F> for OnlineExactDot<F>
where
    F: TwoProduct,
    OnlineExactSum<F>: SumAccumulator<F>,
{
    fn zero() -> Self {
        OnlineExactDot::from(F::zero())
    }

    #[inline]
    fn dot(self) -> F {
        self.s.sum()
    }
}

impl<F> Add<(F, F)> for OnlineExactDot<F>
where
    F: TwoProduct,
    OnlineExactSum<F>: SumAccumulator<F>,
{
    type Output = Self;

    #[inline]
    fn add(mut self, (a, b): (F, F)) -> Self::Output {
        let (h, r1) = two_product(a, b);
        self.s = (self.s + h) + r1;
        self
    }
}

impl<F> From<F> for OnlineExactDot<F>
where
    OnlineExactSum<F>: SumAccumulator<F>,
{
    fn from(x: F) -> Self {
        OnlineExactDot {
            s: OnlineExactSum::from(x),
        }
    }
}

impl<F> Add for OnlineExactDot<F>
where
    OnlineExactSum<F>: Add<Output = OnlineExactSum<F>>,
{
    type Output = Self;

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

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