numo-linalg-randsvd 0.1.0 → 0.2.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: 89e0d0d104cdc34ed4b941a4a6dd3460e72a477ea73c9cf3028a5074bd5234b1
-  data.tar.gz: 621681cad42d247f99993a7db1a4c636f87b1902d2ea4d9d12dc349eb0013c16
+  metadata.gz: 71312ec59a5898236280d11b100f6b43c7e42c4851e69198dc105d71745aef67
+  data.tar.gz: d95416e7779f463e7017d119e701fd232280d2da98767b4c58940b0582174b9f
 SHA512:
-  metadata.gz: 90fb17e4b7865c757e826ca2d8712f58f8e372a97acedd3eb45a03aa1775f4e23e20597218424669a8952350dc4e10504e2459cfdbf86f3d0a02d82d6a2687b0
-  data.tar.gz: c7202a45e21ac88a52502b1b2af2b616466e452af922472db686b49e01ea1827114312c64a482568c2a9406986096edfc9acfac2a38e32f69ccaa1b77231b348
+  metadata.gz: 2c1b847a52683c315bbea34e4b685d03332719a1b9b198318c41d55990cb7502b67fcec1e33ac4d8ab13fbf34abf92721e0633e3fd06182603658bf39fd60b0e
+  data.tar.gz: be658a95443b24be18f82f2895ee49b07133c5ef12355d0d3fcffe49f5b5501af2bead44dde05701df51989e0cd01d41cb427921aec87cd53960b01ba23ac7c6

data/CHANGELOG.md

@@ -1,3 +1,13 @@
+## [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.
+
 ## [0.1.0] - 2022-10-09
 
-- Initial release
+- Initial release.

data/README.md

@@ -1,5 +1,10 @@
 # Numo::Linalg.randsvd
 
+[![Build Status](https://github.com/yoshoku/numo-linalg-randsvd/actions/workflows/main.yml/badge.svg?branch=main)](https://github.com/yoshoku/numo-linalg-randsvd/actions/workflows/main.yml)
+[![Gem Version](https://badge.fury.io/rb/numo-linalg-randsvd.svg)](https://badge.fury.io/rb/numo-linalg-randsvd)
+[![BSD 3-Clause License](https://img.shields.io/badge/License-BSD%203--Clause-orange.svg)](https://github.com/yoshoku/numo-linalg-randsvd/blob/main/LICENSE.txt)
+[![Documentation](https://img.shields.io/badge/api-reference-blue.svg)](https://yoshoku.github.io/numo-linalg-randsvd/doc/)
+
 Numo::Linalg.randsvd is a module function on Numo::Linalg for truncated singular value decomposition with randomized algorithm.
 This gem re-implements [RandSVD](https://github.com/yoshoku/randsvd) using [Numo::NArray](https://github.com/ruby-numo/numo-narray) and
 [Numo::Linalg](https://github.com/ruby-numo/numo-linalg) instead of [NMatrix](https://github.com/SciRuby/nmatrix).

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.0
+  version: 0.2.0
 platform: ruby
 authors:
 - yoshoku
 autorequire:
 bindir: exe
 cert_chain: []
-date: 2022-10-09 00:00:00.000000000 Z
+date: 2022-11-12 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.0
+  type: :runtime
+  prerelease: false
+  version_requirements: !ruby/object:Gem::Requirement
+    requirements:
+    - - "~>"
+      - !ruby/object:Gem::Version
+        version: 0.4.0
 description: |
   Numo::Linalg.randsvd is a module function on Numo::Linalg for
   truncated singular value decomposition with randomized algorithm.
@@ -59,7 +73,7 @@ metadata:
   homepage_uri: https://github.com/yoshoku/numo-linalg-randsvd
   source_code_uri: https://github.com/yoshoku/numo-linalg-randsvd
   changelog_uri: https://github.com/yoshoku/numo-linalg-randsvd/blob/main/CHANGELOG.md
-  documentation_uri: https://rubydoc.info/gems/numo-linalg-randsvd
+  documentation_uri: https://yoshoku.github.io/numo-linalg-randsvd/doc/
   rubygems_mfa_required: 'true'
 post_install_message:
 rdoc_options: []
@@ -76,7 +90,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.2.33
+rubygems_version: 3.3.7
 signing_key:
 specification_version: 4
 summary: Numo::Linalg.randsvd is a module function on Numo::Linalg for truncated singular