ruby-lapack 1.4.1a → 1.5

data/ext/slatrz.c

@@ -42,6 +42,7 @@ rblapack_slatrz(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  l = NUM2INT(rblapack_l);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -52,7 +53,6 @@ rblapack_slatrz(int argc, VALUE *argv, VALUE self){
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
   m = lda;
-  l = NUM2INT(rblapack_l);
   {
     int shape[1];
     shape[0] = m;

data/ext/slatzm.c

@@ -61,7 +61,6 @@ rblapack_slatzm(int argc, VALUE *argv, VALUE self){
   side = StringValueCStr(rblapack_side)[0];
   n = NUM2INT(rblapack_n);
   incv = NUM2INT(rblapack_incv);
-  m = NUM2INT(rblapack_m);
   if (!NA_IsNArray(rblapack_c2))
     rb_raise(rb_eArgError, "c2 (8th argument) must be NArray");
   if (NA_RANK(rblapack_c2) != 2)
@@ -72,6 +71,8 @@ rblapack_slatzm(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_c2) != NA_SFLOAT)
     rblapack_c2 = na_change_type(rblapack_c2, NA_SFLOAT);
   c2 = NA_PTR_TYPE(rblapack_c2, real*);
+  m = NUM2INT(rblapack_m);
+  tau = (real)NUM2DBL(rblapack_tau);
   if (!NA_IsNArray(rblapack_v))
     rb_raise(rb_eArgError, "v (4th argument) must be NArray");
   if (NA_RANK(rblapack_v) != 1)
@@ -81,7 +82,6 @@ rblapack_slatzm(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_v) != NA_SFLOAT)
     rblapack_v = na_change_type(rblapack_v, NA_SFLOAT);
   v = NA_PTR_TYPE(rblapack_v, real*);
-  tau = (real)NUM2DBL(rblapack_tau);
   if (!NA_IsNArray(rblapack_c1))
     rb_raise(rb_eArgError, "c1 (7th argument) must be NArray");
   if (NA_RANK(rblapack_c1) != 2)

data/ext/slauu2.c

@@ -40,6 +40,7 @@ rblapack_slauu2(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -49,7 +50,6 @@ rblapack_slauu2(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
   {
     int shape[2];
     shape[0] = lda;

data/ext/slauum.c

@@ -40,6 +40,7 @@ rblapack_slauum(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -49,7 +50,6 @@ rblapack_slauum(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
   {
     int shape[2];
     shape[0] = lda;

data/ext/sopgtr.c

@@ -46,14 +46,6 @@ rblapack_sopgtr(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  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_SFLOAT)
-    rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
-  ap = NA_PTR_TYPE(rblapack_ap, real*);
   uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_tau))
     rb_raise(rb_eArgError, "tau (3th argument) must be NArray");
@@ -64,6 +56,14 @@ rblapack_sopgtr(int argc, VALUE *argv, VALUE self){
     rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
   tau = NA_PTR_TYPE(rblapack_tau, real*);
   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_SFLOAT)
+    rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
+  ap = NA_PTR_TYPE(rblapack_ap, real*);
   ldq = MAX(1,n);
   {
     int shape[2];

data/ext/sopmtr.c

@@ -57,6 +57,7 @@ rblapack_sopmtr(int argc, VALUE *argv, VALUE self){
   }
 
   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)
@@ -66,7 +67,7 @@ rblapack_sopmtr(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_c) != NA_SFLOAT)
     rblapack_c = na_change_type(rblapack_c, NA_SFLOAT);
   c = NA_PTR_TYPE(rblapack_c, real*);
-  trans = StringValueCStr(rblapack_trans)[0];
+  uplo = StringValueCStr(rblapack_uplo)[0];
   m = NUM2INT(rblapack_m);
   if (!NA_IsNArray(rblapack_tau))
     rb_raise(rb_eArgError, "tau (6th argument) must be NArray");
@@ -86,7 +87,6 @@ rblapack_sopmtr(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_ap) != NA_SFLOAT)
     rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
   ap = NA_PTR_TYPE(rblapack_ap, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
   {
     int shape[2];
     shape[0] = ldc;

data/ext/sorbdb.c

@@ -83,6 +83,9 @@ rblapack_sorbdb(int argc, VALUE *argv, VALUE self){
     rblapack_lwork = Qnil;
   }
 
+  trans = StringValueCStr(rblapack_trans)[0];
+  m = NUM2INT(rblapack_m);
+  signs = StringValueCStr(rblapack_signs)[0];
   if (!NA_IsNArray(rblapack_x11))
     rb_raise(rb_eArgError, "x11 (4th argument) must be NArray");
   if (NA_RANK(rblapack_x11) != 2)
@@ -92,12 +95,8 @@ rblapack_sorbdb(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x11) != NA_SFLOAT)
     rblapack_x11 = na_change_type(rblapack_x11, NA_SFLOAT);
   x11 = NA_PTR_TYPE(rblapack_x11, real*);
-  signs = StringValueCStr(rblapack_signs)[0];
-  trans = StringValueCStr(rblapack_trans)[0];
-  m = NUM2INT(rblapack_m);
   p = ldx11;
   ldx21 = p;
-  ldx22 = p;
   if (!NA_IsNArray(rblapack_x21))
     rb_raise(rb_eArgError, "x21 (6th argument) must be NArray");
   if (NA_RANK(rblapack_x21) != 2)
@@ -114,6 +113,18 @@ rblapack_sorbdb(int argc, VALUE *argv, VALUE self){
   else {
     lwork = NUM2INT(rblapack_lwork);
   }
+  ldx22 = p;
+  if (!NA_IsNArray(rblapack_x22))
+    rb_raise(rb_eArgError, "x22 (7th argument) must be NArray");
+  if (NA_RANK(rblapack_x22) != 2)
+    rb_raise(rb_eArgError, "rank of x22 (7th argument) must be %d", 2);
+  if (NA_SHAPE0(rblapack_x22) != ldx22)
+    rb_raise(rb_eRuntimeError, "shape 0 of x22 must be p");
+  if (NA_SHAPE1(rblapack_x22) != (m-q))
+    rb_raise(rb_eRuntimeError, "shape 1 of x22 must be %d", m-q);
+  if (NA_TYPE(rblapack_x22) != NA_SFLOAT)
+    rblapack_x22 = na_change_type(rblapack_x22, NA_SFLOAT);
+  x22 = NA_PTR_TYPE(rblapack_x22, real*);
   ldx12 = p;
   if (!NA_IsNArray(rblapack_x12))
     rb_raise(rb_eArgError, "x12 (5th argument) must be NArray");
@@ -126,17 +137,6 @@ rblapack_sorbdb(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x12) != NA_SFLOAT)
     rblapack_x12 = na_change_type(rblapack_x12, NA_SFLOAT);
   x12 = NA_PTR_TYPE(rblapack_x12, real*);
-  if (!NA_IsNArray(rblapack_x22))
-    rb_raise(rb_eArgError, "x22 (7th argument) must be NArray");
-  if (NA_RANK(rblapack_x22) != 2)
-    rb_raise(rb_eArgError, "rank of x22 (7th argument) must be %d", 2);
-  if (NA_SHAPE0(rblapack_x22) != ldx22)
-    rb_raise(rb_eRuntimeError, "shape 0 of x22 must be p");
-  if (NA_SHAPE1(rblapack_x22) != (m-q))
-    rb_raise(rb_eRuntimeError, "shape 1 of x22 must be %d", m-q);
-  if (NA_TYPE(rblapack_x22) != NA_SFLOAT)
-    rblapack_x22 = na_change_type(rblapack_x22, NA_SFLOAT);
-  x22 = NA_PTR_TYPE(rblapack_x22, real*);
   {
     int shape[1];
     shape[0] = q;

data/ext/sorcsd.c

@@ -89,6 +89,13 @@ rblapack_sorcsd(int argc, VALUE *argv, VALUE self){
   }
 
   jobu1 = StringValueCStr(rblapack_jobu1)[0];
+  jobv1t = StringValueCStr(rblapack_jobv1t)[0];
+  trans = StringValueCStr(rblapack_trans)[0];
+  m = NUM2INT(rblapack_m);
+  lwork = NUM2INT(rblapack_lwork);
+  jobu2 = StringValueCStr(rblapack_jobu2)[0];
+  signs = StringValueCStr(rblapack_signs)[0];
+  jobv2t = StringValueCStr(rblapack_jobv2t)[0];
   if (!NA_IsNArray(rblapack_x11))
     rb_raise(rb_eArgError, "x11 (8th argument) must be NArray");
   if (NA_RANK(rblapack_x11) != 2)
@@ -98,19 +105,8 @@ rblapack_sorcsd(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x11) != NA_SFLOAT)
     rblapack_x11 = na_change_type(rblapack_x11, NA_SFLOAT);
   x11 = NA_PTR_TYPE(rblapack_x11, real*);
-  jobu2 = StringValueCStr(rblapack_jobu2)[0];
-  signs = StringValueCStr(rblapack_signs)[0];
-  lwork = NUM2INT(rblapack_lwork);
-  jobv1t = StringValueCStr(rblapack_jobv1t)[0];
-  trans = StringValueCStr(rblapack_trans)[0];
-  m = NUM2INT(rblapack_m);
   p = ldx11;
   ldx21 = p;
-  ldx22 = p;
-  ldx12 = p;
-  ldu1 = lsame_(&jobu1,"Y") ? MAX(1,p) : 0;
-  ldu2 = lsame_(&jobu2,"Y") ? MAX(1,m-p) : 0;
-  jobv2t = StringValueCStr(rblapack_jobv2t)[0];
   if (!NA_IsNArray(rblapack_x21))
     rb_raise(rb_eArgError, "x21 (10th argument) must be NArray");
   if (NA_RANK(rblapack_x21) != 2)
@@ -122,6 +118,9 @@ rblapack_sorcsd(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x21) != NA_SFLOAT)
     rblapack_x21 = na_change_type(rblapack_x21, NA_SFLOAT);
   x21 = NA_PTR_TYPE(rblapack_x21, real*);
+  ldv2t = lsame_(&jobv2t,"Y") ? MAX(1,m-q) : 0;
+  ldu1 = lsame_(&jobu1,"Y") ? MAX(1,p) : 0;
+  ldx12 = p;
   if (!NA_IsNArray(rblapack_x12))
     rb_raise(rb_eArgError, "x12 (9th argument) must be NArray");
   if (NA_RANK(rblapack_x12) != 2)
@@ -133,7 +132,8 @@ rblapack_sorcsd(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x12) != NA_SFLOAT)
     rblapack_x12 = na_change_type(rblapack_x12, NA_SFLOAT);
   x12 = NA_PTR_TYPE(rblapack_x12, real*);
-  ldv2t = lsame_(&jobv2t,"Y") ? MAX(1,m-q) : 0;
+  ldv1t = lsame_(&jobv1t,"Y") ? MAX(1,q) : 0;
+  ldx22 = p;
   if (!NA_IsNArray(rblapack_x22))
     rb_raise(rb_eArgError, "x22 (11th argument) must be NArray");
   if (NA_RANK(rblapack_x22) != 2)
@@ -145,7 +145,7 @@ rblapack_sorcsd(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x22) != NA_SFLOAT)
     rblapack_x22 = na_change_type(rblapack_x22, NA_SFLOAT);
   x22 = NA_PTR_TYPE(rblapack_x22, real*);
-  ldv1t = lsame_(&jobv1t,"Y") ? MAX(1,q) : 0;
+  ldu2 = lsame_(&jobu2,"Y") ? MAX(1,m-p) : 0;
   {
     int shape[1];
     shape[0] = MIN(MIN(MIN(p,m-p),q),m-q);

data/ext/sorg2l.c

@@ -45,6 +45,15 @@ rblapack_sorg2l(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  m = NUM2INT(rblapack_m);
+  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);
+  k = NA_SHAPE0(rblapack_tau);
+  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
+    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
+  tau = NA_PTR_TYPE(rblapack_tau, real*);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -54,15 +63,6 @@ rblapack_sorg2l(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
-  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);
-  k = NA_SHAPE0(rblapack_tau);
-  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
-    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
-  tau = NA_PTR_TYPE(rblapack_tau, real*);
-  m = NUM2INT(rblapack_m);
   {
     int shape[2];
     shape[0] = lda;

data/ext/sorg2r.c

@@ -45,6 +45,15 @@ rblapack_sorg2r(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  m = NUM2INT(rblapack_m);
+  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);
+  k = NA_SHAPE0(rblapack_tau);
+  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
+    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
+  tau = NA_PTR_TYPE(rblapack_tau, real*);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -54,15 +63,6 @@ rblapack_sorg2r(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
-  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);
-  k = NA_SHAPE0(rblapack_tau);
-  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
-    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
-  tau = NA_PTR_TYPE(rblapack_tau, real*);
-  m = NUM2INT(rblapack_m);
   {
     int shape[2];
     shape[0] = lda;

data/ext/sorgbr.c

@@ -56,6 +56,7 @@ rblapack_sorgbr(int argc, VALUE *argv, VALUE self){
     rblapack_lwork = Qnil;
   }
 
+  vect = StringValueCStr(rblapack_vect)[0];
   k = NUM2INT(rblapack_k);
   m = NUM2INT(rblapack_m);
   if (!NA_IsNArray(rblapack_tau))
@@ -81,7 +82,6 @@ rblapack_sorgbr(int argc, VALUE *argv, VALUE self){
   else {
     lwork = NUM2INT(rblapack_lwork);
   }
-  vect = StringValueCStr(rblapack_vect)[0];
   {
     int shape[1];
     shape[0] = MAX(1,lwork);

data/ext/sorghr.c

@@ -53,6 +53,7 @@ rblapack_sorghr(int argc, VALUE *argv, VALUE self){
     rblapack_lwork = Qnil;
   }
 
+  ilo = NUM2INT(rblapack_ilo);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (3th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -62,6 +63,12 @@ rblapack_sorghr(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
+  ihi = NUM2INT(rblapack_ihi);
+  if (rblapack_lwork == Qnil)
+    lwork = ihi-ilo;
+  else {
+    lwork = NUM2INT(rblapack_lwork);
+  }
   if (!NA_IsNArray(rblapack_tau))
     rb_raise(rb_eArgError, "tau (4th argument) must be NArray");
   if (NA_RANK(rblapack_tau) != 1)
@@ -71,13 +78,6 @@ rblapack_sorghr(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
     rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
   tau = NA_PTR_TYPE(rblapack_tau, real*);
-  ilo = NUM2INT(rblapack_ilo);
-  ihi = NUM2INT(rblapack_ihi);
-  if (rblapack_lwork == Qnil)
-    lwork = ihi-ilo;
-  else {
-    lwork = NUM2INT(rblapack_lwork);
-  }
   {
     int shape[1];
     shape[0] = MAX(1,lwork);

data/ext/sorgl2.c

@@ -52,6 +52,7 @@ rblapack_sorgl2(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
+  m = lda;
   if (!NA_IsNArray(rblapack_tau))
     rb_raise(rb_eArgError, "tau (2th argument) must be NArray");
   if (NA_RANK(rblapack_tau) != 1)
@@ -60,7 +61,6 @@ rblapack_sorgl2(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
     rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
   tau = NA_PTR_TYPE(rblapack_tau, real*);
-  m = lda;
   {
     int shape[2];
     shape[0] = lda;

data/ext/sorglq.c

@@ -51,6 +51,15 @@ rblapack_sorglq(int argc, VALUE *argv, VALUE self){
     rblapack_lwork = Qnil;
   }
 
+  m = NUM2INT(rblapack_m);
+  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);
+  k = NA_SHAPE0(rblapack_tau);
+  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
+    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
+  tau = NA_PTR_TYPE(rblapack_tau, real*);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -60,15 +69,6 @@ rblapack_sorglq(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_a) != NA_SFLOAT)
     rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
   a = NA_PTR_TYPE(rblapack_a, real*);
-  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);
-  k = NA_SHAPE0(rblapack_tau);
-  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
-    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
-  tau = NA_PTR_TYPE(rblapack_tau, real*);
-  m = NUM2INT(rblapack_m);
   if (rblapack_lwork == Qnil)
     lwork = m;
   else {

data/ext/sorgql.c

@@ -51,6 +51,15 @@ rblapack_sorgql(int argc, VALUE *argv, VALUE self){
     rblapack_lwork = Qnil;
   }
 
+  m = NUM2INT(rblapack_m);
+  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);
+  k = NA_SHAPE0(rblapack_tau);
+  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
+    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
+  tau = NA_PTR_TYPE(rblapack_tau, real*);
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (2th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -65,15 +74,6 @@ rblapack_sorgql(int argc, VALUE *argv, VALUE self){
   else {
     lwork = NUM2INT(rblapack_lwork);
   }
-  m = NUM2INT(rblapack_m);
-  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);
-  k = NA_SHAPE0(rblapack_tau);
-  if (NA_TYPE(rblapack_tau) != NA_SFLOAT)
-    rblapack_tau = na_change_type(rblapack_tau, NA_SFLOAT);
-  tau = NA_PTR_TYPE(rblapack_tau, real*);
   {
     int shape[1];
     shape[0] = MAX(1,lwork);