narr 0.0.5

Numeric Calc
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

use std::ops::{Sub, SubAssign};

use crate::arr::{
    core::{adaptative_result_shape, dim_expansion_for_shape, plain_index},
    shape_iterator::ShapeIterator,
    vec::core::ArrVec,
};

impl<T> ArrVec<T>
where
    T: Default + Clone + Copy + Sub<Output = T>,
{
    fn sub_to_left(lhs: &mut Self, rhs: &Self) {
        let out_shape = adaptative_result_shape(&lhs.shape, &rhs.shape);
        assert_eq!(lhs.shape, out_shape);
        let rhs_expanded_shape = dim_expansion_for_shape(&rhs.shape, out_shape.len());
        for index in ShapeIterator::new(&out_shape) {
            let out_i = plain_index(&out_shape, &out_shape, &index);
            let lhs_i = plain_index(&out_shape, &out_shape, &index);
            let rhs_i = plain_index(&rhs_expanded_shape, &out_shape, &index);
            lhs.buff[out_i] = lhs.buff[lhs_i] - rhs.buff[rhs_i];
        }
    }

    fn sub_to_new(lhs: &Self, rhs: &Self) -> Self {
        let out_shape = adaptative_result_shape(&lhs.shape, &rhs.shape);
        let mut out = ArrVec::<T>::new(T::default(), out_shape);
        Self::sub_to_out(lhs, rhs, &mut out);
        out
    }

    fn sub_to_out(lhs: &Self, rhs: &Self, out: &mut Self) {
        let lhs_expanded_shape = dim_expansion_for_shape(&lhs.shape, out.shape.len());
        let rhs_expanded_shape = dim_expansion_for_shape(&rhs.shape, out.shape.len());
        for index in ShapeIterator::new(&out.shape) {
            let out_i = plain_index(&out.shape, &out.shape, &index);
            let lhs_i = plain_index(&lhs_expanded_shape, &out.shape, &index);
            let rhs_i = plain_index(&rhs_expanded_shape, &out.shape, &index);
            out.buff[out_i] = lhs.buff[lhs_i] - rhs.buff[rhs_i];
        }
    }
}

impl<T> Sub for &ArrVec<T>
where
    T: Default + Clone + Copy + Sub<Output = T>,
{
    type Output = ArrVec<T>;

    fn sub(self, rhs: Self) -> Self::Output {
        ArrVec::<T>::sub_to_new(self, rhs)
    }
}

impl<T> Sub for ArrVec<T>
where
    T: Default + Clone + Copy + Sub<Output = T>,
{
    type Output = ArrVec<T>;

    fn sub(self, rhs: Self) -> Self::Output {
        ArrVec::<T>::sub_to_new(&self, &rhs)
    }
}

impl<T> SubAssign for ArrVec<T>
where
    T: Default + Clone + Copy + Sub<Output = T>,
{
    fn sub_assign(&mut self, rhs: Self) {
        ArrVec::<T>::sub_to_left(self, &rhs);
    }
}

impl<T> SubAssign<&ArrVec<T>> for ArrVec<T>
where
    T: Default + Clone + Copy + Sub<Output = T>,
{
    fn sub_assign(&mut self, rhs: &Self) {
        ArrVec::<T>::sub_to_left(self, rhs);
    }
}