faiss 0.2.6 → 0.2.7

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 5a12656e91180940d60fdd412b6f179cd0faa8a6b473aaadfa0bc4a371131464
-  data.tar.gz: 1076b01164ff1702abc8a9184410e0f8417199056698628ff2b0e5bf641cf168
+  metadata.gz: 8d0e67c5b23ef30f3632f8f878bba1ed8be7cbf8212c21a799daca8a65c23fcb
+  data.tar.gz: d3bcbf4453ece0531cb3cfa017bf60327fffa47401a27451ccb501c11431911e
 SHA512:
-  metadata.gz: 1298e0ed4b4455c478ee94384dbd02b4a029ea96409d66cf0b0e4292e51acde782d6bd875470c093a07198b99e222a366878d7019b6dc499221ee8b126a95cd7
-  data.tar.gz: 4b61b0fa059d240b568e176abdffff5b56a06e16598c185795e4d437bb0581a1cafb36794dd799168512b010ba71858837a30f3197be9b781510e50d3d4a268b
+  metadata.gz: 50f4c6ef3b637b708e0860f524eabb64b708868ce8210153f0868ff846a67ec24c4f229cf9146ab1dea0113b033e2287866ced15e750297ae79822cf9f131fe2
+  data.tar.gz: a3b0df6f2377c7381fdaab54f56e97a80ebaf5e752b92ba04b5d7a256c493a2bcbb14d60e4eeeb26a0d044f3159b4e416699ab6ecf774eb92f6e6014958d20b7

data/CHANGELOG.md

@@ -1,3 +1,7 @@
+## 0.2.7 (2023-04-20)
+
+- Updated Faiss to 1.7.4
+
 ## 0.2.6 (2023-04-11)
 
 - Added `add_with_ids` method

data/ext/faiss/extconf.rb

