numo-tiny_linalg 0.0.4 → 0.1.0

data/ext/numo/tiny_linalg/lapack/hegv.hpp

@@ -24,7 +24,7 @@ struct CHeGv {
   }
 };
 
-template <int nary_dtype_id, int nary_rtype_id, typename DType, typename RType, typename FncType>
+template <int nary_dtype_id, int nary_rtype_id, typename dtype, typename rtype, class LapackFn>
 class HeGv {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -40,15 +40,15 @@ private:
   };
 
   static void iter_hegv(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* b = (DType*)NDL_PTR(lp, 1);
-    RType* w = (RType*)NDL_PTR(lp, 2);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* b = (dtype*)NDL_PTR(lp, 1);
+    rtype* w = (rtype*)NDL_PTR(lp, 2);
     int* info = (int*)NDL_PTR(lp, 3);
     hegv_opt* opt = (hegv_opt*)(lp->opt_ptr);
     const lapack_int n = NDL_SHAPE(lp, 0)[1];
     const lapack_int lda = NDL_SHAPE(lp, 0)[0];
     const lapack_int ldb = NDL_SHAPE(lp, 1)[0];
-    const lapack_int i = FncType().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
+    const lapack_int i = LapackFn().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
     *info = static_cast<int>(i);
   }
 
@@ -63,10 +63,10 @@ private:
     ID kw_table[4] = { rb_intern("itype"), rb_intern("jobz"), rb_intern("uplo"), rb_intern("order") };
     VALUE kw_values[4] = { Qundef, Qundef, Qundef, Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 4, kw_values);
-    const lapack_int itype = kw_values[0] != Qundef ? get_itype(kw_values[0]) : 1;
-    const char jobz = kw_values[1] != Qundef ? get_jobz(kw_values[1]) : 'V';
-    const char uplo = kw_values[2] != Qundef ? get_uplo(kw_values[2]) : 'U';
-    const int matrix_layout = kw_values[3] != Qundef ? get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
+    const lapack_int itype = kw_values[0] != Qundef ? Util().get_itype(kw_values[0]) : 1;
+    const char jobz = kw_values[1] != Qundef ? Util().get_jobz(kw_values[1]) : 'V';
+    const char uplo = kw_values[2] != Qundef ? Util().get_uplo(kw_values[2]) : 'U';
+    const int matrix_layout = kw_values[3] != Qundef ? Util().get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -116,52 +116,6 @@ private:
 
     return ret;
   }
-
-  static 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;
-  }
-
-  static char get_jobz(VALUE val) {
-    const char jobz = NUM2CHR(val);
-
-    if (jobz != 'n' && jobz != 'N' && jobz != 'v' && jobz != 'V') {
-      rb_raise(rb_eArgError, "jobz must be 'N' or 'V'");
-    }
-
-    return jobz;
-  }
-
-  static char get_uplo(VALUE val) {
-    const char uplo = NUM2CHR(val);
-
-    if (uplo != 'u' && uplo != 'U' && uplo != 'l' && uplo != 'L') {
-      rb_raise(rb_eArgError, "uplo must be 'U' or 'L'");
-    }
-
-    return uplo;
-  }
-
-  static 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::TinyLinalg::Lapack.sygv does not support column major.");
-      break;
-    }
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg

data/ext/numo/tiny_linalg/lapack/hegvd.hpp

@@ -24,7 +24,7 @@ struct CHeGvd {
   }
 };
 
-template <int nary_dtype_id, int nary_rtype_id, typename DType, typename RType, typename FncType>
+template <int nary_dtype_id, int nary_rtype_id, typename dtype, typename rtype, class LapackFn>
 class HeGvd {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -40,15 +40,15 @@ private:
   };
 
   static void iter_hegvd(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* b = (DType*)NDL_PTR(lp, 1);
-    RType* w = (RType*)NDL_PTR(lp, 2);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* b = (dtype*)NDL_PTR(lp, 1);
+    rtype* w = (rtype*)NDL_PTR(lp, 2);
     int* info = (int*)NDL_PTR(lp, 3);
     hegvd_opt* opt = (hegvd_opt*)(lp->opt_ptr);
     const lapack_int n = NDL_SHAPE(lp, 0)[1];
     const lapack_int lda = NDL_SHAPE(lp, 0)[0];
     const lapack_int ldb = NDL_SHAPE(lp, 1)[0];
-    const lapack_int i = FncType().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
+    const lapack_int i = LapackFn().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
     *info = static_cast<int>(i);
   }
 
