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
|
// Copyright Yahoo. Licensed under the terms of the Apache 2.0 license. See LICENSE in the project root.
#pragma once
#include <cstring>
#include <vespa/vespalib/util/arrayref.h>
#include <vespa/vespalib/hwaccelrated/iaccelrated.h>
template <typename T, size_t VLEN>
static double hw_l2_sq_dist(const T * af, const T * bf, size_t sz)
{
constexpr const size_t OpsPerV = VLEN/sizeof(T);
typedef T V __attribute__ ((vector_size (VLEN), aligned(VLEN)));
const V * a = reinterpret_cast<const V *>(af);
const V * b = reinterpret_cast<const V *>(bf);
V tmp_diff;
V tmp_squa;
V tmp_sum;
memset(&tmp_sum, 0, sizeof(tmp_sum));
const size_t numOps = sz/OpsPerV;
for (size_t i = 0; i < numOps; ++i) {
tmp_diff = a[i] - b[i];
tmp_squa = tmp_diff * tmp_diff;
tmp_sum += tmp_squa;
}
double sum = 0;
for (size_t i = 0; i < OpsPerV; ++i) {
sum += tmp_sum[i];
}
return sum;
}
template <typename FltType = float>
struct L2DistCalc {
const vespalib::hwaccelrated::IAccelrated & _hw;
L2DistCalc() : _hw(vespalib::hwaccelrated::IAccelrated::getAccelerator()) {}
using Arr = vespalib::ArrayRef<FltType>;
using ConstArr = vespalib::ConstArrayRef<FltType>;
double product(const FltType *v1, const FltType *v2, size_t sz) {
return _hw.dotProduct(v1, v2, sz);
}
double product(ConstArr v1, ConstArr v2) {
const FltType *p1 = v1.begin();
const FltType *p2 = v2.begin();
return _hw.dotProduct(p1, p2, v1.size());
}
double l2sq(ConstArr vector) {
const FltType *v = vector.begin();
return _hw.dotProduct(v, v, vector.size());
}
double l2sq_dist(ConstArr v1, ConstArr v2, Arr tmp) {
for (size_t i = 0; i < v1.size(); ++i) {
tmp[i] = (v1[i] - v2[i]);
}
return l2sq(tmp);
}
double l2sq_dist(ConstArr v1, ConstArr v2) {
return hw_l2_sq_dist<FltType, 32>(v1.cbegin(), v2.cbegin(), v1.size());
}
};
static L2DistCalc l2distCalc;
|
