[go: up one dir, main page]

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
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106

//! Algorithms for summation

pub mod traits;

mod cascaded;
mod ifastsum;
mod kahan;
mod naive;
mod onlineexactsum;
mod sumk;

pub use self::cascaded::{Klein, Neumaier};
pub use self::ifastsum::i_fast_sum_in_place;
pub use self::kahan::Kahan;
pub use self::naive::NaiveSum;
pub use self::onlineexactsum::OnlineExactSum;
pub use self::sumk::{Sum2, Sum3, Sum4, Sum5, Sum6, Sum7, Sum8, Sum9, SumK};

#[cfg(feature = "parallel")]
use num_traits::Zero;

#[cfg(feature = "parallel")]
use rayon::iter::plumbing::{Consumer, Folder, UnindexedConsumer};

#[cfg(feature = "parallel")]
use self::traits::ParallelSumAccumulator;
#[cfg(feature = "parallel")]
use self::traits::SumAccumulator;
#[cfg(feature = "parallel")]
use crate::util::AddReducer;

/// Adapts a `SumAccumulator` into a `Folder`
#[cfg(feature = "parallel")]
#[derive(Copy, Clone, Debug)]
pub struct SumFolder<Acc>(Acc);

#[cfg(feature = "parallel")]
impl<Acc, F> Folder<F> for SumFolder<Acc>
where
    Acc: SumAccumulator<F>,
{
    type Result = Acc;

    #[inline]
    fn consume(self, item: F) -> Self {
        SumFolder(self.0 + item)
    }

    #[inline]
    fn complete(self) -> Self::Result {
        self.0
    }

    #[inline]
    fn full(&self) -> bool {
        false
    }
}

/// Adapts a `ParallelSumAccumulator` into a `Consumer`
#[cfg(feature = "parallel")]
#[derive(Copy, Clone, Debug)]
pub struct SumConsumer<Acc>(Acc);

#[cfg(feature = "parallel")]
impl<Acc, F> Consumer<F> for SumConsumer<Acc>
where
    Acc: ParallelSumAccumulator<F>,
    F: Zero + Send,
{
    type Folder = SumFolder<Acc>;
    type Reducer = AddReducer;
    type Result = Acc;

    #[inline]
    fn split_at(self, _index: usize) -> (Self, Self, Self::Reducer) {
        (self, Acc::zero().into_consumer(), AddReducer)
    }

    #[inline]
    fn into_folder(self) -> Self::Folder {
        SumFolder(self.0)
    }

    #[inline]
    fn full(&self) -> bool {
        false
    }
}

#[cfg(feature = "parallel")]
impl<Acc, F> UnindexedConsumer<F> for SumConsumer<Acc>
where
    Acc: ParallelSumAccumulator<F>,
    F: Zero + Send,
{
    #[inline]
    fn split_off_left(&self) -> Self {
        Acc::zero().into_consumer()
    }

    #[inline]
    fn to_reducer(&self) -> Self::Reducer {
        AddReducer
    }
}