nmatrix-lapacke 0.2.1 → 0.2.3

data/spec/plugins/lapacke/lapacke_spec.rb

@@ -298,6 +298,193 @@ describe "NMatrix::LAPACK functions implemented with LAPACKE interface" do
         expect(vr).to be_within(err).of(vr_true)
         expect(vl).to be_within(err).of(vl_true)
       end
+      
+      it "exposes lapacke_geqrf" do
+        a = NMatrix.new(3, [12.0, -51.0,   4.0, 
+                             6.0, 167.0, -68.0, 
+                            -4.0,  24.0, -41.0] , dtype: dtype)
+
+        b = NMatrix.new([3,1], 0, dtype: dtype)
+
+        NMatrix::LAPACK::lapacke_geqrf(:row, a.shape[0], a.shape[1], a, a.shape[1], b)
+
+        x = NMatrix.new([3,1], TAU_SOLUTION_ARRAY, dtype: dtype)
+     
+        y = NMatrix.new([3,3], GEQRF_SOLUTION_ARRAY, dtype: dtype)
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+        
+        expect(b).to be_within(err).of(x)
+        expect(a).to be_within(err).of(y)      
+      end
+
+      it "calculates QR decomposition in a compressed format using geqrf!" do
+        a = NMatrix.new(3, [12.0, -51.0,   4.0, 
+                             6.0, 167.0, -68.0, 
+                            -4.0,  24.0, -41.0] , dtype: dtype)
+
+        tau = a.geqrf!
+    
+        x = NMatrix.new([3,1], TAU_SOLUTION_ARRAY, dtype: dtype)
+     
+        y = NMatrix.new([3,3], GEQRF_SOLUTION_ARRAY, dtype: dtype)
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+        
+        expect(tau).to be_within(err).of(x)
+        expect(a).to be_within(err).of(y)      
+      end
+
+      it "exposes lapacke_ormqr and lapacke_unmqr" do
+        a = NMatrix.new([4,2], [34.0,  21.0, 
+                                23.0,  53.0, 
+                                26.0, 346.0, 
+                                23.0, 121.0] , dtype: dtype)
+
+        tau = NMatrix.new([2,1], dtype: dtype)
+        result = NMatrix.identity(4, dtype: dtype)
+        
+        # get tau from geqrf, use for ormqr  
+        NMatrix::LAPACK::lapacke_geqrf(:row, a.shape[0], a.shape[1], a, a.shape[1], tau)
+
+        #Q is stored in result 
+        a.complex_dtype? ?
+          NMatrix::LAPACK::lapacke_unmqr(:row, :left, false, result.shape[0], result.shape[1], tau.shape[0], 
+                                                                a, a.shape[1], tau, result, result.shape[1])
+          :
+
+          NMatrix::LAPACK::lapacke_ormqr(:row, :left, false, result.shape[0], result.shape[1], tau.shape[0], 
+                                                                a, a.shape[1], tau, result, result.shape[1])
+
+        x = NMatrix.new([4,4], Q_SOLUTION_ARRAY_1, dtype: dtype)
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+
+        expect(result).to be_within(err).of(x)      
+      end
+
+      it "calculates the product of the orthogonal matrix with an arbitrary matrix" do
+        a = N.new([2,2], [34.0, 21, 23, 53] , dtype: dtype)
+
+        tau = NMatrix.new([2,1], dtype: dtype)
+        
+        #Result is the multiplicand that gets overriden : result = Q * result
+        result   = NMatrix.new([2,2], [2,0,0,2], dtype: dtype)
+        
+        # get tau from geqrf, use for ormqr  
+        NMatrix::LAPACK::lapacke_geqrf(:row, a.shape[0], a.shape[1], a, a.shape[1], tau)
+
+        #Q is stored in result 
+        a.complex_dtype? ?
+          NMatrix::LAPACK::lapacke_unmqr(:row, :left, false, result.shape[0], result.shape[1], tau.shape[0], 
+                                                                a, a.shape[1], tau, result, result.shape[1])
+          :
+
+          NMatrix::LAPACK::lapacke_ormqr(:row, :left, false, result.shape[0], result.shape[1], tau.shape[0], 
+                                                                a, a.shape[1], tau, result, result.shape[1])
+
+        x = NMatrix.new([2,2], [-1.6565668262559257 , -1.1206187354084205, 
+                                -1.1206187354084205 , 1.6565668262559263], dtype: dtype)
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+
+        expect(result).to be_within(err).of(x)      
+      end
+      
+      it "calculates the orthogonal matrix Q using ormqr/unmqr after geqrf!" do
+        a = NMatrix.new([4,2], [34.0,  21.0, 
+                                23.0,  53.0, 
+                                26.0, 346.0, 
+                                23.0, 121.0] , dtype: dtype)
+        
+        # get tau from geqrf, use for ormqr  
+        tau = a.geqrf!
+
+        #Q is stored in result 
+        result = a.complex_dtype? ? a.unmqr(tau) : a.ormqr(tau)
+          
+
+        x = NMatrix.new([4,4], Q_SOLUTION_ARRAY_1, dtype: dtype)
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+
+        expect(result).to be_within(err).of(x)      
+      end
+    end
+    
+    it "calculates the transpose of Q using ormqr/unmqr after geqrf!" do
+        a = NMatrix.new([4,2], [34.0,  21.0, 
+                                23.0,  53.0, 
+                                26.0, 346.0, 
+                                23.0, 121.0] , dtype: dtype)
+        
+        # get tau from geqrf, use for ormqr  
+        tau = a.geqrf!
+
+        #Q is stored in result 
+        result = a.complex_dtype? ? a.unmqr(tau, :left, :complex_conjugate) : a.ormqr(tau, :left, :transpose)
+          
+
+        x = NMatrix.new([4,4], Q_SOLUTION_ARRAY_1, dtype: dtype)
+        x = x.transpose
+
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+
+        expect(result).to be_within(err).of(x)      
+    end
+
+    it "calculates the multiplication c * Q using ormqr/unmqr after geqrf!" do
+        a = NMatrix.new(3, [12.0, -51.0,   4.0, 
+                             6.0, 167.0, -68.0, 
+                            -4.0,  24.0, -41.0] , dtype: dtype)
+        
+        # get tau from geqrf, use for ormqr  
+        tau = a.geqrf!
+        c = NMatrix.new([2,3], [1,0,1,0,0,1], dtype: dtype)
+
+        #Q is stored in result 
+        result = a.complex_dtype? ? a.unmqr(tau, :right, false, c) : a.ormqr(tau, :right, false, c)
+          
+        solution = NMatrix.new([2,3], [-0.5714285714285714,   0.2228571428571429, 1.2742857142857142,
+                                        0.28571428571428575, -0.1714285714285714, 0.9428571428571428] , dtype: dtype)
+        err = case dtype
+                when :float32, :complex64
+                  1e-4
+                when :float64, :complex128
+                  1e-14
+              end
+
+        expect(result).to be_within(err).of(solution)      
     end
   end
 end

