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
use crate::vector::ops::{Asum, Axpy, Copy, Dot, Iamax, Nrm2, Scal};
use num::traits::NumCast;
use num_complex::{Complex32, Complex64};
pub mod ll;
pub mod ops;
pub trait Vector<T> {
fn inc(&self) -> u32 {
1
}
fn len(&self) -> u32;
fn as_ptr(&self) -> *const T;
fn as_mut_ptr(&mut self) -> *mut T;
fn is_empty(&self) -> bool {
self.len() == 0u32
}
}
impl<'a, T> Into<Vec<T>> for &'a dyn 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,
{
fn update(&mut self, alpha: &T, x: &dyn Vector<T>) -> &mut Self {
Axpy::axpy(alpha, x, self);
self
}
fn scale(&mut self, alpha: &T) -> &mut Self {
Scal::scal(alpha, self);
self
}
fn dot(&self, x: &dyn Vector<T>) -> T {
Dot::dot(self, x)
}
fn abs_sum(&self) -> T {
Asum::asum(self)
}
fn norm(&self) -> T {
Nrm2::nrm2(self)
}
fn max_index(&self) -> usize {
Iamax::iamax(self)
}
}
impl<T> Vector<T> for Vec<T> {
fn len(&self) -> u32 {
let l: Option<u32> = NumCast::from(Vec::len(self));
match l {
Some(l) => l,
None => panic!(),
}
}
fn as_ptr(&self) -> *const T {
self[..].as_ptr()
}
fn as_mut_ptr(&mut self) -> *mut T {
(&mut self[..]).as_mut_ptr()
}
}
impl<T> Vector<T> for [T] {
fn len(&self) -> u32 {
let l: Option<u32> = NumCast::from(<[T]>::len(self));
match l {
Some(l) => l,
None => panic!(),
}
}
fn as_ptr(&self) -> *const T {
<[T]>::as_ptr(self)
}
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);