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
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
use std::mem;
use std::rc::Rc;
use {
ArrayBase,
Dimension,
ViewRepr,
};
pub unsafe trait Data : Sized {
type Elem;
#[doc(hidden)]
fn _data_slice(&self) -> &[Self::Elem];
}
pub unsafe trait DataMut : Data {
#[doc(hidden)]
#[inline]
fn ensure_unique<D>(&mut ArrayBase<Self, D>)
where Self: Sized,
D: Dimension
{ }
#[doc(hidden)]
#[inline]
fn is_unique(&mut self) -> bool {
true
}
}
pub unsafe trait DataClone : Data {
#[doc(hidden)]
unsafe fn clone_with_ptr(&self, ptr: *mut Self::Elem) -> (Self, *mut Self::Elem);
}
unsafe impl<A> Data for Rc<Vec<A>> {
type Elem = A;
fn _data_slice(&self) -> &[A] {
self
}
}
unsafe impl<A> DataMut for Rc<Vec<A>>
where A: Clone
{
fn ensure_unique<D>(self_: &mut ArrayBase<Self, D>)
where Self: Sized,
D: Dimension
{
if Rc::get_mut(&mut self_.data).is_some() {
return;
}
if self_.dim.size() <= self_.data.len() / 2 {
unsafe {
*self_ = ArrayBase::from_shape_vec_unchecked(self_.dim.clone(),
self_.iter()
.cloned()
.collect());
}
return;
}
let our_off = (self_.ptr as isize - self_.data.as_ptr() as isize) /
mem::size_of::<A>() as isize;
let rvec = Rc::make_mut(&mut self_.data);
unsafe {
self_.ptr = rvec.as_mut_ptr().offset(our_off);
}
}
fn is_unique(&mut self) -> bool {
Rc::get_mut(self).is_some()
}
}
unsafe impl<A> DataClone for Rc<Vec<A>> {
unsafe fn clone_with_ptr(&self, ptr: *mut Self::Elem) -> (Self, *mut Self::Elem) {
(self.clone(), ptr)
}
}
unsafe impl<A> Data for Vec<A> {
type Elem = A;
fn _data_slice(&self) -> &[A] {
self
}
}
unsafe impl<A> DataMut for Vec<A> { }
unsafe impl<A> DataClone for Vec<A>
where A: Clone
{
unsafe fn clone_with_ptr(&self, ptr: *mut Self::Elem) -> (Self, *mut Self::Elem) {
let mut u = self.clone();
let our_off = (self.as_ptr() as isize - ptr as isize) /
mem::size_of::<A>() as isize;
let new_ptr = u.as_mut_ptr().offset(our_off);
(u, new_ptr)
}
}
unsafe impl<'a, A> Data for ViewRepr<&'a A> {
type Elem = A;
fn _data_slice(&self) -> &[A] {
&[]
}
}
unsafe impl<'a, A> DataClone for ViewRepr<&'a A> {
unsafe fn clone_with_ptr(&self, ptr: *mut Self::Elem) -> (Self, *mut Self::Elem) {
(*self, ptr)
}
}
unsafe impl<'a, A> Data for ViewRepr<&'a mut A> {
type Elem = A;
fn _data_slice(&self) -> &[A] {
&[]
}
}
unsafe impl<'a, A> DataMut for ViewRepr<&'a mut A> { }
pub unsafe trait DataOwned : Data {
#[doc(hidden)]
fn new(elements: Vec<Self::Elem>) -> Self;
#[doc(hidden)]
fn into_shared(self) -> Rc<Vec<Self::Elem>>;
}
pub unsafe trait DataShared : Clone + DataClone { }
unsafe impl<A> DataShared for Rc<Vec<A>> {}
unsafe impl<'a, A> DataShared for ViewRepr<&'a A> {}
unsafe impl<A> DataOwned for Vec<A> {
fn new(elements: Vec<A>) -> Self {
elements
}
fn into_shared(self) -> Rc<Vec<A>> {
Rc::new(self)
}
}
unsafe impl<A> DataOwned for Rc<Vec<A>> {
fn new(elements: Vec<A>) -> Self {
Rc::new(elements)
}
fn into_shared(self) -> Rc<Vec<A>> {
self
}
}