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
use num::traits::NumCast;
use num::complex::{Complex32, Complex64};
use vector::ops::{Copy, Axpy, Scal, Dot, Nrm2, Asum, Iamax};
pub mod ll;
pub mod ops;
pub trait Vector<T> {
fn inc(&self) -> i32 { 1 }
fn len(&self) -> i32;
unsafe fn as_ptr(&self) -> *const T;
unsafe fn as_mut_ptr(&mut self) -> *mut T;
}
impl<'a, T> Into<Vec<T>> for &'a Vector<T>
where T: Copy {
fn into(self) -> Vec<T> {
let n = self.len() as usize;
let mut x = Vec::with_capacity(n);
unsafe { x.set_len(n); }
Copy::copy(self, &mut x);
x
}
}
pub trait VectorOperations<T>: Sized + Vector<T>
where T: Copy + Axpy + Scal + Dot + Nrm2 + Asum + Iamax {
#[inline]
fn update(&mut self, alpha: &T, x: &Vector<T>) -> &mut Self {
Axpy::axpy(alpha, x, self);
self
}
#[inline]
fn scale(&mut self, alpha: &T) -> &mut Self {
Scal::scal(alpha, self);
self
}
#[inline]
fn dot(&self, x: &Vector<T>) -> T {
Dot::dot(self, x)
}
#[inline]
fn abs_sum(&self) -> T {
Asum::asum(self)
}
#[inline]
fn norm(&self) -> T {
Nrm2::nrm2(self)
}
#[inline]
fn max_index(&self) -> usize {
Iamax::iamax(self)
}
}
impl<T> Vector<T> for Vec<T> {
#[inline]
fn len(&self) -> i32 {
let l: Option<i32> = NumCast::from(Vec::len(self));
match l {
Some(l) => l,
None => panic!(),
}
}
#[inline]
unsafe fn as_ptr(&self) -> *const T { self[..].as_ptr() }
#[inline]
unsafe fn as_mut_ptr(&mut self) -> *mut T { (&mut self[..]).as_mut_ptr() }
}
impl<T> Vector<T> for [T] {
#[inline]
fn len(&self) -> i32 {
let l: Option<i32> = NumCast::from(<[T]>::len(self));
match l {
Some(l) => l,
None => panic!(),
}
}
#[inline]
unsafe fn as_ptr(&self) -> *const T { <[T]>::as_ptr(self) }
#[inline]
unsafe fn as_mut_ptr(&mut self) -> *mut T { <[T]>::as_mut_ptr(self) }
}
macro_rules! operations_impl(
($v: ident, $($t: ty), +) => (
$( impl VectorOperations<$t> for $v<$t> {} )+
)
);
operations_impl!(Vec, f32, f64, Complex32, Complex64);