linalg 1.0.2

data/lib/linalg/exception.rb

@@ -0,0 +1,32 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+module Linalg
+   module Exception
+
+      #
+      # The algorithm failed to converge.
+      #
+      class Diverged < RuntimeError
+      end
+
+      #
+      # One or more parameters failed to meet the
+      # prerequisite dimensions.
+      #
+      class DimensionError < RuntimeError
+      end
+
+      #
+      # A singular matrix was encountered where a nonsingular one
+      # was expected.
+      #
+      class SingularMatrix < RuntimeError
+      end
+
+   end
+end
+

data/lib/linalg/iterators.rb

@@ -0,0 +1,221 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+
+module Linalg
+   
+   #
+   # <tt>Enumerable</tt>s and <tt>Enumerable</tt>-like
+   # methods for matrices.
+   #
+   # For blocks which yield <tt>i,j</tt> indexes, the visitation is
+   # always in column-major order -- that is, +i+ changes faster than +j+.
+   #
+   # To visit in row-major order, one would say
+   #   m.rows.each_with_index { |row, i|
+   #      row.elems.each_with_index { |e, j| 
+   #         # e == m[i,j]
+   #      }
+   #   }
+   #
+   module Iterators
+
+      #
+      # Common +Enumerable+ for matrices
+      #
+      class MatrixEnum
+         include Enumerable
+         def initialize(m)
+            @m = m
+         end
+      end
+
+      #
+      # This class doesn't do much except +yield+ matrix elements.
+      #
+      class DiagEnum < MatrixEnum
+         def each
+            min = @m.hsize < @m.vsize ? @m.hsize : @m.vsize
+            min.times { |i|
+               yield @m[i,i]
+            }
+            @m
+         end
+      end
+
+      #
+      # This class doesn't do much except +yield+ matrix elements.
+      #
+      class ElemEnum < MatrixEnum
+         def each
+            (0...@m.hsize).each { |j|
+               (0...@m.vsize).each { |i|
+                  yield @m[i,j]
+               }
+            }
+            @m
+         end
+      end
+
+      #
+      # This class doesn't do much except +yield+ columns.
+      #
+      class ColumnEnum < MatrixEnum
+         def each
+            @m.hsize.times { |i|
+               yield @m.column(i)
+            }
+            @m
+         end
+      end
+
+      #
+      # This class doesn't do much except +yield+ rows.
+      #
+      class RowEnum < MatrixEnum
+         def each
+            @m.vsize.times { |i|
+               yield @m.row(i)
+            }
+            @m
+         end
+      end
+
+      #
+      # Returns an +Enumerable+ of the rows
+      #
+      def rows
+         RowEnum.new(self)
+      end
+
+      #
+      # Returns an +Enumerable+ of the columns
+      #
+      def columns
+         ColumnEnum.new(self)
+      end
+
+      #
+      # Returns an +Enumerable+ of the matrix elements
+      #
+      def elems
+         ElemEnum.new(self)
+      end
+
+      #
+      # Returns an +Enumerable+ of the diagonal elements
+      #
+      def diags
+         DiagEnum.new(self)
+      end
+
+      #
+      # Like <tt>Enumerable#each_with_index</tt>, but with an index pair
+      #
+      def each_with_index # :yields: e, i, j
+         (0...hsize).each { |j|
+            (0...vsize).each { |i|
+               yield self[i,j], i, j
+            }
+         }
+         self
+      end
+
+      #
+      # Like <tt>Array#each_index</tt>, but with an index pair
+      #
+      def each_index
+         (0...hsize).each { |j|
+            (0...vsize).each { |i|
+               yield i, j
+            }
+         }
+         self
+      end
+      
+      #
+      # Like <tt>Enumerable#map</tt>, but the resultant is a another matrix
+      #
+      def map # :yields: e
+         m = self.class.reserve(self.vsize, self.hsize)
+         each_with_index { |e, i, j|
+            m[i,j] = yield e
+         }
+         m
+      end
+
+      #
+      # +map+ with an index pair
+      #
+      def map_with_index # :yields: e, i, j
+         m = self.class.reserve(self.vsize, self.hsize)
+         each_with_index { |e, i, j|
+            m[i,j] = yield e, i, j
+         }
+         m
+      end
+
+      #
+      # In-place +map+
+      #
+      def map! # :yields: e
+         each_with_index { |e, i, j|
+            self[i,j] = yield e
+         }
+         self
+      end
+
+      #
+      # In-place +map_with_index+
+      #
+      def map_with_index! # :yields: e, i, j
+         each_with_index { |e, i, j|
+            self[i,j] = yield e, i, j
+         }
+         self
+      end
+
+      #
+      # Visit below-diagonal elements, with index pair
+      #
+      def each_lower_with_index # :yields: e, i, j
+         (0...hsize).each { |j|
+            ((j+1)...vsize).each { |i|
+               yield self[i,j], i, j
+            }
+         }
+         self
+      end
+
+      #
+      # Visit above-diagonal elements, with index pair
+      #
+      def each_upper_with_index # :yields: e, i, j
+         (1...hsize).each { |j|
+            (0...(j > vsize ? vsize : j)).each { |i|
+               yield self[i,j], i, j
+            }
+         }
+         self
+      end
+
+      private
+
+      # this seems overkill.
+      #
+      # Like <tt>Enumerable#inject</tt>, but with an index pair
+      #
+      #def inject_with_index(initial) # :yields: memo, e, i, j
+      #   res = initial
+      #   each_with_index { |e, i, j|
+      #      res = yield res, e, i, j
+      #   }
+      #   res
+      #end
+
+   end
+end
+   

data/lib/linalg/math.rb

