classifier 2.0.0 → 2.1.0

data/lib/classifier/lsi.rb

@@ -6,32 +6,55 @@
 
 module Classifier
   class LSI
-    # @rbs @gsl_available: bool
-    @gsl_available = false
+    # Backend options: :native, :ruby
+    # @rbs @backend: Symbol
+    @backend = :ruby
 
     class << self
-      # @rbs @gsl_available: bool
-      attr_accessor :gsl_available
+      # @rbs @backend: Symbol
+      attr_accessor :backend
+
+      # Check if using native C extension
+      # @rbs () -> bool
+      def native_available?
+        backend == :native
+      end
+
+      # Get the Vector class for the current backend
+      # @rbs () -> Class
+      def vector_class
+        backend == :native ? Classifier::Linalg::Vector : ::Vector
+      end
+
+      # Get the Matrix class for the current backend
+      # @rbs () -> Class
+      def matrix_class
+        backend == :native ? Classifier::Linalg::Matrix : ::Matrix
+      end
     end
   end
 end
 
+# Backend detection: native extension > pure Ruby
+# Set NATIVE_VECTOR=true to force pure Ruby implementation
+
 begin
-  # to test the native vector class, try `rake test NATIVE_VECTOR=true`
   raise LoadError if ENV['NATIVE_VECTOR'] == 'true'
-  raise LoadError unless Gem::Specification.find_all_by_name('gsl').any?
 
-  require 'gsl'
-  require 'classifier/extensions/vector_serialize'
-  Classifier::LSI.gsl_available = true
+  require 'classifier/classifier_ext'
+  Classifier::LSI.backend = :native
 rescue LoadError
-  unless ENV['SUPPRESS_GSL_WARNING'] == 'true'
-    warn 'Notice: for 10x faster LSI, run `gem install gsl`. Set SUPPRESS_GSL_WARNING=true to hide this.'
+  # Fall back to pure Ruby implementation
+  unless ENV['SUPPRESS_LSI_WARNING'] == 'true'
+    warn 'Notice: for 5-10x faster LSI, install the classifier gem with native extensions. ' \
+         'Set SUPPRESS_LSI_WARNING=true to hide this.'
   end
-  Classifier::LSI.gsl_available = false
+  Classifier::LSI.backend = :ruby
   require 'classifier/extensions/vector'
 end
 
+require 'json'
+require 'mutex_m'
 require 'classifier/lsi/word_list'
 require 'classifier/lsi/content_node'
 require 'classifier/lsi/summary'
@@ -41,14 +64,19 @@ module Classifier
   # data based on underlying semantic relations. For more information on the algorithms used,
   # please consult Wikipedia[http://en.wikipedia.org/wiki/Latent_Semantic_Indexing].
   class LSI
+    include Mutex_m
+
     # @rbs @auto_rebuild: bool
     # @rbs @word_list: WordList
     # @rbs @items: Hash[untyped, ContentNode]
     # @rbs @version: Integer
     # @rbs @built_at_version: Integer
+    # @rbs @singular_values: Array[Float]?
+    # @rbs @dirty: bool
+    # @rbs @storage: Storage::Base?
 
-    attr_reader :word_list
-    attr_accessor :auto_rebuild
+    attr_reader :word_list, :singular_values
+    attr_accessor :auto_rebuild, :storage
 
     # Create a fresh index.
     # If you want to call #build_index manually, use
@@ -56,11 +84,14 @@ module Classifier
     #
     # @rbs (?Hash[Symbol, untyped]) -> void
     def initialize(options = {})
+      super()
       @auto_rebuild = true unless options[:auto_rebuild] == false
       @word_list = WordList.new
       @items = {}
       @version = 0
       @built_at_version = -1
+      @dirty = false
+      @storage = nil
     end
 
     # Returns true if the index needs to be rebuilt.  The index needs
@@ -69,7 +100,26 @@ module Classifier
     #
     # @rbs () -> bool
     def needs_rebuild?
-      (@items.keys.size > 1) && (@version != @built_at_version)
+      synchronize { (@items.keys.size > 1) && (@version != @built_at_version) }
+    end
+
+    # @rbs () -> Array[Hash[Symbol, untyped]]?
+    def singular_value_spectrum
+      return nil unless @singular_values
+
+      total = @singular_values.sum
+      return nil if total.zero?
+
+      cumulative = 0.0
+      @singular_values.map.with_index do |value, i|
+        cumulative += value
+        {
+          dimension: i,
+          value: value,
+          percentage: value / total,
+          cumulative_percentage: cumulative / total
+        }
+      end
     end
 
     # Adds an item to the index. item is assumed to be a string, but
