ruby-lapack 1.7 → 1.7.1

data/ext/zunmtr.c

@@ -0,0 +1,129 @@
+#include "rb_lapack.h"
+
+extern VOID zunmtr_(char* side, char* uplo, char* trans, integer* m, integer* n, doublecomplex* a, integer* lda, doublecomplex* tau, doublecomplex* c, integer* ldc, doublecomplex* work, integer* lwork, integer* info);
+
+
+static VALUE
+rblapack_zunmtr(int argc, VALUE *argv, VALUE self){
+  VALUE rblapack_side;
+  char side; 
+  VALUE rblapack_uplo;
+  char uplo; 
+  VALUE rblapack_trans;
+  char trans; 
+  VALUE rblapack_a;
+  doublecomplex *a; 
+  VALUE rblapack_tau;
+  doublecomplex *tau; 
+  VALUE rblapack_c;
+  doublecomplex *c; 
+  VALUE rblapack_lwork;
+  integer lwork; 
+  VALUE rblapack_work;
+  doublecomplex *work; 
+  VALUE rblapack_info;
+  integer info; 
+  VALUE rblapack_c_out__;
+  doublecomplex *c_out__;
+
+  integer lda;
+  integer m;
+  integer ldc;
+  integer n;
+
+  VALUE rblapack_options;
+  if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
+    argc--;
+    rblapack_options = argv[argc];
+    if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
+      printf("%s\n", "USAGE:\n  work, info, c = NumRu::Lapack.zunmtr( side, uplo, trans, a, tau, c, [:lwork => lwork, :usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n      SUBROUTINE ZUNMTR( SIDE, UPLO, TRANS, M, N, A, LDA, TAU, C, LDC, WORK, LWORK, INFO )\n\n*  Purpose\n*  =======\n*\n*  ZUNMTR overwrites the general complex M-by-N matrix C with\n*\n*                  SIDE = 'L'     SIDE = 'R'\n*  TRANS = 'N':      Q * C          C * Q\n*  TRANS = 'C':      Q**H * C       C * Q**H\n*\n*  where Q is a complex unitary matrix of order nq, with nq = m if\n*  SIDE = 'L' and nq = n if SIDE = 'R'. Q is defined as the product of\n*  nq-1 elementary reflectors, as returned by ZHETRD:\n*\n*  if UPLO = 'U', Q = H(nq-1) . . . H(2) H(1);\n*\n*  if UPLO = 'L', Q = H(1) H(2) . . . H(nq-1).\n*\n\n*  Arguments\n*  =========\n*\n*  SIDE    (input) CHARACTER*1\n*          = 'L': apply Q or Q**H from the Left;\n*          = 'R': apply Q or Q**H from the Right.\n*\n*  UPLO    (input) CHARACTER*1\n*          = 'U': Upper triangle of A contains elementary reflectors\n*                 from ZHETRD;\n*          = 'L': Lower triangle of A contains elementary reflectors\n*                 from ZHETRD.\n*\n*  TRANS   (input) CHARACTER*1\n*          = 'N':  No transpose, apply Q;\n*          = 'C':  Conjugate transpose, apply Q**H.\n*\n*  M       (input) INTEGER\n*          The number of rows of the matrix C. M >= 0.\n*\n*  N       (input) INTEGER\n*          The number of columns of the matrix C. N >= 0.\n*\n*  A       (input) COMPLEX*16 array, dimension\n*                               (LDA,M) if SIDE = 'L'\n*                               (LDA,N) if SIDE = 'R'\n*          The vectors which define the elementary reflectors, as\n*          returned by ZHETRD.\n*\n*  LDA     (input) INTEGER\n*          The leading dimension of the array A.\n*          LDA >= max(1,M) if SIDE = 'L'; LDA >= max(1,N) if SIDE = 'R'.\n*\n*  TAU     (input) COMPLEX*16 array, dimension\n*                               (M-1) if SIDE = 'L'\n*                               (N-1) if SIDE = 'R'\n*          TAU(i) must contain the scalar factor of the elementary\n*          reflector H(i), as returned by ZHETRD.\n*\n*  C       (input/output) COMPLEX*16 array, dimension (LDC,N)\n*          On entry, the M-by-N matrix C.\n*          On exit, C is overwritten by Q*C or Q**H*C or C*Q**H or C*Q.\n*\n*  LDC     (input) INTEGER\n*          The leading dimension of the array C. LDC >= max(1,M).\n*\n*  WORK    (workspace/output) COMPLEX*16 array, dimension (MAX(1,LWORK))\n*          On exit, if INFO = 0, WORK(1) returns the optimal LWORK.\n*\n*  LWORK   (input) INTEGER\n*          The dimension of the array WORK.\n*          If SIDE = 'L', LWORK >= max(1,N);\n*          if SIDE = 'R', LWORK >= max(1,M).\n*          For optimum performance LWORK >= N*NB if SIDE = 'L', and\n*          LWORK >=M*NB if SIDE = 'R', where NB is the optimal\n*          blocksize.\n*\n*          If LWORK = -1, then a workspace query is assumed; the routine\n*          only calculates the optimal size of the WORK array, returns\n*          this value as the first entry of the WORK array, and no error\n*          message related to LWORK is issued by XERBLA.\n*\n*  INFO    (output) INTEGER\n*          = 0:  successful exit\n*          < 0:  if INFO = -i, the i-th argument had an illegal value\n*\n\n*  =====================================================================\n*\n*     .. Local Scalars ..\n      LOGICAL            LEFT, LQUERY, UPPER\n      INTEGER            I1, I2, IINFO, LWKOPT, MI, NB, NI, NQ, NW\n*     ..\n*     .. External Functions ..\n      LOGICAL            LSAME\n      INTEGER            ILAENV\n      EXTERNAL           LSAME, ILAENV\n*     ..\n*     .. External Subroutines ..\n      EXTERNAL           XERBLA, ZUNMQL, ZUNMQR\n*     ..\n*     .. Intrinsic Functions ..\n      INTRINSIC          MAX\n*     ..\n\n");
+      return Qnil;
+    }
+    if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
+      printf("%s\n", "USAGE:\n  work, info, c = NumRu::Lapack.zunmtr( side, uplo, trans, a, tau, c, [:lwork => lwork, :usage => usage, :help => help])\n");
+      return Qnil;
+    } 
+  } else
+    rblapack_options = Qnil;
+  if (argc != 6 && argc != 7)
+    rb_raise(rb_eArgError,"wrong number of arguments (%d for 6)", argc);
+  rblapack_side = argv[0];
+  rblapack_uplo = argv[1];
+  rblapack_trans = argv[2];
+  rblapack_a = argv[3];
+  rblapack_tau = argv[4];
+  rblapack_c = argv[5];
+  if (argc == 7) {
+    rblapack_lwork = argv[6];
+  } else if (rblapack_options != Qnil) {
+    rblapack_lwork = rb_hash_aref(rblapack_options, ID2SYM(rb_intern("lwork")));
+  } else {
+    rblapack_lwork = Qnil;
+  }
+
+  side = StringValueCStr(rblapack_side)[0];
+  trans = StringValueCStr(rblapack_trans)[0];
+  if (!NA_IsNArray(rblapack_c))
+    rb_raise(rb_eArgError, "c (6th argument) must be NArray");
+  if (NA_RANK(rblapack_c) != 2)
+    rb_raise(rb_eArgError, "rank of c (6th argument) must be %d", 2);
+  ldc = NA_SHAPE0(rblapack_c);
+  n = NA_SHAPE1(rblapack_c);
+  if (NA_TYPE(rblapack_c) != NA_DCOMPLEX)
+    rblapack_c = na_change_type(rblapack_c, NA_DCOMPLEX);
+  c = NA_PTR_TYPE(rblapack_c, doublecomplex*);
+  uplo = StringValueCStr(rblapack_uplo)[0];
+  if (!NA_IsNArray(rblapack_a))
+    rb_raise(rb_eArgError, "a (4th argument) must be NArray");
+  if (NA_RANK(rblapack_a) != 2)
+    rb_raise(rb_eArgError, "rank of a (4th argument) must be %d", 2);
+  lda = NA_SHAPE0(rblapack_a);
+  m = NA_SHAPE1(rblapack_a);
+  if (NA_TYPE(rblapack_a) != NA_DCOMPLEX)
+    rblapack_a = na_change_type(rblapack_a, NA_DCOMPLEX);
+  a = NA_PTR_TYPE(rblapack_a, doublecomplex*);
+  if (rblapack_lwork == Qnil)
+    lwork = lsame_(&side,"L") ? n : lsame_(&side,"R") ? m : 0;
+  else {
+    lwork = NUM2INT(rblapack_lwork);
+  }
+  if (!NA_IsNArray(rblapack_tau))
+    rb_raise(rb_eArgError, "tau (5th argument) must be NArray");
+  if (NA_RANK(rblapack_tau) != 1)
+    rb_raise(rb_eArgError, "rank of tau (5th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_tau) != (m-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of tau must be %d", m-1);
+  if (NA_TYPE(rblapack_tau) != NA_DCOMPLEX)
+    rblapack_tau = na_change_type(rblapack_tau, NA_DCOMPLEX);
+  tau = NA_PTR_TYPE(rblapack_tau, doublecomplex*);
+  {
+    na_shape_t shape[1];
+    shape[0] = MAX(1,lwork);
+    rblapack_work = na_make_object(NA_DCOMPLEX, 1, shape, cNArray);
+  }
+  work = NA_PTR_TYPE(rblapack_work, doublecomplex*);
+  {
+    na_shape_t shape[2];
+    shape[0] = ldc;
+    shape[1] = n;
+    rblapack_c_out__ = na_make_object(NA_DCOMPLEX, 2, shape, cNArray);
+  }
+  c_out__ = NA_PTR_TYPE(rblapack_c_out__, doublecomplex*);
+  MEMCPY(c_out__, c, doublecomplex, NA_TOTAL(rblapack_c));
+  rblapack_c = rblapack_c_out__;
+  c = c_out__;
+
+  zunmtr_(&side, &uplo, &trans, &m, &n, a, &lda, tau, c, &ldc, work, &lwork, &info);
+
+  rblapack_info = INT2NUM(info);
+  return rb_ary_new3(3, rblapack_work, rblapack_info, rblapack_c);
+}
+
+void
+init_lapack_zunmtr(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
+  sHelp = sH;
+  sUsage = sU;
+  rblapack_ZERO = zero;
+
+  rb_define_module_function(mLapack, "zunmtr", rblapack_zunmtr, -1);
+}

