kmat 0.0.3

data/lib/kmat/random.rb

@@ -0,0 +1,106 @@
+class Random
+	private def _kmat_rand_d(arg=nil)
+		if arg
+			self.rand(arg).to_f
+		else
+			self.rand()
+		end
+	end
+	private def _kmat_rand_z(arg=nil)
+		if arg
+			Complex(self.rand(arg), self.rand(arg))
+		else
+			Complex(self.rand(), self.rand())
+		end
+	end
+	private def _kmat_rand_i(arg=nil)
+		case arg
+		when NilClass
+			self.rand(2)
+		when Integer, Range
+			self.rand(arg)
+		when Float
+			if 0.0 <= arg && arg <= 1.0
+				( self.rand() < arg ) ? 1 : 0
+			else
+				raise ArgumentError, 'Float argument must be in [0, 1] for vt:int'
+			end
+		else
+			raise ArgumentError, 'for vt:int, arg must be Integer(Uniform(0...arg)), Range(Uniform(arg)), Float(1 with probability arg, 0 with probability 1-arg) or omitted'
+		end
+	end
+	private def _kmat_rand_b(arg=nil)
+		if arg.nil?
+			if self.rand(2) == 0
+				false
+			else
+				true
+			end
+		elsif arg.kind_of?(Numeric) && 0.0 <= arg && arg <= 1.0
+			( self.rand() < arg )
+		else
+			raise ArgumentError, 'for vt:bool, arg must be Float(true with probability arg, false with probability 1-arg) or omitted'
+		end
+	end
+	private def _kmat_rand_v(arg=nil)
+		if arg
+			self.rand(arg)
+		else
+			self.rand()
+		end
+	end
+	def kmat_vt_rand(vt, arg=nil)
+		case Mat.vt_refine(vt)
+		when :float
+			_kmat_rand_d(arg)
+		when :complex
+			_kmat_rand_z(arg)
+		when :int
+			_kmat_rand_i(arg)
+		when :bool
+			_kmat_rand_b(arg)
+		when :object
+			_kmat_rand_v(arg)
+		else
+			raise Mat::InternalError, 'unknown value type'
+		end
+	end
+end
+
+class << Random
+	using(Module.new do
+		refine Random.singleton_class do
+			alias :kmat_srand_org :srand
+		end
+	end)
+	def srand(*args)
+		@kmat_stored = false
+		kmat_srand_org(*args)
+	end
+	def randn(*args)
+		Random::DEFAULT.randn(*args)
+	end
+end
+
+module Kernel
+	module_function def randn(*args)
+		Random::DEFAULT.randn(*args)
+	end
+end
+
+class Mat
+	def rand(arg=nil, random: $MatRandom)
+		_rand0(random, arg)
+	end
+	def randn(random: $MatRandom)
+		_randn0(random)
+	end
+end
+class << Mat
+	def rand(m=1, n=1, random: $MatRandom)
+		_rand0(m, n, random)
+	end
+	def randn(m=1, n=1, random: $MatRandom)
+		_randn0(m, n, random)
+	end
+end

data/lib/kmat/statistics.rb

@@ -0,0 +1,98 @@
+class Mat
+	def sum(arg=:all)
+		case arg
+		when Mat
+			_sum(arg)
+		when :all, :a
+			_sum(Mat.new(1, 1, self.vtype))
+		when :col, :c
+			_sum(Mat.new(1, col_size, self.vtype))
+		when :row, :r
+			_sum(Mat.new(row_size, 1, self.vtype))
+		when :none, :n
+			self.dup
+		end
+	end
+	
+	def mean(arg=:all)
+		foo = sum(arg)
+		if foo.kind_of?(Mat)
+			if foo.row_size != self.row_size
+				foo.e_div(self.row_size)
+			else
+				foo.e_div(self.col_size)
+			end
+		elsif arg.kind_of?(Mat)
+			arg.e_div(self.length)
+			arg[0, 0]
+		else
+			foo.quo(self.length)
+		end
+	end
+	
+	private def _mean(arg, ddof)
+		foo = sum(arg)
+		if foo.kind_of?(Mat)
+			if foo.row_size != self.row_size
+				foo.e_div(self.row_size-ddof)
+			else
+				foo.e_div(self.col_size-ddof)
+			end
+		elsif arg.kind_of?(Mat)
+			arg.e_div(self.length-ddof)
+			arg[0, 0]
+		else
+			foo.quo(self.length)
+		end
+	end
+	def var(arg=:all, ddof: 0)
+		foo = mean(arg)
+		if foo.kind_of?(Mat)
+			bar = foo.repmat(self.row_size.div(foo.col_size), self.col_size.div(foo.col_size))
+			bar.sub!(self).e_mul!(bar)
+		else
+			unless arg.kind_of?(Mat)
+				arg = Mat.new(1, 1, self.vtype) do
+					foo
+				end
+			end
+			bar = Mat.new(self.row_size, self.col_size, self.vtype)
+			bar.fill(foo)
+		end
+		bar.sub!(self)
+		bar.e_mul!(bar)
+		bar.__send__(:_mean, arg, ddof)
+		if foo.kind_of?(Mat)
+			bar
+		else
+			bar[0, 0]
+		end
+	end
+	alias :variance :var
+	
+	def std(arg=:all, ddof: 0)
+		foo = var(arg, ddof: ddof)
+		if foo.kind_of?(Mat)
+			foo.sqrt!
+		elsif arg.kind_of?(Mat)
+			arg.sqrt!
+			arg[0, 0]
+		else
+			Math.sqrt(arg)
+		end
+	end
+	alias :stddev :std
+	alias :standard_deviation :std
+	
+	def geo_mean(arg=:all)
+		foo = self.log.mean(arg)
+		if foo.kind_of?(Mat)
+			foo.exp!
+		elsif arg.kind_of?(Mat)
+			arg.exp!
+			arg[0, 0]
+		else
+			Math.exp(foo)
+		end
+	end
+end