@@ -88,8 +138,11 @@ module Classifier
     # @rbs (String, *String | Symbol) ?{ (String) -> String } -> void
     def add_item(item, *categories, &block)
       clean_word_hash = block ? block.call(item).clean_word_hash : item.to_s.clean_word_hash
-      @items[item] = ContentNode.new(clean_word_hash, *categories)
-      @version += 1
+      synchronize do
+        @items[item] = ContentNode.new(clean_word_hash, *categories)
+        @version += 1
+        @dirty = true
+      end
       build_index if @auto_rebuild
     end
 
@@ -107,25 +160,32 @@ module Classifier
     #
     # @rbs (String) -> Array[String | Symbol]
     def categories_for(item)
-      return [] unless @items[item]
+      synchronize do
+        return [] unless @items[item]
 
-      @items[item].categories
+        @items[item].categories
+      end
     end
 
     # Removes an item from the database, if it is indexed.
     #
     # @rbs (String) -> void
     def remove_item(item)
-      return unless @items.key?(item)
+      removed = synchronize do
+        next false unless @items.key?(item)
 
-      @items.delete(item)
-      @version += 1
+        @items.delete(item)
+        @version += 1
+        @dirty = true
+        true
+      end
+      build_index if removed && @auto_rebuild
     end
 
     # Returns an array of items that are indexed.
     # @rbs () -> Array[untyped]
     def items
-      @items.keys
+      synchronize { @items.keys }
     end
 
     # This function rebuilds the index if needs_rebuild? returns true.
@@ -145,38 +205,30 @@ module Classifier
     #
     # @rbs (?Float) -> void
     def build_index(cutoff = 0.75)
-      return unless needs_rebuild?
-
-      make_word_list
-
-      doc_list = @items.values
-      tda = doc_list.collect { |node| node.raw_vector_with(@word_list) }
-
-      if self.class.gsl_available
-        tdm = GSL::Matrix.alloc(*tda).trans
-        ntdm = build_reduced_matrix(tdm, cutoff)
-
-        ntdm.size[1].times do |col|
-          vec = GSL::Vector.alloc(ntdm.column(col)).row
-          doc_list[col].lsi_vector = vec
-          doc_list[col].lsi_norm = vec.normalize
+      validate_cutoff!(cutoff)
+
+      synchronize do
+        return unless needs_rebuild_unlocked?
+
+        make_word_list
+
+        doc_list = @items.values
+        tda = doc_list.collect { |node| node.raw_vector_with(@word_list) }
+
+        if self.class.native_available?
+          # Convert vectors to arrays for matrix construction
+          tda_arrays = tda.map { |v| v.respond_to?(:to_a) ? v.to_a : v }
+          tdm = self.class.matrix_class.alloc(*tda_arrays).trans
+          ntdm = build_reduced_matrix(tdm, cutoff)
+          assign_native_ext_lsi_vectors(ntdm, doc_list)
+        else
+          tdm = Matrix.rows(tda).trans
+          ntdm = build_reduced_matrix(tdm, cutoff)
+          assign_ruby_lsi_vectors(ntdm, doc_list)
         end
-      else
-        tdm = Matrix.rows(tda).trans
-        ntdm = build_reduced_matrix(tdm, cutoff)
-
-        ntdm.column_size.times do |col|
-          next unless doc_list[col]
-
-          column = ntdm.column(col)
-          next unless column
 
-          doc_list[col].lsi_vector = column
-          doc_list[col].lsi_norm = column.normalize
-        end
+        @built_at_version = @version
       end
-
-      @built_at_version = @version
     end
 
     # This method returns max_chunks entries, ordered by their average semantic rating.
@@ -190,12 +242,14 @@ module Classifier
     #
     # @rbs (?Integer) -> Array[String]
     def highest_relative_content(max_chunks = 10)
-      return [] if needs_rebuild?
+      synchronize do
+        return [] if needs_rebuild_unlocked?
 
-      avg_density = {}
-      @items.each_key { |x| avg_density[x] = proximity_array_for_content(x).sum { |pair| pair[1] } }
+        avg_density = {}
+        @items.each_key { |x| avg_density[x] = proximity_array_for_content_unlocked(x).sum { |pair| pair[1] } }
 
-      avg_density.keys.sort_by { |x| avg_density[x] }.reverse[0..(max_chunks - 1)].map
+        avg_density.keys.sort_by { |x| avg_density[x] }.reverse[0..(max_chunks - 1)].map
+      end
     end
 
     # This function is the primitive that find_related and classify