@@ -0,0 +1,23 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+module Linalg
+   class Math
+      extend ::Math
+
+      def self.min(a, b)
+         a < b ? a : b
+      end
+
+      def self.max(a, b)
+         a > b ? a : b
+      end
+
+      def self.log2(x)
+         ::Math.log(x)/::Math.log(2.0)
+      end
+   end
+end

data/lib/linalg/scomplex.rb

@@ -0,0 +1,15 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+require 'linalg'
+
+module Linalg
+   class SComplex
+      def to_s
+         "#{real} + #{imag}i"
+      end
+   end
+end

data/lib/linalg/version.rb

@@ -0,0 +1,3 @@
+module Linalg
+  VERSION = "1.0.2"
+end

data/lib/linalg/xdata.rb

@@ -0,0 +1,123 @@
+#
+# Copyright (c) 2004-2008 by James M. Lawrence
+#
+# See LICENSE
+#
+
+require 'linalg'
+
+module Linalg
+   module XData
+      
+      class NullClass
+         def to_doublecomplex_ptr
+            0
+         end
+
+         def to_complex_ptr
+            0
+         end
+         
+         def to_doublereal_ptr
+            0
+         end
+
+         def to_real_ptr
+            0
+         end
+
+         def to_integer_ptr
+            0
+         end
+         
+         def to_logical_ptr
+            0
+         end
+
+         def to_char_ptr
+            0
+         end
+         
+         def to_L_fp
+            0
+         end
+      end
+
+      NULL = NullClass.new
+
+      class DReal < DData
+         def initialize(value = 0.0)
+            super(1)
+            self[0] = value
+         end
+
+         def value
+            self[0]
+         end
+         
+         def value=(a)
+            self[0] = a
+         end
+      end
+
+      class SReal < SData
+         def initialize(value = 0.0)
+            super(1)
+            self[0] = value
+         end
+
+         def value
+            self[0]
+         end
+         
+         def value=(a)
+            self[0] = a
+         end
+      end
+
+      class XInteger < IData
+         def initialize(value = 0)
+            super(1)
+            self[0] = value
+         end
+
+         def value
+            self[0]
+         end
+         
+         def value=(a)
+            self[0] = a
+         end
+      end
+
+      class Logical < LData
+         def initialize(value = 0)
+            super(1)
+            self[0] = value
+         end
+
+         def value
+            self[0]
+         end
+         
+         def value=(a)
+            self[0] = a
+         end
+      end
+
+      class Char < CharData
+         def initialize(value = "\0")
+            super(1)
+            self[0] = value
+         end
+
+         def value
+            self[0]
+         end
+         
+         def value=(a)
+            self[0] = a
+         end
+      end
+   end
+end

metadata

@@ -0,0 +1,94 @@
+--- !ruby/object:Gem::Specification
+name: linalg
+version: !ruby/object:Gem::Version
+  version: 1.0.2
+platform: ruby
+authors:
+- James M. Lawrence
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2014-01-04 00:00:00.000000000 Z
+dependencies: []
+description: A Fortran-based linear algebra package
+email:
+- quixoticsycophant@gmail.com
+executables: []
+extensions:
+- ext/lapack/extconf.rb
+- ext/linalg/extconf.rb
+extra_rdoc_files: []
+files:
+- LICENSE
+- ext/g2c_typedefs.h
+- ext/lapack/extconf.rb
+- ext/lapack/include/BLAS.h
+- ext/lapack/include/LAPACK.h
+- ext/lapack/main.c
+- ext/lapack/rb_lapack.h
+- ext/lapack/rb_lapack_c.c
+- ext/lapack/rb_lapack_d.c
+- ext/lapack/rb_lapack_s.c
+- ext/lapack/rb_lapack_x.c
+- ext/lapack/rb_lapack_z.c
+- ext/linalg/dcomplex.c
+- ext/linalg/dcomplex.h
+- ext/linalg/ddata.c
+- ext/linalg/extconf.rb
+- ext/linalg/linalg.c
+- ext/linalg/linalg.h
+- ext/linalg/xdata.c
+- ext/linalg/xdata.h
+- ext/linalg/xmatrix.c.tmpl
+- ext/linalg/xmatrix.h.tmpl
+- ext/linalg/xmatrixc.c.tmpl
+- ext/linalg/xmatrixr.c.tmpl
+- lib/lapack.rb
+- lib/linalg.rb
+- lib/linalg/dcomplex.rb
+- lib/linalg/dmatrix.rb
+- lib/linalg/dmatrix/alias.rb
+- lib/linalg/dmatrix/cholesky.rb
+- lib/linalg/dmatrix/cond.rb
+- lib/linalg/dmatrix/det.rb
+- lib/linalg/dmatrix/eigen.rb
+- lib/linalg/dmatrix/fit.rb
+- lib/linalg/dmatrix/inverse.rb
+- lib/linalg/dmatrix/lu.rb
+- lib/linalg/dmatrix/main.rb
+- lib/linalg/dmatrix/norms.rb
+- lib/linalg/dmatrix/nullspace.rb
+- lib/linalg/dmatrix/qr.rb
+- lib/linalg/dmatrix/schur.rb
+- lib/linalg/dmatrix/solve.rb
+- lib/linalg/dmatrix/svd.rb
+- lib/linalg/exception.rb
+- lib/linalg/iterators.rb
+- lib/linalg/math.rb
+- lib/linalg/scomplex.rb
+- lib/linalg/version.rb
+- lib/linalg/xdata.rb
+homepage: http://github.com/quix/linalg
+licenses: []
+metadata: {}
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 1.3.6
+requirements: []
+rubyforge_project: linalg
+rubygems_version: 2.2.0
+signing_key: 
+specification_version: 4
+summary: Ruby Linear Algebra Library
+test_files: []