data/spec/shortcuts_spec.rb

@@ -50,10 +50,45 @@ describe NMatrix do
     expect(m).to eq identity3
   end
 
+  it "hilbert() creates an hilbert matrix" do
+    m = NMatrix.hilbert(8)
+    expect(m[4, 0]).to be_within(0.000001).of(0.2)
+    expect(m[4, 1]).to be_within(0.000001).of(0.16666666666666666)
+    expect(m[4, 2]).to be_within(0.000001).of(0.14285714285714285)
+    expect(m[4, 3]).to be_within(0.000001).of(0.125)
+
+    m = NMatrix.hilbert(3)
+    hilbert3 = NMatrix.new([3, 3], [1.0, 0.5, 0.3333333333333333,\
+     0.5, 0.3333333333333333, 0.25, 0.3333333333333333, 0.25, 0.2])
+    expect(m).to eq hilbert3
+    0.upto(2) do |i|
+      0.upto(2) do |j|
+        expect(m[i, j]).to be_within(0.000001).of(hilbert3[i,j])
+      end
+    end
+  end
+
+  it "inv_hilbert() creates an inverse hilbert matrix" do
+    m = NMatrix.inv_hilbert(6)
+    inv_hilbert6 = [3360.0,  -88200.0,   564480.0, -1411200.0]
+    expect(m[2,0]).to be_within(0.000001).of(inv_hilbert6[0])
+    expect(m[2,1]).to be_within(0.000001).of(inv_hilbert6[1])
+    expect(m[2,2]).to be_within(0.000001).of(inv_hilbert6[2])
+    expect(m[2,3]).to be_within(0.000001).of(inv_hilbert6[3])
+
+    m = NMatrix.inv_hilbert(3)
+    inv_hilbert3 = NMatrix.new([3, 3], [  9.0,  -36.0,   30.0, -36.0,  192.0, -180.0, 30.0, -180.0,  180.0] )
+    0.upto(2) do |i|
+      0.upto(2) do |j|
+        expect(m[i, j]).to be_within(0.000001).of(inv_hilbert3[i,j])
+      end
+    end
+  end
+
   it "diag() creates a matrix with pre-supplied diagonal" do
     arr = [1,2,3,4]
     m = NMatrix.diag(arr)
