-
Notifications
You must be signed in to change notification settings - Fork 237
Expand file tree
/
Copy pathnumpy_backend.py
More file actions
290 lines (239 loc) · 10.6 KB
/
numpy_backend.py
File metadata and controls
290 lines (239 loc) · 10.6 KB
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
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
import warnings
from typing import Any, List, Optional, Tuple
import numpy as np
from docarray.computation import AbstractComputationalBackend
from docarray.computation.abstract_numpy_based_backend import AbstractNumpyBasedBackend
def _expand_if_single_axis(*matrices: np.ndarray) -> List[np.ndarray]:
"""Expands arrays that only have one axis, at dim 0.
This ensures that all outputs can be treated as matrices, not vectors.
:param matrices: Matrices to be expanded
:return: List of the input matrices,
where single axis matrices are expanded at dim 0.
"""
expanded = []
for m in matrices:
if len(m.shape) == 1:
expanded.append(np.expand_dims(m, axis=0))
else:
expanded.append(m)
return expanded
def _expand_if_scalar(arr: np.ndarray) -> np.ndarray:
if len(arr.shape) == 0: # avoid scalar output
arr = np.expand_dims(arr, axis=0)
return arr
def identity(array: np.ndarray) -> np.ndarray:
return array
class NumpyCompBackend(AbstractNumpyBasedBackend):
"""
Computational backend for Numpy.
"""
_module = np
_cast_output = identity
_get_tensor = identity
@classmethod
def to_device(cls, tensor: 'np.ndarray', device: str) -> 'np.ndarray':
"""Move the tensor to the specified device."""
raise NotImplementedError('Numpy does not support devices (GPU).')
@classmethod
def device(cls, tensor: 'np.ndarray') -> Optional[str]:
"""Return device on which the tensor is allocated."""
return None
@classmethod
def to_numpy(cls, array: 'np.ndarray') -> 'np.ndarray':
return array
@classmethod
def none_value(cls) -> Any:
"""Provide a compatible value that represents None in numpy."""
return None
@classmethod
def detach(cls, tensor: 'np.ndarray') -> 'np.ndarray':
"""
Returns the tensor detached from its current graph.
:param tensor: tensor to be detached
:return: a detached tensor with the same data.
"""
return tensor
@classmethod
def dtype(cls, tensor: 'np.ndarray') -> np.dtype:
"""Get the data type of the tensor."""
return tensor.dtype
@classmethod
def minmax_normalize(
cls,
tensor: 'np.ndarray',
t_range: Tuple = (0, 1),
x_range: Optional[Tuple] = None,
eps: float = 1e-7,
) -> 'np.ndarray':
"""
Normalize values in `tensor` into `t_range`.
`tensor` can be a 1D array or a 2D array. When `tensor` is a 2D array, then
normalization is row-based.
!!! note
- with `t_range=(0, 1)` will normalize the min-value of data to 0, max to 1;
- with `t_range=(1, 0)` will normalize the min-value of data to 1, max value
of the data to 0.
:param tensor: the data to be normalized
:param t_range: a tuple represents the target range.
:param x_range: a tuple represents tensors range.
:param eps: a small jitter to avoid divide by zero
:return: normalized data in `t_range`
"""
a, b = t_range
min_d = x_range[0] if x_range else np.min(tensor, axis=-1, keepdims=True)
max_d = x_range[1] if x_range else np.max(tensor, axis=-1, keepdims=True)
r = (b - a) * (tensor - min_d) / (max_d - min_d + eps) + a
return np.clip(r, *((a, b) if a < b else (b, a)))
@classmethod
def equal(cls, tensor1: 'np.ndarray', tensor2: 'np.ndarray') -> bool:
"""
Check if two tensors are equal.
:param tensor1: the first array
:param tensor2: the second array
:return: True if two arrays are equal, False otherwise.
If one or more of the inputs is not an ndarray, return False.
"""
are_np_arrays = isinstance(tensor1, np.ndarray) and isinstance(
tensor2, np.ndarray
)
return are_np_arrays and np.array_equal(tensor1, tensor2)
class Retrieval(AbstractComputationalBackend.Retrieval[np.ndarray]):
"""
Abstract class for retrieval and ranking functionalities
"""
@staticmethod
def top_k(
values: 'np.ndarray',
k: int,
descending: bool = False,
device: Optional[str] = None,
) -> Tuple['np.ndarray', 'np.ndarray']:
"""
Retrieves the top k smallest values in `values`,
and returns them alongside their indices in the input `values`.
Can also be used to retrieve the top k largest values,
by setting the `descending` flag.
:param values: Torch tensor of values to rank.
Should be of shape (n_queries, n_values_per_query).
Inputs of shape (n_values_per_query,) will be expanded
to (1, n_values_per_query).
:param k: number of values to retrieve
:param descending: retrieve largest values instead of smallest values
:param device: Not supported for this backend
:return: Tuple containing the retrieved values, and their indices.
Both ar of shape (n_queries, k)
"""
if device is not None:
warnings.warn('`device` is not supported for numpy operations')
if len(values.shape) == 1:
values = np.expand_dims(values, axis=0)
if descending:
values = -values
if k >= values.shape[1]:
idx = values.argsort(axis=1)[:, :k]
values = np.take_along_axis(values, idx, axis=1)
else:
idx_ps = values.argpartition(kth=k, axis=1)[:, :k]
values = np.take_along_axis(values, idx_ps, axis=1)
idx_fs = values.argsort(axis=1)
idx = np.take_along_axis(idx_ps, idx_fs, axis=1)
values = np.take_along_axis(values, idx_fs, axis=1)
if descending:
values = -values
return values, idx
class Metrics(AbstractComputationalBackend.Metrics[np.ndarray]):
"""
Abstract base class for metrics (distances and similarities).
"""
@staticmethod
def cosine_sim(
x_mat: np.ndarray,
y_mat: np.ndarray,
eps: float = 1e-7,
device: Optional[str] = None,
) -> np.ndarray:
"""Pairwise cosine similarities between all vectors in x_mat and y_mat.
:param x_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param y_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param eps: a small jitter to avoid divde by zero
:param device: Not supported for this backend
:return: np.ndarray of shape (n_vectors, n_vectors) containing all
pairwise cosine distances.
The index [i_x, i_y] contains the cosine distance between
x_mat[i_x] and y_mat[i_y].
"""
if device is not None:
warnings.warn('`device` is not supported for numpy operations')
x_mat, y_mat = _expand_if_single_axis(x_mat, y_mat)
sims = np.clip(
(np.dot(x_mat, y_mat.T) + eps)
/ (
np.outer(
np.linalg.norm(x_mat, axis=1), np.linalg.norm(y_mat, axis=1)
)
+ eps
),
-1,
1,
).squeeze()
return _expand_if_scalar(sims)
@classmethod
def euclidean_dist(
cls, x_mat: np.ndarray, y_mat: np.ndarray, device: Optional[str] = None
) -> np.ndarray:
"""Pairwise Euclidian distances between all vectors in x_mat and y_mat.
:param x_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param y_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param eps: a small jitter to avoid divde by zero
:param device: Not supported for this backend
:return: np.ndarray of shape (n_vectors, n_vectors) containing all
pairwise euclidian distances.
The index [i_x, i_y] contains the euclidian distance between
x_mat[i_x] and y_mat[i_y].
"""
if device is not None:
warnings.warn('`device` is not supported for numpy operations')
x_mat, y_mat = _expand_if_single_axis(x_mat, y_mat)
return _expand_if_scalar(
np.sqrt(cls.sqeuclidean_dist(x_mat, y_mat)).squeeze()
)
@staticmethod
def sqeuclidean_dist(
x_mat: np.ndarray,
y_mat: np.ndarray,
device: Optional[str] = None,
) -> np.ndarray:
"""Pairwise Squared Euclidian distances between all vectors in
x_mat and y_mat.
:param x_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param y_mat: np.ndarray of shape (n_vectors, n_dim), where n_vectors is
the number of vectors and n_dim is the number of dimensions of each
example.
:param device: Not supported for this backend
:return: np.ndarray of shape (n_vectors, n_vectors) containing all
pairwise Squared Euclidian distances.
The index [i_x, i_y] contains the cosine Squared Euclidian between
x_mat[i_x] and y_mat[i_y].
"""
eps: float = 1e-7 # avoid problems with numerical inaccuracies
if device is not None:
warnings.warn('`device` is not supported for numpy operations')
x_mat, y_mat = _expand_if_single_axis(x_mat, y_mat)
dists = (
np.sum(y_mat**2, axis=1)
+ np.sum(x_mat**2, axis=1)[:, np.newaxis]
- 2 * np.dot(x_mat, y_mat.T)
).squeeze()
# remove numerical artifacts
dists = np.where(np.logical_and(dists < 0, dists > -eps), 0, dists)
return _expand_if_scalar(dists)