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};
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)
}
}
}
#[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,
}
}