@@ -212,20 +266,8 @@ module Classifier
     # text data. See add_item for examples of how this works.
     #
     # @rbs (String) ?{ (String) -> String } -> Array[[String, Float]]
-    def proximity_array_for_content(doc, &)
-      return [] if needs_rebuild?
-
-      content_node = node_for_content(doc, &)
-      result =
-        @items.keys.collect do |item|
-          val = if self.class.gsl_available
-                  content_node.search_vector * @items[item].search_vector.col
-                else
-                  (Matrix[content_node.search_vector] * @items[item].search_vector)[0]
-                end
-          [item, val]
-        end
-      result.sort_by { |x| x[1] }.reverse
+    def proximity_array_for_content(doc, &block)
+      synchronize { proximity_array_for_content_unlocked(doc, &block) }
     end
 
     # Similar to proximity_array_for_content, this function takes similar
@@ -235,20 +277,8 @@ module Classifier
     # the text you're working with. search uses this primitive.
     #
     # @rbs (String) ?{ (String) -> String } -> Array[[String, Float]]
-    def proximity_norms_for_content(doc, &)
-      return [] if needs_rebuild?
-
-      content_node = node_for_content(doc, &)
-      result =
-        @items.keys.collect do |item|
-          val = if self.class.gsl_available
-                  content_node.search_norm * @items[item].search_norm.col
-                else
-                  (Matrix[content_node.search_norm] * @items[item].search_norm)[0]
-                end
-          [item, val]
-        end
-      result.sort_by { |x| x[1] }.reverse
+    def proximity_norms_for_content(doc, &block)
+      synchronize { proximity_norms_for_content_unlocked(doc, &block) }
     end
 
     # This function allows for text-based search of your index. Unlike other functions
@@ -261,11 +291,13 @@ module Classifier
     #
     # @rbs (String, ?Integer) -> Array[String]
     def search(string, max_nearest = 3)
-      return [] if needs_rebuild?
+      synchronize do
+        return [] if needs_rebuild_unlocked?
 
-      carry = proximity_norms_for_content(string)
-      result = carry.collect { |x| x[0] }
-      result[0..(max_nearest - 1)]
+        carry = proximity_norms_for_content_unlocked(string)
+        result = carry.collect { |x| x[0] }
+        result[0..(max_nearest - 1)]
+      end
     end
 
     # This function takes content and finds other documents
@@ -280,10 +312,12 @@ module Classifier
     #
     # @rbs (String, ?Integer) ?{ (String) -> String } -> Array[String]
     def find_related(doc, max_nearest = 3, &block)
-      carry =
-        proximity_array_for_content(doc, &block).reject { |pair| pair[0] == doc }
-      result = carry.collect { |x| x[0] }
-      result[0..(max_nearest - 1)]
+      synchronize do
+        carry =
+          proximity_array_for_content_unlocked(doc, &block).reject { |pair| pair[0] == doc }
+        result = carry.collect { |x| x[0] }
+        result[0..(max_nearest - 1)]
+      end
     end
 
     # This function uses a voting system to categorize documents, based on
@@ -291,32 +325,23 @@ module Classifier
     # find_related function to find related documents, then returns the
     # most obvious category from this list.
     #
-    # cutoff signifies the number of documents to consider when clasifying
-    # text. A cutoff of 1 means that every document in the index votes on
-    # what category the document is in. This may not always make sense.
-    #
     # @rbs (String, ?Float) ?{ (String) -> String } -> String | Symbol
-    def classify(doc, cutoff = 0.30, &)
-      votes = vote(doc, cutoff, &)
+    def classify(doc, cutoff = 0.30, &block)
+      validate_cutoff!(cutoff)
+
+      synchronize do
+        votes = vote_unlocked(doc, cutoff, &block)
 
-      ranking = votes.keys.sort_by { |x| votes[x] }
-      ranking[-1]
+        ranking = votes.keys.sort_by { |x| votes[x] }
+        ranking[-1]
+      end
     end
 
     # @rbs (String, ?Float) ?{ (String) -> String } -> Hash[String | Symbol, Float]
