tblis-src 0.2.5

Crate for TBLIS Source Code
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

#[cfg(test)]
mod tests {
    extern crate tblis_src;
    use tblis_ffi::tblis::*;

    #[test]
    fn test_workable() {
        let mut data_a = [0.0f64; 10 * 9 * 2 * 5];
        let mut data_b = [0.0f64; 7 * 5 * 9 * 8];
        let mut data_c = [0.0f64; 7 * 2 * 10 * 8];

        let mut len_a = [10, 9, 2, 5];
        let mut len_b = [7, 5, 9, 8];
        let mut len_c = [7, 2, 10, 8];
        let mut stride_a = [1, 10, 90, 180];
        let mut stride_b = [1, 7, 35, 315];
        let mut stride_c = [1, 7, 14, 140];

        let a = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 0.0 }, type_: TYPE_DOUBLE },
            data: data_a.as_mut_ptr() as *mut _,
            ndim: 4,
            len: len_a.as_mut_ptr(),
            stride: stride_a.as_mut_ptr(),
        };

        let b = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 0.0 }, type_: TYPE_DOUBLE },
            data: data_b.as_mut_ptr() as *mut _,
            ndim: 4,
            len: len_b.as_mut_ptr(),
            stride: stride_b.as_mut_ptr(),
        };

        let mut c = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 0.0 }, type_: TYPE_DOUBLE },
            data: data_c.as_mut_ptr() as *mut _,
            ndim: 4,
            len: len_c.as_mut_ptr(),
            stride: stride_c.as_mut_ptr(),
        };

        unsafe {
            tblis_tensor_mult(
                std::ptr::null_mut(),
                std::ptr::null_mut(),
                &a,
                c"cebf".as_ptr(),
                &b,
                c"afed".as_ptr(),
                &mut c,
                c"abcd".as_ptr(),
            );
        }

        /* Test case of the following C code.

        ```C
        double data_A[10*9*2*5];
        tblis_tensor A;
        tblis_init_tensor_d(&A, 4, (len_type[]){10, 9, 2, 5},
                            data_A, (stride_type[]){1, 10, 90, 180});

        double data_B[7*5*9*8];
        tblis_tensor B;
        tblis_init_tensor_d(&B, 4, (len_type[]){7, 5, 9, 8},
                            data_B, (stride_type[]){1, 7, 35, 315});

        double data_C[7*2*10*8];
        tblis_tensor C;
        tblis_init_tensor_d(&C, 4, (len_type[]){7, 2, 10, 8},
                            data_C, (stride_type[]){1, 7, 14, 140});

        // initialize data_A and data_B...

        // this computes C[abcd] += A[cebf] B[afed]
        tblis_tensor_mult(NULL, NULL, &A, "cebf", &B, "afed", &C, "abcd");
        ```

        */
    }

    #[test]
    fn test_matmul() {
        let data_a: Vec<f64> = vec![1.0, 2.0, 3.0, 4.0];
        let data_b: Vec<f64> = vec![5.0, 6.0, 7.0, 8.0];
        let mut data_c: Vec<f64> = vec![0.0; 4];

        let mut len = [2, 2];
        let mut stride = [2, 1];

        let a = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 1.0 }, type_: TYPE_DOUBLE },
            data: data_a.as_ptr() as *mut _,
            ndim: 2,
            len: len.as_mut_ptr(),
            stride: stride.as_mut_ptr(),
        };

        let b = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 1.0 }, type_: TYPE_DOUBLE },
            data: data_b.as_ptr() as *mut _,
            ndim: 2,
            len: len.as_mut_ptr(),
            stride: stride.as_mut_ptr(),
        };

        let mut c = tblis_tensor {
            type_: TYPE_DOUBLE,
            conj: 0,
            scalar: tblis_scalar { data: tblis_scalar_scalar { d: 1.0 }, type_: TYPE_DOUBLE },
            data: data_c.as_mut_ptr() as *mut _,
            ndim: 2,
            len: len.as_mut_ptr(),
            stride: stride.as_mut_ptr(),
        };

        unsafe {
            tblis_tensor_mult(
                std::ptr::null_mut(),
                std::ptr::null_mut(),
                &a,
                c"ij".as_ptr(),
                &b,
                c"jk".as_ptr(),
                &mut c,
                c"ik".as_ptr(),
            );
        }

        assert_eq!(data_c, vec![19.0, 22.0, 43.0, 50.0]);
    }
}