ha-ndarray 0.5.0

A hardware-accelerated n-dimensional array
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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154

use std::fmt;

use crate::Error;

use super::{CLElement, OpenCL, TILE_SIZE, WG_SIZE};

pub mod constructors;
pub mod elementwise;
pub mod gather;
pub mod linalg;
pub mod reduce;
pub mod slice;
pub mod view;

#[derive(Clone, Eq, PartialEq, Hash, fmt::Debug)]
pub(super) enum CLExpr {
    Static(&'static str),
    String(String),
}

impl From<&'static str> for CLExpr {
    fn from(s: &'static str) -> Self {
        CLExpr::Static(s)
    }
}

impl From<String> for CLExpr {
    fn from(s: String) -> Self {
        CLExpr::String(s)
    }
}

impl fmt::Display for CLExpr {
    fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result {
        match self {
            Self::Static(s) => f.write_str(s),
            Self::String(s) => f.write_str(s),
        }
    }
}

pub trait Builder {
    fn build(self) -> String;
}

#[derive(Clone, Eq, PartialEq, Hash, fmt::Debug)]
pub struct ElementDual {
    pub(super) i_type: &'static str,
    pub(super) o_type: &'static str,
    pub(super) name: &'static str,
    pub(super) op: CLExpr,
}

impl ElementDual {
    pub(super) fn new<I, O, Op>(name: &'static str, op: Op) -> Self
    where
        I: CLElement,
        O: CLElement,
        Op: Into<CLExpr>,
    {
        Self {
            i_type: I::TYPE,
            o_type: O::TYPE,
            name,
            op: op.into(),
        }
    }
}

impl Builder for ElementDual {
    fn build(self) -> String {
        format!(
            r#"
            inline {o_type} {name}(const {i_type} lhs, const {i_type} rhs) {{
                {op}
            }}
            "#,
            i_type = self.i_type,
            o_type = self.o_type,
            name = self.name,
            op = self.op
        )
    }
}

#[derive(Clone, Eq, PartialEq, Hash, fmt::Debug)]
pub struct ElementUnary {
    pub(super) i_type: &'static str,
    pub(super) o_type: &'static str,
    pub(super) name: &'static str,
    pub(super) op: CLExpr,
}

impl ElementUnary {
    pub(super) fn new<I, O, Op>(name: &'static str, op: Op) -> Self
    where
        I: CLElement,
        O: CLElement,
        Op: Into<CLExpr>,
    {
        Self {
            i_type: I::TYPE,
            o_type: O::TYPE,
            name,
            op: op.into(),
        }
    }
}

impl Builder for ElementUnary {
    fn build(self) -> String {
        format!(
            "
            inline {o_type} {name}(const {i_type} n) {{
                {op}
            }}
            ",
            i_type = self.i_type,
            o_type = self.o_type,
            name = self.name,
            op = self.op,
        )
    }
}

struct ArrayFormat<'a, T> {
    arr: &'a [T],
}

impl<'a, T> From<&'a [T]> for ArrayFormat<'a, T> {
    fn from(arr: &'a [T]) -> Self {
        Self { arr }
    }
}

impl<'a, T: fmt::Display> fmt::Display for ArrayFormat<'a, T> {
    fn fmt(&self, f: &mut fmt::Formatter) -> fmt::Result {
        f.write_str("{ ")?;

        for item in self.arr {
            write!(f, "{item}, ")?;
        }

        f.write_str(" }")
    }
}

#[inline]
fn build(src: &str) -> Result<ocl::Program, Error> {
    ocl::Program::builder()
        .source(src)
        .build(OpenCL::context())
        .map_err(Error::from)
}