-    expect(m.is_a?(NMatrix)).to be_true
+    expect(m.is_a?(NMatrix)).to be true
   end
 
   it "diagonals() contains the seeded values on the diagonal" do
@@ -123,6 +158,112 @@ describe NMatrix do
       expect { NMatrix.random("not an array or integer") }.to raise_error
     end
   end
+  
+  context "::magic" do
+    
+    ALL_DTYPES.each do |dtype|
+      context dtype do
+        it "creates a matrix with numbers from 1 to n^n(n squared)" do
+          a = NMatrix.magic(3, dtype: dtype)
+          magic3 = NMatrix.new([3,3], [4, 9, 2, 3, 5, 7, 8, 1, 6], dtype: dtype)	
+          expect(a).to eq magic3
+      
+          b = NMatrix.magic(4, dtype: dtype)
+          magic4 = NMatrix.new([4,4], [1, 15, 14,  4, 12,  6,  7, 9, 8, 10, 11, 5, 13, 3, 2, 16], dtype: dtype)
+          expect(b).to eq magic4
+      
+          c = NMatrix.magic(6, dtype: dtype)
+          magic6 = NMatrix.new([6,6], [31, 9, 2, 22, 27, 20, 3, 32, 7, 21, 23, 25, 35, 1, 6, 26, 19, 24, 4, 36, 29, 13, 18, 11, 30, 5, 34, 12, 14, 16, 8, 28, 33, 17, 10, 15], dtype: dtype)
+          expect(c).to eq magic6  
+        end
+      end
+    end
+    
+    it "shape of two is not allowed" do
+      expect { NMatrix.magic(2) }.to raise_error(ArgumentError)
+    end
+    
+    it "Only accepts an integer as dimension" do 
+      expect { NMatrix.magic(3.0) }.to raise_error(ArgumentError)
+    end
+  end
+  
+  context "::linspace" do
+    it "creates a row vector when given only one shape parameter" do
+      v = NMatrix.linspace(1, 10, 4)
+      #Expect a row vector only
+      expect(v.shape.length).to eq(1)
+      
+      ans = [1.0,4.0,7.0,10.0]
+
+      expect(v[0]).to be_within(0.000001).of(ans[0])
+      expect(v[1]).to be_within(0.000001).of(ans[1])
+      expect(v[2]).to be_within(0.000001).of(ans[2])
+      expect(v[3]).to be_within(0.000001).of(ans[3])
+    end
+    
+    it "creates a matrix of input shape with each entry linearly spaced in row major order" do
+      v = NMatrix.linspace(1, Math::PI, [2,2])
+      expect(v.dtype).to eq(:float64)
+
+      ans = [1.0, 1.7138642072677612, 2.4277284145355225, 3.1415927410125732]
+      
+      expect(v[0,0]).to be_within(0.000001).of(ans[0])
+      expect(v[0,1]).to be_within(0.000001).of(ans[1])
+      expect(v[1,0]).to be_within(0.000001).of(ans[2])
+      expect(v[1,1]).to be_within(0.000001).of(ans[3])
+    end
+  end
+
+  context "::logspace" do
+    it "creates a logarithmically spaced vector" do
+      v = NMatrix.logspace(1, 2, 10)
+      
+      expect(v.shape.length).to eq(1)
+      
+      #Unit test taken from Matlab R2015b output of logspace(1,2,10)
+      ans = [10.0000, 12.9155, 16.6810, 21.5443, 27.8256, 35.9381, 46.4159, 59.9484, 77.4264, 100.0000]
+      
+      expect(v[0].round(4)).to be_within(0.000001).of(ans[0])
+      expect(v[1].round(4)).to be_within(0.000001).of(ans[1])
+      expect(v[2].round(4)).to be_within(0.000001).of(ans[2])
+      expect(v[3].round(4)).to be_within(0.000001).of(ans[3])
+      expect(v[4].round(4)).to be_within(0.000001).of(ans[4])
+      expect(v[5].round(4)).to be_within(0.000001).of(ans[5])
+      expect(v[6].round(4)).to be_within(0.000001).of(ans[6])
+      expect(v[7].round(4)).to be_within(0.000001).of(ans[7])
+      expect(v[8].round(4)).to be_within(0.000001).of(ans[8])
+      expect(v[9].round(4)).to be_within(0.000001).of(ans[9])
+    end
+    
+    it "creates a logarithmically spaced vector bounded by Math::PI if :pi is pre-supplied" do
+      v = NMatrix.logspace(1, :pi, 7)
+                 
+      #Unit test taken from Matlab R2015b output of logspace(1,pi,10)
+      ans = [10.0000, 8.2450, 6.7980, 5.6050, 4.6213, 3.8103, 3.1416]
+      
+      expect(v[0].round(4)).to be_within(0.000001).of(ans[0])
+      expect(v[1].round(4)).to be_within(0.000001).of(ans[1])
+      expect(v[2].round(4)).to be_within(0.000001).of(ans[2])
+      expect(v[3].round(4)).to be_within(0.000001).of(ans[3])
+      expect(v[4].round(4)).to be_within(0.000001).of(ans[4])
+      expect(v[5].round(4)).to be_within(0.000001).of(ans[5])
+      expect(v[6].round(4)).to be_within(0.000001).of(ans[6])
+    end    
+
+    it "creates a matrix of input shape with each entry logarithmically spaced in row major order" do
+      v = NMatrix.logspace(1, 2, [3,2])
+      
+      ans = [10.0, 15.8489, 25.1189, 39.8107, 63.0957, 100.0]
+      
+      expect(v[0,0].round(4)).to be_within(0.000001).of(ans[0])
+      expect(v[0,1].round(4)).to be_within(0.000001).of(ans[1])
+      expect(v[1,0].round(4)).to be_within(0.000001).of(ans[2])
+      expect(v[1,1].round(4)).to be_within(0.000001).of(ans[3])
+      expect(v[2,0].round(4)).to be_within(0.000001).of(ans[4])
+      expect(v[2,1].round(4)).to be_within(0.000001).of(ans[5])
+    end
+  end
 
   it "seq() creates a matrix of integers, sequentially" do
     m = NMatrix.seq(2) # 2x2 matrix.