data/ext/zupgtr.c

@@ -0,0 +1,91 @@
+#include "rb_lapack.h"
+
+extern VOID zupgtr_(char* uplo, integer* n, doublecomplex* ap, doublecomplex* tau, doublecomplex* q, integer* ldq, doublecomplex* work, integer* info);
+
+
+static VALUE
+rblapack_zupgtr(int argc, VALUE *argv, VALUE self){
+  VALUE rblapack_uplo;
+  char uplo; 
+  VALUE rblapack_ap;
+  doublecomplex *ap; 
+  VALUE rblapack_tau;
+  doublecomplex *tau; 
+  VALUE rblapack_q;
+  doublecomplex *q; 
+  VALUE rblapack_info;
+  integer info; 
+  doublecomplex *work;
+
+  integer ldap;
+  integer ldtau;
+  integer ldq;
+  integer n;
+
+  VALUE rblapack_options;
+  if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
+    argc--;
+    rblapack_options = argv[argc];
+    if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
+      printf("%s\n", "USAGE:\n  q, info = NumRu::Lapack.zupgtr( uplo, ap, tau, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n      SUBROUTINE ZUPGTR( UPLO, N, AP, TAU, Q, LDQ, WORK, INFO )\n\n*  Purpose\n*  =======\n*\n*  ZUPGTR generates a complex unitary matrix Q which is defined as the\n*  product of n-1 elementary reflectors H(i) of order n, as returned by\n*  ZHPTRD using packed storage:\n*\n*  if UPLO = 'U', Q = H(n-1) . . . H(2) H(1),\n*\n*  if UPLO = 'L', Q = H(1) H(2) . . . H(n-1).\n*\n\n*  Arguments\n*  =========\n*\n*  UPLO    (input) CHARACTER*1\n*          = 'U': Upper triangular packed storage used in previous\n*                 call to ZHPTRD;\n*          = 'L': Lower triangular packed storage used in previous\n*                 call to ZHPTRD.\n*\n*  N       (input) INTEGER\n*          The order of the matrix Q. N >= 0.\n*\n*  AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)\n*          The vectors which define the elementary reflectors, as\n*          returned by ZHPTRD.\n*\n*  TAU     (input) COMPLEX*16 array, dimension (N-1)\n*          TAU(i) must contain the scalar factor of the elementary\n*          reflector H(i), as returned by ZHPTRD.\n*\n*  Q       (output) COMPLEX*16 array, dimension (LDQ,N)\n*          The N-by-N unitary matrix Q.\n*\n*  LDQ     (input) INTEGER\n*          The leading dimension of the array Q. LDQ >= max(1,N).\n*\n*  WORK    (workspace) COMPLEX*16 array, dimension (N-1)\n*\n*  INFO    (output) INTEGER\n*          = 0:  successful exit\n*          < 0:  if INFO = -i, the i-th argument had an illegal value\n*\n\n*  =====================================================================\n*\n\n");
+      return Qnil;
+    }
+    if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
+      printf("%s\n", "USAGE:\n  q, info = NumRu::Lapack.zupgtr( uplo, ap, tau, [:usage => usage, :help => help])\n");
+      return Qnil;
+    } 
+  } else
+    rblapack_options = Qnil;
+  if (argc != 3 && argc != 3)
+    rb_raise(rb_eArgError,"wrong number of arguments (%d for 3)", argc);
+  rblapack_uplo = argv[0];
+  rblapack_ap = argv[1];
+  rblapack_tau = argv[2];
+  if (argc == 3) {
+  } else if (rblapack_options != Qnil) {
+  } else {
+  }
+
+  uplo = StringValueCStr(rblapack_uplo)[0];
+  if (!NA_IsNArray(rblapack_tau))
+    rb_raise(rb_eArgError, "tau (3th argument) must be NArray");
+  if (NA_RANK(rblapack_tau) != 1)
+    rb_raise(rb_eArgError, "rank of tau (3th argument) must be %d", 1);
+  ldtau = NA_SHAPE0(rblapack_tau);
+  if (NA_TYPE(rblapack_tau) != NA_DCOMPLEX)
+    rblapack_tau = na_change_type(rblapack_tau, NA_DCOMPLEX);
+  tau = NA_PTR_TYPE(rblapack_tau, doublecomplex*);
+  n = ldtau+1;
+  if (!NA_IsNArray(rblapack_ap))
+    rb_raise(rb_eArgError, "ap (2th argument) must be NArray");
+  if (NA_RANK(rblapack_ap) != 1)
+    rb_raise(rb_eArgError, "rank of ap (2th argument) must be %d", 1);
+  ldap = NA_SHAPE0(rblapack_ap);
+  if (NA_TYPE(rblapack_ap) != NA_DCOMPLEX)
+    rblapack_ap = na_change_type(rblapack_ap, NA_DCOMPLEX);
+  ap = NA_PTR_TYPE(rblapack_ap, doublecomplex*);
+  ldq = MAX(1,n);
+  {
+    na_shape_t shape[2];
+    shape[0] = ldq;
+    shape[1] = n;
+    rblapack_q = na_make_object(NA_DCOMPLEX, 2, shape, cNArray);
+  }
+  q = NA_PTR_TYPE(rblapack_q, doublecomplex*);
+  work = ALLOC_N(doublecomplex, (n-1));
+
+  zupgtr_(&uplo, &n, ap, tau, q, &ldq, work, &info);
+
+  free(work);
+  rblapack_info = INT2NUM(info);
+  return rb_ary_new3(2, rblapack_q, rblapack_info);
+}
+
+void
+init_lapack_zupgtr(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
+  sHelp = sH;
+  sUsage = sU;
+  rblapack_ZERO = zero;
+
+  rb_define_module_function(mLapack, "zupgtr", rblapack_zupgtr, -1);
+}

data/ext/zupmtr.c

