convolver 0.3.1 → 0.4.0

data/spec/convolver_fftw3_spec.rb

@@ -0,0 +1,166 @@
+# frozen_string_literal: true
+
+require 'helpers'
+
+describe Convolver do
+  describe '#convolve_fftw3' do
+    it 'works like the example in the README' do
+      a = NArray[0.3, 0.4, 0.5]
+      b = NArray[1.3, -0.5]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[0.19, 0.27]
+    end
+
+    it 'convolves 1D arrays with a variety of signal and kernel lengths' do
+      a = NArray[0.3]
+      b = NArray[-0.7]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[-0.21]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2]
+      b = NArray[-0.7]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[-0.21, -0.28, -0.35, -0.14]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2]
+      b = NArray[1.1, -0.7]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[0.05, 0.09, 0.41]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2]
+      b = NArray[1.1, -0.7, -0.2]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[-0.05, 0.05]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2, 0.6]
+      b = NArray[1.1, -0.7]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[0.05, 0.09, 0.41, -0.2]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2, 0.6]
+      b = NArray[1.1, -0.7, 2.1]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[1.1, 0.51, 1.67]
+
+      a = NArray[0.3, 0.4, 0.5, 0.2, 0.6]
+      b = NArray[0.6, -0.5, -0.4, 0.7]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[-0.08, 0.33]
+    end
+
+    it 'calculates a 2D convolution' do
+      a = NArray[[0.3, 0.4, 0.5], [0.6, 0.8, 0.2], [0.9, 1.0, 0.1]]
+      b = NArray[[1.2, -0.5], [0.5, -1.3]]
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[[-0.58, 0.37], [-0.53, 1.23]]
+    end
+
+    it 'calculates a 3D convolution' do
+      # 5x4x3
+      a = NArray[
+        [[1.0, 0.6, 1.1, 0.2, 0.9], [1.0, 0.7, 0.8, 1.0, 1.0], [0.2, 0.6, 0.1, 0.2, 0.5],
+         [0.5, 0.9, 0.2, 0.1, 0.6]],
+        [[0.4, 0.9, 0.4, 0.0, 0.6], [0.2, 1.1, 0.2, 0.4, 0.1], [0.4, 0.2, 0.5, 0.8, 0.7],
+         [0.1, 0.9, 0.7, 0.1, 0.3]],
+        [[0.8, 0.6, 1.0, 0.1, 0.4], [0.3, 0.8, 0.6, 0.7, 1.1], [0.9, 1.0, 0.3, 0.4, 0.6],
+         [0.2, 0.5, 0.4, 0.7, 0.2]]
+      ]
+
+      # 3x3x3
+      b = NArray[
+        [[-0.9, 1.2, 0.8], [0.9, 0.1, -0.5], [1.1, 0.1, -1.1]],
+        [[-0.2, -1.0, 1.4], [-1.4, 0.0, 1.3], [0.3, 1.0, -0.5]],
+        [[0.6, 0.0, 0.7],   [-0.7, 1.1, 1.2], [1.3, 0.7, 0.0]]
+      ]
+
+      # Should be 3x2x1
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[[[5.51, 3.04, 4.3], [3.04, 6.31, 3.87]]]
+    end
+
+    it 'calculates a 4D convolution' do
+      # 3x4x5x3
+      a = NArray[
+        [[[0.5, 0.4, 0.9], [0.1, 0.9, 0.8], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]]],
+        [[[0.5, 0.4, 0.9], [0.1, 0.9, 0.8], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]]],
+        [[[0.5, 0.4, 0.9], [0.1, 0.9, 0.8], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]],
+         [[0.0, 0.4, 0.0], [0.2, 0.3, 0.8], [0.6, 0.3, 0.2], [0.7, 0.4, 0.3]],
+         [[0.3, 0.3, 0.1], [0.6, 0.9, 0.4], [0.4, 0.0, 0.1], [0.8, 0.3, 0.4]]] ]
+
+      # 2x3x3x2
+      b = NArray[ [
+        [[1.1, 0.6], [1.2, 0.6], [0.8, 0.1]], [[-0.4, 0.8], [0.5, 0.4], [1.2, 0.2]],
+        [[0.8, 0.2], [0.5, 0.0], [1.4, 1.3]]
+      ],
+                  [[[1.1, 0.6], [1.2, 0.6], [0.8, 0.1]], [[-0.4, 0.8], [0.5, 0.4], [1.2, 0.2]],
+                   [[0.8, 0.2], [0.5, 0.0], [1.4, 1.3]]] ]
+
+      # Should be 2x2x3x2
+      c = described_class.convolve_fftw3(a, b)
+      expect(c).to be_narray_like NArray[
+        [[[8.5, 8.2], [11.34, 9.68]], [[7.68, 6.56], [11.24, 7.16]], [[9.14, 6.54], [12.44, 9.2]]],
+        [[[8.5, 8.2], [11.34, 9.68]], [[7.68, 6.56], [11.24, 7.16]], [[9.14, 6.54], [12.44, 9.2]]]
+      ]
+    end
+
+    describe 'compared with #convolve' do
+      it 'produces same results for 1D arrays' do
+        (1..30).each do |signal_length|
+          (1..signal_length).each do |kernel_length|
+            signal = NArray.sfloat(signal_length).random
+            kernel = NArray.sfloat(kernel_length).random
+            expect_result = described_class.convolve_basic(signal, kernel)
+            got_result = described_class.convolve_fftw3(signal, kernel)
+            expect(got_result).to be_narray_like expect_result
+          end
+        end
+      end
+
+      it 'produces same results for 2D arrays' do
+        (3..10).each do |signal_x|
+          ((signal_x - 2)..(signal_x + 2)).each do |signal_y|
+            (1..signal_x).each do |kernel_x|
+              (1..signal_y).each do |kernel_y|
+                signal = NArray.sfloat(signal_x, signal_y).random
+                kernel = NArray.sfloat(kernel_x, kernel_y).random
+                expect_result = described_class.convolve_basic(signal, kernel)
+                got_result = described_class.convolve_fftw3(signal, kernel)
+                expect(got_result).to be_narray_like expect_result
+              end
+            end
+          end
+        end
+      end
+
+      it 'produces same results for 3D arrays' do
+        (3..5).each do |signal_x|
+          ((signal_x - 2)..(signal_x + 2)).each do |signal_y|
+            ((signal_x - 2)..(signal_x + 2)).each do |signal_z|
+              (1..signal_x).each do |kernel_x|
+                (1..signal_y).each do |kernel_y|
+                  (1..signal_z).each do |kernel_z|
+                    signal = NArray.sfloat(signal_x, signal_y, signal_z).random
+                    kernel = NArray.sfloat(kernel_x, kernel_y, kernel_z).random
+                    expect_result = described_class.convolve_basic(signal, kernel)
+                    got_result = described_class.convolve_fftw3(signal, kernel)
+                    expect(got_result).to be_narray_like expect_result
+                  end
+                end
+              end
+            end
+          end
+        end
+      end
+    end
+  end
+end

