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

use super::misc::ToFloat;
use std::ops::{Add, Sub, Div};

/// An iterator of a sequence of evenly spaced floats.
///
/// Iterator element type is `F`.
pub struct Linspace<F> {
    start: F,
    step: F,
    len: usize,
}

impl<F> Iterator for Linspace<F> where
    F: Copy + Add<Output=F>,
{
    type Item = F;

    #[inline]
    fn next(&mut self) -> Option<F>
    {
        if self.len == 0 {
            None
        } else {
            self.len -= 1;
            let elt = self.start;
            self.start = self.start + self.step;
            Some(elt)
        }
    }
}

/// Return an iterator with `n` elements, where the first
/// element is `a` and the last element is `b`.
///
/// Iterator element type is `F`.
///
/// ```
/// use itertools as it;
/// let mut xs = it::linspace::<f32>(0., 1., 5);
/// assert_eq!(xs.collect::<Vec<_>>(),
///            vec![0., 0.25, 0.5, 0.75, 1.0]);
/// ```
#[inline]
pub fn linspace<F>(a: F, b: F, n: usize) -> Linspace<F> where
    F: Copy + Sub<Output=F> + Div<Output=F>,
    usize: ToFloat<F>,
{
    let nf: F = n.to_float();
    let step = (b - a)/(nf - 1usize.to_float());
    Linspace {
        start: a,
        step: step,
        len: n,
    }
}