matrix_vs_nmatrix.rb

require 'benchmark'
require 'matrix'
require 'nmatrix'

class Matrix
  public :"[]=", :set_element, :set_component
end

m = Matrix.zero 1000 ;
n = NMatrix.new [1000, 1000], [0]*10**6

Benchmark.bm do |x|
  x.report "Matrix getter" do
    (10**7).times { m[555, 555] }
  end

  x.report "NMatrix getter" do
    (10**7).times { n[555, 555] }
  end

  x.report "Matrix setter" do
    (10**7).times { m[555, 555] = 100 }
  end

  x.report "NMatrix setter" do
    (10**7).times { n[555, 555] = 100 }
  end
end