data/spec/convolver_spec.rb

@@ -0,0 +1,53 @@
+# frozen_string_literal: true
+
+require 'helpers'
+
+describe Convolver do
+  describe '#convolve' do
+    it 'works like the example in the README' do
+      a = NArray[0.3, 0.4, 0.5]
+      b = NArray[1.3, -0.5]
+      c = described_class.convolve(a, b)
+      expect(c).to be_narray_like NArray[0.19, 0.27]
+    end
+
+    it 'processes convolutions of different sizes' do
+      # The variety here is to ensure all branches of optimisation algorithm
+      # are covered
+      [10, 30, 60, 90, 100, 120, 130, 150, 175, 200].each do |asize|
+        [5, 10, 12, 15, 20, 30, 40, 50].each do |bsize|
+          next unless bsize < asize
+
+          a = NArray.sfloat(asize, asize).random
+          b = NArray.sfloat(bsize, bsize).random
+          c = described_class.convolve(a, b)
+
+          # We should always match output of convolve_basic irrespective
+          # of what the optimal choice of algorithm is (larger error allowed here due to rounding)
+          expect_result = described_class.convolve_basic(a, b)
+          expect(c).to be_narray_like(expect_result, 1e-6)
+        end
+      end
+    end
+
+    it 'chooses #convolve_basic for small inputs' do
+      a = NArray.sfloat(50, 50).random
+      b = NArray.sfloat(10, 10).random
+      allow(described_class).to receive(:convolve_basic)
+      allow(described_class).to receive(:convolve_fftw3)
+      described_class.convolve(a, b)
+      expect(described_class).to have_received(:convolve_basic).once
+      expect(described_class).not_to have_received(:convolve_fftw3)
+    end
+
+    it 'chooses #convolve_fftw3 for large inputs' do
+      a = NArray.sfloat(500, 500).random
+      b = NArray.sfloat(100, 100).random
+      allow(described_class).to receive(:convolve_fftw3)
+      allow(described_class).to receive(:convolve_basic)
+      described_class.convolve(a, b)
+      expect(described_class).to have_received(:convolve_fftw3).once
+      expect(described_class).not_to have_received(:convolve_basic)
+    end
+  end
+end

