numo-linalg-randsvd 0.1.1 → 0.2.1

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 13231de1b8ebda4466ffc207054984fb294ee73653dcd7c954089e2739d75ea0
-  data.tar.gz: f09bf78c3de2598053be56ce84083fdd293f3d4609c277ac6642c18a0ff3843f
+  metadata.gz: 6124cfa5d92e207e25176f5a76c10edf7efe094c7467b38bedaef3221ebf1125
+  data.tar.gz: c38604617b8dc031414205bd2359ecf54ceeac3d496cdd6c28b1e35b049baf88
 SHA512:
-  metadata.gz: 5349bfc846c6ac5a52139747db1c8261a2d0d3ba2ac1b2e25b958deaf8cc3283f2786401ba9def1aa8d2210fb39483bc5d1d0d2a90ef6cfd2b1f3e36d58efb87
-  data.tar.gz: 5add7d6b98a45739ab2c35f2c57e6d92bdcb210fa43327c913b59a3761816a0fca12090adc9fbea795fbad151fb58434c2a80a7b2674ac9ba4c99959efd787d9
+  metadata.gz: '04126899d6e0fdf582778a435d027380eaaff1eb6826614220dce3533f8e055ed518e3128e7bf822fe52e766858ea45ab08badb57d3f6ae83d92911e55abd056'
+  data.tar.gz: be04847c850344d151508bcd5a6051da2bd373480fdc8a2c91fd9e4403b41e8bc2ede31932aefa29f8271c7c3b1f4786f4a164593d04e773d63e041d1fce2aa9

data/CHANGELOG.md

@@ -1,3 +1,13 @@
+## [0.2.1] - 2022-11-2022
+
+- Fix dependent version of numo-random.
+
+## [0.2.0] - 2022-11-12
+
+**Breaking Changes**
+
+- Change to use numo-random gem for generating random matrix.
+
 ## [0.1.1] - 2022-10-23
 
 - Fix documentation uri.

data/lib/numo/linalg/randsvd.rb

@@ -1,6 +1,7 @@
 # frozen_string_literal: true
 
 require 'numo/narray'
+require 'numo/random'
 
 # Ruby/Numo (NUmerical MOdules)
 module Numo
@@ -15,24 +16,25 @@ module Numo
     # @param t [Integer] The number of iterations for orthogonalization.
     # @param driver [String] The driver parameter of Numo::Linalg.svd.
     # @param job [String] The job parameter of Numo::Linalg.svd.
-    def randsvd(a, k, t = 0, driver: 'svd', job: 'A')
+    # @param seed [Integer] The seed of random number generator.
+    def randsvd(a, k, t = 0, driver: 'svd', job: 'A', seed: nil)
       n = a.shape[1]
-      q = _orthonormal_mat(a, [k + 10, n].min, t)
+      q = _orthonormal_mat(a, [k + 10, n].min, t, seed)
       b = a.dot(q)
       s, u, vt = Numo::Linalg.svd(b, driver: driver, job: job)
       vtqt = vt.dot(q.transpose)
       _truncated_mat(s, u, vtqt, k)
     end
 
-    def _rand_normal(shape, dtype = Numo::DFloat, mu = 0.0, sigma = 1.0)
-      a = dtype.new(shape).rand
-      b = dtype.new(shape).rand
-      ((Numo::NMath.sqrt(Numo::NMath.log(a) * -2.0) * Numo::NMath.sin(b * 2.0 * Math::PI)) * sigma) + mu
+    def _rand_normal(shape, dtype = Numo::DFloat, seed = nil, mu = 0.0, sigma = 1.0)
+      Numo::Random::Generator.new(seed: seed).normal(
+        shape: shape, dtype: dtype.name.split('::')[-1].downcase.to_sym, loc: mu, scale: sigma
+      )
     end
 
-    def _orthonormal_mat(a, l, t)
+    def _orthonormal_mat(a, l, t, seed)
       m = a.shape[0]
-      r = _rand_normal([m, l], a.class)
+      r = _rand_normal([m, l], a.class, seed)
       q, = Numo::Linalg.qr(a.transpose.dot(r), mode: 'economic')
       t.times do
         r, = Numo::Linalg.qr(a.dot(q), mode: 'economic')

metadata

@@ -1,14 +1,14 @@
 --- !ruby/object:Gem::Specification
 name: numo-linalg-randsvd
 version: !ruby/object:Gem::Version
-  version: 0.1.1
+  version: 0.2.1
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2022-10-23 00:00:00.000000000 Z
+date: 2022-11-20 00:00:00.000000000 Z
 dependencies:
 - !ruby/object:Gem::Dependency
   name: numo-linalg
@@ -38,6 +38,20 @@ dependencies:
     - - ">="
       - !ruby/object:Gem::Version
         version: 0.9.1
+- !ruby/object:Gem::Dependency
+  name: numo-random
+  requirement: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: '0.4'
 description: |
   Numo::Linalg.randsvd is a module function on Numo::Linalg for
   truncated singular value decomposition with randomized algorithm.