-    def vote(doc, cutoff = 0.30, &)
-      icutoff = (@items.size * cutoff).round
-      carry = proximity_array_for_content(doc, &)
-      carry = carry[0..(icutoff - 1)]
-      votes = {}
-      carry.each do |pair|
-        categories = @items[pair[0]].categories
-        categories.each do |category|
-          votes[category] ||= 0.0
-          votes[category] += pair[1]
-        end
-      end
-      votes
+    def vote(doc, cutoff = 0.30, &block)
+      validate_cutoff!(cutoff)
+
+      synchronize { vote_unlocked(doc, cutoff, &block) }
     end
 
     # Returns the same category as classify() but also returns
@@ -331,15 +356,19 @@ module Classifier
     #
     # See classify() for argument docs
     # @rbs (String, ?Float) ?{ (String) -> String } -> [String | Symbol | nil, Float?]
-    def classify_with_confidence(doc, cutoff = 0.30, &)
-      votes = vote(doc, cutoff, &)
-      votes_sum = votes.values.sum
-      return [nil, nil] if votes_sum.zero?
-
-      ranking = votes.keys.sort_by { |x| votes[x] }
-      winner = ranking[-1]
-      vote_share = votes[winner] / votes_sum.to_f
-      [winner, vote_share]
+    def classify_with_confidence(doc, cutoff = 0.30, &block)
+      validate_cutoff!(cutoff)
+
+      synchronize do
+        votes = vote_unlocked(doc, cutoff, &block)
+        votes_sum = votes.values.sum
+        return [nil, nil] if votes_sum.zero?
+
+        ranking = votes.keys.sort_by { |x| votes[x] }
+        winner = ranking[-1]
+        vote_share = votes[winner] / votes_sum.to_f
+        [winner, vote_share]
+      end
     end
 
     # Prototype, only works on indexed documents.
@@ -347,45 +376,314 @@ module Classifier
     # it's supposed to.
     # @rbs (String, ?Integer) -> Array[Symbol]
     def highest_ranked_stems(doc, count = 3)
-      raise 'Requested stem ranking on non-indexed content!' unless @items[doc]
+      synchronize do
+        raise 'Requested stem ranking on non-indexed content!' unless @items[doc]
+
+        arr = node_for_content_unlocked(doc).lsi_vector.to_a
+        top_n = arr.sort.reverse[0..(count - 1)]
+        top_n.collect { |x| @word_list.word_for_index(arr.index(x)) }
+      end
+    end
+
+    # Custom marshal serialization to exclude mutex state
+    # @rbs () -> Array[untyped]
+    def marshal_dump
+      [@auto_rebuild, @word_list, @items, @version, @built_at_version, @dirty]
+    end
+
+    # Custom marshal deserialization to recreate mutex
+    # @rbs (Array[untyped]) -> void
+    def marshal_load(data)
+      mu_initialize
+      @auto_rebuild, @word_list, @items, @version, @built_at_version, @dirty = data
+      @storage = nil
+    end
+
+    # Returns a hash representation of the LSI index.
+    # Only source data (word_hash, categories) is included, not computed vectors.
+    # This can be converted to JSON or used directly.
+    #
+    # @rbs () -> untyped
+    def as_json(*)
+      items_data = @items.transform_values do |node|
+        {
+          word_hash: node.word_hash.transform_keys(&:to_s),
+          categories: node.categories.map(&:to_s)
+        }
+      end
+
+      {
+        version: 1,
+        type: 'lsi',
+        auto_rebuild: @auto_rebuild,
+        items: items_data
+      }
+    end
+
+    # Serializes the LSI index to a JSON string.
+    # Only source data (word_hash, categories) is serialized, not computed vectors.
+    # On load, the index will be rebuilt automatically.
+    #
+    # @rbs () -> String
+    def to_json(*)
+      as_json.to_json
+    end
+
+    # Loads an LSI index from a JSON string or Hash created by #to_json or #as_json.
+    # The index will be rebuilt after loading.
+    #
+    # @rbs (String | Hash[String, untyped]) -> LSI
+    def self.from_json(json)
+      data = json.is_a?(String) ? JSON.parse(json) : json
+      raise ArgumentError, "Invalid classifier type: #{data['type']}" unless data['type'] == 'lsi'
+
+      # Create instance with auto_rebuild disabled during loading
+      instance = new(auto_rebuild: false)
+
+      # Restore items (categories stay as strings, matching original storage)
+      data['items'].each do |item_key, item_data|
+        word_hash = item_data['word_hash'].transform_keys(&:to_sym)
+        categories = item_data['categories']
+        instance.instance_variable_get(:@items)[item_key] = ContentNode.new(word_hash, *categories)
+        instance.instance_variable_set(:@version, instance.instance_variable_get(:@version) + 1)
+      end
+
+      # Restore auto_rebuild setting and rebuild index
+      instance.auto_rebuild = data['auto_rebuild']
+      instance.build_index
+      instance
+    end
+
+    # Saves the LSI index to the configured storage.
+    # Raises ArgumentError if no storage is configured.
+    #
+    # @rbs () -> void
+    def save
+      raise ArgumentError, 'No storage configured. Use save_to_file(path) or set storage=' unless storage
+
+      storage.write(to_json)
+      @dirty = false
+    end
+
+    # Saves the LSI index to a file (legacy API).
+    #
+    # @rbs (String) -> Integer
+    def save_to_file(path)
+      result = File.write(path, to_json)
+      @dirty = false
+      result
+    end
+
+    # Reloads the LSI index from the configured storage.
+    # Raises UnsavedChangesError if there are unsaved changes.
+    # Use reload! to force reload and discard changes.
+    #
+    # @rbs () -> self
+    def reload
+      raise ArgumentError, 'No storage configured' unless storage
+      raise UnsavedChangesError, 'Unsaved changes would be lost. Call save first or use reload!' if @dirty
 