data/spec/helpers.rb

@@ -1,43 +1,43 @@
+# frozen_string_literal: true
+
 # convolver/spec/helpers.rb
+require 'simplecov'
+SimpleCov.start
+
 require 'convolver'
-require 'mocha/api'
 
 # Matcher compares NArrays numerically
-RSpec::Matchers.define :be_narray_like do |expected_narray|
+RSpec::Matchers.define :be_narray_like do |expected_narray, mse = 1e-9|
   match do |given|
     @error = nil
-    if ! given.is_a?(NArray)
-      @error = "Wrong class."
+    if !given.is_a?(NArray)
+      @error = 'Wrong class.'
     elsif given.shape != expected_narray.shape
-      @error = "Shapes are different."
+      @error = 'Shapes are different.'
     else
       d = given - expected_narray
-      difference =  ( d * d ).sum / d.size
-      if difference > 1e-9
-        @error = "Numerical difference with mean square error #{difference}"
-      end
+      difference = (d * d).sum / d.size
+      @error = "Numerical difference with mean square error #{difference}" if difference > mse
     end
     @given = given.clone
 
-    if @error
-      @expected = expected_narray.clone
-    end
+    @expected = expected_narray.clone if @error
 
-    ! @error
+    !@error
   end
 
-  failure_message_for_should do
+  failure_message do
     "NArray does not match supplied example. #{@error}
     Expected: #{@expected.inspect}
     Got: #{@given.inspect}"
   end
 
-  failure_message_for_should_not do
+  failure_message_when_negated do
     "NArray is too close to unwanted example.
     Unwanted: #{@given.inspect}"
   end
 
-  description do |given, expected|
-    "numerically very close to example"
+  description do |_given, _expected|
+    'numerically very close to example'
   end
 end

metadata

@@ -1,127 +1,42 @@
 --- !ruby/object:Gem::Specification
 name: convolver
 version: !ruby/object:Gem::Version
-  version: 0.3.1
+  version: 0.4.0
 platform: ruby
 authors:
 - Neil Slater
-autorequire: 
 bindir: bin
 cert_chain: []
-date: 2013-11-14 00:00:00.000000000 Z
+date: 1980-01-02 00:00:00.000000000 Z
 dependencies:
-- !ruby/object:Gem::Dependency
-  name: narray
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.6.0.8
-  type: :runtime
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.6.0.8
 - !ruby/object:Gem::Dependency
   name: fftw3
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0.3'
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
         version: '0.3'
 - !ruby/object:Gem::Dependency
-  name: yard
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.8.7.2
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.8.7.2
-- !ruby/object:Gem::Dependency
-  name: bundler
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: '1.3'
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: '1.3'
-- !ruby/object:Gem::Dependency
-  name: rspec
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 2.13.0
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 2.13.0
-- !ruby/object:Gem::Dependency
-  name: mocha
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.14.0
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 0.14.0
-- !ruby/object:Gem::Dependency
-  name: rake
-  requirement: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 1.9.1
-  type: :development
-  prerelease: false
-  version_requirements: !ruby/object:Gem::Requirement
-    requirements:
-    - - '>='
-      - !ruby/object:Gem::Version
-        version: 1.9.1
-- !ruby/object:Gem::Dependency
-  name: rake-compiler
+  name: narray
   requirement: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
-        version: 0.8.3
-  type: :development
+        version: 0.6.0.8
+  type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
-    - - '>='
+    - - ">="
       - !ruby/object:Gem::Version
