rlst 0.10.0

A Rust native linear algebra library.
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

//! Container representing application of a unary operator

use crate::{
    ContainerType,
    dense::array::{Array, Shape, UnsafeRandomAccessByValue},
    traits::{accessors::UnsafeRandom1DAccessByValue, base_operations::BaseItem},
};

/// Application of a unitary Operator
pub struct ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, const NDIM: usize>
where
    Op: Fn(OpItem) -> OpTarget,
{
    arr: Array<ArrayImpl, NDIM>,
    op: Op,
    _marker: std::marker::PhantomData<(OpItem, OpTarget)>,
}

impl<OpItem, OpTarget, ArrayImpl, Op: Fn(OpItem) -> OpTarget, const NDIM: usize>
    ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
{
    /// Create new
    pub fn new(arr: Array<ArrayImpl, NDIM>, op: Op) -> Self {
        Self {
            arr,
            op,
            _marker: std::marker::PhantomData,
        }
    }
}

impl<OpItem, OpTarget, ArrayImpl, Op: Fn(OpItem) -> OpTarget, const NDIM: usize> ContainerType
    for ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
where
    ArrayImpl: ContainerType,
{
    type Type = ArrayImpl::Type;
}

impl<OpItem, OpTarget, ArrayImpl, Op, const NDIM: usize> BaseItem
    for ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
where
    ArrayImpl: BaseItem<Item = OpItem>,
    OpTarget: Copy + Default,
    Op: Fn(OpItem) -> OpTarget,
{
    type Item = OpTarget;
}

impl<OpItem, OpTarget, ArrayImpl, Op, const NDIM: usize> UnsafeRandomAccessByValue<NDIM>
    for ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
where
    ArrayImpl: UnsafeRandomAccessByValue<NDIM, Item = OpItem>,
    OpTarget: Copy + Default,
    Op: Fn(OpItem) -> OpTarget,
{
    #[inline(always)]
    unsafe fn get_value_unchecked(&self, multi_index: [usize; NDIM]) -> Self::Item {
        unsafe { (self.op)(self.arr.get_value_unchecked(multi_index)) }
    }
}

impl<OpItem, OpTarget, ArrayImpl, Op, const NDIM: usize> UnsafeRandom1DAccessByValue
    for ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
where
    ArrayImpl: UnsafeRandom1DAccessByValue<Item = OpItem>,
    OpTarget: Copy + Default,
    Op: Fn(OpItem) -> OpTarget,
{
    #[inline(always)]
    unsafe fn get_value_1d_unchecked(&self, index: usize) -> Self::Item {
        unsafe { (self.op)(self.arr.imp().get_value_1d_unchecked(index)) }
    }
}

impl<OpItem, OpTarget, ArrayImpl, Op: Fn(OpItem) -> OpTarget, const NDIM: usize> Shape<NDIM>
    for ArrayUnaryOperator<OpItem, OpTarget, ArrayImpl, Op, NDIM>
where
    ArrayImpl: Shape<NDIM>,
{
    #[inline(always)]
    fn shape(&self) -> [usize; NDIM] {
        self.arr.shape()
    }
}