@@ -63,10 +63,10 @@ private:
     ID kw_table[4] = { rb_intern("itype"), rb_intern("jobz"), rb_intern("uplo"), rb_intern("order") };
     VALUE kw_values[4] = { Qundef, Qundef, Qundef, Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 4, kw_values);
-    const lapack_int itype = kw_values[0] != Qundef ? get_itype(kw_values[0]) : 1;
-    const char jobz = kw_values[1] != Qundef ? get_jobz(kw_values[1]) : 'V';
-    const char uplo = kw_values[2] != Qundef ? get_uplo(kw_values[2]) : 'U';
-    const int matrix_layout = kw_values[3] != Qundef ? get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
+    const lapack_int itype = kw_values[0] != Qundef ? Util().get_itype(kw_values[0]) : 1;
+    const char jobz = kw_values[1] != Qundef ? Util().get_jobz(kw_values[1]) : 'V';
+    const char uplo = kw_values[2] != Qundef ? Util().get_uplo(kw_values[2]) : 'U';
+    const int matrix_layout = kw_values[3] != Qundef ? Util().get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -116,52 +116,6 @@ private:
 
     return ret;
   }
-
-  static 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;
-  }
-
-  static char get_jobz(VALUE val) {
-    const char jobz = NUM2CHR(val);
-
-    if (jobz != 'n' && jobz != 'N' && jobz != 'v' && jobz != 'V') {
-      rb_raise(rb_eArgError, "jobz must be 'N' or 'V'");
-    }
-
-    return jobz;
-  }
-
-  static char get_uplo(VALUE val) {
-    const char uplo = NUM2CHR(val);
-
-    if (uplo != 'u' && uplo != 'U' && uplo != 'l' && uplo != 'L') {
-      rb_raise(rb_eArgError, "uplo must be 'U' or 'L'");
-    }
-
-    return uplo;
-  }
-
-  static 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::TinyLinalg::Lapack.sygvd does not support column major.");
-      break;
-    }
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg

data/ext/numo/tiny_linalg/lapack/hegvx.hpp

@@ -18,7 +18,7 @@ struct CHeGvx {
   }
 };
 