-        version: 0.8.3
+        version: 0.6.0.8
 description: Convolution for NArray
 email:
 - slobo777@gmail.com
@@ -130,8 +45,9 @@ extensions:
 - ext/convolver/extconf.rb
 extra_rdoc_files: []
 files:
-- .gitignore
-- .travis.yml
+- ".github/workflows/ci.yml"
+- ".gitignore"
+- ".rubocop.yml"
 - Gemfile
 - LICENSE.txt
 - README.md
@@ -145,37 +61,30 @@ files:
 - ext/convolver/narray_shared.h
 - lib/convolver.rb
 - lib/convolver/version.rb
-- spec/convolve_basic_spec.rb
-- spec/convolve_fftw3_spec.rb
-- spec/convolve_spec.rb
+- spec/convolver_basic_spec.rb
+- spec/convolver_fftw3_spec.rb
+- spec/convolver_spec.rb
 - spec/helpers.rb
 homepage: http://github.com/neilslater/convolver
 licenses:
 - MIT
-metadata: {}
-post_install_message: 
+metadata:
+  rubygems_mfa_required: 'true'
 rdoc_options: []
 require_paths:
 - lib
 required_ruby_version: !ruby/object:Gem::Requirement
   requirements:
-  - - '>='
+  - - ">="
     - !ruby/object:Gem::Version
-      version: '0'
+      version: '3.2'
 required_rubygems_version: !ruby/object:Gem::Requirement
   requirements:
-  - - '>='
+  - - ">="
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubyforge_project: 
-rubygems_version: 2.1.8
-signing_key: 
+rubygems_version: 4.0.1
 specification_version: 4
 summary: Convolution for NArray
-test_files:
-- spec/convolve_basic_spec.rb
-- spec/convolve_fftw3_spec.rb
-- spec/convolve_spec.rb
-- spec/helpers.rb
-has_rdoc: 
+test_files: []

data/.travis.yml

@@ -1,11 +0,0 @@
-before_install:
- - sudo apt-get update -qq
- - sudo apt-get install -qq libfftw3-dev
-language: ruby
-rvm:
-  - "1.8.7"
-  - "1.9.3"
-  - "2.0.0"
-  - rbx-18mode
-  - rbx-19mode
-  - ree

data/spec/convolve_basic_spec.rb