@@ -26,7 +26,7 @@ $CXXFLAGS << " -std=c++17 $(optflags) -DFINTEGER=int " << with_config("optflags"
 ext = File.expand_path(".", __dir__)
 vendor = File.expand_path("../../vendor/faiss", __dir__)
 
-$srcs = Dir["{#{ext},#{vendor}/faiss,#{vendor}/faiss/{impl,invlists,utils}}/*.{cpp}"]
+$srcs = Dir["{#{ext},#{vendor}/faiss,#{vendor}/faiss/{impl,invlists,utils}/**}/*.{cpp}"]
 $objs = $srcs.map { |v| v.sub(/cpp\z/, "o") }
 abort "Faiss not found" unless find_header("faiss/Index.h", vendor)
 $VPATH << vendor

data/lib/faiss/version.rb

@@ -1,3 +1,3 @@
 module Faiss
-  VERSION = "0.2.6"
+  VERSION = "0.2.7"
 end

data/lib/faiss.rb

@@ -2,7 +2,7 @@
 require "numo/narray"
 
 # ext
-require "faiss/ext"
+require_relative "faiss/ext"
 
 # modules
-require "faiss/version"
+require_relative "faiss/version"

data/vendor/faiss/faiss/AutoTune.cpp

@@ -32,6 +32,7 @@
 #include <faiss/IndexPreTransform.h>
 #include <faiss/IndexRefine.h>
 #include <faiss/IndexScalarQuantizer.h>
+#include <faiss/IndexShardsIVF.h>
 #include <faiss/MetaIndexes.h>
 #include <faiss/VectorTransform.h>
 
@@ -354,7 +355,7 @@ void ParameterSpace::initialize(const Index* index) {
         index = ix->index;
     }
 
-    if (DC(IndexIVF)) {
+    if (DC(IndexIVFInterface)) {
         {
             ParameterRange& pr = add_range("nprobe");
             for (int i = 0; i < 13; i++) {
@@ -461,6 +462,16 @@ void ParameterSpace::set_index_parameter(
         set_index_parameter(ix->index, name, val);
         return;
     }
+    if (DC(IndexShardsIVF)) {
+        // special handling because the nprobe is set at the sub-class level
+        // but other params are set on the class itself
+        if (name.find("quantizer_") == 0 && name != "nprobe" &&
+            name != "quantizer_nprobe") {
+            std::string sub_name = name.substr(strlen("quantizer_"));
+            set_index_parameter(ix->quantizer, sub_name, val);
+            return;
+        }
+    }
     if (DC(ThreadedIndex<Index>)) {
         // call on all sub-indexes
         auto fn = [this, name, val](int /* no */, Index* subIndex) {
@@ -608,7 +619,7 @@ void ParameterSpace::explore(
     if (n_experiments == 0) {
         for (size_t cno = 0; cno < n_comb; cno++) {
             set_index_parameters(index, cno);
-            std::vector<Index::idx_t> I(nq * crit.nnn);
+            std::vector<idx_t> I(nq * crit.nnn);
             std::vector<float> D(nq * crit.nnn);
 
             double t0 = getmillisecs();
@@ -677,7 +688,7 @@ void ParameterSpace::explore(
         }
 
         set_index_parameters(index, cno);
-        std::vector<Index::idx_t> I(nq * crit.nnn);
+        std::vector<idx_t> I(nq * crit.nnn);
         std::vector<float> D(nq * crit.nnn);
 
         double t0 = getmillisecs();
@@ -688,7 +699,7 @@ void ParameterSpace::explore(
         do {
             if (thread_over_batches) {
 #pragma omp parallel for
-                for (Index::idx_t q0 = 0; q0 < nq; q0 += batchsize) {
+                for (idx_t q0 = 0; q0 < nq; q0 += batchsize) {
                     size_t q1 = q0 + batchsize;
                     if (q1 > nq)
                         q1 = nq;

data/vendor/faiss/faiss/AutoTune.h

@@ -24,7 +24,6 @@ namespace faiss {
  * higher is better.
  */
 struct AutoTuneCriterion {
-    typedef Index::idx_t idx_t;
     idx_t nq;     ///< nb of queries this criterion is evaluated on
     idx_t nnn;    ///< nb of NNs that the query should request
     idx_t gt_nnn; ///< nb of GT NNs required to evaluate criterion

data/vendor/faiss/faiss/Clustering.cpp

@@ -88,8 +88,6 @@ void Clustering::train(
 
 namespace {
 
-using idx_t = Clustering::idx_t;
-
 idx_t subsample_training_set(
         const Clustering& clus,
         idx_t nx,
@@ -374,7 +372,7 @@ void Clustering::train_encoded(
     std::unique_ptr<float[]> dis(new float[nx]);
 
     // remember best iteration for redo
-    bool lower_is_better = index.metric_type != METRIC_INNER_PRODUCT;
+    bool lower_is_better = !is_similarity_metric(index.metric_type);
     float best_obj = lower_is_better ? HUGE_VALF : -HUGE_VALF;
     std::vector<ClusteringIterationStats> best_iteration_stats;
     std::vector<float> best_centroids;
@@ -624,8 +622,6 @@ ProgressiveDimClustering::ProgressiveDimClustering(
 
 namespace {
 
-using idx_t = Index::idx_t;
-
 void copy_columns(idx_t n, idx_t d1, const float* src, idx_t d2, float* dest) {
     idx_t d = std::min(d1, d2);
     for (idx_t i = 0; i < n; i++) {

data/vendor/faiss/faiss/Clustering.h

@@ -61,7 +61,6 @@ struct ClusteringIterationStats {
  *
  */
 struct Clustering : ClusteringParameters {
-    typedef Index::idx_t idx_t;
     size_t d; ///< dimension of the vectors
     size_t k; ///< nb of centroids
 
@@ -154,7 +153,6 @@ struct ProgressiveDimIndexFactory {
  * https://arxiv.org/abs/1509.05195
  */
 struct ProgressiveDimClustering : ProgressiveDimClusteringParameters {
-    using idx_t = Index::idx_t;
     size_t d; ///< dimension of the vectors
     size_t k; ///< nb of centroids
 

data/vendor/faiss/faiss/IVFlib.h

@@ -48,8 +48,6 @@ IndexIVF* try_extract_index_ivf(Index* index);
  */
 void merge_into(Index* index0, Index* index1, bool shift_ids);
 
-typedef Index::idx_t idx_t;
-
 /* Returns the cluster the embeddings belong to.
  *
  * @param index      Index, which should be an IVF index

data/vendor/faiss/faiss/Index.h

@@ -18,7 +18,7 @@
 
 #define FAISS_VERSION_MAJOR 1
 #define FAISS_VERSION_MINOR 7
-#define FAISS_VERSION_PATCH 3
+#define FAISS_VERSION_PATCH 4
 
 /**
  * @namespace faiss
@@ -62,7 +62,6 @@ struct SearchParameters {
  * although the internal representation may vary.
  */
 struct Index {
-    using idx_t = int64_t; ///< all indices are this type
     using component_t = float;
     using distance_t = float;
 

data/vendor/faiss/faiss/IndexAdditiveQuantizer.cpp

@@ -484,7 +484,7 @@ ResidualCoarseQuantizer::ResidualCoarseQuantizer(
         int d,        ///< dimensionality of the input vectors
         const std::vector<size_t>& nbits,
         MetricType metric)
-        : AdditiveCoarseQuantizer(d, &rq, metric), rq(d, nbits), beam_factor(4.0) {
+        : AdditiveCoarseQuantizer(d, &rq, metric), rq(d, nbits) {
     FAISS_THROW_IF_NOT(rq.tot_bits <= 63);
     is_trained = false;
 }
@@ -520,10 +520,15 @@ void ResidualCoarseQuantizer::search(
         idx_t k,
         float* distances,
         idx_t* labels,
-        const SearchParameters * params
+        const SearchParameters * params_in
         ) const {
 
-    FAISS_THROW_IF_NOT_MSG(!params, "search params not supported for this index");
+    float beam_factor = this->beam_factor;
+    if (params_in) {
+        auto params = dynamic_cast<const SearchParametersResidualCoarseQuantizer*>(params_in);
+        FAISS_THROW_IF_NOT_MSG(params, "need SearchParametersResidualCoarseQuantizer parameters");
+        beam_factor = params->beam_factor;
+    }
 
     if (beam_factor < 0) {
         AdditiveCoarseQuantizer::search(n, x, k, distances, labels);
@@ -585,6 +590,15 @@ void ResidualCoarseQuantizer::search(
     }
 }
 
+void ResidualCoarseQuantizer::initialize_from(const ResidualCoarseQuantizer &other) {
+    FAISS_THROW_IF_NOT(rq.M <= other.rq.M);
+    rq.initialize_from(other.rq);
+    set_beam_factor(other.beam_factor);
+    is_trained = other.is_trained;
+    ntotal = (idx_t)1 << aq->tot_bits;
+}
+
+
 /**************************************************************************************
  * LocalSearchCoarseQuantizer
  **************************************************************************************/

data/vendor/faiss/faiss/IndexAdditiveQuantizer.h

@@ -187,6 +187,11 @@ struct AdditiveCoarseQuantizer : Index {
     void reset() override;
 };
 
+struct SearchParametersResidualCoarseQuantizer : SearchParameters {
+    float beam_factor = 4.0f;
+    ~SearchParametersResidualCoarseQuantizer() {}
+};
+
 /** The ResidualCoarseQuantizer is a bit specialized compared to the
  * default AdditiveCoarseQuantizer because it can use a beam search
  * at search time (slow but may be useful for very large vocabularies) */
@@ -196,7 +201,7 @@ struct ResidualCoarseQuantizer : AdditiveCoarseQuantizer {
 
     /// factor between the beam size and the search k
     /// if negative, use exact search-to-centroid
-    float beam_factor;
+    float beam_factor = 4.0f;
 
     /// computes centroid norms if required
     void set_beam_factor(float new_beam_factor);
@@ -226,6 +231,10 @@ struct ResidualCoarseQuantizer : AdditiveCoarseQuantizer {
             idx_t* labels,
             const SearchParameters* params = nullptr) const override;
 
+    /** Copy the M first codebook levels from other. Useful to crop a
+     * ResidualQuantizer to its first M quantizers. */
+    void initialize_from(const ResidualCoarseQuantizer& other);
+
     ResidualCoarseQuantizer();
 };
 

data/vendor/faiss/faiss/IndexBinary.h

@@ -32,7 +32,6 @@ struct RangeSearchResult;
  * vectors.
  */
 struct IndexBinary {
-    using idx_t = Index::idx_t; ///< all indices are this type
     using component_t = uint8_t;
     using distance_t = int32_t;
 

data/vendor/faiss/faiss/IndexBinaryFlat.cpp

@@ -60,7 +60,8 @@ void IndexBinaryFlat::search(
                     xb.data(),
                     ntotal,
                     code_size,
-                    /* ordered = */ true);
+                    /* ordered = */ true,
+                    approx_topk_mode);
         } else {
             hammings_knn_mc(
                     x + s * code_size,

data/vendor/faiss/faiss/IndexBinaryFlat.h

@@ -14,6 +14,8 @@
 
 #include <faiss/IndexBinary.h>
 
+#include <faiss/utils/approx_topk/mode.h>
+
 namespace faiss {
 
 /** Index that stores the full vectors and performs exhaustive search. */
@@ -28,6 +30,8 @@ struct IndexBinaryFlat : IndexBinary {
 
     size_t query_batch_size = 32;
 
+    ApproxTopK_mode_t approx_topk_mode = ApproxTopK_mode_t::EXACT_TOPK;
+
     explicit IndexBinaryFlat(idx_t d);
 
     void add(idx_t n, const uint8_t* x) override;

data/vendor/faiss/faiss/IndexBinaryHash.cpp

@@ -108,8 +108,6 @@ struct FlipEnumerator {
     }
 };
 
-using idx_t = Index::idx_t;
-
 struct RangeSearchResults {
     int radius;
     RangeQueryResult& qres;
@@ -353,7 +351,7 @@ template <class HammingComputer, class SearchResults>
 static void verify_shortlist(
         const IndexBinaryFlat& index,
         const uint8_t* q,
-        const std::unordered_set<Index::idx_t>& shortlist,
+        const std::unordered_set<idx_t>& shortlist,
         SearchResults& res) {
     size_t code_size = index.code_size;
     size_t nlist = 0, ndis = 0, n0 = 0;