@@ -136,7 +277,6 @@ describe NMatrix do
     end
   end
 
-
   it "indgen() creates a matrix of integers as well as seq()" do
     m = NMatrix.indgen(2) # 2x2 matrix.
     value = 0
@@ -189,19 +329,19 @@ describe NMatrix do
   it "column() returns a NMatrix" do
     m = NMatrix.random(3)
 
-    expect(m.column(2).is_a?(NMatrix)).to be_true
+    expect(m.column(2).is_a?(NMatrix)).to be true
   end
 
   it "row() returns a NMatrix" do
     m = NMatrix.random(3)
 
-    expect(m.row(2).is_a?(NMatrix)).to be_true
+    expect(m.row(2).is_a?(NMatrix)).to be true
   end
 
   it "diagonals() creates an NMatrix" do
     arr = [1,2,3,4]
     m = NMatrix.diagonals(arr)
-    expect(m.is_a?(NMatrix)).to be_true
+    expect(m.is_a?(NMatrix)).to be true
   end
 
   it "diagonals() contains the seeded values on the diagonal" do

data/spec/spec_helper.rb

@@ -32,7 +32,7 @@ require "./lib/nmatrix/rspec"
 
 ALL_DTYPES = [:byte,:int8,:int16,:int32,:int64, :float32,:float64, :object,
   :complex64, :complex128]