@@ -1,84 +0,0 @@
-require 'helpers'
-
-describe Convolver do
-  describe "#convolve_basic" do
-
-    it "should work like the example in the README" do
-      a = NArray[ 0.3, 0.4, 0.5 ]
-      b = NArray[ 1.3, -0.5 ]
-      c = Convolver.convolve_basic( a, b )
-      c.should be_narray_like NArray[ 0.19, 0.27 ]
-    end
-
-    it "should calculate a 2D convolution" do
-      a = NArray[ [ 0.3, 0.4, 0.5 ], [ 0.6, 0.8, 0.2 ], [ 0.9, 1.0, 0.1 ] ]
-      b = NArray[ [ 1.2, -0.5 ], [ 0.5, -1.3 ] ]
-      c = Convolver.convolve_basic( a, b )
-      c.should be_narray_like NArray[ [ -0.58, 0.37 ], [ -0.53, 1.23 ] ]
-    end
-
-    it "should calculate a 2D convolution with rectangular arrays" do
-      a = NArray[ [ 0.3, 0.4, 0.5, 0.3, 0.4 ], [ 0.6, 0.8, 0.2, 0.8, 0.2 ],
-        [ 0.9, 1.0, 0.1, 0.9, 1.0 ], [ 0.5, 0.9, 0.3, 0.2, 0.8 ], [ 0.7, 0.1, 0.3, 0.0, 0.1 ],
-        [ 0.4, 0.5, 0.6, 0.7, 0.8 ], [ 0.5, 0.4, 0.3, 0.2, 0.1 ] ]
-      b = NArray[ [ 1.2, -0.5, 0.2 ], [ 1.8, 0.5, -1.3 ] ]
-      c = Convolver.convolve_basic( a, b )
-      c.should be_narray_like NArray[ [ 1.48, 0.79, 1.03 ], [ 2.35, 1.7, -0.79 ], [ 1.56, 2.84, -0.53 ],
-        [ 1.13, 1.3, 0.83 ], [ 1.04, 0.26, 0.77 ], [ 1.06, 1.05, 1.04 ] ]
-    end
-
-    it "should calculate a 3D convolution" do
-      # 5x4x3
-      a = NArray[
-        [ [ 1.0, 0.6, 1.1, 0.2, 0.9 ], [ 1.0, 0.7, 0.8, 1.0, 1.0 ], [ 0.2, 0.6, 0.1, 0.2, 0.5 ], [ 0.5, 0.9, 0.2, 0.1, 0.6 ] ],
-        [ [ 0.4, 0.9, 0.4, 0.0, 0.6 ], [ 0.2, 1.1, 0.2, 0.4, 0.1 ], [ 0.4, 0.2, 0.5, 0.8, 0.7 ], [ 0.1, 0.9, 0.7, 0.1, 0.3 ] ],
-        [ [ 0.8, 0.6, 1.0, 0.1, 0.4 ], [ 0.3, 0.8, 0.6, 0.7, 1.1 ], [ 0.9, 1.0, 0.3, 0.4, 0.6 ], [ 0.2, 0.5, 0.4, 0.7, 0.2 ] ]
-      ]
-
-      # 3x3x3
-      b = NArray[
-        [ [ -0.9, 1.2, 0.8  ], [ 0.9, 0.1, -0.5 ], [ 1.1, 0.1, -1.1 ] ],
-        [ [ -0.2, -1.0, 1.4 ], [ -1.4, 0.0, 1.3 ], [ 0.3, 1.0, -0.5 ] ],
-        [ [ 0.6, 0.0, 0.7 ],   [ -0.7, 1.1, 1.2 ], [ 1.3, 0.7, 0.0  ] ]
-      ]
-
-      # Should be 3x2x1
-      c = Convolver.convolve_basic( a, b )
-      c.should be_narray_like NArray[ [ [ 5.51, 3.04, 4.3 ], [ 3.04, 6.31, 3.87 ] ] ]
-    end
-
-    it "should calculate a 4D convolution" do
-      # 3x4x5x3
-      a = NArray[
-        [ [ [ 0.5, 0.4, 0.9 ], [ 0.1, 0.9, 0.8 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ] ],
-        [ [ [ 0.5, 0.4, 0.9 ], [ 0.1, 0.9, 0.8 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ] ],
-        [ [ [ 0.5, 0.4, 0.9 ], [ 0.1, 0.9, 0.8 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ],
-          [ [ 0.0, 0.4, 0.0 ], [ 0.2, 0.3, 0.8 ], [ 0.6, 0.3, 0.2 ], [ 0.7, 0.4, 0.3 ] ],
-          [ [ 0.3, 0.3, 0.1 ], [ 0.6, 0.9, 0.4 ], [ 0.4, 0.0, 0.1 ], [ 0.8, 0.3, 0.4 ] ] ] ]
-
-      # 2x3x3x2
-      b = NArray[ [
-        [ [ 1.1, 0.6 ], [ 1.2, 0.6 ], [ 0.8, 0.1 ] ], [ [ -0.4, 0.8 ], [ 0.5, 0.4 ], [ 1.2, 0.2 ] ],
-        [ [ 0.8, 0.2 ], [ 0.5, 0.0 ], [ 1.4, 1.3 ] ] ],
-        [ [ [ 1.1, 0.6 ], [ 1.2, 0.6 ], [ 0.8, 0.1 ] ], [ [ -0.4, 0.8 ], [ 0.5, 0.4 ], [ 1.2, 0.2 ] ],
-        [ [ 0.8, 0.2 ], [ 0.5, 0.0 ], [ 1.4, 1.3 ] ] ] ]
-
-      # Should be 2x2x3x2
-      c = Convolver.convolve_basic( a, b )
-      c.should be_narray_like NArray[
-        [ [ [ 8.5, 8.2 ], [ 11.34, 9.68 ] ], [ [ 7.68, 6.56 ], [ 11.24, 7.16 ] ], [ [ 9.14, 6.54 ], [ 12.44, 9.2 ] ] ],
-        [ [ [ 8.5, 8.2 ], [ 11.34, 9.68 ] ], [ [ 7.68, 6.56 ], [ 11.24, 7.16 ] ], [ [ 9.14, 6.54 ], [ 12.44, 9.2 ] ] ]
-      ]
-    end
-  end
-end