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

use array::Array;

/// Returns evenly spaced elements from given interval
///
/// # Arguments
///
/// * `start` - start of interval
/// * `stop` - end of interval
/// * `num` - number of elements
///
/// # Examples
///
/// ```
///  use numas::array::Array;
///
///  let array = numas::factory::numeric::linspace(1.0, 5.0, 9);
///  let data = array.collect();
///
///  assert_eq!(data, vec![1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0]);
/// ```
pub fn linspace(start: f64, stop: f64, num: usize) -> Array<f64> {
    let mut data: Vec<f64> = Vec::with_capacity(num);

    let step: f64 = (stop - start) / (num - 1) as f64;
    let mut tmp = start;

    for _ in 0..num {
        data.push(tmp);
        tmp += step;
    }

    return Array::new(data, vec![num as i32]);
}

/// Returns evenly spaced interval within a given interval
///
/// # Arguments
///
/// * `start` - start of interval
/// * `stop` - end of interval
/// * `step` - step of interval
///
/// # Examples
///
/// ```
///  use numas::array::Array;
///
///  let array = numas::factory::numeric::arange(1.0, 5.0, 1.5);
///  let data = array.collect();
///
///  assert_eq!(data, vec![1.0, 2.5, 4.0]);
/// ```
pub fn arange(start: f64, stop: f64, step: f64) -> Array<f64> {
    let num = f64::ceil((stop - start) / step) as i32;
    let mut data: Vec<f64> = Vec::with_capacity(num as usize);
    let mut tmp = start;

    for _ in 0..num {
        data.push(tmp);
        tmp = tmp + step;
    }

    return Array::new(data, vec![num]);
}

/// Returns evenly spaced elements on log space from given interval
///
/// # Arguments
///
/// * `start` - start of interval
/// * `stop` - end of interval
/// * `base` - base of log
/// * `num` - number of elements
pub fn logspace(start: f64, stop: f64, base: f64, num: usize) -> Array<f64> {
    let mut data: Vec<f64> = Vec::with_capacity(num);

    let step: f64 = (stop - start) / (num - 1) as f64;
    let mut tmp = start;

    for _ in 0..num {
        data.push(base.powf(tmp));
        tmp += step;
    }

    return Array::new(data, vec![num as i32]);
}