data/lib/kmat/version.rb

@@ -0,0 +1,3 @@
+class Mat
+  VERSION = "0.0.3"
+end

metadata

@@ -0,0 +1,156 @@
+--- !ruby/object:Gem::Specification
+name: kmat
+version: !ruby/object:Gem::Version
+  version: 0.0.3
+platform: ruby
+authors:
+- KAZOON
+autorequire: 
+bindir: exe
+cert_chain: []
+date: 2019-12-03 00:00:00.000000000 Z
+dependencies:
+- !ruby/object:Gem::Dependency
+  name: bundler
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '2.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '2.0'
+- !ruby/object:Gem::Dependency
+  name: rake
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '10.0'
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: '10.0'
+- !ruby/object:Gem::Dependency
+  name: rake-compiler
+  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: pry
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.12.2
+  type: :development
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - ">="
+      - !ruby/object:Gem::Version
+        version: 0.12.2
+description: 
+email:
+- cycloawaodorin+gem@gmail.com
+executables: []
+extensions:
+- ext/kmat/extconf.rb
+extra_rdoc_files: []
+files:
+- ".gitattributes"
+- ".gitignore"
+- CHANGELOG.md
+- Gemfile
+- LICENSE.md
+- README.md
+- Rakefile
+- bin/console
+- bin/setup
+- ext/kmat/arith/binary.c
+- ext/kmat/arith/logical.c
+- ext/kmat/arith/math.c
+- ext/kmat/arith/statistics.c
+- ext/kmat/arith/unary.c
+- ext/kmat/auto_collect.rb
+- ext/kmat/elementwise_function.rb
+- ext/kmat/extconf.rb
+- ext/kmat/id.txt
+- ext/kmat/id_sym.rb
+- ext/kmat/km_util.h
+- ext/kmat/kmat.h
+- ext/kmat/lapack_headers/blas.h
+- ext/kmat/lapack_headers/lapacke.h
+- ext/kmat/lapack_headers/lapacke_config.h
+- ext/kmat/lapack_headers/lapacke_mangling.h
+- ext/kmat/lapack_headers/lapacke_utils.h
+- ext/kmat/linalg/dla.c
+- ext/kmat/linalg/linalg.c
+- ext/kmat/linalg/norm.c
+- ext/kmat/linalg/vla.c
+- ext/kmat/linalg/working.c
+- ext/kmat/main.c
+- ext/kmat/smat/accessor.c
+- ext/kmat/smat/array.c
+- ext/kmat/smat/boxmuller.c
+- ext/kmat/smat/constructer.c
+- ext/kmat/smat/convert.c
+- ext/kmat/smat/elem.c
+- ext/kmat/smat/fund.c
+- ext/kmat/smat/share.c
+- ext/kmat/smat/smat.c
+- ext/kmat/smat/sort.c
+- ext/kmat/sym.txt
+- kmat.gemspec
+- lib/kmat.rb
+- lib/kmat/accessor.rb
+- lib/kmat/arith.rb
+- lib/kmat/linalg.rb
+- lib/kmat/logical.rb
+- lib/kmat/misc.rb
+- lib/kmat/random.rb
+- lib/kmat/statistics.rb
+- lib/kmat/version.rb
+homepage: https://github.com/cycloawaodorin/
+licenses:
+- GPL-3.0
+metadata:
+  homepage_uri: https://github.com/cycloawaodorin/
+  source_code_uri: https://github.com/cycloawaodorin/kmat
+  changelog_uri: https://github.com/cycloawaodorin/kmat/blob/master/CHANGELOG.md
+post_install_message: 
+rdoc_options: []
+require_paths:
+- lib
+required_ruby_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: 2.6.0
+required_rubygems_version: !ruby/object:Gem::Requirement
+  requirements:
+  - - ">="
+    - !ruby/object:Gem::Version
+      version: '0'
+requirements: []
+rubygems_version: 3.0.6
+signing_key: 
+specification_version: 4
+summary: Kmat is a Ruby gem for matrix operations. Kmat uses BLAS/LAPACK as back-end.
+test_files: []