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
use libnum::{Zero, One, Float};
use std::any::Any;
use std::fmt;
use std::ops::{Add, Sub, Mul, Div};
#[cfg(feature="assign_ops")]
use std::ops::{
AddAssign,
SubAssign,
MulAssign,
DivAssign,
RemAssign,
};
use ScalarOperand;
#[cfg(feature="rblas")]
use std::any::TypeId;
#[cfg(feature="rblas")]
use ShapeError;
#[cfg(feature="rblas")]
use error::{from_kind, ErrorKind};
#[cfg(feature="rblas")]
use blas::{AsBlas, BlasArrayView};
#[cfg(feature="rblas")]
use imp_prelude::*;
pub trait LinalgScalar :
Any +
Copy +
Zero + One +
Add<Output=Self> +
Sub<Output=Self> +
Mul<Output=Self> +
Div<Output=Self>
{ }
impl<T> LinalgScalar for T
where T:
Any +
Copy +
Zero + One +
Add<Output=T> +
Sub<Output=T> +
Mul<Output=T> +
Div<Output=T>
{ }
#[cfg(not(feature="assign_ops"))]
pub trait NdFloat :
Float +
fmt::Display + fmt::Debug + fmt::LowerExp + fmt::UpperExp +
ScalarOperand + LinalgScalar
{ }
#[cfg(feature="assign_ops")]
pub trait NdFloat :
Float +
AddAssign + SubAssign + MulAssign + DivAssign + RemAssign +
fmt::Display + fmt::Debug + fmt::LowerExp + fmt::UpperExp +
ScalarOperand + LinalgScalar
{ }
impl NdFloat for f32 { }
impl NdFloat for f64 { }
#[cfg(feature = "rblas")]
pub trait AsBlasAny<A, S, D> : AsBlas<A, S, D> {
fn blas_view_as_type<T: Any>(&self) -> Result<BlasArrayView<T, D>, ShapeError>
where S: Data;
}
#[cfg(feature = "rblas")]
impl<A, S, D> AsBlasAny<A, S, D> for ArrayBase<S, D>
where S: Data<Elem=A>,
D: Dimension,
A: Any,
{
fn blas_view_as_type<T: Any>(&self) -> Result<BlasArrayView<T, D>, ShapeError>
where S: Data
{
if TypeId::of::<A>() == TypeId::of::<T>() {
unsafe {
let v = self.view();
let u = ArrayView::new_(v.ptr as *const T, v.dim, v.strides);
Priv(u).into_blas_view()
}
} else {
Err(from_kind(ErrorKind::IncompatibleLayout))
}
}
}