@@ -0,0 +1,116 @@
+#include "rb_lapack.h"
+
+extern VOID zupmtr_(char* side, char* uplo, char* trans, integer* m, integer* n, doublecomplex* ap, doublecomplex* tau, doublecomplex* c, integer* ldc, doublecomplex* work, integer* info);
+
+
+static VALUE
+rblapack_zupmtr(int argc, VALUE *argv, VALUE self){
+  VALUE rblapack_side;
+  char side; 
+  VALUE rblapack_uplo;
+  char uplo; 
+  VALUE rblapack_trans;
+  char trans; 
+  VALUE rblapack_m;
+  integer m; 
+  VALUE rblapack_ap;
+  doublecomplex *ap; 
+  VALUE rblapack_tau;
+  doublecomplex *tau; 
+  VALUE rblapack_c;
+  doublecomplex *c; 
+  VALUE rblapack_info;
+  integer info; 
+  VALUE rblapack_c_out__;
+  doublecomplex *c_out__;
+  doublecomplex *work;
+
+  integer ldc;
+  integer n;
+
+  VALUE rblapack_options;
+  if (argc > 0 && TYPE(argv[argc-1]) == T_HASH) {
+    argc--;
+    rblapack_options = argv[argc];
+    if (rb_hash_aref(rblapack_options, sHelp) == Qtrue) {
+      printf("%s\n", "USAGE:\n  info, c = NumRu::Lapack.zupmtr( side, uplo, trans, m, ap, tau, c, [:usage => usage, :help => help])\n\n\nFORTRAN MANUAL\n      SUBROUTINE ZUPMTR( SIDE, UPLO, TRANS, M, N, AP, TAU, C, LDC, WORK, INFO )\n\n*  Purpose\n*  =======\n*\n*  ZUPMTR overwrites the general complex M-by-N matrix C with\n*\n*                  SIDE = 'L'     SIDE = 'R'\n*  TRANS = 'N':      Q * C          C * Q\n*  TRANS = 'C':      Q**H * C       C * Q**H\n*\n*  where Q is a complex unitary matrix of order nq, with nq = m if\n*  SIDE = 'L' and nq = n if SIDE = 'R'. Q is defined as the product of\n*  nq-1 elementary reflectors, as returned by ZHPTRD using packed\n*  storage:\n*\n*  if UPLO = 'U', Q = H(nq-1) . . . H(2) H(1);\n*\n*  if UPLO = 'L', Q = H(1) H(2) . . . H(nq-1).\n*\n\n*  Arguments\n*  =========\n*\n*  SIDE    (input) CHARACTER*1\n*          = 'L': apply Q or Q**H from the Left;\n*          = 'R': apply Q or Q**H from the Right.\n*\n*  UPLO    (input) CHARACTER*1\n*          = 'U': Upper triangular packed storage used in previous\n*                 call to ZHPTRD;\n*          = 'L': Lower triangular packed storage used in previous\n*                 call to ZHPTRD.\n*\n*  TRANS   (input) CHARACTER*1\n*          = 'N':  No transpose, apply Q;\n*          = 'C':  Conjugate transpose, apply Q**H.\n*\n*  M       (input) INTEGER\n*          The number of rows of the matrix C. M >= 0.\n*\n*  N       (input) INTEGER\n*          The number of columns of the matrix C. N >= 0.\n*\n*  AP      (input) COMPLEX*16 array, dimension\n*                               (M*(M+1)/2) if SIDE = 'L'\n*                               (N*(N+1)/2) if SIDE = 'R'\n*          The vectors which define the elementary reflectors, as\n*          returned by ZHPTRD.  AP is modified by the routine but\n*          restored on exit.\n*\n*  TAU     (input) COMPLEX*16 array, dimension (M-1) if SIDE = 'L'\n*                                     or (N-1) if SIDE = 'R'\n*          TAU(i) must contain the scalar factor of the elementary\n*          reflector H(i), as returned by ZHPTRD.\n*\n*  C       (input/output) COMPLEX*16 array, dimension (LDC,N)\n*          On entry, the M-by-N matrix C.\n*          On exit, C is overwritten by Q*C or Q**H*C or C*Q**H or C*Q.\n*\n*  LDC     (input) INTEGER\n*          The leading dimension of the array C. LDC >= max(1,M).\n*\n*  WORK    (workspace) COMPLEX*16 array, dimension\n*                                   (N) if SIDE = 'L'\n*                                   (M) if SIDE = 'R'\n*\n*  INFO    (output) INTEGER\n*          = 0:  successful exit\n*          < 0:  if INFO = -i, the i-th argument had an illegal value\n*\n\n*  =====================================================================\n*\n\n");
+      return Qnil;
+    }
+    if (rb_hash_aref(rblapack_options, sUsage) == Qtrue) {
+      printf("%s\n", "USAGE:\n  info, c = NumRu::Lapack.zupmtr( side, uplo, trans, m, ap, tau, c, [:usage => usage, :help => help])\n");
+      return Qnil;
+    } 
+  } else
+    rblapack_options = Qnil;
+  if (argc != 7 && argc != 7)
+    rb_raise(rb_eArgError,"wrong number of arguments (%d for 7)", argc);
+  rblapack_side = argv[0];
+  rblapack_uplo = argv[1];
+  rblapack_trans = argv[2];
+  rblapack_m = argv[3];
+  rblapack_ap = argv[4];
+  rblapack_tau = argv[5];
+  rblapack_c = argv[6];
+  if (argc == 7) {
+  } else if (rblapack_options != Qnil) {
+  } else {
+  }
+
+  side = StringValueCStr(rblapack_side)[0];
+  trans = StringValueCStr(rblapack_trans)[0];
+  if (!NA_IsNArray(rblapack_c))
+    rb_raise(rb_eArgError, "c (7th argument) must be NArray");
+  if (NA_RANK(rblapack_c) != 2)
+    rb_raise(rb_eArgError, "rank of c (7th argument) must be %d", 2);
+  ldc = NA_SHAPE0(rblapack_c);
+  n = NA_SHAPE1(rblapack_c);
+  if (NA_TYPE(rblapack_c) != NA_DCOMPLEX)
+    rblapack_c = na_change_type(rblapack_c, NA_DCOMPLEX);
+  c = NA_PTR_TYPE(rblapack_c, doublecomplex*);
+  uplo = StringValueCStr(rblapack_uplo)[0];
+  m = NUM2INT(rblapack_m);
+  if (!NA_IsNArray(rblapack_tau))
+    rb_raise(rb_eArgError, "tau (6th argument) must be NArray");
+  if (NA_RANK(rblapack_tau) != 1)
+    rb_raise(rb_eArgError, "rank of tau (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_tau) != (m-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of tau must be %d", m-1);
+  if (NA_TYPE(rblapack_tau) != NA_DCOMPLEX)
+    rblapack_tau = na_change_type(rblapack_tau, NA_DCOMPLEX);
+  tau = NA_PTR_TYPE(rblapack_tau, doublecomplex*);
+  if (!NA_IsNArray(rblapack_ap))
+    rb_raise(rb_eArgError, "ap (5th argument) must be NArray");
+  if (NA_RANK(rblapack_ap) != 1)
+    rb_raise(rb_eArgError, "rank of ap (5th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_ap) != (m*(m+1)/2))
+    rb_raise(rb_eRuntimeError, "shape 0 of ap must be %d", m*(m+1)/2);
+  if (NA_TYPE(rblapack_ap) != NA_DCOMPLEX)
+    rblapack_ap = na_change_type(rblapack_ap, NA_DCOMPLEX);
+  ap = NA_PTR_TYPE(rblapack_ap, doublecomplex*);
+  {
+    na_shape_t shape[2];
+    shape[0] = ldc;
+    shape[1] = n;
+    rblapack_c_out__ = na_make_object(NA_DCOMPLEX, 2, shape, cNArray);
+  }
+  c_out__ = NA_PTR_TYPE(rblapack_c_out__, doublecomplex*);
+  MEMCPY(c_out__, c, doublecomplex, NA_TOTAL(rblapack_c));
+  rblapack_c = rblapack_c_out__;
+  c = c_out__;
+  work = ALLOC_N(doublecomplex, (lsame_(&side,"L") ? n : lsame_(&side,"R") ? m : 0));
+
+  zupmtr_(&side, &uplo, &trans, &m, &n, ap, tau, c, &ldc, work, &info);
+
+  free(work);
+  rblapack_info = INT2NUM(info);
+  return rb_ary_new3(2, rblapack_info, rblapack_c);
+}
+
+void
+init_lapack_zupmtr(VALUE mLapack, VALUE sH, VALUE sU, VALUE zero){
+  sHelp = sH;
+  sUsage = sU;
+  rblapack_ZERO = zero;
+
+  rb_define_module_function(mLapack, "zupmtr", rblapack_zupmtr, -1);
+}

metadata

@@ -1,7 +1,7 @@
 --- !ruby/object:Gem::Specification
 name: ruby-lapack
 version: !ruby/object:Gem::Version
-  version: '1.7'
+  version: 1.7.1
 platform: ruby
 authors:
 - Seiya Nishizawa
@@ -1670,9 +1670,1640 @@ files:
 - dev/make_csrc.rb
 - dev/mkdoc.rb
 - dev/parse.rb
+- ext/cbbcsd.c
+- ext/cbdsqr.c
+- ext/cgbbrd.c
+- ext/cgbcon.c
+- ext/cgbequ.c
+- ext/cgbequb.c
+- ext/cgbrfs.c
+- ext/cgbrfsx.c
+- ext/cgbsv.c
+- ext/cgbsvx.c
+- ext/cgbsvxx.c
+- ext/cgbtf2.c
+- ext/cgbtrf.c
+- ext/cgbtrs.c
+- ext/cgebak.c
+- ext/cgebal.c
+- ext/cgebd2.c
+- ext/cgebrd.c
+- ext/cgecon.c
+- ext/cgeequ.c
+- ext/cgeequb.c
+- ext/cgees.c
+- ext/cgeesx.c
+- ext/cgeev.c
+- ext/cgeevx.c
+- ext/cgegs.c
+- ext/cgegv.c
+- ext/cgehd2.c
+- ext/cgehrd.c
+- ext/cgelq2.c
+- ext/cgelqf.c
+- ext/cgels.c
+- ext/cgelsd.c
+- ext/cgelss.c
+- ext/cgelsx.c
+- ext/cgelsy.c
+- ext/cgeql2.c
+- ext/cgeqlf.c
+- ext/cgeqp3.c
+- ext/cgeqpf.c
+- ext/cgeqr2.c
+- ext/cgeqr2p.c
+- ext/cgeqrf.c
+- ext/cgeqrfp.c
+- ext/cgerfs.c
+- ext/cgerfsx.c
+- ext/cgerq2.c
+- ext/cgerqf.c
+- ext/cgesc2.c
+- ext/cgesdd.c
+- ext/cgesv.c
+- ext/cgesvd.c
+- ext/cgesvx.c
+- ext/cgesvxx.c
+- ext/cgetc2.c
+- ext/cgetf2.c
+- ext/cgetrf.c
+- ext/cgetri.c
+- ext/cgetrs.c
+- ext/cggbak.c
+- ext/cggbal.c
+- ext/cgges.c
+- ext/cggesx.c
+- ext/cggev.c
+- ext/cggevx.c
+- ext/cggglm.c
+- ext/cgghrd.c
+- ext/cgglse.c
+- ext/cggqrf.c
+- ext/cggrqf.c
+- ext/cggsvd.c
+- ext/cggsvp.c
+- ext/cgtcon.c
+- ext/cgtrfs.c
+- ext/cgtsv.c
+- ext/cgtsvx.c
+- ext/cgttrf.c
+- ext/cgttrs.c
+- ext/cgtts2.c
+- ext/chbev.c
+- ext/chbevd.c
+- ext/chbevx.c
+- ext/chbgst.c
+- ext/chbgv.c
+- ext/chbgvd.c
+- ext/chbgvx.c
+- ext/chbtrd.c
+- ext/checon.c
+- ext/cheequb.c
+- ext/cheev.c
+- ext/cheevd.c
+- ext/cheevr.c
+- ext/cheevx.c
+- ext/chegs2.c
+- ext/chegst.c
+- ext/chegv.c
+- ext/chegvd.c
+- ext/chegvx.c
+- ext/cherfs.c
+- ext/cherfsx.c
+- ext/chesv.c
+- ext/chesvx.c
+- ext/chesvxx.c
+- ext/chetd2.c
+- ext/chetf2.c
+- ext/chetrd.c
+- ext/chetrf.c
+- ext/chetri.c
+- ext/chetrs.c
+- ext/chetrs2.c
+- ext/chfrk.c
+- ext/chgeqz.c
+- ext/chla_transtype.c
+- ext/chpcon.c
+- ext/chpev.c
+- ext/chpevd.c
+- ext/chpevx.c
+- ext/chpgst.c
+- ext/chpgv.c
+- ext/chpgvd.c
+- ext/chpgvx.c
+- ext/chprfs.c
+- ext/chpsv.c
+- ext/chpsvx.c
+- ext/chptrd.c
+- ext/chptrf.c
+- ext/chptri.c
+- ext/chptrs.c
+- ext/chsein.c
+- ext/chseqr.c
+- ext/cla_gbamv.c
+- ext/cla_gbrcond_c.c
+- ext/cla_gbrcond_x.c
+- ext/cla_gbrfsx_extended.c
+- ext/cla_gbrpvgrw.c
+- ext/cla_geamv.c
+- ext/cla_gercond_c.c
+- ext/cla_gercond_x.c
+- ext/cla_gerfsx_extended.c
+- ext/cla_heamv.c
+- ext/cla_hercond_c.c
+- ext/cla_hercond_x.c
+- ext/cla_herfsx_extended.c
+- ext/cla_herpvgrw.c
+- ext/cla_lin_berr.c
+- ext/cla_porcond_c.c
+- ext/cla_porcond_x.c
+- ext/cla_porfsx_extended.c
+- ext/cla_porpvgrw.c
+- ext/cla_rpvgrw.c
+- ext/cla_syamv.c
+- ext/cla_syrcond_c.c
+- ext/cla_syrcond_x.c
+- ext/cla_syrfsx_extended.c
+- ext/cla_syrpvgrw.c
+- ext/cla_wwaddw.c
+- ext/clabrd.c
+- ext/clacgv.c
+- ext/clacn2.c
+- ext/clacon.c
+- ext/clacp2.c
+- ext/clacpy.c
+- ext/clacrm.c
+- ext/clacrt.c
+- ext/cladiv.c
+- ext/claed0.c
+- ext/claed7.c
+- ext/claed8.c
+- ext/claein.c
+- ext/claesy.c
+- ext/claev2.c
+- ext/clag2z.c
+- ext/clags2.c
+- ext/clagtm.c
+- ext/clahef.c
+- ext/clahqr.c
+- ext/clahr2.c
+- ext/clahrd.c
+- ext/claic1.c
+- ext/clals0.c
+- ext/clalsa.c
+- ext/clalsd.c
+- ext/clangb.c
+- ext/clange.c
+- ext/clangt.c
+- ext/clanhb.c
+- ext/clanhe.c
+- ext/clanhf.c
+- ext/clanhp.c
+- ext/clanhs.c
+- ext/clanht.c
+- ext/clansb.c
+- ext/clansp.c
+- ext/clansy.c
+- ext/clantb.c
+- ext/clantp.c
+- ext/clantr.c
+- ext/clapll.c
+- ext/clapmr.c
+- ext/clapmt.c
+- ext/claqgb.c
+- ext/claqge.c
+- ext/claqhb.c
+- ext/claqhe.c
+- ext/claqhp.c
+- ext/claqp2.c
+- ext/claqps.c
+- ext/claqr0.c
+- ext/claqr1.c
+- ext/claqr2.c
+- ext/claqr3.c
+- ext/claqr4.c
+- ext/claqr5.c
+- ext/claqsb.c
+- ext/claqsp.c
+- ext/claqsy.c
+- ext/clar1v.c
+- ext/clar2v.c
+- ext/clarcm.c
+- ext/clarf.c
+- ext/clarfb.c
+- ext/clarfg.c
+- ext/clarfgp.c
+- ext/clarft.c
+- ext/clarfx.c
+- ext/clargv.c
+- ext/clarnv.c
+- ext/clarrv.c
+- ext/clarscl2.c
+- ext/clartg.c
+- ext/clartv.c
+- ext/clarz.c
+- ext/clarzb.c
+- ext/clarzt.c
+- ext/clascl.c
+- ext/clascl2.c
+- ext/claset.c
+- ext/clasr.c
+- ext/classq.c
+- ext/claswp.c
+- ext/clasyf.c
+- ext/clatbs.c
+- ext/clatdf.c
+- ext/clatps.c
+- ext/clatrd.c
+- ext/clatrs.c
+- ext/clatrz.c
+- ext/clatzm.c
+- ext/clauu2.c
+- ext/clauum.c
+- ext/cpbcon.c
+- ext/cpbequ.c
+- ext/cpbrfs.c
+- ext/cpbstf.c
+- ext/cpbsv.c
+- ext/cpbsvx.c
+- ext/cpbtf2.c
+- ext/cpbtrf.c
+- ext/cpbtrs.c
+- ext/cpftrf.c
+- ext/cpftri.c
+- ext/cpftrs.c
+- ext/cpocon.c
+- ext/cpoequ.c
+- ext/cpoequb.c
+- ext/cporfs.c
+- ext/cporfsx.c
+- ext/cposv.c
+- ext/cposvx.c
+- ext/cposvxx.c
+- ext/cpotf2.c
+- ext/cpotrf.c
+- ext/cpotri.c
+- ext/cpotrs.c
+- ext/cppcon.c
+- ext/cppequ.c
+- ext/cpprfs.c
+- ext/cppsv.c
+- ext/cppsvx.c
+- ext/cpptrf.c
+- ext/cpptri.c
+- ext/cpptrs.c
+- ext/cpstf2.c
+- ext/cpstrf.c
+- ext/cptcon.c
+- ext/cpteqr.c
+- ext/cptrfs.c
+- ext/cptsv.c
+- ext/cptsvx.c
+- ext/cpttrf.c
+- ext/cpttrs.c
+- ext/cptts2.c
+- ext/crot.c
+- ext/cspcon.c
+- ext/cspmv.c
+- ext/cspr.c
+- ext/csprfs.c
+- ext/cspsv.c
+- ext/cspsvx.c
+- ext/csptrf.c
+- ext/csptri.c
+- ext/csptrs.c
+- ext/csrscl.c
+- ext/cstedc.c
+- ext/cstegr.c
+- ext/cstein.c
+- ext/cstemr.c
+- ext/csteqr.c
+- ext/csycon.c
+- ext/csyconv.c
+- ext/csyequb.c
+- ext/csymv.c
+- ext/csyr.c
+- ext/csyrfs.c
+- ext/csyrfsx.c
+- ext/csysv.c
+- ext/csysvx.c
+- ext/csysvxx.c
+- ext/csyswapr.c
+- ext/csytf2.c
+- ext/csytrf.c
+- ext/csytri.c
+- ext/csytri2.c
+- ext/csytri2x.c
+- ext/csytrs.c
+- ext/csytrs2.c
+- ext/ctbcon.c
+- ext/ctbrfs.c
+- ext/ctbtrs.c
+- ext/ctfsm.c
+- ext/ctftri.c
+- ext/ctfttp.c
+- ext/ctfttr.c
+- ext/ctgevc.c
+- ext/ctgex2.c
+- ext/ctgexc.c
+- ext/ctgsen.c
+- ext/ctgsja.c
+- ext/ctgsna.c
+- ext/ctgsy2.c
+- ext/ctgsyl.c
+- ext/ctpcon.c
+- ext/ctprfs.c
+- ext/ctptri.c
+- ext/ctptrs.c
+- ext/ctpttf.c
+- ext/ctpttr.c
+- ext/ctrcon.c
+- ext/ctrevc.c
+- ext/ctrexc.c
+- ext/ctrrfs.c
+- ext/ctrsen.c
+- ext/ctrsna.c
+- ext/ctrsyl.c
+- ext/ctrti2.c
+- ext/ctrtri.c
+- ext/ctrtrs.c
+- ext/ctrttf.c
+- ext/ctrttp.c
+- ext/ctzrqf.c
+- ext/ctzrzf.c
+- ext/cunbdb.c
+- ext/cuncsd.c
+- ext/cung2l.c
+- ext/cung2r.c
+- ext/cungbr.c
+- ext/cunghr.c
+- ext/cungl2.c
+- ext/cunglq.c
+- ext/cungql.c
+- ext/cungqr.c
+- ext/cungr2.c
+- ext/cungrq.c
+- ext/cungtr.c
+- ext/cunm2l.c
+- ext/cunm2r.c
+- ext/cunmbr.c
+- ext/cunmhr.c
+- ext/cunml2.c
+- ext/cunmlq.c
+- ext/cunmql.c
+- ext/cunmqr.c
+- ext/cunmr2.c
+- ext/cunmr3.c
+- ext/cunmrq.c
+- ext/cunmrz.c
+- ext/cunmtr.c
+- ext/cupgtr.c
+- ext/cupmtr.c
+- ext/dbbcsd.c
+- ext/dbdsdc.c
+- ext/dbdsqr.c
+- ext/ddisna.c
+- ext/dgbbrd.c
+- ext/dgbcon.c
+- ext/dgbequ.c
+- ext/dgbequb.c
+- ext/dgbrfs.c
+- ext/dgbrfsx.c
+- ext/dgbsv.c
+- ext/dgbsvx.c
+- ext/dgbsvxx.c
+- ext/dgbtf2.c
+- ext/dgbtrf.c
+- ext/dgbtrs.c
+- ext/dgebak.c
+- ext/dgebal.c
+- ext/dgebd2.c
+- ext/dgebrd.c
+- ext/dgecon.c
+- ext/dgeequ.c
+- ext/dgeequb.c
+- ext/dgees.c
+- ext/dgeesx.c
+- ext/dgeev.c
+- ext/dgeevx.c
+- ext/dgegs.c
+- ext/dgegv.c
+- ext/dgehd2.c
+- ext/dgehrd.c
+- ext/dgejsv.c
+- ext/dgelq2.c
+- ext/dgelqf.c
+- ext/dgels.c
+- ext/dgelsd.c
+- ext/dgelss.c
+- ext/dgelsx.c
+- ext/dgelsy.c
+- ext/dgeql2.c
+- ext/dgeqlf.c
+- ext/dgeqp3.c
+- ext/dgeqpf.c
+- ext/dgeqr2.c
+- ext/dgeqr2p.c
+- ext/dgeqrf.c
+- ext/dgeqrfp.c
+- ext/dgerfs.c
+- ext/dgerfsx.c
+- ext/dgerq2.c
+- ext/dgerqf.c
+- ext/dgesc2.c
+- ext/dgesdd.c
+- ext/dgesv.c
+- ext/dgesvd.c
+- ext/dgesvj.c
+- ext/dgesvx.c
+- ext/dgesvxx.c
+- ext/dgetc2.c
+- ext/dgetf2.c
+- ext/dgetrf.c
+- ext/dgetri.c
+- ext/dgetrs.c
+- ext/dggbak.c
+- ext/dggbal.c
+- ext/dgges.c
+- ext/dggesx.c
+- ext/dggev.c
+- ext/dggevx.c
+- ext/dggglm.c
+- ext/dgghrd.c
+- ext/dgglse.c
+- ext/dggqrf.c
+- ext/dggrqf.c
+- ext/dggsvd.c
+- ext/dggsvp.c
+- ext/dgsvj0.c
+- ext/dgsvj1.c
+- ext/dgtcon.c
+- ext/dgtrfs.c
+- ext/dgtsv.c
+- ext/dgtsvx.c
+- ext/dgttrf.c
+- ext/dgttrs.c
+- ext/dgtts2.c
+- ext/dhgeqz.c
+- ext/dhsein.c
+- ext/dhseqr.c
+- ext/disnan.c
+- ext/dla_gbamv.c
+- ext/dla_gbrcond.c
+- ext/dla_gbrfsx_extended.c
+- ext/dla_gbrpvgrw.c
+- ext/dla_geamv.c
+- ext/dla_gercond.c
+- ext/dla_gerfsx_extended.c
+- ext/dla_lin_berr.c
+- ext/dla_porcond.c
+- ext/dla_porfsx_extended.c
+- ext/dla_porpvgrw.c
+- ext/dla_rpvgrw.c
+- ext/dla_syamv.c
+- ext/dla_syrcond.c
+- ext/dla_syrfsx_extended.c
+- ext/dla_syrpvgrw.c
+- ext/dla_wwaddw.c
+- ext/dlabad.c
+- ext/dlabrd.c
+- ext/dlacn2.c
+- ext/dlacon.c
+- ext/dlacpy.c
+- ext/dladiv.c
+- ext/dlae2.c
+- ext/dlaebz.c
+- ext/dlaed0.c
+- ext/dlaed1.c
+- ext/dlaed2.c
+- ext/dlaed3.c
+- ext/dlaed4.c
+- ext/dlaed5.c
+- ext/dlaed6.c
+- ext/dlaed7.c
+- ext/dlaed8.c
+- ext/dlaed9.c
+- ext/dlaeda.c
+- ext/dlaein.c
+- ext/dlaev2.c
+- ext/dlaexc.c
+- ext/dlag2.c
+- ext/dlag2s.c
+- ext/dlags2.c
+- ext/dlagtf.c
+- ext/dlagtm.c
+- ext/dlagts.c
+- ext/dlagv2.c
+- ext/dlahqr.c
+- ext/dlahr2.c
+- ext/dlahrd.c
+- ext/dlaic1.c
+- ext/dlaln2.c
+- ext/dlals0.c
+- ext/dlalsa.c
+- ext/dlalsd.c
+- ext/dlamrg.c
+- ext/dlaneg.c
+- ext/dlangb.c
+- ext/dlange.c
+- ext/dlangt.c
+- ext/dlanhs.c
+- ext/dlansb.c
+- ext/dlansf.c
+- ext/dlansp.c
+- ext/dlanst.c
+- ext/dlansy.c
+- ext/dlantb.c
+- ext/dlantp.c
+- ext/dlantr.c
+- ext/dlanv2.c
+- ext/dlapll.c
+- ext/dlapmr.c
+- ext/dlapmt.c
+- ext/dlapy2.c
+- ext/dlapy3.c
+- ext/dlaqgb.c
+- ext/dlaqge.c
+- ext/dlaqp2.c
+- ext/dlaqps.c
+- ext/dlaqr0.c
+- ext/dlaqr1.c
+- ext/dlaqr2.c
+- ext/dlaqr3.c
+- ext/dlaqr4.c
+- ext/dlaqr5.c
+- ext/dlaqsb.c
+- ext/dlaqsp.c
+- ext/dlaqsy.c
+- ext/dlaqtr.c
+- ext/dlar1v.c
+- ext/dlar2v.c
+- ext/dlarf.c
+- ext/dlarfb.c
+- ext/dlarfg.c
+- ext/dlarfgp.c
+- ext/dlarft.c
+- ext/dlarfx.c
+- ext/dlargv.c
+- ext/dlarnv.c
+- ext/dlarra.c
+- ext/dlarrb.c
+- ext/dlarrc.c
+- ext/dlarrd.c
+- ext/dlarre.c
+- ext/dlarrf.c
+- ext/dlarrj.c
+- ext/dlarrk.c
+- ext/dlarrr.c
+- ext/dlarrv.c
+- ext/dlarscl2.c
+- ext/dlartg.c
+- ext/dlartgp.c
+- ext/dlartgs.c
+- ext/dlartv.c
+- ext/dlaruv.c
+- ext/dlarz.c
+- ext/dlarzb.c
+- ext/dlarzt.c
+- ext/dlas2.c
+- ext/dlascl.c
+- ext/dlascl2.c
+- ext/dlasd0.c
+- ext/dlasd1.c
+- ext/dlasd2.c
+- ext/dlasd3.c
+- ext/dlasd4.c
+- ext/dlasd5.c
+- ext/dlasd6.c
+- ext/dlasd7.c
+- ext/dlasd8.c
+- ext/dlasda.c
+- ext/dlasdq.c
+- ext/dlasdt.c
+- ext/dlaset.c
+- ext/dlasq1.c
+- ext/dlasq2.c
+- ext/dlasq3.c
+- ext/dlasq4.c
+- ext/dlasq5.c
+- ext/dlasq6.c
+- ext/dlasr.c
+- ext/dlasrt.c
+- ext/dlassq.c
+- ext/dlasv2.c
+- ext/dlaswp.c
+- ext/dlasy2.c
+- ext/dlasyf.c
+- ext/dlat2s.c
+- ext/dlatbs.c
+- ext/dlatdf.c
+- ext/dlatps.c
+- ext/dlatrd.c
+- ext/dlatrs.c
+- ext/dlatrz.c
+- ext/dlatzm.c
+- ext/dlauu2.c
+- ext/dlauum.c
+- ext/dopgtr.c
+- ext/dopmtr.c
+- ext/dorbdb.c
+- ext/dorcsd.c
+- ext/dorg2l.c
+- ext/dorg2r.c
+- ext/dorgbr.c
+- ext/dorghr.c
+- ext/dorgl2.c
+- ext/dorglq.c
+- ext/dorgql.c
+- ext/dorgqr.c
+- ext/dorgr2.c
+- ext/dorgrq.c
+- ext/dorgtr.c
+- ext/dorm2l.c
+- ext/dorm2r.c
+- ext/dormbr.c
+- ext/dormhr.c
+- ext/dorml2.c
+- ext/dormlq.c
+- ext/dormql.c
+- ext/dormqr.c
+- ext/dormr2.c
+- ext/dormr3.c
+- ext/dormrq.c
+- ext/dormrz.c
+- ext/dormtr.c
+- ext/dpbcon.c
+- ext/dpbequ.c
+- ext/dpbrfs.c
+- ext/dpbstf.c
+- ext/dpbsv.c
+- ext/dpbsvx.c
+- ext/dpbtf2.c
+- ext/dpbtrf.c
+- ext/dpbtrs.c
+- ext/dpftrf.c
+- ext/dpftri.c
+- ext/dpftrs.c
+- ext/dpocon.c
+- ext/dpoequ.c
+- ext/dpoequb.c
+- ext/dporfs.c
+- ext/dporfsx.c
+- ext/dposv.c
+- ext/dposvx.c
+- ext/dposvxx.c
+- ext/dpotf2.c
+- ext/dpotrf.c
+- ext/dpotri.c
+- ext/dpotrs.c
+- ext/dppcon.c
+- ext/dppequ.c
+- ext/dpprfs.c
+- ext/dppsv.c
+- ext/dppsvx.c
+- ext/dpptrf.c
+- ext/dpptri.c
+- ext/dpptrs.c
+- ext/dpstf2.c
+- ext/dpstrf.c
+- ext/dptcon.c
+- ext/dpteqr.c
+- ext/dptrfs.c
+- ext/dptsv.c
+- ext/dptsvx.c
+- ext/dpttrf.c
+- ext/dpttrs.c
+- ext/dptts2.c
+- ext/drscl.c
+- ext/dsbev.c
+- ext/dsbevd.c
+- ext/dsbevx.c
+- ext/dsbgst.c
+- ext/dsbgv.c
+- ext/dsbgvd.c
+- ext/dsbgvx.c
+- ext/dsbtrd.c
+- ext/dsfrk.c
+- ext/dsgesv.c
+- ext/dspcon.c
+- ext/dspev.c
+- ext/dspevd.c
+- ext/dspevx.c
+- ext/dspgst.c
+- ext/dspgv.c
+- ext/dspgvd.c
+- ext/dspgvx.c
+- ext/dsposv.c
+- ext/dsprfs.c
+- ext/dspsv.c
+- ext/dspsvx.c
+- ext/dsptrd.c
+- ext/dsptrf.c
+- ext/dsptri.c
+- ext/dsptrs.c
+- ext/dstebz.c
+- ext/dstedc.c
+- ext/dstegr.c
+- ext/dstein.c
+- ext/dstemr.c
+- ext/dsteqr.c
+- ext/dsterf.c
+- ext/dstev.c
+- ext/dstevd.c
+- ext/dstevr.c
+- ext/dstevx.c
+- ext/dsycon.c
+- ext/dsyconv.c
+- ext/dsyequb.c
+- ext/dsyev.c
+- ext/dsyevd.c
+- ext/dsyevr.c
+- ext/dsyevx.c
+- ext/dsygs2.c
+- ext/dsygst.c
+- ext/dsygv.c
+- ext/dsygvd.c
+- ext/dsygvx.c
+- ext/dsyrfs.c
+- ext/dsyrfsx.c
+- ext/dsysv.c
+- ext/dsysvx.c
+- ext/dsysvxx.c
+- ext/dsyswapr.c
+- ext/dsytd2.c
+- ext/dsytf2.c
+- ext/dsytrd.c
+- ext/dsytrf.c
+- ext/dsytri.c
+- ext/dsytri2.c
+- ext/dsytri2x.c
+- ext/dsytrs.c
+- ext/dsytrs2.c
+- ext/dtbcon.c
+- ext/dtbrfs.c
+- ext/dtbtrs.c
+- ext/dtfsm.c
+- ext/dtftri.c
+- ext/dtfttp.c
+- ext/dtfttr.c
+- ext/dtgevc.c
+- ext/dtgex2.c
+- ext/dtgexc.c
+- ext/dtgsen.c
+- ext/dtgsja.c
+- ext/dtgsna.c
+- ext/dtgsy2.c
+- ext/dtgsyl.c
+- ext/dtpcon.c
+- ext/dtprfs.c
+- ext/dtptri.c
+- ext/dtptrs.c
+- ext/dtpttf.c
+- ext/dtpttr.c
+- ext/dtrcon.c
+- ext/dtrevc.c
+- ext/dtrexc.c
+- ext/dtrrfs.c
+- ext/dtrsen.c
+- ext/dtrsna.c
+- ext/dtrsyl.c
+- ext/dtrti2.c
+- ext/dtrtri.c
+- ext/dtrtrs.c
+- ext/dtrttf.c
+- ext/dtrttp.c
+- ext/dtzrqf.c
+- ext/dtzrzf.c
+- ext/dzsum1.c
 - ext/extconf.rb
 - ext/f2c_minimal.h
+- ext/icmax1.c
+- ext/ieeeck.c
+- ext/ilaclc.c
+- ext/ilaclr.c
+- ext/iladiag.c
+- ext/iladlc.c
+- ext/iladlr.c
+- ext/ilaenv.c
+- ext/ilaprec.c
+- ext/ilaslc.c
+- ext/ilaslr.c
+- ext/ilatrans.c
+- ext/ilauplo.c
+- ext/ilaver.c
+- ext/ilazlc.c
+- ext/ilazlr.c
+- ext/iparmq.c
+- ext/izmax1.c
+- ext/lsamen.c
+- ext/rb_lapack.c
 - ext/rb_lapack.h
+- ext/sbbcsd.c
+- ext/sbdsdc.c
+- ext/sbdsqr.c
+- ext/scsum1.c
+- ext/sdisna.c
+- ext/sgbbrd.c
+- ext/sgbcon.c
+- ext/sgbequ.c
+- ext/sgbequb.c
+- ext/sgbrfs.c
+- ext/sgbrfsx.c
+- ext/sgbsv.c
+- ext/sgbsvx.c
+- ext/sgbsvxx.c
+- ext/sgbtf2.c
+- ext/sgbtrf.c
+- ext/sgbtrs.c
+- ext/sgebak.c
+- ext/sgebal.c
+- ext/sgebd2.c
+- ext/sgebrd.c
+- ext/sgecon.c
+- ext/sgeequ.c
+- ext/sgeequb.c
+- ext/sgees.c
+- ext/sgeesx.c
+- ext/sgeev.c
+- ext/sgeevx.c
+- ext/sgegs.c
+- ext/sgegv.c
+- ext/sgehd2.c
+- ext/sgehrd.c
+- ext/sgejsv.c
+- ext/sgelq2.c
+- ext/sgelqf.c
+- ext/sgels.c
+- ext/sgelsd.c
+- ext/sgelss.c
+- ext/sgelsx.c
+- ext/sgelsy.c
+- ext/sgeql2.c
+- ext/sgeqlf.c
+- ext/sgeqp3.c
+- ext/sgeqpf.c
+- ext/sgeqr2.c
+- ext/sgeqr2p.c
+- ext/sgeqrf.c
+- ext/sgeqrfp.c
+- ext/sgerfs.c
+- ext/sgerfsx.c
+- ext/sgerq2.c
+- ext/sgerqf.c
+- ext/sgesc2.c
+- ext/sgesdd.c
+- ext/sgesv.c
+- ext/sgesvd.c
+- ext/sgesvj.c
+- ext/sgesvx.c
+- ext/sgesvxx.c
+- ext/sgetc2.c
+- ext/sgetf2.c
+- ext/sgetrf.c
+- ext/sgetri.c
+- ext/sgetrs.c
+- ext/sggbak.c
+- ext/sggbal.c
+- ext/sgges.c
+- ext/sggesx.c
+- ext/sggev.c
+- ext/sggevx.c
+- ext/sggglm.c
+- ext/sgghrd.c
+- ext/sgglse.c
+- ext/sggqrf.c
+- ext/sggrqf.c
+- ext/sggsvd.c
+- ext/sggsvp.c
+- ext/sgsvj0.c
+- ext/sgsvj1.c
+- ext/sgtcon.c
+- ext/sgtrfs.c
+- ext/sgtsv.c
+- ext/sgtsvx.c
+- ext/sgttrf.c
+- ext/sgttrs.c
+- ext/sgtts2.c
+- ext/shgeqz.c
+- ext/shsein.c
+- ext/shseqr.c
+- ext/sisnan.c
+- ext/sla_gbamv.c
+- ext/sla_gbrcond.c
+- ext/sla_gbrfsx_extended.c
+- ext/sla_gbrpvgrw.c
+- ext/sla_geamv.c
+- ext/sla_gercond.c
+- ext/sla_gerfsx_extended.c
+- ext/sla_lin_berr.c
+- ext/sla_porcond.c
+- ext/sla_porfsx_extended.c
+- ext/sla_porpvgrw.c
+- ext/sla_rpvgrw.c
+- ext/sla_syamv.c
+- ext/sla_syrcond.c
+- ext/sla_syrfsx_extended.c
+- ext/sla_syrpvgrw.c
+- ext/sla_wwaddw.c
+- ext/slabad.c
+- ext/slabrd.c
+- ext/slacn2.c
+- ext/slacon.c
+- ext/slacpy.c
+- ext/sladiv.c
+- ext/slae2.c
+- ext/slaebz.c
+- ext/slaed0.c
+- ext/slaed1.c
+- ext/slaed2.c
+- ext/slaed3.c
+- ext/slaed4.c
+- ext/slaed5.c
+- ext/slaed6.c
+- ext/slaed7.c
+- ext/slaed8.c
+- ext/slaed9.c
+- ext/slaeda.c
+- ext/slaein.c
+- ext/slaev2.c
+- ext/slaexc.c
+- ext/slag2.c
+- ext/slag2d.c
+- ext/slags2.c
+- ext/slagtf.c
+- ext/slagtm.c
+- ext/slagts.c
+- ext/slagv2.c
+- ext/slahqr.c
+- ext/slahr2.c
+- ext/slahrd.c
+- ext/slaic1.c
+- ext/slaln2.c
+- ext/slals0.c
+- ext/slalsa.c
+- ext/slalsd.c
+- ext/slamrg.c
+- ext/slaneg.c
+- ext/slangb.c
+- ext/slange.c
+- ext/slangt.c
+- ext/slanhs.c
+- ext/slansb.c
+- ext/slansf.c
+- ext/slansp.c
+- ext/slanst.c
+- ext/slansy.c
+- ext/slantb.c
+- ext/slantp.c
+- ext/slantr.c
+- ext/slanv2.c
+- ext/slapll.c
+- ext/slapmr.c
+- ext/slapmt.c
+- ext/slapy2.c
+- ext/slapy3.c
+- ext/slaqgb.c
+- ext/slaqge.c
+- ext/slaqp2.c
+- ext/slaqps.c
+- ext/slaqr0.c
+- ext/slaqr1.c
+- ext/slaqr2.c
+- ext/slaqr3.c
+- ext/slaqr4.c
+- ext/slaqr5.c
+- ext/slaqsb.c
+- ext/slaqsp.c
+- ext/slaqsy.c
+- ext/slaqtr.c
+- ext/slar1v.c
+- ext/slar2v.c
+- ext/slarf.c
+- ext/slarfb.c
+- ext/slarfg.c
+- ext/slarfgp.c
+- ext/slarft.c
+- ext/slarfx.c
+- ext/slargv.c
+- ext/slarnv.c
+- ext/slarra.c
+- ext/slarrb.c
+- ext/slarrc.c
+- ext/slarrd.c
+- ext/slarre.c
+- ext/slarrf.c
+- ext/slarrj.c
+- ext/slarrk.c
+- ext/slarrr.c
+- ext/slarrv.c
+- ext/slarscl2.c
+- ext/slartg.c
+- ext/slartgp.c
+- ext/slartgs.c
+- ext/slartv.c
+- ext/slaruv.c
+- ext/slarz.c
+- ext/slarzb.c
+- ext/slarzt.c
+- ext/slas2.c
+- ext/slascl.c
+- ext/slascl2.c
+- ext/slasd0.c
+- ext/slasd1.c
+- ext/slasd2.c
+- ext/slasd3.c
+- ext/slasd4.c
+- ext/slasd5.c
+- ext/slasd6.c
+- ext/slasd7.c
+- ext/slasd8.c
+- ext/slasda.c
+- ext/slasdq.c
+- ext/slasdt.c
+- ext/slaset.c
+- ext/slasq1.c
+- ext/slasq2.c
+- ext/slasq3.c
+- ext/slasq4.c
+- ext/slasq5.c
+- ext/slasq6.c
+- ext/slasr.c
+- ext/slasrt.c
+- ext/slassq.c
+- ext/slasv2.c
+- ext/slaswp.c
+- ext/slasy2.c
+- ext/slasyf.c
+- ext/slatbs.c
+- ext/slatdf.c
+- ext/slatps.c
+- ext/slatrd.c
+- ext/slatrs.c
+- ext/slatrz.c
+- ext/slatzm.c
+- ext/slauu2.c
+- ext/slauum.c
+- ext/sopgtr.c
+- ext/sopmtr.c
+- ext/sorbdb.c
+- ext/sorcsd.c
+- ext/sorg2l.c
+- ext/sorg2r.c
+- ext/sorgbr.c
+- ext/sorghr.c
+- ext/sorgl2.c
+- ext/sorglq.c
+- ext/sorgql.c
+- ext/sorgqr.c
+- ext/sorgr2.c
+- ext/sorgrq.c
+- ext/sorgtr.c
+- ext/sorm2l.c
+- ext/sorm2r.c
+- ext/sormbr.c
+- ext/sormhr.c
+- ext/sorml2.c
+- ext/sormlq.c
+- ext/sormql.c
+- ext/sormqr.c
+- ext/sormr2.c
+- ext/sormr3.c
+- ext/sormrq.c
+- ext/sormrz.c
+- ext/sormtr.c
+- ext/spbcon.c
+- ext/spbequ.c
+- ext/spbrfs.c
+- ext/spbstf.c
+- ext/spbsv.c
+- ext/spbsvx.c
+- ext/spbtf2.c
+- ext/spbtrf.c
+- ext/spbtrs.c
+- ext/spftrf.c
+- ext/spftri.c
+- ext/spftrs.c
+- ext/spocon.c
+- ext/spoequ.c
+- ext/spoequb.c
+- ext/sporfs.c
+- ext/sporfsx.c
+- ext/sposv.c
+- ext/sposvx.c
+- ext/sposvxx.c
+- ext/spotf2.c
+- ext/spotrf.c
+- ext/spotri.c
+- ext/spotrs.c
+- ext/sppcon.c
+- ext/sppequ.c
+- ext/spprfs.c
+- ext/sppsv.c
+- ext/sppsvx.c
+- ext/spptrf.c
+- ext/spptri.c
+- ext/spptrs.c
+- ext/spstf2.c
+- ext/spstrf.c
+- ext/sptcon.c
+- ext/spteqr.c
+- ext/sptrfs.c
+- ext/sptsv.c
+- ext/sptsvx.c
+- ext/spttrf.c
+- ext/spttrs.c
+- ext/sptts2.c
+- ext/srscl.c
+- ext/ssbev.c
+- ext/ssbevd.c
+- ext/ssbevx.c
+- ext/ssbgst.c
+- ext/ssbgv.c
+- ext/ssbgvd.c
+- ext/ssbgvx.c
+- ext/ssbtrd.c
+- ext/ssfrk.c
+- ext/sspcon.c
+- ext/sspev.c
+- ext/sspevd.c
+- ext/sspevx.c
+- ext/sspgst.c
+- ext/sspgv.c
+- ext/sspgvd.c
+- ext/sspgvx.c
+- ext/ssprfs.c
+- ext/sspsv.c
+- ext/sspsvx.c
+- ext/ssptrd.c
+- ext/ssptrf.c
+- ext/ssptri.c
+- ext/ssptrs.c
+- ext/sstebz.c
+- ext/sstedc.c
+- ext/sstegr.c
+- ext/sstein.c
+- ext/sstemr.c
+- ext/ssteqr.c
+- ext/ssterf.c
+- ext/sstev.c
+- ext/sstevd.c
+- ext/sstevr.c
+- ext/sstevx.c
+- ext/ssycon.c
+- ext/ssyconv.c
+- ext/ssyequb.c
+- ext/ssyev.c
+- ext/ssyevd.c
+- ext/ssyevr.c
+- ext/ssyevx.c
+- ext/ssygs2.c
+- ext/ssygst.c
+- ext/ssygv.c
+- ext/ssygvd.c
+- ext/ssygvx.c
+- ext/ssyrfs.c
+- ext/ssyrfsx.c
+- ext/ssysv.c
+- ext/ssysvx.c
+- ext/ssysvxx.c
+- ext/ssyswapr.c
+- ext/ssytd2.c
+- ext/ssytf2.c
+- ext/ssytrd.c
+- ext/ssytrf.c
+- ext/ssytri.c
+- ext/ssytri2.c
+- ext/ssytri2x.c
+- ext/ssytrs.c
+- ext/ssytrs2.c
+- ext/stbcon.c
+- ext/stbrfs.c
+- ext/stbtrs.c
+- ext/stfsm.c
+- ext/stftri.c
+- ext/stfttp.c
+- ext/stfttr.c
+- ext/stgevc.c
+- ext/stgex2.c
+- ext/stgexc.c
+- ext/stgsen.c
+- ext/stgsja.c
+- ext/stgsna.c
+- ext/stgsy2.c
+- ext/stgsyl.c
+- ext/stpcon.c
+- ext/stprfs.c
+- ext/stptri.c
+- ext/stptrs.c
+- ext/stpttf.c
+- ext/stpttr.c
+- ext/strcon.c
+- ext/strevc.c
+- ext/strexc.c
+- ext/strrfs.c
+- ext/strsen.c
+- ext/strsna.c
+- ext/strsyl.c
+- ext/strti2.c
+- ext/strtri.c
+- ext/strtrs.c
+- ext/strttf.c
+- ext/strttp.c
+- ext/stzrqf.c
+- ext/stzrzf.c
+- ext/xerbla.c
+- ext/xerbla_array.c
+- ext/zbbcsd.c
+- ext/zbdsqr.c
+- ext/zcgesv.c
+- ext/zcposv.c
+- ext/zdrscl.c
+- ext/zgbbrd.c
+- ext/zgbcon.c
+- ext/zgbequ.c
+- ext/zgbequb.c
+- ext/zgbrfs.c
+- ext/zgbrfsx.c
+- ext/zgbsv.c
+- ext/zgbsvx.c
+- ext/zgbsvxx.c
+- ext/zgbtf2.c
+- ext/zgbtrf.c
+- ext/zgbtrs.c
+- ext/zgebak.c
+- ext/zgebal.c
+- ext/zgebd2.c
+- ext/zgebrd.c
+- ext/zgecon.c
+- ext/zgeequ.c
+- ext/zgeequb.c
+- ext/zgees.c
+- ext/zgeesx.c
+- ext/zgeev.c
+- ext/zgeevx.c
+- ext/zgegs.c
+- ext/zgegv.c
+- ext/zgehd2.c
+- ext/zgehrd.c
+- ext/zgelq2.c
+- ext/zgelqf.c
+- ext/zgels.c
+- ext/zgelsd.c
+- ext/zgelss.c
+- ext/zgelsx.c
+- ext/zgelsy.c
+- ext/zgeql2.c
+- ext/zgeqlf.c
+- ext/zgeqp3.c
+- ext/zgeqpf.c
+- ext/zgeqr2.c
+- ext/zgeqr2p.c
+- ext/zgeqrf.c
+- ext/zgeqrfp.c
+- ext/zgerfs.c
+- ext/zgerfsx.c
+- ext/zgerq2.c
+- ext/zgerqf.c
+- ext/zgesc2.c
+- ext/zgesdd.c
+- ext/zgesv.c
+- ext/zgesvd.c
+- ext/zgesvx.c
+- ext/zgesvxx.c
+- ext/zgetc2.c
+- ext/zgetf2.c
+- ext/zgetrf.c
+- ext/zgetri.c
+- ext/zgetrs.c
+- ext/zggbak.c
+- ext/zggbal.c
+- ext/zgges.c
+- ext/zggesx.c
+- ext/zggev.c
+- ext/zggevx.c
+- ext/zggglm.c
+- ext/zgghrd.c
+- ext/zgglse.c
+- ext/zggqrf.c
+- ext/zggrqf.c
+- ext/zggsvd.c
+- ext/zggsvp.c
+- ext/zgtcon.c
+- ext/zgtrfs.c
+- ext/zgtsv.c
+- ext/zgtsvx.c
+- ext/zgttrf.c
+- ext/zgttrs.c
+- ext/zgtts2.c
+- ext/zhbev.c
+- ext/zhbevd.c
+- ext/zhbevx.c
+- ext/zhbgst.c
+- ext/zhbgv.c
+- ext/zhbgvd.c
+- ext/zhbgvx.c
+- ext/zhbtrd.c
+- ext/zhecon.c
+- ext/zheequb.c
+- ext/zheev.c
+- ext/zheevd.c
+- ext/zheevr.c
+- ext/zheevx.c
+- ext/zhegs2.c
+- ext/zhegst.c
+- ext/zhegv.c
+- ext/zhegvd.c
+- ext/zhegvx.c
+- ext/zherfs.c
+- ext/zherfsx.c
+- ext/zhesv.c
+- ext/zhesvx.c
+- ext/zhesvxx.c
+- ext/zhetd2.c
+- ext/zhetf2.c
+- ext/zhetrd.c
+- ext/zhetrf.c
+- ext/zhetri.c
+- ext/zhetrs.c
+- ext/zhetrs2.c
+- ext/zhfrk.c
+- ext/zhgeqz.c
+- ext/zhpcon.c
+- ext/zhpev.c
+- ext/zhpevd.c
+- ext/zhpevx.c
+- ext/zhpgst.c
+- ext/zhpgv.c
+- ext/zhpgvd.c
+- ext/zhpgvx.c
+- ext/zhprfs.c
+- ext/zhpsv.c
+- ext/zhpsvx.c
+- ext/zhptrd.c
+- ext/zhptrf.c
+- ext/zhptri.c
+- ext/zhptrs.c
+- ext/zhsein.c
+- ext/zhseqr.c
+- ext/zla_gbamv.c
+- ext/zla_gbrcond_c.c
+- ext/zla_gbrcond_x.c
+- ext/zla_gbrfsx_extended.c
+- ext/zla_gbrpvgrw.c
+- ext/zla_geamv.c
+- ext/zla_gercond_c.c
+- ext/zla_gercond_x.c
+- ext/zla_gerfsx_extended.c
+- ext/zla_heamv.c
+- ext/zla_hercond_c.c
+- ext/zla_hercond_x.c
+- ext/zla_herfsx_extended.c
+- ext/zla_herpvgrw.c
+- ext/zla_lin_berr.c
+- ext/zla_porcond_c.c
+- ext/zla_porcond_x.c
+- ext/zla_porfsx_extended.c
+- ext/zla_porpvgrw.c
+- ext/zla_rpvgrw.c
+- ext/zla_syamv.c
+- ext/zla_syrcond_c.c
+- ext/zla_syrcond_x.c
+- ext/zla_syrfsx_extended.c
+- ext/zla_syrpvgrw.c
+- ext/zla_wwaddw.c
+- ext/zlabrd.c
+- ext/zlacgv.c
+- ext/zlacn2.c
+- ext/zlacon.c
+- ext/zlacp2.c
+- ext/zlacpy.c
+- ext/zlacrm.c
+- ext/zlacrt.c
+- ext/zladiv.c
+- ext/zlaed0.c
+- ext/zlaed7.c
+- ext/zlaed8.c
+- ext/zlaein.c
+- ext/zlaesy.c
+- ext/zlaev2.c
+- ext/zlag2c.c
+- ext/zlags2.c
+- ext/zlagtm.c
+- ext/zlahef.c
+- ext/zlahqr.c
+- ext/zlahr2.c
+- ext/zlahrd.c
+- ext/zlaic1.c
+- ext/zlals0.c
+- ext/zlalsa.c
+- ext/zlalsd.c
+- ext/zlangb.c
+- ext/zlange.c
+- ext/zlangt.c
+- ext/zlanhb.c
+- ext/zlanhe.c
+- ext/zlanhf.c
+- ext/zlanhp.c
+- ext/zlanhs.c
+- ext/zlanht.c
+- ext/zlansb.c
+- ext/zlansp.c
+- ext/zlansy.c
+- ext/zlantb.c
+- ext/zlantp.c
+- ext/zlantr.c
+- ext/zlapll.c
+- ext/zlapmr.c
+- ext/zlapmt.c
+- ext/zlaqgb.c
+- ext/zlaqge.c
+- ext/zlaqhb.c
+- ext/zlaqhe.c
+- ext/zlaqhp.c
+- ext/zlaqp2.c
+- ext/zlaqps.c
+- ext/zlaqr0.c
+- ext/zlaqr1.c
+- ext/zlaqr2.c
+- ext/zlaqr3.c
+- ext/zlaqr4.c
+- ext/zlaqr5.c
+- ext/zlaqsb.c
+- ext/zlaqsp.c
+- ext/zlaqsy.c
+- ext/zlar1v.c
+- ext/zlar2v.c
+- ext/zlarcm.c
+- ext/zlarf.c
+- ext/zlarfb.c
+- ext/zlarfg.c
+- ext/zlarfgp.c
+- ext/zlarft.c
+- ext/zlarfx.c
+- ext/zlargv.c
+- ext/zlarnv.c
+- ext/zlarrv.c
+- ext/zlarscl2.c
+- ext/zlartg.c
+- ext/zlartv.c
+- ext/zlarz.c
+- ext/zlarzb.c
+- ext/zlarzt.c
+- ext/zlascl.c
+- ext/zlascl2.c
+- ext/zlaset.c
+- ext/zlasr.c
+- ext/zlassq.c
+- ext/zlaswp.c
+- ext/zlasyf.c
+- ext/zlat2c.c
+- ext/zlatbs.c
+- ext/zlatdf.c
+- ext/zlatps.c
+- ext/zlatrd.c
+- ext/zlatrs.c
+- ext/zlatrz.c
+- ext/zlatzm.c
+- ext/zlauu2.c
+- ext/zlauum.c
+- ext/zpbcon.c
+- ext/zpbequ.c
+- ext/zpbrfs.c
+- ext/zpbstf.c
+- ext/zpbsv.c
+- ext/zpbsvx.c
+- ext/zpbtf2.c
+- ext/zpbtrf.c
+- ext/zpbtrs.c
+- ext/zpftrf.c
+- ext/zpftri.c
+- ext/zpftrs.c
+- ext/zpocon.c
+- ext/zpoequ.c
+- ext/zpoequb.c
+- ext/zporfs.c
+- ext/zporfsx.c
+- ext/zposv.c
+- ext/zposvx.c
+- ext/zposvxx.c
+- ext/zpotf2.c
+- ext/zpotrf.c
+- ext/zpotri.c
+- ext/zpotrs.c
+- ext/zppcon.c
+- ext/zppequ.c
+- ext/zpprfs.c
+- ext/zppsv.c
+- ext/zppsvx.c
+- ext/zpptrf.c
+- ext/zpptri.c
+- ext/zpptrs.c
+- ext/zpstf2.c
+- ext/zpstrf.c
+- ext/zptcon.c
+- ext/zpteqr.c
+- ext/zptrfs.c
+- ext/zptsv.c
+- ext/zptsvx.c
+- ext/zpttrf.c
+- ext/zpttrs.c
+- ext/zptts2.c
+- ext/zrot.c
+- ext/zspcon.c
+- ext/zspmv.c
+- ext/zspr.c
+- ext/zsprfs.c
+- ext/zspsv.c
+- ext/zspsvx.c
+- ext/zsptrf.c
+- ext/zsptri.c
+- ext/zsptrs.c
+- ext/zstedc.c
+- ext/zstegr.c
+- ext/zstein.c
+- ext/zstemr.c
+- ext/zsteqr.c
+- ext/zsycon.c
+- ext/zsyconv.c
+- ext/zsyequb.c
+- ext/zsymv.c
+- ext/zsyr.c
+- ext/zsyrfs.c
+- ext/zsyrfsx.c
+- ext/zsysv.c
+- ext/zsysvx.c
+- ext/zsysvxx.c
+- ext/zsyswapr.c
+- ext/zsytf2.c
+- ext/zsytrf.c
+- ext/zsytri.c
+- ext/zsytri2.c
+- ext/zsytri2x.c
+- ext/zsytrs.c
+- ext/zsytrs2.c
+- ext/ztbcon.c
+- ext/ztbrfs.c
+- ext/ztbtrs.c
+- ext/ztfsm.c
+- ext/ztftri.c
+- ext/ztfttp.c
+- ext/ztfttr.c
+- ext/ztgevc.c
+- ext/ztgex2.c
+- ext/ztgexc.c
+- ext/ztgsen.c
+- ext/ztgsja.c
+- ext/ztgsna.c
+- ext/ztgsy2.c
+- ext/ztgsyl.c
+- ext/ztpcon.c
+- ext/ztprfs.c
+- ext/ztptri.c
+- ext/ztptrs.c
+- ext/ztpttf.c
+- ext/ztpttr.c
+- ext/ztrcon.c
+- ext/ztrevc.c
+- ext/ztrexc.c
+- ext/ztrrfs.c
+- ext/ztrsen.c
+- ext/ztrsna.c
+- ext/ztrsyl.c
+- ext/ztrti2.c
+- ext/ztrtri.c
+- ext/ztrtrs.c
+- ext/ztrttf.c
+- ext/ztrttp.c
+- ext/ztzrqf.c
+- ext/ztzrzf.c
+- ext/zunbdb.c
+- ext/zuncsd.c
+- ext/zung2l.c
+- ext/zung2r.c
+- ext/zungbr.c
+- ext/zunghr.c
+- ext/zungl2.c
+- ext/zunglq.c
+- ext/zungql.c
+- ext/zungqr.c
+- ext/zungr2.c
+- ext/zungrq.c
+- ext/zungtr.c
+- ext/zunm2l.c
+- ext/zunm2r.c
+- ext/zunmbr.c
+- ext/zunmhr.c
+- ext/zunml2.c
+- ext/zunmlq.c
+- ext/zunmql.c
+- ext/zunmqr.c
+- ext/zunmr2.c
+- ext/zunmr3.c
+- ext/zunmrq.c
+- ext/zunmrz.c
+- ext/zunmtr.c
+- ext/zupgtr.c
+- ext/zupmtr.c
 - lib/numru/lapack.rb
 - samples/dpotrf.rb
 - samples/dsyevr.rb