-template <int nary_dtype_id, int nary_rtype_id, typename DType, typename RType, typename FncType>
+template <int nary_dtype_id, int nary_rtype_id, typename dtype, typename rtype, class LapackFn>
 class HeGvx {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -32,18 +32,18 @@ private:
     char jobz;
     char range;
     char uplo;
-    RType vl;
-    RType vu;
+    rtype vl;
+    rtype vu;
     lapack_int il;
     lapack_int iu;
   };
 
   static void iter_hegvx(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* b = (DType*)NDL_PTR(lp, 1);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* b = (dtype*)NDL_PTR(lp, 1);
     int* m = (int*)NDL_PTR(lp, 2);
-    RType* w = (RType*)NDL_PTR(lp, 3);
-    DType* z = (DType*)NDL_PTR(lp, 4);
+    rtype* w = (rtype*)NDL_PTR(lp, 3);
+    dtype* z = (dtype*)NDL_PTR(lp, 4);
     int* ifail = (int*)NDL_PTR(lp, 5);
     int* info = (int*)NDL_PTR(lp, 6);
     hegvx_opt* opt = (hegvx_opt*)(lp->opt_ptr);
@@ -51,8 +51,8 @@ private:
     const lapack_int lda = NDL_SHAPE(lp, 0)[0];
     const lapack_int ldb = NDL_SHAPE(lp, 1)[0];
     const lapack_int ldz = opt->range != 'I' ? n : opt->iu - opt->il + 1;
-    const RType abstol = 0.0;
-    const lapack_int i = FncType().call(
+    const rtype abstol = 0.0;
+    const lapack_int i = LapackFn().call(
       opt->matrix_layout, opt->itype, opt->jobz, opt->range, opt->uplo, n, a, lda, b, ldb,
       opt->vl, opt->vu, opt->il, opt->iu, abstol, m, w, z, ldz, ifail);
     *info = static_cast<int>(i);
@@ -70,15 +70,15 @@ private:
                        rb_intern("vl"), rb_intern("vu"), rb_intern("il"), rb_intern("iu"), rb_intern("order") };
     VALUE kw_values[9] = { Qundef, Qundef, Qundef, Qundef, Qundef, Qundef, Qundef, Qundef, Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 9, kw_values);
-    const lapack_int itype = kw_values[0] != Qundef ? get_itype(kw_values[0]) : 1;
-    const char jobz = kw_values[1] != Qundef ? get_jobz(kw_values[1]) : 'V';
-    const char range = kw_values[2] != Qundef ? get_range(kw_values[2]) : 'A';
-    const char uplo = kw_values[3] != Qundef ? get_uplo(kw_values[3]) : 'U';
-    const RType vl = kw_values[4] != Qundef ? NUM2DBL(kw_values[4]) : 0.0;
-    const RType vu = kw_values[5] != Qundef ? NUM2DBL(kw_values[5]) : 0.0;
+    const lapack_int itype = kw_values[0] != Qundef ? Util().get_itype(kw_values[0]) : 1;
+    const char jobz = kw_values[1] != Qundef ? Util().get_jobz(kw_values[1]) : 'V';
+    const char range = kw_values[2] != Qundef ? Util().get_range(kw_values[2]) : 'A';
+    const char uplo = kw_values[3] != Qundef ? Util().get_uplo(kw_values[3]) : 'U';
+    const rtype vl = kw_values[4] != Qundef ? NUM2DBL(kw_values[4]) : 0.0;
+    const rtype vu = kw_values[5] != Qundef ? NUM2DBL(kw_values[5]) : 0.0;
     const lapack_int il = kw_values[6] != Qundef ? NUM2INT(kw_values[6]) : 0;
     const lapack_int iu = kw_values[7] != Qundef ? NUM2INT(kw_values[7]) : 0;
-    const int matrix_layout = kw_values[8] != Qundef ? get_matrix_layout(kw_values[8]) : LAPACK_ROW_MAJOR;
+    const int matrix_layout = kw_values[8] != Qundef ? Util().get_matrix_layout(kw_values[8]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -132,62 +132,6 @@ private:
 
     return ret;
   }
-
-  static 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;
-  }
-
-  static 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;
-  }
-
-  static 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;
-  }
-
-  static 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;
-  }
-
-  static 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::TinyLinalg::Lapack.hegvx does not support column major.");
-      break;
-    }
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg

data/ext/numo/tiny_linalg/lapack/orgqr.hpp

@@ -14,7 +14,7 @@ struct SOrgQr {
   }
 };
 
-template <int nary_dtype_id, typename DType, typename FncType>
+template <int nary_dtype_id, typename dtype, class LapackFn>
 class OrgQr {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -27,15 +27,15 @@ private:
   };
 
   static void iter_orgqr(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* tau = (DType*)NDL_PTR(lp, 1);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* tau = (dtype*)NDL_PTR(lp, 1);
     int* info = (int*)NDL_PTR(lp, 2);
     orgqr_opt* opt = (orgqr_opt*)(lp->opt_ptr);
     const lapack_int m = NDL_SHAPE(lp, 0)[0];
     const lapack_int n = NDL_SHAPE(lp, 0)[1];
     const lapack_int k = NDL_SHAPE(lp, 1)[0];
     const lapack_int lda = n;
-    const lapack_int i = FncType().call(opt->matrix_layout, m, n, k, a, lda, tau);
+    const lapack_int i = LapackFn().call(opt->matrix_layout, m, n, k, a, lda, tau);
     *info = static_cast<int>(i);
   }
 
