numo-linalg-alt 0.5.0 → 0.7.0

data/ext/numo/linalg/lapack/sytrf.h

@@ -8,7 +8,7 @@
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "../util.h"
+#include "lapack_util.h"
 
 void define_linalg_lapack_sytrf(VALUE mLapack);
 

data/ext/numo/linalg/lapack/trtrs.h

@@ -8,7 +8,7 @@
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "../util.h"
+#include "lapack_util.h"
 
 void define_linalg_lapack_trtrs(VALUE mLapack);
 

data/ext/numo/linalg/lapack/unghr.h

@@ -8,7 +8,7 @@
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "../util.h"
+#include "lapack_util.h"
 
 void define_linalg_lapack_unghr(VALUE mLapack);
 

data/ext/numo/linalg/lapack/ungqr.h

@@ -8,7 +8,7 @@
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "../util.h"
+#include "lapack_util.h"
 
 void define_linalg_lapack_ungqr(VALUE mLapack);
 

data/ext/numo/linalg/lapack/ungrq.h

@@ -8,7 +8,7 @@
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "../util.h"
+#include "lapack_util.h"
 
 void define_linalg_lapack_ungrq(VALUE mLapack);
 

data/ext/numo/linalg/linalg.c

@@ -224,8 +224,10 @@ void Init_linalg(void) {
    */
   rb_mLinalgLapack = rb_define_module_under(rb_mLinalg, "Lapack");
 
+#ifdef HAVE_OPENBLAS_CONFIG_H
   /* The version of OpenBLAS used in background library. */
   rb_define_const(rb_mLinalg, "OPENBLAS_VERSION", rb_str_new_cstr(OPENBLAS_VERSION));
+#endif
 
   /**
    * Returns BLAS char ([sdcz]) defined by data-type of arguments.

data/ext/numo/linalg/linalg.h

@@ -35,15 +35,23 @@
 
 #include <ruby.h>
 
-#include <cblas.h>
-#include <lapacke.h>
-#include <openblas_config.h>
-
 #include <numo/narray.h>
 #include <numo/template.h>
 
-#include "converter.h"
-#include "util.h"
+#ifndef _DEFINED_SCOMPLEX
+#define _DEFINED_SCOMPLEX 1
+#endif
+#ifndef _DEFINED_DCOMPLEX
+#define _DEFINED_DCOMPLEX 1
+#endif
+
+#include <cblas.h>
+#include <lapacke.h>
+
+#include "extconf.h"
+#ifdef HAVE_OPENBLAS_CONFIG_H
+#include <openblas_config.h>
+#endif
 
 #include "blas/dot.h"
 #include "blas/dot_sub.h"

data/lib/numo/linalg/version.rb

@@ -5,6 +5,6 @@ module Numo
   # Numo::Linalg Alternative (numo-linalg-alt) is an alternative to Numo::Linalg.
   module Linalg
     # The version of numo-linalg-alt you install.
-    VERSION = '0.5.0'
+    VERSION = '0.7.0'
   end
 end

data/lib/numo/linalg.rb

@@ -8,6 +8,9 @@ require_relative 'linalg/linalg'
 module Numo
   # Numo::Linalg Alternative (numo-linalg-alt) is an alternative to Numo::Linalg.
   module Linalg # rubocop:disable Metrics/ModuleLength
+    # Exception class for errors occurred in LAPACK functions.
+    class LapackError < StandardError; end
+
     module_function
 
     # Computes the eigenvalues and eigenvectors of a symmetric / Hermitian matrix
@@ -340,7 +343,9 @@ module Numo
       raise ArgumentError, "invalid array type: #{a.class}" if bchr == 'n'
 
       fnc = :"#{bchr}potrs"
-      x, _info = Numo::Linalg::Lapack.send(fnc, a, b.dup, uplo: uplo)
+      x, info = Numo::Linalg::Lapack.send(fnc, a, b.dup, uplo: uplo)
+      raise LapackError, "the #{-info}-th argument of potrs had illegal value" if info.negative?
+
       x
     end
 
@@ -364,17 +369,17 @@ module Numo
 
       getrf = :"#{bchr}getrf"
       lu, piv, info = Numo::Linalg::Lapack.send(getrf, a.dup)
+      raise LapackError, "the #{-info}-th argument of getrf had illegal value" if info.negative?
 
-      if info.zero?
-        det_l = 1
-        det_u = lu.diagonal.prod
-        det_p = piv.map_with_index { |v, i| v == i + 1 ? 1 : -1 }.prod
-        det_l * det_u * det_p
-      elsif info.positive?
-        raise 'the factor U is singular, and the inverse matrix could not be computed.'
-      else
-        raise "the #{-info}-th argument of getrf had illegal value"
-      end
+      # info > 0 means the factor U has a zero diagonal element and is singular.
+      # In this case, the determinant is zero. The method should simply return 0.0.
+      # Therefore, the error is not raised here.
+      # raise 'the factor U is singular, ...' if info.positive?
+
+      det_l = 1
+      det_u = lu.diagonal.prod
+      det_p = piv.map_with_index { |v, i| v == i + 1 ? 1 : -1 }.prod
+      det_l * det_u * det_p
     end
 
     # Computes the inverse matrix of a square matrix.
@@ -407,13 +412,13 @@ module Numo
       getri = :"#{bchr}getri"
 
       lu, piv, info = Numo::Linalg::Lapack.send(getrf, a.dup)
-      if info.zero?
-        Numo::Linalg::Lapack.send(getri, lu, piv)[0]
-      elsif info.positive?
-        raise 'the factor U is singular, and the inverse matrix could not be computed.'
-      else
-        raise "the #{-info}-th argument of getrf had illegal value"
-      end
+      raise LapackError, "the #{-info}-th argument of getrf had illegal value" if info.negative?
+
+      a_inv, info = Numo::Linalg::Lapack.send(getri, lu, piv)
+      raise LapackError, "the #{-info}-th argument of getrf had illegal value" if info.negative?
+      raise LapackError, 'The matrix is singular, and the inverse matrix could not be computed.' if info.positive?
+
+      a_inv
     end
 
     # Computes the (Moore-Penrose) pseudo-inverse of a matrix using singular value decomposition.
@@ -441,7 +446,7 @@ module Numo
       rank = s.gt(rcond * s[0]).count
 
       u = u[true, 0...rank] / s[0...rank]
-      u.dot(vh[0...rank, true]).conj.transpose
+      u.dot(vh[0...rank, true]).conj.transpose.dup
     end
 
     # Computes the polar decomposition of a matrix.
@@ -612,7 +617,7 @@ module Numo
 
       fnc = :"#{bchr}gerqf"
       rq, tau, info = Numo::Linalg::Lapack.send(fnc, a.dup)
-      raise "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
 
       m, n = rq.shape
       r = rq.triu(n - m).dup
@@ -630,7 +635,7 @@ module Numo
             end
 
       q, info = Numo::Linalg::Lapack.send(fnc, tmp, tau)
-      raise "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
 
       [r, q]
     end
@@ -675,8 +680,21 @@ module Numo
         aa, bb, _alpha, _beta, q, z, _sdim, info = Numo::Linalg::Lapack.send(fnc, a.dup, b.dup)
       end
 
-      raise "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the QZ algorithm failed.' if info.positive?
+      n = a.shape[0]
+      raise LapackError, "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'something other than QZ iteration failed.' if info == n + 1
+      raise LapackError, "reordering failed in #{bchr}tgsen" if info == n + 3
+
+      if info == n + 2
+        raise LapackError, 'after reordering, roundoff changed values of some eigenvalues ' \
+                           'so that leading eigenvalues in the Generalized Schur form no ' \
+                           'longer satisfy the sorting condition.'
+      end
+
+      if info.positive? && info <= n
+        warn('the QZ iteration failed. (a, b) are not in Schur form, ' \
+             "but alpha[i] and beta[i] for i = #{info},...,n are correct.")
+      end
 
       [aa, bb, q, z]
     end
@@ -732,10 +750,15 @@ module Numo
       end
 
       n = a.shape[0]
-      raise "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the QR algorithm failed to compute all the eigenvalues.' if info.positive? && info <= n
-      raise 'the eigenvalues could not be reordered.' if info == n + 1
-      raise 'after reordering, roundoff changed values of some complex eigenvalues.' if info == n + 2
+      raise LapackError, "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'the QR algorithm failed to compute all the eigenvalues.' if info.positive? && info <= n
+      raise LapackError, 'the eigenvalues could not be reordered.' if info == n + 1
+
+      if info == n + 2
+        raise LapackError, 'after reordering, roundoff changed values of some eigenvalues ' \
+                           'so that leading eigenvalues in the Schur form no longer satisfy ' \
+                           'the sorting condition.'
+      end
 
       [b, v, sdim]
     end
@@ -779,12 +802,12 @@ module Numo
       func = :"#{bchr}gebal"
       b, ilo, ihi, _, info = Numo::Linalg::Lapack.send(func, a.dup)
 
-      raise "the #{-info}-th argument of #{func} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{func} had illegal value" if info.negative?
 
       func = :"#{bchr}gehrd"
       hq, tau, info = Numo::Linalg::Lapack.send(func, b, ilo: ilo, ihi: ihi)
 
-      raise "the #{-info}-th argument of #{func} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{func} had illegal value" if info.negative?
 
       h = hq.triu(-1)
       return h unless calc_q
@@ -792,7 +815,7 @@ module Numo
       func = %w[d s].include?(bchr) ? :"#{bchr}orghr" : :"#{bchr}unghr"
       q, info = Numo::Linalg::Lapack.send(func, hq, tau, ilo: ilo, ihi: ihi)
 
-      raise "the #{-info}-th argument of #{func} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{func} had illegal value" if info.negative?
 
       [h, q]
     end
@@ -830,7 +853,15 @@ module Numo
       raise ArgumentError, "invalid array type: #{a.class}, #{b.class}" if bchr == 'n'
 
       gesv = :"#{bchr}gesv"
-      Numo::Linalg::Lapack.send(gesv, a.dup, b.dup)[1]
+      _lu, x, _ipiv, info = Numo::Linalg::Lapack.send(gesv, a.dup, b.dup)
+      raise LapackError, "the #{-info}-th argument of getrf had illegal value" if info.negative?
+
+      if info.positive?
+        warn('the factorization has been completed, but the factor is singular, ' \
+             'so the solution could not be computed.')
+      end
+
+      x
     end
 
     # Solves linear equation `A * x = b` or `A * X = B` for `x` assuming `A` is a triangular matrix.
@@ -867,7 +898,7 @@ module Numo
       trtrs = :"#{bchr}trtrs"
       uplo = lower ? 'L' : 'U'
       x, info = Numo::Linalg::Lapack.send(trtrs, a, b.dup, uplo: uplo)
-      raise "wrong value is given to the #{info}-th argument of #{trtrs} used internally" if info.negative?
+      raise LapackError, "wrong value is given to the #{info}-th argument of #{trtrs} used internally" if info.negative?
 
       x
     end
@@ -923,9 +954,9 @@ module Numo
         raise ArgumentError, "invalid driver: #{driver}"
       end
 
-      raise "the #{info.abs}-th argument had illegal value" if info.negative?
-      raise 'input array has a NAN entry' if info == -4
-      raise 'svd did not converge' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument had illegal value" if info.negative?
+      raise LapackError, 'the input array has a NAN entry' if info == -4
+      raise LapackError, 'the did not converge' if info.positive?
 
       [s, u, vt]
     end
@@ -1012,9 +1043,9 @@ module Numo
         raise ArgumentError, "invalid driver: #{driver}"
       end
 
-      raise "the #{info.abs}-th argument had illegal value" if info.negative?
-      raise 'input array has a NAN entry' if info == -4
-      raise 'svd did not converge' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument had illegal value" if info.negative?
+      raise LapackError, 'the input array has a NAN entry' if info == -4
+      raise LapackError, 'the decomposition did not converge' if info.positive?
 
       s
     end
@@ -1140,8 +1171,12 @@ module Numo
       getrf = :"#{bchr}getrf"
       lu, piv, info = Numo::Linalg::Lapack.send(getrf, a.dup)
 
-      raise "the #{info.abs}-th argument of getrf had illegal value" if info.negative?
-      raise "the U(#{info}, #{info}) is exactly zero. The factorization has been completed." if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of getrf had illegal value" if info.negative?
+
+      if info.positive?
+        warn("the factorization has been completed, but the factor U[#{info - 1}, #{info - 1}] is " \
+             'exactly zero, indicating that the matrix is singular.')
+      end
 
       [lu, piv]
     end
@@ -1179,7 +1214,7 @@ module Numo
       getrs = :"#{bchr}getrs"
       x, info = Numo::Linalg::Lapack.send(getrs, lu, ipiv, b.dup)
 
-      raise "the #{info.abs}-th argument of getrs had illegal value" if info.negative?
+      raise LapackError, "the #{info.abs}-th argument of getrs had illegal value" if info.negative?
 
       x
     end
@@ -1200,8 +1235,12 @@ module Numo
       fnc = :"#{bchr}potrf"
       c, info = Numo::Linalg::Lapack.send(fnc, a.dup, uplo: uplo)
 
-      raise "the #{info}-th leading minor of the array is not positive definite, and the factorization could not be completed." if info.positive?
-      raise "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, "the #{-info}-th argument of #{fnc} had illegal value" if info.negative?
+
+      if info.positive?
+        raise LapackError, "the leading principal minor of order #{info} is not positive, " \
+                           'and the factorization could not be completed.'
+      end
 
       c
     end
@@ -1297,7 +1336,7 @@ module Numo
       fnc = :"#{bchr}gebal"
       b, lo, hi, prm_scl, info = Numo::Linalg::Lapack.send(fnc, a.dup, job: job)
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
 
       # convert from Fortran style index to Ruby style index.
       lo -= 1
@@ -1374,8 +1413,8 @@ module Numo
         wr, wi, vl, vr, info = Numo::Linalg::Lapack.send(fnc, a.dup, jobvl: jobvl, jobvr: jobvr)
       end
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the QR algorithm failed to compute all the eigenvalues.' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'the QR algorithm failed to compute all the eigenvalues.' if info.positive?
 
       if %w[d s].include?(bchr)
         w = wr + (wi * 1.0i)
@@ -1419,8 +1458,8 @@ module Numo
         w = wr + (wi * 1.0i)
       end
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the QR algorithm failed to compute all the eigenvalues.' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'the QR algorithm failed to compute all the eigenvalues.' if info.positive?
 
       w
     end
@@ -1470,8 +1509,12 @@ module Numo
       lud = a.dup
       ipiv, info = Numo::Linalg::Lapack.send(fnc, lud, uplo: uplo)
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the factorization has been completed' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+
+      if info.positive?
+        warn("the factorization has been completed, but the D[#{info - 1}, #{info - 1}] is " \
+             'exactly zero, indicating that the block diagonal matrix is singular.')
+      end
 
       _lud_permutation(lud, ipiv, uplo: uplo, hermitian: hermitian)
     end
@@ -1559,8 +1602,8 @@ module Numo
       fnc = :"#{bchr}gelsd"
       s, rank, info = Numo::Linalg::Lapack.send(fnc, a.dup, x, rcond: rcond)
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the algorithm for computing the SVD failed to converge' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'the algorithm for computing the SVD failed to converge' if info.positive?
 
       resids = x.class[]
       if m > n
@@ -1616,6 +1659,22 @@ module Numo
       a_expm
     end
 
+    # Computes the matrix logarithm using its eigenvalue decomposition.
+    #
+    # @param a [Numo::NArray] The n-by-n square matrix.
+    # @return [Numo::NArray] The matrix logarithm of `a`.
+    def logm(a)
+      raise Numo::NArray::ShapeError, 'input array a must be 2-dimensional' if a.ndim != 2
+      raise Numo::NArray::ShapeError, 'input array a must be square' if a.shape[0] != a.shape[1]
+
+      ev, vl, = eig(a, left: true, right: false)
+      v = vl.transpose.conj
+      inv_v = Numo::Linalg.inv(v)
+      log_ev = Numo::NMath.log(ev)
+
+      inv_v.dot(log_ev.diag).dot(v)
+    end
+
     # Computes the matrix sine using the matrix exponential.
     #
     # @param a [Numo::NArray] The n-by-n square matrix.
@@ -1705,6 +1764,35 @@ module Numo
       a_sinh.dot(Numo::Linalg.inv(a_cosh))
     end
 
+    # Computes the square root of a matrix using its eigenvalue decomposition.
+    #
+    # @param a [Numo::NArray] The n-by-n square matrix.
+    # @return [Numo::NArray] The matrix square root of `a`.
+    def sqrtm(a)
+      raise Numo::NArray::ShapeError, 'input array a must be 2-dimensional' if a.ndim != 2
+      raise Numo::NArray::ShapeError, 'input array a must be square' if a.shape[0] != a.shape[1]
+
+      ev, vl, = eig(a, left: true, right: false)
+      v = vl.transpose.conj
+      inv_v = Numo::Linalg.inv(v)
+      sqrt_ev = Numo::NMath.sqrt(ev)
+
+      inv_v.dot(sqrt_ev.diag).dot(v)
+    end
+
+    # Computes the matrix sign function using its inverse and square root matrices.
+    #
+    # @param a [Numo::NArray] The n-by-n square matrix.
+    # @return [Numo::NArray] The matrix sign function of `a`.
+    def signm(a)
+      raise Numo::NArray::ShapeError, 'input array a must be 2-dimensional' if a.ndim != 2
+      raise Numo::NArray::ShapeError, 'input array a must be square' if a.shape[0] != a.shape[1]
+
+      a_sqrt = sqrtm(a.dot(a))
+      a_inv = Numo::Linalg.inv(a)
+      a_inv.dot(a_sqrt)
+    end
+
     # Computes the inverse of a matrix using its LU decomposition.
     #
     # @param lu [Numo::NArray] The LU decomposition of the n-by-n matrix `A`.
@@ -1717,8 +1805,8 @@ module Numo
       fnc = :"#{bchr}getri"
       inv, info = Numo::Linalg::Lapack.send(fnc, lu.dup, ipiv)
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise 'the matrix is singular and its inverse could not be computed' if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+      raise LapackError, 'the matrix is singular and its inverse could not be computed' if info.positive?
 
       inv
     end
@@ -1748,8 +1836,12 @@ module Numo
       fnc = :"#{bchr}potri"
       inv, info = Numo::Linalg::Lapack.send(fnc, a.dup, uplo: uplo)
 
-      raise "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
-      raise "the (#{info}, #info)-th element of the factor U or L is zero, and the inverse could not be computed." if info.positive?
+      raise LapackError, "the #{info.abs}-th argument of #{fnc} had illegal value" if info.negative?
+
+      if info.positive?
+        raise LapackError, "the (#{info - 1}, #{info - 1})-th element of the factor U or L is zero, " \
+                           'and the inverse could not be computed.'
+      end
 
       inv
     end

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: numo-linalg-alt
 version: !ruby/object:Gem::Version
-  version: 0.5.0
+  version: 0.7.0
 platform: ruby
 authors:
 - yoshoku
@@ -37,6 +37,11 @@ files:
 - CODE_OF_CONDUCT.md
 - LICENSE.txt
 - README.md
+- ext/numo/linalg/blas/blas_common.h
+- ext/numo/linalg/blas/blas_util.c
+- ext/numo/linalg/blas/blas_util.h
+- ext/numo/linalg/blas/converter.c
+- ext/numo/linalg/blas/converter.h
 - ext/numo/linalg/blas/dot.c
 - ext/numo/linalg/blas/dot.h
 - ext/numo/linalg/blas/dot_sub.c
@@ -47,8 +52,6 @@ files:
 - ext/numo/linalg/blas/gemv.h
 - ext/numo/linalg/blas/nrm2.c
 - ext/numo/linalg/blas/nrm2.h
-- ext/numo/linalg/converter.c
-- ext/numo/linalg/converter.h
 - ext/numo/linalg/extconf.rb
 - ext/numo/linalg/lapack/gebal.c
 - ext/numo/linalg/lapack/gebal.h
@@ -94,6 +97,8 @@ files:
 - ext/numo/linalg/lapack/hetrf.h
 - ext/numo/linalg/lapack/lange.c
 - ext/numo/linalg/lapack/lange.h
+- ext/numo/linalg/lapack/lapack_util.c
+- ext/numo/linalg/lapack/lapack_util.h
 - ext/numo/linalg/lapack/orghr.c
 - ext/numo/linalg/lapack/orghr.h
 - ext/numo/linalg/lapack/orgqr.c
@@ -130,8 +135,6 @@ files:
 - ext/numo/linalg/lapack/ungrq.h
 - ext/numo/linalg/linalg.c
 - ext/numo/linalg/linalg.h
-- ext/numo/linalg/util.c
-- ext/numo/linalg/util.h
 - lib/numo/linalg.rb
 - lib/numo/linalg/version.rb
 - vendor/tmp/.gitkeep
@@ -142,7 +145,7 @@ metadata:
   homepage_uri: https://github.com/yoshoku/numo-linalg-alt
   source_code_uri: https://github.com/yoshoku/numo-linalg-alt
   changelog_uri: https://github.com/yoshoku/numo-linalg-alt/blob/main/CHANGELOG.md
-  documentation_uri: https://gemdocs.org/gems/numo-linalg-alt/0.5.0/
+  documentation_uri: https://gemdocs.org/gems/numo-linalg-alt/0.7.0/
   rubygems_mfa_required: 'true'
 rdoc_options: []
 require_paths:
@@ -158,7 +161,7 @@ required_rubygems_version: !ruby/object:Gem::Requirement
     - !ruby/object:Gem::Version
       version: '0'
 requirements: []
-rubygems_version: 3.6.9
+rubygems_version: 3.7.2
 specification_version: 4
 summary: Numo::Linalg Alternative (numo-linalg-alt) is an alternative to Numo::Linalg.
 test_files: []

data/ext/numo/linalg/util.c

@@ -1,103 +0,0 @@
-#include "util.h"
-
-lapack_int get_itype(VALUE val) {
-  const lapack_int itype = NUM2INT(val);
-
-  if (itype != 1 && itype != 2 && itype != 3) {
-    rb_raise(rb_eArgError, "itype must be 1, 2 or 3");
-  }
-
-  return itype;
-}
-
-char get_jobz(VALUE val) {
-  const char jobz = NUM2CHR(val);
-
-  if (jobz != 'N' && jobz != 'V') {
-    rb_raise(rb_eArgError, "jobz must be 'N' or 'V'");
-  }
-
-  return jobz;
-}
-
-char get_jobvs(VALUE val) {
-  const char jobvs = NUM2CHR(val);
-  if (jobvs != 'N' && jobvs != 'V') {
-    rb_raise(rb_eArgError, "jobvs must be 'N' or 'V'");
-  }
-  return jobvs;
-}
-
-char get_range(VALUE val) {
-  const char range = NUM2CHR(val);
-
-  if (range != 'A' && range != 'V' && range != 'I') {
-    rb_raise(rb_eArgError, "range must be 'A', 'V' or 'I'");
-  }
-
-  return range;
-}
-
-char get_uplo(VALUE val) {
-  const char uplo = NUM2CHR(val);
-
-  if (uplo != 'U' && uplo != 'L') {
-    rb_raise(rb_eArgError, "uplo must be 'U' or 'L'");
-  }
-
-  return uplo;
-}
-
-int get_matrix_layout(VALUE val) {
-  const char option = NUM2CHR(val);
-
-  switch (option) {
-  case 'r':
-  case 'R':
-    break;
-  case 'c':
-  case 'C':
-    rb_warn("Numo::Linalg does not support column major.");
-    break;
-  }
-
-  return LAPACK_ROW_MAJOR;
-}
-
-enum CBLAS_TRANSPOSE get_cblas_trans(VALUE val) {
-  const char option = NUM2CHR(val);
-  enum CBLAS_TRANSPOSE res = CblasNoTrans;
-
-  switch (option) {
-  case 'n':
-  case 'N':
-    res = CblasNoTrans;
-    break;
-  case 't':
-  case 'T':
-    res = CblasTrans;
-    break;
-  case 'c':
-  case 'C':
-    res = CblasConjTrans;
-    break;
-  }
-
-  return res;
-}
-
-enum CBLAS_ORDER get_cblas_order(VALUE val) {
-  const char option = NUM2CHR(val);
-
-  switch (option) {
-  case 'r':
-  case 'R':
-    break;
-  case 'c':
-  case 'C':
-    rb_warn("Numo::Linalg does not support column major.");
-    break;
-  }
-
-  return CblasRowMajor;
-}

data/ext/numo/linalg/util.h

@@ -1,18 +0,0 @@
-#ifndef NUMO_LINALG_ALT_UTIL_H
-#define NUMO_LINALG_ALT_UTIL_H 1
-
-#include <ruby.h>
-
-#include <cblas.h>
-#include <lapacke.h>
-
-lapack_int get_itype(VALUE val);
-char get_jobz(VALUE val);
-char get_jobvs(VALUE val);
-char get_range(VALUE val);
-char get_uplo(VALUE val);
-int get_matrix_layout(VALUE val);
-enum CBLAS_TRANSPOSE get_cblas_trans(VALUE val);
-enum CBLAS_ORDER get_cblas_order(VALUE val);
-
-#endif // NUMO_LINALG_ALT_UTIL_H