faiss 0.1.7 → 0.2.3
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/CHANGELOG.md +18 -0
- data/README.md +7 -7
- data/ext/faiss/ext.cpp +1 -1
- data/ext/faiss/extconf.rb +8 -2
- data/ext/faiss/index.cpp +102 -69
- data/ext/faiss/index_binary.cpp +24 -30
- data/ext/faiss/kmeans.cpp +20 -16
- data/ext/faiss/numo.hpp +867 -0
- data/ext/faiss/pca_matrix.cpp +13 -14
- data/ext/faiss/product_quantizer.cpp +23 -24
- data/ext/faiss/utils.cpp +10 -37
- data/ext/faiss/utils.h +2 -13
- data/lib/faiss/version.rb +1 -1
- data/lib/faiss.rb +0 -5
- data/vendor/faiss/faiss/AutoTune.cpp +292 -291
- data/vendor/faiss/faiss/AutoTune.h +55 -56
- data/vendor/faiss/faiss/Clustering.cpp +334 -195
- data/vendor/faiss/faiss/Clustering.h +88 -35
- data/vendor/faiss/faiss/IVFlib.cpp +171 -195
- data/vendor/faiss/faiss/IVFlib.h +48 -51
- data/vendor/faiss/faiss/Index.cpp +85 -103
- data/vendor/faiss/faiss/Index.h +54 -48
- data/vendor/faiss/faiss/Index2Layer.cpp +139 -164
- data/vendor/faiss/faiss/Index2Layer.h +22 -22
- data/vendor/faiss/faiss/IndexBinary.cpp +45 -37
- data/vendor/faiss/faiss/IndexBinary.h +140 -132
- data/vendor/faiss/faiss/IndexBinaryFlat.cpp +73 -53
- data/vendor/faiss/faiss/IndexBinaryFlat.h +29 -24
- data/vendor/faiss/faiss/IndexBinaryFromFloat.cpp +46 -43
- data/vendor/faiss/faiss/IndexBinaryFromFloat.h +16 -15
- data/vendor/faiss/faiss/IndexBinaryHNSW.cpp +215 -232
- data/vendor/faiss/faiss/IndexBinaryHNSW.h +25 -24
- data/vendor/faiss/faiss/IndexBinaryHash.cpp +182 -177
- data/vendor/faiss/faiss/IndexBinaryHash.h +41 -34
- data/vendor/faiss/faiss/IndexBinaryIVF.cpp +489 -461
- data/vendor/faiss/faiss/IndexBinaryIVF.h +97 -68
- data/vendor/faiss/faiss/IndexFlat.cpp +116 -147
- data/vendor/faiss/faiss/IndexFlat.h +35 -46
- data/vendor/faiss/faiss/IndexHNSW.cpp +372 -348
- data/vendor/faiss/faiss/IndexHNSW.h +57 -41
- data/vendor/faiss/faiss/IndexIVF.cpp +474 -454
- data/vendor/faiss/faiss/IndexIVF.h +146 -113
- data/vendor/faiss/faiss/IndexIVFFlat.cpp +248 -250
- data/vendor/faiss/faiss/IndexIVFFlat.h +48 -51
- data/vendor/faiss/faiss/IndexIVFPQ.cpp +457 -516
- data/vendor/faiss/faiss/IndexIVFPQ.h +74 -66
- data/vendor/faiss/faiss/IndexIVFPQFastScan.cpp +406 -372
- data/vendor/faiss/faiss/IndexIVFPQFastScan.h +82 -57
- data/vendor/faiss/faiss/IndexIVFPQR.cpp +104 -102
- data/vendor/faiss/faiss/IndexIVFPQR.h +33 -28
- data/vendor/faiss/faiss/IndexIVFSpectralHash.cpp +125 -133
- data/vendor/faiss/faiss/IndexIVFSpectralHash.h +19 -21
- data/vendor/faiss/faiss/IndexLSH.cpp +75 -96
- data/vendor/faiss/faiss/IndexLSH.h +21 -26
- data/vendor/faiss/faiss/IndexLattice.cpp +42 -56
- data/vendor/faiss/faiss/IndexLattice.h +11 -16
- data/vendor/faiss/faiss/IndexNNDescent.cpp +231 -0
- data/vendor/faiss/faiss/IndexNNDescent.h +72 -0
- data/vendor/faiss/faiss/IndexNSG.cpp +303 -0
- data/vendor/faiss/faiss/IndexNSG.h +85 -0
- data/vendor/faiss/faiss/IndexPQ.cpp +405 -464
- data/vendor/faiss/faiss/IndexPQ.h +64 -67
- data/vendor/faiss/faiss/IndexPQFastScan.cpp +143 -170
- data/vendor/faiss/faiss/IndexPQFastScan.h +46 -32
- data/vendor/faiss/faiss/IndexPreTransform.cpp +120 -150
- data/vendor/faiss/faiss/IndexPreTransform.h +33 -36
- data/vendor/faiss/faiss/IndexRefine.cpp +115 -131
- data/vendor/faiss/faiss/IndexRefine.h +22 -23
- data/vendor/faiss/faiss/IndexReplicas.cpp +147 -153
- data/vendor/faiss/faiss/IndexReplicas.h +62 -56
- data/vendor/faiss/faiss/IndexResidual.cpp +291 -0
- data/vendor/faiss/faiss/IndexResidual.h +152 -0
- data/vendor/faiss/faiss/IndexScalarQuantizer.cpp +120 -155
- data/vendor/faiss/faiss/IndexScalarQuantizer.h +41 -45
- data/vendor/faiss/faiss/IndexShards.cpp +256 -240
- data/vendor/faiss/faiss/IndexShards.h +85 -73
- data/vendor/faiss/faiss/MatrixStats.cpp +112 -97
- data/vendor/faiss/faiss/MatrixStats.h +7 -10
- data/vendor/faiss/faiss/MetaIndexes.cpp +135 -157
- data/vendor/faiss/faiss/MetaIndexes.h +40 -34
- data/vendor/faiss/faiss/MetricType.h +7 -7
- data/vendor/faiss/faiss/VectorTransform.cpp +652 -474
- data/vendor/faiss/faiss/VectorTransform.h +61 -89
- data/vendor/faiss/faiss/clone_index.cpp +77 -73
- data/vendor/faiss/faiss/clone_index.h +4 -9
- data/vendor/faiss/faiss/gpu/GpuAutoTune.cpp +33 -38
- data/vendor/faiss/faiss/gpu/GpuAutoTune.h +11 -9
- data/vendor/faiss/faiss/gpu/GpuCloner.cpp +197 -170
- data/vendor/faiss/faiss/gpu/GpuCloner.h +53 -35
- data/vendor/faiss/faiss/gpu/GpuClonerOptions.cpp +12 -14
- data/vendor/faiss/faiss/gpu/GpuClonerOptions.h +27 -25
- data/vendor/faiss/faiss/gpu/GpuDistance.h +116 -112
- data/vendor/faiss/faiss/gpu/GpuFaissAssert.h +1 -2
- data/vendor/faiss/faiss/gpu/GpuIndex.h +134 -137
- data/vendor/faiss/faiss/gpu/GpuIndexBinaryFlat.h +76 -73
- data/vendor/faiss/faiss/gpu/GpuIndexFlat.h +173 -162
- data/vendor/faiss/faiss/gpu/GpuIndexIVF.h +67 -64
- data/vendor/faiss/faiss/gpu/GpuIndexIVFFlat.h +89 -86
- data/vendor/faiss/faiss/gpu/GpuIndexIVFPQ.h +150 -141
- data/vendor/faiss/faiss/gpu/GpuIndexIVFScalarQuantizer.h +101 -103
- data/vendor/faiss/faiss/gpu/GpuIndicesOptions.h +17 -16
- data/vendor/faiss/faiss/gpu/GpuResources.cpp +116 -128
- data/vendor/faiss/faiss/gpu/GpuResources.h +182 -186
- data/vendor/faiss/faiss/gpu/StandardGpuResources.cpp +433 -422
- data/vendor/faiss/faiss/gpu/StandardGpuResources.h +131 -130
- data/vendor/faiss/faiss/gpu/impl/InterleavedCodes.cpp +468 -456
- data/vendor/faiss/faiss/gpu/impl/InterleavedCodes.h +25 -19
- data/vendor/faiss/faiss/gpu/impl/RemapIndices.cpp +22 -20
- data/vendor/faiss/faiss/gpu/impl/RemapIndices.h +9 -8
- data/vendor/faiss/faiss/gpu/perf/IndexWrapper-inl.h +39 -44
- data/vendor/faiss/faiss/gpu/perf/IndexWrapper.h +16 -14
- data/vendor/faiss/faiss/gpu/perf/PerfClustering.cpp +77 -71
- data/vendor/faiss/faiss/gpu/perf/PerfIVFPQAdd.cpp +109 -88
- data/vendor/faiss/faiss/gpu/perf/WriteIndex.cpp +75 -64
- data/vendor/faiss/faiss/gpu/test/TestCodePacking.cpp +230 -215
- data/vendor/faiss/faiss/gpu/test/TestGpuIndexBinaryFlat.cpp +80 -86
- data/vendor/faiss/faiss/gpu/test/TestGpuIndexFlat.cpp +284 -277
- data/vendor/faiss/faiss/gpu/test/TestGpuIndexIVFFlat.cpp +416 -416
- data/vendor/faiss/faiss/gpu/test/TestGpuIndexIVFPQ.cpp +611 -517
- data/vendor/faiss/faiss/gpu/test/TestGpuIndexIVFScalarQuantizer.cpp +166 -164
- data/vendor/faiss/faiss/gpu/test/TestGpuMemoryException.cpp +61 -53
- data/vendor/faiss/faiss/gpu/test/TestUtils.cpp +274 -238
- data/vendor/faiss/faiss/gpu/test/TestUtils.h +73 -57
- data/vendor/faiss/faiss/gpu/test/demo_ivfpq_indexing_gpu.cpp +47 -50
- data/vendor/faiss/faiss/gpu/utils/DeviceUtils.h +79 -72
- data/vendor/faiss/faiss/gpu/utils/StackDeviceMemory.cpp +140 -146
- data/vendor/faiss/faiss/gpu/utils/StackDeviceMemory.h +69 -71
- data/vendor/faiss/faiss/gpu/utils/StaticUtils.h +21 -16
- data/vendor/faiss/faiss/gpu/utils/Timer.cpp +25 -29
- data/vendor/faiss/faiss/gpu/utils/Timer.h +30 -29
- data/vendor/faiss/faiss/impl/AdditiveQuantizer.cpp +270 -0
- data/vendor/faiss/faiss/impl/AdditiveQuantizer.h +115 -0
- data/vendor/faiss/faiss/impl/AuxIndexStructures.cpp +90 -120
- data/vendor/faiss/faiss/impl/AuxIndexStructures.h +81 -65
- data/vendor/faiss/faiss/impl/FaissAssert.h +73 -58
- data/vendor/faiss/faiss/impl/FaissException.cpp +56 -48
- data/vendor/faiss/faiss/impl/FaissException.h +41 -29
- data/vendor/faiss/faiss/impl/HNSW.cpp +595 -611
- data/vendor/faiss/faiss/impl/HNSW.h +179 -200
- data/vendor/faiss/faiss/impl/LocalSearchQuantizer.cpp +672 -0
- data/vendor/faiss/faiss/impl/LocalSearchQuantizer.h +172 -0
- data/vendor/faiss/faiss/impl/NNDescent.cpp +487 -0
- data/vendor/faiss/faiss/impl/NNDescent.h +154 -0
- data/vendor/faiss/faiss/impl/NSG.cpp +682 -0
- data/vendor/faiss/faiss/impl/NSG.h +199 -0
- data/vendor/faiss/faiss/impl/PolysemousTraining.cpp +484 -454
- data/vendor/faiss/faiss/impl/PolysemousTraining.h +52 -55
- data/vendor/faiss/faiss/impl/ProductQuantizer-inl.h +26 -47
- data/vendor/faiss/faiss/impl/ProductQuantizer.cpp +469 -459
- data/vendor/faiss/faiss/impl/ProductQuantizer.h +76 -87
- data/vendor/faiss/faiss/impl/ResidualQuantizer.cpp +448 -0
- data/vendor/faiss/faiss/impl/ResidualQuantizer.h +130 -0
- data/vendor/faiss/faiss/impl/ResultHandler.h +96 -132
- data/vendor/faiss/faiss/impl/ScalarQuantizer.cpp +648 -701
- data/vendor/faiss/faiss/impl/ScalarQuantizer.h +48 -46
- data/vendor/faiss/faiss/impl/ThreadedIndex-inl.h +129 -131
- data/vendor/faiss/faiss/impl/ThreadedIndex.h +61 -55
- data/vendor/faiss/faiss/impl/index_read.cpp +547 -479
- data/vendor/faiss/faiss/impl/index_write.cpp +497 -407
- data/vendor/faiss/faiss/impl/io.cpp +75 -94
- data/vendor/faiss/faiss/impl/io.h +31 -41
- data/vendor/faiss/faiss/impl/io_macros.h +40 -29
- data/vendor/faiss/faiss/impl/lattice_Zn.cpp +137 -186
- data/vendor/faiss/faiss/impl/lattice_Zn.h +40 -51
- data/vendor/faiss/faiss/impl/platform_macros.h +29 -8
- data/vendor/faiss/faiss/impl/pq4_fast_scan.cpp +77 -124
- data/vendor/faiss/faiss/impl/pq4_fast_scan.h +39 -48
- data/vendor/faiss/faiss/impl/pq4_fast_scan_search_1.cpp +41 -52
- data/vendor/faiss/faiss/impl/pq4_fast_scan_search_qbs.cpp +80 -117
- data/vendor/faiss/faiss/impl/simd_result_handlers.h +109 -137
- data/vendor/faiss/faiss/index_factory.cpp +269 -218
- data/vendor/faiss/faiss/index_factory.h +6 -7
- data/vendor/faiss/faiss/index_io.h +23 -26
- data/vendor/faiss/faiss/invlists/BlockInvertedLists.cpp +67 -75
- data/vendor/faiss/faiss/invlists/BlockInvertedLists.h +22 -24
- data/vendor/faiss/faiss/invlists/DirectMap.cpp +96 -112
- data/vendor/faiss/faiss/invlists/DirectMap.h +29 -33
- data/vendor/faiss/faiss/invlists/InvertedLists.cpp +307 -364
- data/vendor/faiss/faiss/invlists/InvertedLists.h +151 -151
- data/vendor/faiss/faiss/invlists/InvertedListsIOHook.cpp +29 -34
- data/vendor/faiss/faiss/invlists/InvertedListsIOHook.h +17 -18
- data/vendor/faiss/faiss/invlists/OnDiskInvertedLists.cpp +257 -293
- data/vendor/faiss/faiss/invlists/OnDiskInvertedLists.h +50 -45
- data/vendor/faiss/faiss/python/python_callbacks.cpp +23 -26
- data/vendor/faiss/faiss/python/python_callbacks.h +9 -16
- data/vendor/faiss/faiss/utils/AlignedTable.h +79 -44
- data/vendor/faiss/faiss/utils/Heap.cpp +40 -48
- data/vendor/faiss/faiss/utils/Heap.h +186 -209
- data/vendor/faiss/faiss/utils/WorkerThread.cpp +67 -76
- data/vendor/faiss/faiss/utils/WorkerThread.h +32 -33
- data/vendor/faiss/faiss/utils/distances.cpp +301 -310
- data/vendor/faiss/faiss/utils/distances.h +133 -118
- data/vendor/faiss/faiss/utils/distances_simd.cpp +456 -516
- data/vendor/faiss/faiss/utils/extra_distances-inl.h +117 -0
- data/vendor/faiss/faiss/utils/extra_distances.cpp +113 -232
- data/vendor/faiss/faiss/utils/extra_distances.h +30 -29
- data/vendor/faiss/faiss/utils/hamming-inl.h +260 -209
- data/vendor/faiss/faiss/utils/hamming.cpp +375 -469
- data/vendor/faiss/faiss/utils/hamming.h +62 -85
- data/vendor/faiss/faiss/utils/ordered_key_value.h +16 -18
- data/vendor/faiss/faiss/utils/partitioning.cpp +393 -318
- data/vendor/faiss/faiss/utils/partitioning.h +26 -21
- data/vendor/faiss/faiss/utils/quantize_lut.cpp +78 -66
- data/vendor/faiss/faiss/utils/quantize_lut.h +22 -20
- data/vendor/faiss/faiss/utils/random.cpp +39 -63
- data/vendor/faiss/faiss/utils/random.h +13 -16
- data/vendor/faiss/faiss/utils/simdlib.h +4 -2
- data/vendor/faiss/faiss/utils/simdlib_avx2.h +88 -85
- data/vendor/faiss/faiss/utils/simdlib_emulated.h +226 -165
- data/vendor/faiss/faiss/utils/simdlib_neon.h +832 -0
- data/vendor/faiss/faiss/utils/utils.cpp +304 -287
- data/vendor/faiss/faiss/utils/utils.h +53 -48
- metadata +26 -12
- data/lib/faiss/index.rb +0 -20
- data/lib/faiss/index_binary.rb +0 -20
- data/lib/faiss/kmeans.rb +0 -15
- data/lib/faiss/pca_matrix.rb +0 -15
- data/lib/faiss/product_quantizer.rb +0 -22
@@ -9,108 +9,96 @@
|
|
9
9
|
|
10
10
|
#include <faiss/IndexFlat.h>
|
11
11
|
|
12
|
-
#include <
|
12
|
+
#include <faiss/impl/AuxIndexStructures.h>
|
13
|
+
#include <faiss/impl/FaissAssert.h>
|
14
|
+
#include <faiss/utils/Heap.h>
|
13
15
|
#include <faiss/utils/distances.h>
|
14
16
|
#include <faiss/utils/extra_distances.h>
|
15
17
|
#include <faiss/utils/utils.h>
|
16
|
-
#include <
|
17
|
-
#include <faiss/impl/FaissAssert.h>
|
18
|
-
#include <faiss/impl/AuxIndexStructures.h>
|
19
|
-
|
18
|
+
#include <cstring>
|
20
19
|
|
21
20
|
namespace faiss {
|
22
21
|
|
23
|
-
IndexFlat::IndexFlat
|
24
|
-
Index(d, metric)
|
25
|
-
{
|
26
|
-
}
|
27
|
-
|
22
|
+
IndexFlat::IndexFlat(idx_t d, MetricType metric) : Index(d, metric) {}
|
28
23
|
|
29
|
-
|
30
|
-
void IndexFlat::add (idx_t n, const float *x) {
|
24
|
+
void IndexFlat::add(idx_t n, const float* x) {
|
31
25
|
xb.insert(xb.end(), x, x + n * d);
|
32
26
|
ntotal += n;
|
33
27
|
}
|
34
28
|
|
35
|
-
|
36
29
|
void IndexFlat::reset() {
|
37
30
|
xb.clear();
|
38
31
|
ntotal = 0;
|
39
32
|
}
|
40
33
|
|
34
|
+
void IndexFlat::search(
|
35
|
+
idx_t n,
|
36
|
+
const float* x,
|
37
|
+
idx_t k,
|
38
|
+
float* distances,
|
39
|
+
idx_t* labels) const {
|
40
|
+
FAISS_THROW_IF_NOT(k > 0);
|
41
41
|
|
42
|
-
void IndexFlat::search (idx_t n, const float *x, idx_t k,
|
43
|
-
float *distances, idx_t *labels) const
|
44
|
-
{
|
45
42
|
// we see the distances and labels as heaps
|
46
43
|
|
47
44
|
if (metric_type == METRIC_INNER_PRODUCT) {
|
48
|
-
float_minheap_array_t res = {
|
49
|
-
|
50
|
-
knn_inner_product (x, xb.data(), d, n, ntotal, &res);
|
45
|
+
float_minheap_array_t res = {size_t(n), size_t(k), labels, distances};
|
46
|
+
knn_inner_product(x, xb.data(), d, n, ntotal, &res);
|
51
47
|
} else if (metric_type == METRIC_L2) {
|
52
|
-
float_maxheap_array_t res = {
|
53
|
-
|
54
|
-
knn_L2sqr (x, xb.data(), d, n, ntotal, &res);
|
48
|
+
float_maxheap_array_t res = {size_t(n), size_t(k), labels, distances};
|
49
|
+
knn_L2sqr(x, xb.data(), d, n, ntotal, &res);
|
55
50
|
} else {
|
56
|
-
float_maxheap_array_t res = {
|
57
|
-
|
58
|
-
|
59
|
-
metric_type, metric_arg,
|
60
|
-
&res);
|
51
|
+
float_maxheap_array_t res = {size_t(n), size_t(k), labels, distances};
|
52
|
+
knn_extra_metrics(
|
53
|
+
x, xb.data(), d, n, ntotal, metric_type, metric_arg, &res);
|
61
54
|
}
|
62
55
|
}
|
63
56
|
|
64
|
-
void IndexFlat::range_search
|
65
|
-
|
66
|
-
|
57
|
+
void IndexFlat::range_search(
|
58
|
+
idx_t n,
|
59
|
+
const float* x,
|
60
|
+
float radius,
|
61
|
+
RangeSearchResult* result) const {
|
67
62
|
switch (metric_type) {
|
68
|
-
|
69
|
-
|
70
|
-
|
71
|
-
|
72
|
-
|
73
|
-
|
74
|
-
|
75
|
-
|
76
|
-
|
63
|
+
case METRIC_INNER_PRODUCT:
|
64
|
+
range_search_inner_product(
|
65
|
+
x, xb.data(), d, n, ntotal, radius, result);
|
66
|
+
break;
|
67
|
+
case METRIC_L2:
|
68
|
+
range_search_L2sqr(x, xb.data(), d, n, ntotal, radius, result);
|
69
|
+
break;
|
70
|
+
default:
|
71
|
+
FAISS_THROW_MSG("metric type not supported");
|
77
72
|
}
|
78
73
|
}
|
79
74
|
|
80
|
-
|
81
|
-
|
82
|
-
|
83
|
-
|
84
|
-
|
85
|
-
|
86
|
-
const idx_t *labels) const
|
87
|
-
{
|
75
|
+
void IndexFlat::compute_distance_subset(
|
76
|
+
idx_t n,
|
77
|
+
const float* x,
|
78
|
+
idx_t k,
|
79
|
+
float* distances,
|
80
|
+
const idx_t* labels) const {
|
88
81
|
switch (metric_type) {
|
89
82
|
case METRIC_INNER_PRODUCT:
|
90
|
-
fvec_inner_products_by_idx
|
91
|
-
|
92
|
-
x, xb.data(), labels, d, n, k);
|
83
|
+
fvec_inner_products_by_idx(
|
84
|
+
distances, x, xb.data(), labels, d, n, k);
|
93
85
|
break;
|
94
86
|
case METRIC_L2:
|
95
|
-
fvec_L2sqr_by_idx (
|
96
|
-
distances,
|
97
|
-
x, xb.data(), labels, d, n, k);
|
87
|
+
fvec_L2sqr_by_idx(distances, x, xb.data(), labels, d, n, k);
|
98
88
|
break;
|
99
89
|
default:
|
100
90
|
FAISS_THROW_MSG("metric type not supported");
|
101
91
|
}
|
102
|
-
|
103
92
|
}
|
104
93
|
|
105
|
-
size_t IndexFlat::remove_ids
|
106
|
-
{
|
94
|
+
size_t IndexFlat::remove_ids(const IDSelector& sel) {
|
107
95
|
idx_t j = 0;
|
108
96
|
for (idx_t i = 0; i < ntotal; i++) {
|
109
|
-
if (sel.is_member
|
97
|
+
if (sel.is_member(i)) {
|
110
98
|
// should be removed
|
111
99
|
} else {
|
112
100
|
if (i > j) {
|
113
|
-
memmove
|
101
|
+
memmove(&xb[d * j], &xb[d * i], sizeof(xb[0]) * d);
|
114
102
|
}
|
115
103
|
j++;
|
116
104
|
}
|
@@ -118,23 +106,21 @@ size_t IndexFlat::remove_ids (const IDSelector & sel)
|
|
118
106
|
size_t nremove = ntotal - j;
|
119
107
|
if (nremove > 0) {
|
120
108
|
ntotal = j;
|
121
|
-
xb.resize
|
109
|
+
xb.resize(ntotal * d);
|
122
110
|
}
|
123
111
|
return nremove;
|
124
112
|
}
|
125
113
|
|
126
|
-
|
127
114
|
namespace {
|
128
115
|
|
129
|
-
|
130
116
|
struct FlatL2Dis : DistanceComputer {
|
131
117
|
size_t d;
|
132
118
|
Index::idx_t nb;
|
133
|
-
const float
|
134
|
-
const float
|
119
|
+
const float* q;
|
120
|
+
const float* b;
|
135
121
|
size_t ndis;
|
136
122
|
|
137
|
-
float operator
|
123
|
+
float operator()(idx_t i) override {
|
138
124
|
ndis++;
|
139
125
|
return fvec_L2sqr(q, b + i * d, d);
|
140
126
|
}
|
@@ -143,14 +129,14 @@ struct FlatL2Dis : DistanceComputer {
|
|
143
129
|
return fvec_L2sqr(b + j * d, b + i * d, d);
|
144
130
|
}
|
145
131
|
|
146
|
-
explicit FlatL2Dis(const IndexFlat& storage, const float
|
147
|
-
|
148
|
-
|
149
|
-
|
150
|
-
|
151
|
-
|
132
|
+
explicit FlatL2Dis(const IndexFlat& storage, const float* q = nullptr)
|
133
|
+
: d(storage.d),
|
134
|
+
nb(storage.ntotal),
|
135
|
+
q(q),
|
136
|
+
b(storage.xb.data()),
|
137
|
+
ndis(0) {}
|
152
138
|
|
153
|
-
void set_query(const float
|
139
|
+
void set_query(const float* x) override {
|
154
140
|
q = x;
|
155
141
|
}
|
156
142
|
};
|
@@ -158,128 +144,106 @@ struct FlatL2Dis : DistanceComputer {
|
|
158
144
|
struct FlatIPDis : DistanceComputer {
|
159
145
|
size_t d;
|
160
146
|
Index::idx_t nb;
|
161
|
-
const float
|
162
|
-
const float
|
147
|
+
const float* q;
|
148
|
+
const float* b;
|
163
149
|
size_t ndis;
|
164
150
|
|
165
|
-
float operator
|
151
|
+
float operator()(idx_t i) override {
|
166
152
|
ndis++;
|
167
|
-
return fvec_inner_product
|
153
|
+
return fvec_inner_product(q, b + i * d, d);
|
168
154
|
}
|
169
155
|
|
170
156
|
float symmetric_dis(idx_t i, idx_t j) override {
|
171
|
-
return fvec_inner_product
|
157
|
+
return fvec_inner_product(b + j * d, b + i * d, d);
|
172
158
|
}
|
173
159
|
|
174
|
-
explicit FlatIPDis(const IndexFlat& storage, const float
|
175
|
-
|
176
|
-
|
177
|
-
|
178
|
-
|
179
|
-
|
160
|
+
explicit FlatIPDis(const IndexFlat& storage, const float* q = nullptr)
|
161
|
+
: d(storage.d),
|
162
|
+
nb(storage.ntotal),
|
163
|
+
q(q),
|
164
|
+
b(storage.xb.data()),
|
165
|
+
ndis(0) {}
|
180
166
|
|
181
|
-
void set_query(const float
|
167
|
+
void set_query(const float* x) override {
|
182
168
|
q = x;
|
183
169
|
}
|
184
170
|
};
|
185
171
|
|
172
|
+
} // namespace
|
186
173
|
|
187
|
-
|
188
|
-
|
189
|
-
} // namespace
|
190
|
-
|
191
|
-
|
192
|
-
DistanceComputer * IndexFlat::get_distance_computer() const {
|
174
|
+
DistanceComputer* IndexFlat::get_distance_computer() const {
|
193
175
|
if (metric_type == METRIC_L2) {
|
194
176
|
return new FlatL2Dis(*this);
|
195
177
|
} else if (metric_type == METRIC_INNER_PRODUCT) {
|
196
178
|
return new FlatIPDis(*this);
|
197
179
|
} else {
|
198
|
-
return get_extra_distance_computer
|
199
|
-
|
180
|
+
return get_extra_distance_computer(
|
181
|
+
d, metric_type, metric_arg, ntotal, xb.data());
|
200
182
|
}
|
201
183
|
}
|
202
184
|
|
203
|
-
|
204
|
-
|
205
|
-
{
|
206
|
-
memcpy (recons, &(xb[key * d]), sizeof(*recons) * d);
|
185
|
+
void IndexFlat::reconstruct(idx_t key, float* recons) const {
|
186
|
+
memcpy(recons, &(xb[key * d]), sizeof(*recons) * d);
|
207
187
|
}
|
208
188
|
|
209
|
-
|
210
189
|
/* The standalone codec interface */
|
211
|
-
size_t IndexFlat::sa_code_size
|
212
|
-
{
|
190
|
+
size_t IndexFlat::sa_code_size() const {
|
213
191
|
return sizeof(float) * d;
|
214
192
|
}
|
215
193
|
|
216
|
-
void IndexFlat::sa_encode
|
217
|
-
|
218
|
-
memcpy (bytes, x, sizeof(float) * d * n);
|
194
|
+
void IndexFlat::sa_encode(idx_t n, const float* x, uint8_t* bytes) const {
|
195
|
+
memcpy(bytes, x, sizeof(float) * d * n);
|
219
196
|
}
|
220
197
|
|
221
|
-
void IndexFlat::sa_decode
|
222
|
-
|
223
|
-
memcpy (x, bytes, sizeof(float) * d * n);
|
198
|
+
void IndexFlat::sa_decode(idx_t n, const uint8_t* bytes, float* x) const {
|
199
|
+
memcpy(x, bytes, sizeof(float) * d * n);
|
224
200
|
}
|
225
201
|
|
226
|
-
|
227
|
-
|
228
|
-
|
229
|
-
|
230
|
-
|
231
202
|
/***************************************************
|
232
203
|
* IndexFlat1D
|
233
204
|
***************************************************/
|
234
205
|
|
235
|
-
|
236
|
-
|
237
|
-
IndexFlatL2 (1),
|
238
|
-
continuous_update (continuous_update)
|
239
|
-
{
|
240
|
-
}
|
206
|
+
IndexFlat1D::IndexFlat1D(bool continuous_update)
|
207
|
+
: IndexFlatL2(1), continuous_update(continuous_update) {}
|
241
208
|
|
242
209
|
/// if not continuous_update, call this between the last add and
|
243
210
|
/// the first search
|
244
|
-
void IndexFlat1D::update_permutation
|
245
|
-
|
246
|
-
perm.resize (ntotal);
|
211
|
+
void IndexFlat1D::update_permutation() {
|
212
|
+
perm.resize(ntotal);
|
247
213
|
if (ntotal < 1000000) {
|
248
|
-
fvec_argsort
|
214
|
+
fvec_argsort(ntotal, xb.data(), (size_t*)perm.data());
|
249
215
|
} else {
|
250
|
-
fvec_argsort_parallel
|
216
|
+
fvec_argsort_parallel(ntotal, xb.data(), (size_t*)perm.data());
|
251
217
|
}
|
252
218
|
}
|
253
219
|
|
254
|
-
void IndexFlat1D::add
|
255
|
-
|
256
|
-
IndexFlatL2::add (n, x);
|
220
|
+
void IndexFlat1D::add(idx_t n, const float* x) {
|
221
|
+
IndexFlatL2::add(n, x);
|
257
222
|
if (continuous_update)
|
258
223
|
update_permutation();
|
259
224
|
}
|
260
225
|
|
261
|
-
void IndexFlat1D::reset()
|
262
|
-
{
|
226
|
+
void IndexFlat1D::reset() {
|
263
227
|
IndexFlatL2::reset();
|
264
228
|
perm.clear();
|
265
229
|
}
|
266
230
|
|
267
|
-
void IndexFlat1D::search
|
268
|
-
|
269
|
-
|
270
|
-
|
271
|
-
|
272
|
-
|
273
|
-
|
274
|
-
|
275
|
-
|
231
|
+
void IndexFlat1D::search(
|
232
|
+
idx_t n,
|
233
|
+
const float* x,
|
234
|
+
idx_t k,
|
235
|
+
float* distances,
|
236
|
+
idx_t* labels) const {
|
237
|
+
FAISS_THROW_IF_NOT(k > 0);
|
238
|
+
|
239
|
+
FAISS_THROW_IF_NOT_MSG(
|
240
|
+
perm.size() == ntotal, "Call update_permutation before search");
|
276
241
|
|
277
242
|
#pragma omp parallel for
|
278
243
|
for (idx_t i = 0; i < n; i++) {
|
279
|
-
|
280
244
|
float q = x[i]; // query
|
281
|
-
float
|
282
|
-
idx_t
|
245
|
+
float* D = distances + i * k;
|
246
|
+
idx_t* I = labels + i * k;
|
283
247
|
|
284
248
|
// binary search
|
285
249
|
idx_t i0 = 0, i1 = ntotal;
|
@@ -297,8 +261,10 @@ void IndexFlat1D::search (
|
|
297
261
|
|
298
262
|
while (i0 + 1 < i1) {
|
299
263
|
idx_t imed = (i0 + i1) / 2;
|
300
|
-
if (xb[perm[imed]] <= q)
|
301
|
-
|
264
|
+
if (xb[perm[imed]] <= q)
|
265
|
+
i0 = imed;
|
266
|
+
else
|
267
|
+
i1 = imed;
|
302
268
|
}
|
303
269
|
|
304
270
|
// query is between xb[perm[i0]] and xb[perm[i1]]
|
@@ -311,13 +277,19 @@ void IndexFlat1D::search (
|
|
311
277
|
if (q - xleft < xright - q) {
|
312
278
|
D[wp] = q - xleft;
|
313
279
|
I[wp] = perm[i0];
|
314
|
-
i0--;
|
315
|
-
|
280
|
+
i0--;
|
281
|
+
wp++;
|
282
|
+
if (i0 < 0) {
|
283
|
+
goto finish_right;
|
284
|
+
}
|
316
285
|
} else {
|
317
286
|
D[wp] = xright - q;
|
318
287
|
I[wp] = perm[i1];
|
319
|
-
i1++;
|
320
|
-
|
288
|
+
i1++;
|
289
|
+
wp++;
|
290
|
+
if (i1 >= ntotal) {
|
291
|
+
goto finish_left;
|
292
|
+
}
|
321
293
|
}
|
322
294
|
}
|
323
295
|
goto done;
|
@@ -350,11 +322,8 @@ void IndexFlat1D::search (
|
|
350
322
|
}
|
351
323
|
wp++;
|
352
324
|
}
|
353
|
-
done
|
325
|
+
done:;
|
354
326
|
}
|
355
|
-
|
356
327
|
}
|
357
328
|
|
358
|
-
|
359
|
-
|
360
329
|
} // namespace faiss
|
@@ -14,33 +14,31 @@
|
|
14
14
|
|
15
15
|
#include <faiss/Index.h>
|
16
16
|
|
17
|
-
|
18
17
|
namespace faiss {
|
19
18
|
|
20
19
|
/** Index that stores the full vectors and performs exhaustive search */
|
21
|
-
struct IndexFlat: Index {
|
22
|
-
|
20
|
+
struct IndexFlat : Index {
|
23
21
|
/// database vectors, size ntotal * d
|
24
22
|
std::vector<float> xb;
|
25
23
|
|
26
|
-
explicit IndexFlat
|
24
|
+
explicit IndexFlat(idx_t d, MetricType metric = METRIC_L2);
|
27
25
|
|
28
26
|
void add(idx_t n, const float* x) override;
|
29
27
|
|
30
28
|
void reset() override;
|
31
29
|
|
32
30
|
void search(
|
33
|
-
|
34
|
-
|
35
|
-
|
36
|
-
|
37
|
-
|
31
|
+
idx_t n,
|
32
|
+
const float* x,
|
33
|
+
idx_t k,
|
34
|
+
float* distances,
|
35
|
+
idx_t* labels) const override;
|
38
36
|
|
39
37
|
void range_search(
|
40
|
-
|
41
|
-
|
42
|
-
|
43
|
-
|
38
|
+
idx_t n,
|
39
|
+
const float* x,
|
40
|
+
float radius,
|
41
|
+
RangeSearchResult* result) const override;
|
44
42
|
|
45
43
|
void reconstruct(idx_t key, float* recons) const override;
|
46
44
|
|
@@ -52,59 +50,51 @@ struct IndexFlat: Index {
|
|
52
50
|
* @param distances
|
53
51
|
* corresponding output distances, size n * k
|
54
52
|
*/
|
55
|
-
void compute_distance_subset
|
53
|
+
void compute_distance_subset(
|
56
54
|
idx_t n,
|
57
|
-
const float
|
55
|
+
const float* x,
|
58
56
|
idx_t k,
|
59
|
-
float
|
60
|
-
const idx_t
|
57
|
+
float* distances,
|
58
|
+
const idx_t* labels) const;
|
61
59
|
|
62
60
|
/** remove some ids. NB that Because of the structure of the
|
63
61
|
* indexing structure, the semantics of this operation are
|
64
62
|
* different from the usual ones: the new ids are shifted */
|
65
63
|
size_t remove_ids(const IDSelector& sel) override;
|
66
64
|
|
67
|
-
IndexFlat
|
65
|
+
IndexFlat() {}
|
68
66
|
|
69
|
-
DistanceComputer
|
67
|
+
DistanceComputer* get_distance_computer() const override;
|
70
68
|
|
71
69
|
/* The stanadlone codec interface (just memcopies in this case) */
|
72
|
-
size_t sa_code_size
|
73
|
-
|
74
|
-
void sa_encode (idx_t n, const float *x,
|
75
|
-
uint8_t *bytes) const override;
|
70
|
+
size_t sa_code_size() const override;
|
76
71
|
|
77
|
-
void
|
78
|
-
float *x) const override;
|
72
|
+
void sa_encode(idx_t n, const float* x, uint8_t* bytes) const override;
|
79
73
|
|
74
|
+
void sa_decode(idx_t n, const uint8_t* bytes, float* x) const override;
|
80
75
|
};
|
81
76
|
|
82
|
-
|
83
|
-
|
84
|
-
|
85
|
-
explicit IndexFlatIP (idx_t d): IndexFlat (d, METRIC_INNER_PRODUCT) {}
|
86
|
-
IndexFlatIP () {}
|
77
|
+
struct IndexFlatIP : IndexFlat {
|
78
|
+
explicit IndexFlatIP(idx_t d) : IndexFlat(d, METRIC_INNER_PRODUCT) {}
|
79
|
+
IndexFlatIP() {}
|
87
80
|
};
|
88
81
|
|
89
|
-
|
90
|
-
|
91
|
-
|
92
|
-
IndexFlatL2 () {}
|
82
|
+
struct IndexFlatL2 : IndexFlat {
|
83
|
+
explicit IndexFlatL2(idx_t d) : IndexFlat(d, METRIC_L2) {}
|
84
|
+
IndexFlatL2() {}
|
93
85
|
};
|
94
86
|
|
95
|
-
|
96
|
-
|
97
87
|
/// optimized version for 1D "vectors".
|
98
|
-
struct IndexFlat1D:IndexFlatL2 {
|
88
|
+
struct IndexFlat1D : IndexFlatL2 {
|
99
89
|
bool continuous_update; ///< is the permutation updated continuously?
|
100
90
|
|
101
91
|
std::vector<idx_t> perm; ///< sorted database indices
|
102
92
|
|
103
|
-
explicit IndexFlat1D
|
93
|
+
explicit IndexFlat1D(bool continuous_update = true);
|
104
94
|
|
105
95
|
/// if not continuous_update, call this between the last add and
|
106
96
|
/// the first search
|
107
|
-
void update_permutation
|
97
|
+
void update_permutation();
|
108
98
|
|
109
99
|
void add(idx_t n, const float* x) override;
|
110
100
|
|
@@ -112,14 +102,13 @@ struct IndexFlat1D:IndexFlatL2 {
|
|
112
102
|
|
113
103
|
/// Warn: the distances returned are L1 not L2
|
114
104
|
void search(
|
115
|
-
|
116
|
-
|
117
|
-
|
118
|
-
|
119
|
-
|
105
|
+
idx_t n,
|
106
|
+
const float* x,
|
107
|
+
idx_t k,
|
108
|
+
float* distances,
|
109
|
+
idx_t* labels) const override;
|
120
110
|
};
|
121
111
|
|
122
|
-
|
123
|
-
}
|
112
|
+
} // namespace faiss
|
124
113
|
|
125
114
|
#endif
|