@@ -49,7 +49,7 @@ private:
     ID kw_table[1] = { rb_intern("order") };
     VALUE kw_values[1] = { Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 1, kw_values);
-    const int matrix_layout = kw_values[0] != Qundef ? get_matrix_layout(kw_values[0]) : LAPACK_ROW_MAJOR;
+    const int matrix_layout = kw_values[0] != Qundef ? Util().get_matrix_layout(kw_values[0]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -90,26 +90,6 @@ private:
 
     return ret;
   }
-
-  static int get_matrix_layout(VALUE val) {
-    const char* option_str = StringValueCStr(val);
-
-    if (std::strlen(option_str) > 0) {
-      switch (option_str[0]) {
-      case 'r':
-      case 'R':
-        break;
-      case 'c':
-      case 'C':
-        rb_warn("Numo::TinyLinalg::Lapack.getrf does not support column major.");
-        break;
-      }
-    }
-
-    RB_GC_GUARD(val);
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg

data/ext/numo/tiny_linalg/lapack/sygv.hpp

@@ -16,7 +16,7 @@ struct SSyGv {
   }
 };
 
-template <int nary_dtype_id, typename DType, typename FncType>
+template <int nary_dtype_id, typename dtype, class LapackFn>
 class SyGv {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -32,15 +32,15 @@ private:
   };
 
   static void iter_sygv(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* b = (DType*)NDL_PTR(lp, 1);
-    DType* w = (DType*)NDL_PTR(lp, 2);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* b = (dtype*)NDL_PTR(lp, 1);
+    dtype* w = (dtype*)NDL_PTR(lp, 2);
     int* info = (int*)NDL_PTR(lp, 3);
     sygv_opt* opt = (sygv_opt*)(lp->opt_ptr);
     const lapack_int n = NDL_SHAPE(lp, 0)[1];
     const lapack_int lda = NDL_SHAPE(lp, 0)[0];
     const lapack_int ldb = NDL_SHAPE(lp, 1)[0];
-    const lapack_int i = FncType().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
+    const lapack_int i = LapackFn().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
     *info = static_cast<int>(i);
   }
 
@@ -54,10 +54,10 @@ private:
     ID kw_table[4] = { rb_intern("itype"), rb_intern("jobz"), rb_intern("uplo"), rb_intern("order") };
     VALUE kw_values[4] = { Qundef, Qundef, Qundef, Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 4, kw_values);
-    const lapack_int itype = kw_values[0] != Qundef ? get_itype(kw_values[0]) : 1;
-    const char jobz = kw_values[1] != Qundef ? get_jobz(kw_values[1]) : 'V';
-    const char uplo = kw_values[2] != Qundef ? get_uplo(kw_values[2]) : 'U';
-    const int matrix_layout = kw_values[3] != Qundef ? get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
+    const lapack_int itype = kw_values[0] != Qundef ? Util().get_itype(kw_values[0]) : 1;
+    const char jobz = kw_values[1] != Qundef ? Util().get_jobz(kw_values[1]) : 'V';
+    const char uplo = kw_values[2] != Qundef ? Util().get_uplo(kw_values[2]) : 'U';
+    const int matrix_layout = kw_values[3] != Qundef ? Util().get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -107,52 +107,6 @@ private:
 
     return ret;
   }
-
-  static 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;
-  }
-
-  static char get_jobz(VALUE val) {
-    const char jobz = NUM2CHR(val);
-
-    if (jobz != 'n' && jobz != 'N' && jobz != 'v' && jobz != 'V') {
-      rb_raise(rb_eArgError, "jobz must be 'N' or 'V'");
-    }
-
-    return jobz;
-  }
-
-  static char get_uplo(VALUE val) {
-    const char uplo = NUM2CHR(val);
-
-    if (uplo != 'u' && uplo != 'U' && uplo != 'l' && uplo != 'L') {
-      rb_raise(rb_eArgError, "uplo must be 'U' or 'L'");
-    }
-
-    return uplo;
-  }
-
-  static 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::TinyLinalg::Lapack.sygv does not support column major.");
-      break;
-    }
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg

data/ext/numo/tiny_linalg/lapack/sygvd.hpp

@@ -16,7 +16,7 @@ struct SSyGvd {
   }
 };
 
-template <int nary_dtype_id, typename DType, typename FncType>
+template <int nary_dtype_id, typename dtype, class LapackFn>
 class SyGvd {
 public:
   static void define_module_function(VALUE mLapack, const char* fnc_name) {
@@ -32,15 +32,15 @@ private:
   };
 
   static void iter_sygvd(na_loop_t* const lp) {
-    DType* a = (DType*)NDL_PTR(lp, 0);
-    DType* b = (DType*)NDL_PTR(lp, 1);
-    DType* w = (DType*)NDL_PTR(lp, 2);
+    dtype* a = (dtype*)NDL_PTR(lp, 0);
+    dtype* b = (dtype*)NDL_PTR(lp, 1);
+    dtype* w = (dtype*)NDL_PTR(lp, 2);
     int* info = (int*)NDL_PTR(lp, 3);
     sygvd_opt* opt = (sygvd_opt*)(lp->opt_ptr);
     const lapack_int n = NDL_SHAPE(lp, 0)[1];
     const lapack_int lda = NDL_SHAPE(lp, 0)[0];
     const lapack_int ldb = NDL_SHAPE(lp, 1)[0];
-    const lapack_int i = FncType().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
+    const lapack_int i = LapackFn().call(opt->matrix_layout, opt->itype, opt->jobz, opt->uplo, n, a, lda, b, ldb, w);
     *info = static_cast<int>(i);
   }
 
@@ -54,10 +54,10 @@ private:
     ID kw_table[4] = { rb_intern("itype"), rb_intern("jobz"), rb_intern("uplo"), rb_intern("order") };
     VALUE kw_values[4] = { Qundef, Qundef, Qundef, Qundef };
     rb_get_kwargs(kw_args, kw_table, 0, 4, kw_values);
-    const lapack_int itype = kw_values[0] != Qundef ? get_itype(kw_values[0]) : 1;
-    const char jobz = kw_values[1] != Qundef ? get_jobz(kw_values[1]) : 'V';
-    const char uplo = kw_values[2] != Qundef ? get_uplo(kw_values[2]) : 'U';
-    const int matrix_layout = kw_values[3] != Qundef ? get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
+    const lapack_int itype = kw_values[0] != Qundef ? Util().get_itype(kw_values[0]) : 1;
+    const char jobz = kw_values[1] != Qundef ? Util().get_jobz(kw_values[1]) : 'V';
+    const char uplo = kw_values[2] != Qundef ? Util().get_uplo(kw_values[2]) : 'U';
+    const int matrix_layout = kw_values[3] != Qundef ? Util().get_matrix_layout(kw_values[3]) : LAPACK_ROW_MAJOR;
 
     if (CLASS_OF(a_vnary) != nary_dtype) {
       a_vnary = rb_funcall(nary_dtype, rb_intern("cast"), 1, a_vnary);
@@ -107,52 +107,6 @@ private:
 
     return ret;
   }
-
-  static 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;
-  }
-
-  static char get_jobz(VALUE val) {
-    const char jobz = NUM2CHR(val);
-
-    if (jobz != 'n' && jobz != 'N' && jobz != 'v' && jobz != 'V') {
-      rb_raise(rb_eArgError, "jobz must be 'N' or 'V'");
-    }
-
-    return jobz;
-  }
-
-  static char get_uplo(VALUE val) {
-    const char uplo = NUM2CHR(val);
-
-    if (uplo != 'u' && uplo != 'U' && uplo != 'l' && uplo != 'L') {
-      rb_raise(rb_eArgError, "uplo must be 'U' or 'L'");
-    }
-
-    return uplo;
-  }
-
-  static 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::TinyLinalg::Lapack.sygvd does not support column major.");
-      break;
-    }
-
-    return LAPACK_ROW_MAJOR;
-  }
 };
 
 } // namespace TinyLinalg