-
+  
 NON_INTEGER_DTYPES = [:float32, :float64, :complex64, :complex128,
   :object]
 
@@ -48,6 +48,29 @@ COMPLEX_MATRIX32A_ARRAY = MATRIX32A_ARRAY.zip(MATRIX32A_ARRAY.reverse).collect {
 GETRF_EXAMPLE_ARRAY = [-1,0,10,4,9,2,3,5,7,8,1,6]
 GETRF_SOLUTION_ARRAY = [9.0, 2.0, 3.0, 5.0, 7.0/9, 58.0/9, -4.0/3, 19.0/9, -1.0/9, 1.0/29, 301.0/29, 130.0/29]
 
+TAU_SOLUTION_ARRAY = [1.8571428571428572,1.9938461538461538, 0.0]
+
+GEQRF_SOLUTION_ARRAY =[                -14.0,                -21.0, 14.000000000000002,
+                         0.23076923076923078,  -175.00000000000003,  70.00000000000001,
+                        -0.15384615384615385, 0.055555555555555546,              -35.0]
+
+R_SOLUTION_ARRAY   = [-159.2388143638353, -41.00131005172065, -56.75123892439876,  -90.75048729628048, 
+                                     0.0, 25.137473501580676,  2.073591725046292,   9.790607357775713, 
+                                     0.0,                0.0, -20.83259700334131, -17.592414929551445]
+
+Q_SOLUTION_ARRAY_1 = [-0.632455532033676, -0.5209522876558295, -0.3984263084135902,  -0.41214704991068,
+                    -0.42783756578748666, -0.20837937347171134, 0.876505919951498, 0.07259770177184455,
+                    -0.48364246567281094, 0.8265854747306287,-0.015758658987033422, -0.2873988222474053,
+                    -0.42783756578748666,  0.044081783789183565, -0.26971376257215296, 0.8615487797670971]
+
+Q_SOLUTION_ARRAY_2 = [-0.8571428571428572,   0.3942857142857143,  0.33142857142857146, 
+                      -0.4285714285714286,  -0.9028571428571428, -0.03428571428571425, 
+                       0.28571428571428575, -0.1714285714285714,   0.9428571428571428]
+
+Q_SOLUTION_ARRAY_3 = [-0.7724247413634004, -0.026670393594597247, -0.6345460653374136, 
+                      -0.5777485870360393,  -0.38541856437557026,  0.7194853024298236,
+                      -0.26375478973384403,   0.9223563413020934, 0.28229805268947933]
+
 def create_matrix(stype) #:nodoc:
   m = NMatrix.new([3,3], 0, dtype: :int32, stype: stype, default: 0)
 

metadata

@@ -1,16 +1,15 @@
 --- !ruby/object:Gem::Specification
 name: nmatrix-lapacke
 version: !ruby/object:Gem::Version
-  version: 0.2.1
+  version: 0.2.3
 platform: ruby
 authors:
+- Will Levine
 - John Woods
-- Chris Wailes
-- Aleksey Timin
 autorequire: 
 bindir: bin
 cert_chain: []
-date: 2016-01-18 00:00:00.000000000 Z
+date: 2016-07-25 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: nmatrix
@@ -18,14 +17,14 @@ dependencies:
     requirements:
     - - '='
       - !ruby/object:Gem::Version
-        version: 0.2.1
+        version: 0.2.3
   type: :runtime
   prerelease: false
   version_requirements: !ruby/object:Gem::Requirement
     requirements:
     - - '='
       - !ruby/object:Gem::Version
-        version: 0.2.1
+        version: 0.2.3
 description: For using linear algebra fuctions provided by LAPACK and BLAS
 email:
 - john.o.woods@gmail.com
@@ -48,6 +47,7 @@ files:
 - ext/nmatrix/math/imax.h
 - ext/nmatrix/math/laswp.h
 - ext/nmatrix/math/long_dtype.h
+- ext/nmatrix/math/magnitude.h
 - ext/nmatrix/math/math.h
 - ext/nmatrix/math/nrm2.h
 - ext/nmatrix/math/rot.h
@@ -84,6 +84,8 @@ files:
 - ext/nmatrix_lapacke/lapacke/include/lapacke_utils.h
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cgeev.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cgeev_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_cgeqrf.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_cgeqrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cgesdd.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cgesdd_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cgesvd.c
@@ -100,8 +102,12 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cpotri_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cpotrs.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_cpotrs_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_cunmqr.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_cunmqr_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgeev.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgeev_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_dgeqrf.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_dgeqrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgesdd.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgesdd_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgesvd.c
@@ -112,6 +118,8 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgetri_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgetrs.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dgetrs_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_dormqr.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_dormqr_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dpotrf.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dpotrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dpotri.c
@@ -120,6 +128,8 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_dpotrs_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgeev.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgeev_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_sgeqrf.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_sgeqrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgesdd.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgesdd_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgesvd.c
@@ -130,6 +140,8 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgetri_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgetrs.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_sgetrs_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_sormqr.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_sormqr_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_spotrf.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_spotrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_spotri.c
@@ -138,6 +150,8 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_spotrs_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zgeev.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zgeev_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_zgeqrf.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_zgeqrf_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zgesdd.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zgesdd_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zgesvd.c
@@ -154,12 +168,16 @@ files:
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zpotri_work.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zpotrs.c
 - ext/nmatrix_lapacke/lapacke/src/lapacke_zpotrs_work.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_zunmqr.c
+- ext/nmatrix_lapacke/lapacke/src/lapacke_zunmqr_work.c
+- ext/nmatrix_lapacke/lapacke/utils/lapacke_c_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_cge_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_cge_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_cpo_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_cpo_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_ctr_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_ctr_trans.c
+- ext/nmatrix_lapacke/lapacke/utils/lapacke_d_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_dge_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_dge_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_dpo_nancheck.c
@@ -167,6 +185,7 @@ files:
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_dtr_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_dtr_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_lsame.c
+- ext/nmatrix_lapacke/lapacke/utils/lapacke_s_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_sge_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_sge_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_spo_nancheck.c
@@ -174,6 +193,7 @@ files:
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_str_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_str_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_xerbla.c
+- ext/nmatrix_lapacke/lapacke/utils/lapacke_z_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_zge_nancheck.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_zge_trans.c
 - ext/nmatrix_lapacke/lapacke/utils/lapacke_zpo_nancheck.c
@@ -219,7 +239,7 @@ files:
 - spec/utm5940.mtx
 homepage: http://sciruby.com
 licenses:
-- BSD 3-clause
+- BSD-3-Clause
 metadata: {}
 post_install_message: 
 rdoc_options: []
@@ -237,7 +257,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
       version: '0'
 requirements: []
 rubyforge_project: 
-rubygems_version: 2.4.5
+rubygems_version: 2.5.1
 signing_key: 
 specification_version: 4
 summary: general LAPACK backend for nmatrix using LAPACKE interface
@@ -246,17 +266,27 @@ test_files:
 - spec/01_enum_spec.rb
 - spec/02_slice_spec.rb
 - spec/03_nmatrix_monkeys_spec.rb
+- spec/2x2_dense_double.mat
+- spec/4x4_sparse.mat
+- spec/4x5_dense.mat
 - spec/blas_spec.rb
 - spec/elementwise_spec.rb
 - spec/homogeneous_spec.rb
 - spec/io/fortran_format_spec.rb
 - spec/io/harwell_boeing_spec.rb
+- spec/io/test.rua
 - spec/io_spec.rb
 - spec/lapack_core_spec.rb
+- spec/leakcheck.rb
 - spec/math_spec.rb
+- spec/nmatrix_yale_resize_test_associations.yaml
 - spec/nmatrix_yale_spec.rb
+- spec/rspec_monkeys.rb
 - spec/rspec_spec.rb
 - spec/shortcuts_spec.rb
 - spec/slice_set_spec.rb
+- spec/spec_helper.rb
 - spec/stat_spec.rb
+- spec/test.pcd
+- spec/utm5940.mtx
 - spec/plugins/lapacke/lapacke_spec.rb