-      arr = node_for_content(doc).lsi_vector.to_a
-      top_n = arr.sort.reverse[0..(count - 1)]
-      top_n.collect { |x| @word_list.word_for_index(arr.index(x)) }
+      data = storage.read
+      raise StorageError, 'No saved state found' unless data
+
+      restore_from_json(data)
+      @dirty = false
+      self
+    end
+
+    # Force reloads the LSI index from storage, discarding any unsaved changes.
+    #
+    # @rbs () -> self
+    def reload!
+      raise ArgumentError, 'No storage configured' unless storage
+
+      data = storage.read
+      raise StorageError, 'No saved state found' unless data
+
+      restore_from_json(data)
+      @dirty = false
+      self
+    end
+
+    # Returns true if there are unsaved changes.
+    #
+    # @rbs () -> bool
+    def dirty?
+      @dirty
+    end
+
+    # Loads an LSI index from the configured storage.
+    # The storage is set on the returned instance.
+    #
+    # @rbs (storage: Storage::Base) -> LSI
+    def self.load(storage:)
+      data = storage.read
+      raise StorageError, 'No saved state found' unless data
+
+      instance = from_json(data)
+      instance.storage = storage
+      instance
+    end
+
+    # Loads an LSI index from a file (legacy API).
+    #
+    # @rbs (String) -> LSI
+    def self.load_from_file(path)
+      from_json(File.read(path))
     end
 
     private
 
