numo-linalg 0.0.1

data/lib/numo/linalg/linalg.rb

@@ -0,0 +1,7 @@
+require "numo/narray"
+require "numo/linalg/blas"
+require "numo/linalg/lapack"
+require "numo/linalg/function"
+require "numo/linalg/site_conf"
+require "numo/linalg/loader"
+require "numo/linalg/version"

data/lib/numo/linalg/loader.rb

@@ -0,0 +1,174 @@
+require "fiddle"
+
+module Numo
+  module Linalg
+    module Loader
+
+      @@libs = nil
+
+      module_function
+
+      def libs
+        @@libs
+      end
+
+      def load_mkl(path=nil,exc:true)
+        a = case path
+            when String; [path]
+            when Array ;  path
+            when NilClass; [MKL_LIBPATH,""].compact
+            else raise ArgumentError,"invalid path"
+            end
+        a.each do |d|
+          if d.empty?
+            f_iomp5 = "libiomp5.so"
+            f_mkl_core = "libmkl_core.so"
+            f_intel_thread = "libmkl_intel_thread.so"
+            f_intel_lp64 = "libmkl_intel_lp64.so"
+          else
+            e = d.sub(/\/mkl\//, "/")
+            f_iomp5 = File.join(e,"libiomp5.so")
+            f_mkl_core = File.join(d,"libmkl_core.so")
+            f_intel_thread = File.join(d,"libmkl_intel_thread.so")
+            f_intel_lp64 = File.join(d,"libmkl_intel_lp64.so")
+          end
+          begin
+            Fiddle.dlopen(f_iomp5)
+            Fiddle.dlopen(f_mkl_core)
+            Fiddle.dlopen(f_intel_thread)
+            Blas.dlopen(f_intel_lp64)
+            Lapack.dlopen(f_intel_lp64)
+            @@libs = [ f_iomp5, f_mkl_core, f_intel_thread, f_intel_lp64 ]
+            if $DEBUG
+              $stderr.puts "Numo::Linalg: use #{f_intel_lp64}"
+            end
+            return true
+          rescue
+          end
+        end
+        if exc
+          raise RuntimeError, "cannot find MKL library"
+        end
+        false
+      end
+
+      def load_openblas(path=nil,exc:true)
+        a = case path
+            when String; [path]
+            when Array ;  path
+            when NilClass; [OPENBLAS_LIBPATH,""].compact
+            else raise ArgumentError,"invalid path"
+            end
+        a.each do |d|
+          %w[libopenblaso.so libopenblasp.so libopenblas.so].each do |f|
+            f = File.join(d,f) if !d.empty?
+            begin
+              Blas.dlopen(f)
+              Lapack.dlopen(f)
+              @@libs = [ f ]
+              if $DEBUG
+                $stderr.puts "Numo::Linalg: use #{f}"
+              end
+              return true
+            rescue
+            end
+          end
+        end
+        if exc
+          raise RuntimeError, "cannot find OpenBLAS library"
+        end
+        false
+      end
+
+      def load_atlas(path=nil,exc:true)
+        a = case path
+            when String; [path]
+            when Array ;  path
+            when NilClass; [ATLAS_LIBPATH,""].compact
+            else raise ArgumentError,"invalid path"
+            end
+        a.each do |d|
+          %w[libtatlas.so libatlas.so libsatlas.so].each do |f|
+            f = File.join(d,f) if !d.empty?
+            begin
+              Fiddle.dlopen(f)
+              Blas.dlopen(f)
+              l = LAPACK_LIBPATH
+              e = "liblapacke.so"
+              e = File.join(l,e) if l && !l.empty?
+              Lapack.dlopen(e)
+              @@libs = [ f, e ]
+              if $DEBUG
+                $stderr.puts "Numo::Linalg: use #{f}"
+              end
+              return true
+            rescue
+            end
+          end
+        end
+        if exc
+          raise RuntimeError, "cannot find ATLAS library"
+        end
+        false
+      end
+
+      def load_lapack(path=nil,exc:true)
+        a = case path
+            when String; [[path,path]]
+            when Array ;  [path]
+            when NilClass
+              b = BLAS_LIBPATH || LAPACK_LIBPATH
+              l = LAPACK_LIBPATH || BLAS_LIBPATH
+              (b ? [[b,l]] : []) + ["",""]
+            else raise ArgumentError,"invalid path"
+            end
+        a.each do |arg|
+          b,l = *arg
+          if b.empty?
+            f_blas = "libblas.so"
+            f_cblas = "libcblas.so"
+            f_lapack = "liblapack.so"
+            f_lapacke = "liblapacke.so"
+          else
+            f_blas = File.join(b,"libblas.so")
+            f_cblas = File.join(b,"libcblas.so")
+            f_lapack = File.join(l,"liblapack.so")
+            f_lapacke = File.join(l,"liblapacke.so")
+          end
+          begin
+            Fiddle.dlopen(f_blas)
+            Fiddle.dlopen(f_lapack)
+            Blas.dlopen(f_cblas)
+            Lapack.dlopen(f_lapacke)
+            @@libs = [ f_blas, f_lapack, f_cblas, f_lapacke ]
+            if $DEBUG
+              $stderr.puts "Numo::Linalg: use #{f_blas} and #{f_lapack}"
+            end
+            return true
+          rescue
+          end
+        end
+        if exc
+          raise RuntimeError, "cannot find BLAS/LAPABK library"
+        end
+        false
+      end
+
+      def load_library
+        case BACKEND
+        when /mkl/i         ; return if load_mkl(exc:false)
+        when /^openblas/i   ; return if load_openblas(exc:false)
+        when /^atlas/i      ; return if load_atlas(exc:false)
+        when /lapack|blas/i ; return if load_lapack(exc:false)
+        else
+          return if load_mkl(exc:false)
+          return if load_openblas(exc:false)
+          return if load_atlas(exc:false)
+          return if load_lapack(exc:false)
+        end
+        raise RuntimeError, "cannot find backend library for Numo::Linalg"
+      end
+
+    end
+  end
+end

data/lib/numo/linalg/use/atlas.rb

@@ -0,0 +1,3 @@
+require "numo/linalg/linalg"
+
+Numo::Linalg::Loader.load_atlas

data/lib/numo/linalg/use/lapack.rb

@@ -0,0 +1,3 @@
+require "numo/linalg/linalg"
+
+Numo::Linalg::Loader.load_lapack

data/lib/numo/linalg/use/mkl.rb

@@ -0,0 +1,3 @@
+require "numo/linalg/linalg"
+
+Numo::Linalg::Loader.load_mkl

data/lib/numo/linalg/use/openblas.rb

@@ -0,0 +1,3 @@
+require "numo/linalg/linalg"
+
+Numo::Linalg::Loader.load_openblas

data/lib/numo/linalg/version.rb

@@ -0,0 +1,5 @@
+module Numo
+  module Linalg
+    VERSION = "0.0.1"
+  end
+end

data/numo-linalg.gemspec

@@ -0,0 +1,26 @@
+lib = File.expand_path('../lib', __FILE__)
+$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
+
+require_relative "lib/numo/linalg/version"
+
+Gem::Specification.new do |spec|
+  spec.name          = "numo-linalg"
+  spec.version       = Numo::Linalg::VERSION
+  spec.authors       = ["Masahiro TANAKA"]
+  spec.email         = ["masa16.tanaka@gmail.com"]
+  spec.description   = %q{Linalg development version.}
+  spec.summary       = %q{Linalg development version}
+  spec.homepage      = "https://github.com/ruby-numo/numo-linalg"
+  spec.license       = "BSD-3-Clause"
+
+  spec.files         = `git ls-files Gemfile README.md Rakefile lib ext numo-linalg.gemspec spec`.split($/)
+  spec.executables   = spec.files.grep(%r{^bin/}) { |f| File.basename(f) }
+  spec.test_files    = spec.files.grep(%r{^(test|spec|features)/})
+  spec.require_paths = ["lib"]
+  spec.extensions    = ["ext/numo/linalg/blas/extconf.rb",
+                        "ext/numo/linalg/lapack/extconf.rb"]
+
+  spec.add_development_dependency "bundler", "~> 1.3"
+  spec.add_development_dependency "rake", "~> 0"
+  spec.add_runtime_dependency "numo-narray", "~> 0.9.0.7"
+end

data/spec/lapack_spec.rb

@@ -0,0 +1,13 @@
+require 'numo/narray'
+require File.join(File.dirname(__FILE__), "../ext/numo/lapack/lapack")
+#Numo::NArray.debug = true
+
+RSpec.configure do |config|
+  config.filter_run :focus
+  config.run_all_when_everything_filtered = true
+end
+#context :focus=>true do ... end
+
+describe 'something' do
+  it{expect(Numo::LAPACK).to be_a(Module)}
+end

metadata

@@ -0,0 +1,172 @@
+--- !ruby/object:Gem::Specification
+name: numo-linalg
+version: !ruby/object:Gem::Version
+  version: 0.0.1
+platform: ruby
+authors:
+- Masahiro TANAKA
+autorequire: 
+bindir: bin
+cert_chain: []
+date: 2017-05-28 00:00:00.000000000 Z
+dependencies:
+- !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: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0'
+- !ruby/object:Gem::Dependency
+  name: numo-narray
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.0.7
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.9.0.7
+description: Linalg development version.
+email:
+- masa16.tanaka@gmail.com
+executables: []
+extensions:
+- ext/numo/linalg/blas/extconf.rb
+- ext/numo/linalg/lapack/extconf.rb
+extra_rdoc_files: []
+files:
+- Gemfile
+- README.md
+- Rakefile
+- ext/numo/linalg/blas/blas.c
+- ext/numo/linalg/blas/cblas.h
+- ext/numo/linalg/blas/cblas_t.h
+- ext/numo/linalg/blas/depend.erb
+- ext/numo/linalg/blas/extconf.rb
+- ext/numo/linalg/blas/gen/cogen.rb
+- ext/numo/linalg/blas/gen/decl.rb
+- ext/numo/linalg/blas/gen/desc.rb
+- ext/numo/linalg/blas/gen/erbpp2.rb
+- ext/numo/linalg/blas/gen/replace_cblas_h.rb
+- ext/numo/linalg/blas/gen/spec.rb
+- ext/numo/linalg/blas/numo_blas.h
+- ext/numo/linalg/blas/tmpl/axpy.c
+- ext/numo/linalg/blas/tmpl/copy.c
+- ext/numo/linalg/blas/tmpl/def_c.c
+- ext/numo/linalg/blas/tmpl/def_d.c
+- ext/numo/linalg/blas/tmpl/def_s.c
+- ext/numo/linalg/blas/tmpl/def_z.c
+- ext/numo/linalg/blas/tmpl/dot.c
+- ext/numo/linalg/blas/tmpl/ger.c
+- ext/numo/linalg/blas/tmpl/init_class.c
+- ext/numo/linalg/blas/tmpl/init_module.c
+- ext/numo/linalg/blas/tmpl/lib.c
+- ext/numo/linalg/blas/tmpl/mm.c
+- ext/numo/linalg/blas/tmpl/module.c
+- ext/numo/linalg/blas/tmpl/mv.c
+- ext/numo/linalg/blas/tmpl/nrm2.c
+- ext/numo/linalg/blas/tmpl/rot.c
+- ext/numo/linalg/blas/tmpl/rotm.c
+- ext/numo/linalg/blas/tmpl/scal.c
+- ext/numo/linalg/blas/tmpl/sdsdot.c
+- ext/numo/linalg/blas/tmpl/set_prefix.c
+- ext/numo/linalg/blas/tmpl/swap.c
+- ext/numo/linalg/blas/tmpl/syr.c
+- ext/numo/linalg/blas/tmpl/syr2.c
+- ext/numo/linalg/blas/tmpl/syr2k.c
+- ext/numo/linalg/blas/tmpl/syrk.c
+- ext/numo/linalg/lapack/depend.erb
+- ext/numo/linalg/lapack/extconf.rb
+- ext/numo/linalg/lapack/gen/cogen.rb
+- ext/numo/linalg/lapack/gen/desc.rb
+- ext/numo/linalg/lapack/gen/replace_lapacke_h.rb
+- ext/numo/linalg/lapack/gen/spec.rb
+- ext/numo/linalg/lapack/lapack.c
+- ext/numo/linalg/lapack/lapacke.h
+- ext/numo/linalg/lapack/lapacke_config.h
+- ext/numo/linalg/lapack/lapacke_mangling.h
+- ext/numo/linalg/lapack/lapacke_t.h
+- ext/numo/linalg/lapack/numo_lapack.h
+- ext/numo/linalg/lapack/tmpl/def_c.c
+- ext/numo/linalg/lapack/tmpl/def_d.c
+- ext/numo/linalg/lapack/tmpl/def_s.c
+- ext/numo/linalg/lapack/tmpl/def_z.c
+- ext/numo/linalg/lapack/tmpl/fact.c
+- ext/numo/linalg/lapack/tmpl/geev.c
+- ext/numo/linalg/lapack/tmpl/gels.c
+- ext/numo/linalg/lapack/tmpl/gesv.c
+- ext/numo/linalg/lapack/tmpl/gesvd.c
+- ext/numo/linalg/lapack/tmpl/ggev.c
+- ext/numo/linalg/lapack/tmpl/gqr.c
+- ext/numo/linalg/lapack/tmpl/init_class.c
+- ext/numo/linalg/lapack/tmpl/init_module.c
+- ext/numo/linalg/lapack/tmpl/lange.c
+- ext/numo/linalg/lapack/tmpl/lib.c
+- ext/numo/linalg/lapack/tmpl/module.c
+- ext/numo/linalg/lapack/tmpl/syev.c
+- ext/numo/linalg/lapack/tmpl/sygv.c
+- ext/numo/linalg/lapack/tmpl/trf.c
+- ext/numo/linalg/numo_linalg.h
+- lib/numo/linalg.rb
+- lib/numo/linalg/function.rb
+- lib/numo/linalg/linalg.rb
+- lib/numo/linalg/loader.rb
+- lib/numo/linalg/use/atlas.rb
+- lib/numo/linalg/use/lapack.rb
+- lib/numo/linalg/use/mkl.rb
+- lib/numo/linalg/use/openblas.rb
+- lib/numo/linalg/version.rb
+- numo-linalg.gemspec
+- spec/lapack_spec.rb
+homepage: https://github.com/ruby-numo/numo-linalg
+licenses:
+- BSD-3-Clause
+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: '0'
+requirements: []
+rubyforge_project: 
+rubygems_version: 2.6.11
+signing_key: 
+specification_version: 4
+summary: Linalg development version
+test_files:
+- spec/lapack_spec.rb