+    # Restores LSI state from a JSON string (used by reload)
+    # @rbs (String) -> void
+    def restore_from_json(json)
+      data = JSON.parse(json)
+      raise ArgumentError, "Invalid classifier type: #{data['type']}" unless data['type'] == 'lsi'
+
+      synchronize do
+        # Recreate the items
+        @items = {}
+        data['items'].each do |item_key, item_data|
+          word_hash = item_data['word_hash'].transform_keys(&:to_sym)
+          categories = item_data['categories']
+          @items[item_key] = ContentNode.new(word_hash, *categories)
+        end
+
+        # Restore settings
+        @auto_rebuild = data['auto_rebuild']
+        @version += 1
+        @built_at_version = -1
+        @word_list = WordList.new
+        @dirty = false
+      end
+
+      # Rebuild the index
+      build_index
+    end
+
+    # @rbs (Float) -> void
+    def validate_cutoff!(cutoff)
+      return if cutoff.positive? && cutoff < 1
+
+      raise ArgumentError, "cutoff must be between 0 and 1 (exclusive), got #{cutoff}"
+    end
+
+    # Assigns LSI vectors using native C extension
+    # @rbs (untyped, Array[ContentNode]) -> void
+    def assign_native_ext_lsi_vectors(ntdm, doc_list)
+      ntdm.size[1].times do |col|
+        vec = self.class.vector_class.alloc(ntdm.column(col).to_a).row
+        doc_list[col].lsi_vector = vec
+        doc_list[col].lsi_norm = vec.normalize
+      end
+    end
+
+    # Assigns LSI vectors using pure Ruby Matrix
+    # @rbs (untyped, Array[ContentNode]) -> void
+    def assign_ruby_lsi_vectors(ntdm, doc_list)
+      ntdm.column_size.times do |col|
+        next unless doc_list[col]
+
+        column = ntdm.column(col)
+        next unless column
+
+        doc_list[col].lsi_vector = column
+        doc_list[col].lsi_norm = column.normalize
+      end
+    end
+
+    # Unlocked version of needs_rebuild? for internal use when lock is already held
+    # @rbs () -> bool
+    def needs_rebuild_unlocked?
+      (@items.keys.size > 1) && (@version != @built_at_version)
+    end
+
+    # Unlocked version of proximity_array_for_content for internal use
+    # @rbs (String) ?{ (String) -> String } -> Array[[String, Float]]
+    def proximity_array_for_content_unlocked(doc, &)
+      return [] if needs_rebuild_unlocked?
+
+      content_node = node_for_content_unlocked(doc, &)
+      result =
+        @items.keys.collect do |item|
+          val = if self.class.native_available?
+                  content_node.search_vector * @items[item].search_vector.col
+                else
+                  (Matrix[content_node.search_vector] * @items[item].search_vector)[0]
+                end
+          [item, val]
+        end
+      result.sort_by { |x| x[1] }.reverse
+    end
+
+    # Unlocked version of proximity_norms_for_content for internal use
+    # @rbs (String) ?{ (String) -> String } -> Array[[String, Float]]
+    def proximity_norms_for_content_unlocked(doc, &)
+      return [] if needs_rebuild_unlocked?
+
+      content_node = node_for_content_unlocked(doc, &)
+      result =
+        @items.keys.collect do |item|
+          val = if self.class.native_available?
+                  content_node.search_norm * @items[item].search_norm.col
+                else
+                  (Matrix[content_node.search_norm] * @items[item].search_norm)[0]
+                end
+          [item, val]
+        end
+      result.sort_by { |x| x[1] }.reverse
+    end
+
+    # Unlocked version of vote for internal use
+    # @rbs (String, ?Float) ?{ (String) -> String } -> Hash[String | Symbol, Float]
+    def vote_unlocked(doc, cutoff = 0.30, &)
+      icutoff = (@items.size * cutoff).round
+      carry = proximity_array_for_content_unlocked(doc, &)
+      carry = carry[0..(icutoff - 1)]
+      votes = {}
+      carry.each do |pair|
+        categories = @items[pair[0]].categories
+        categories.each do |category|
+          votes[category] ||= 0.0
+          votes[category] += pair[1]
+        end
+      end
+      votes
+    end
+
+    # Unlocked version of node_for_content for internal use
+    # @rbs (String) ?{ (String) -> String } -> ContentNode
+    def node_for_content_unlocked(item, &block)
+      return @items[item] if @items[item]
+
+      clean_word_hash = block ? block.call(item).clean_word_hash : item.to_s.clean_word_hash
+      cn = ContentNode.new(clean_word_hash, &block)
+      cn.raw_vector_with(@word_list) unless needs_rebuild_unlocked?
+      cn
+    end
+
     # @rbs (untyped, ?Float) -> untyped
     def build_reduced_matrix(matrix, cutoff = 0.75)
       # TODO: Check that M>=N on these dimensions! Transpose helps assure this
       u, v, s = matrix.SV_decomp
 
-      # TODO: Better than 75% term, please. :\
-      s_cutoff = s.sort.reverse[(s.size * cutoff).round - 1]
+      @singular_values = s.sort.reverse
+
+      s_cutoff_index = [(s.size * cutoff).round - 1, 0].max
+      s_cutoff = @singular_values[s_cutoff_index]
       s.size.times do |ord|
         s[ord] = 0.0 if s[ord] < s_cutoff
       end
       # Reconstruct the term document matrix, only with reduced rank
-      result = u * (self.class.gsl_available ? GSL::Matrix : ::Matrix).diag(s) * v.trans
+      result = u * self.class.matrix_class.diag(s) * v.trans
 
-      # Native Ruby SVD returns transposed dimensions when row_size < column_size
+      # SVD may return transposed dimensions when row_size < column_size
       # Ensure result matches input dimensions
-      result = result.trans if !self.class.gsl_available && result.row_size != matrix.row_size
+      result = result.trans if result.row_size != matrix.row_size
 
       result
     end
 
-    # @rbs (String) ?{ (String) -> String } -> ContentNode
-    def node_for_content(item, &block)
-      return @items[item] if @items[item]
-
-      clean_word_hash = block ? block.call(item).clean_word_hash : item.to_s.clean_word_hash
-      cn = ContentNode.new(clean_word_hash, &block)
-      cn.raw_vector_with(@word_list) unless needs_rebuild?
-      cn
-    end
-
     # @rbs () -> void
     def make_word_list
       @word_list = WordList.new