ruby-lapack 1.4.1a → 1.5

data/ext/sgtcon.c

@@ -56,7 +56,7 @@ rblapack_sgtcon(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  anorm = (real)NUM2DBL(rblapack_anorm);
+  norm = StringValueCStr(rblapack_norm)[0];
   if (!NA_IsNArray(rblapack_d))
     rb_raise(rb_eArgError, "d (3th argument) must be NArray");
   if (NA_RANK(rblapack_d) != 1)
@@ -65,7 +65,6 @@ rblapack_sgtcon(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_d) != NA_SFLOAT)
     rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
   d = NA_PTR_TYPE(rblapack_d, real*);
-  norm = StringValueCStr(rblapack_norm)[0];
   if (!NA_IsNArray(rblapack_ipiv))
     rb_raise(rb_eArgError, "ipiv (6th argument) must be NArray");
   if (NA_RANK(rblapack_ipiv) != 1)
@@ -75,24 +74,6 @@ rblapack_sgtcon(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_ipiv) != NA_LINT)
     rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
   ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
-  if (!NA_IsNArray(rblapack_du))
-    rb_raise(rb_eArgError, "du (4th argument) must be NArray");
-  if (NA_RANK(rblapack_du) != 1)
-    rb_raise(rb_eArgError, "rank of du (4th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_du) != (n-1))
-    rb_raise(rb_eRuntimeError, "shape 0 of du must be %d", n-1);
-  if (NA_TYPE(rblapack_du) != NA_SFLOAT)
-    rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
-  du = NA_PTR_TYPE(rblapack_du, real*);
-  if (!NA_IsNArray(rblapack_du2))
-    rb_raise(rb_eArgError, "du2 (5th argument) must be NArray");
-  if (NA_RANK(rblapack_du2) != 1)
-    rb_raise(rb_eArgError, "rank of du2 (5th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_du2) != (n-2))
-    rb_raise(rb_eRuntimeError, "shape 0 of du2 must be %d", n-2);
-  if (NA_TYPE(rblapack_du2) != NA_SFLOAT)
-    rblapack_du2 = na_change_type(rblapack_du2, NA_SFLOAT);
-  du2 = NA_PTR_TYPE(rblapack_du2, real*);
   if (!NA_IsNArray(rblapack_dl))
     rb_raise(rb_eArgError, "dl (2th argument) must be NArray");
   if (NA_RANK(rblapack_dl) != 1)
@@ -102,6 +83,25 @@ rblapack_sgtcon(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_dl) != NA_SFLOAT)
     rblapack_dl = na_change_type(rblapack_dl, NA_SFLOAT);
   dl = NA_PTR_TYPE(rblapack_dl, real*);
+  if (!NA_IsNArray(rblapack_du2))
+    rb_raise(rb_eArgError, "du2 (5th argument) must be NArray");
+  if (NA_RANK(rblapack_du2) != 1)
+    rb_raise(rb_eArgError, "rank of du2 (5th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_du2) != (n-2))
+    rb_raise(rb_eRuntimeError, "shape 0 of du2 must be %d", n-2);
+  if (NA_TYPE(rblapack_du2) != NA_SFLOAT)
+    rblapack_du2 = na_change_type(rblapack_du2, NA_SFLOAT);
+  du2 = NA_PTR_TYPE(rblapack_du2, real*);
+  if (!NA_IsNArray(rblapack_du))
+    rb_raise(rb_eArgError, "du (4th argument) must be NArray");
+  if (NA_RANK(rblapack_du) != 1)
+    rb_raise(rb_eArgError, "rank of du (4th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_du) != (n-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of du must be %d", n-1);
+  if (NA_TYPE(rblapack_du) != NA_SFLOAT)
+    rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
+  du = NA_PTR_TYPE(rblapack_du, real*);
+  anorm = (real)NUM2DBL(rblapack_anorm);
   work = ALLOC_N(real, (2*n));
   iwork = ALLOC_N(integer, (n));
 

data/ext/sgtrfs.c

@@ -75,11 +75,30 @@ rblapack_sgtrfs(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  trans = StringValueCStr(rblapack_trans)[0];
+  if (!NA_IsNArray(rblapack_d))
+    rb_raise(rb_eArgError, "d (3th argument) must be NArray");
+  if (NA_RANK(rblapack_d) != 1)
+    rb_raise(rb_eArgError, "rank of d (3th argument) must be %d", 1);
+  n = NA_SHAPE0(rblapack_d);
+  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
+    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
+  d = NA_PTR_TYPE(rblapack_d, real*);
+  if (!NA_IsNArray(rblapack_df))
+    rb_raise(rb_eArgError, "df (6th argument) must be NArray");
+  if (NA_RANK(rblapack_df) != 1)
+    rb_raise(rb_eArgError, "rank of df (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_df) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of df must be the same as shape 0 of d");
+  if (NA_TYPE(rblapack_df) != NA_SFLOAT)
+    rblapack_df = na_change_type(rblapack_df, NA_SFLOAT);
+  df = NA_PTR_TYPE(rblapack_df, real*);
   if (!NA_IsNArray(rblapack_ipiv))
     rb_raise(rb_eArgError, "ipiv (9th argument) must be NArray");
   if (NA_RANK(rblapack_ipiv) != 1)
     rb_raise(rb_eArgError, "rank of ipiv (9th argument) must be %d", 1);
-  n = NA_SHAPE0(rblapack_ipiv);
+  if (NA_SHAPE0(rblapack_ipiv) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 0 of d");
   if (NA_TYPE(rblapack_ipiv) != NA_LINT)
     rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
   ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
@@ -92,35 +111,15 @@ rblapack_sgtrfs(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_x) != NA_SFLOAT)
     rblapack_x = na_change_type(rblapack_x, NA_SFLOAT);
   x = NA_PTR_TYPE(rblapack_x, real*);
-  if (!NA_IsNArray(rblapack_df))
-    rb_raise(rb_eArgError, "df (6th argument) must be NArray");
-  if (NA_RANK(rblapack_df) != 1)
-    rb_raise(rb_eArgError, "rank of df (6th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_df) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of df must be the same as shape 0 of ipiv");
-  if (NA_TYPE(rblapack_df) != NA_SFLOAT)
-    rblapack_df = na_change_type(rblapack_df, NA_SFLOAT);
-  df = NA_PTR_TYPE(rblapack_df, real*);
-  trans = StringValueCStr(rblapack_trans)[0];
-  if (!NA_IsNArray(rblapack_b))
-    rb_raise(rb_eArgError, "b (10th argument) must be NArray");
-  if (NA_RANK(rblapack_b) != 2)
-    rb_raise(rb_eArgError, "rank of b (10th argument) must be %d", 2);
-  ldb = NA_SHAPE0(rblapack_b);
-  if (NA_SHAPE1(rblapack_b) != nrhs)
-    rb_raise(rb_eRuntimeError, "shape 1 of b must be the same as shape 1 of x");
-  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
-    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
-  b = NA_PTR_TYPE(rblapack_b, real*);
-  if (!NA_IsNArray(rblapack_d))
-    rb_raise(rb_eArgError, "d (3th argument) must be NArray");
-  if (NA_RANK(rblapack_d) != 1)
-    rb_raise(rb_eArgError, "rank of d (3th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_d) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of d must be the same as shape 0 of ipiv");
-  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
-    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
-  d = NA_PTR_TYPE(rblapack_d, real*);
+  if (!NA_IsNArray(rblapack_dl))
+    rb_raise(rb_eArgError, "dl (2th argument) must be NArray");
+  if (NA_RANK(rblapack_dl) != 1)
+    rb_raise(rb_eArgError, "rank of dl (2th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_dl) != (n-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of dl must be %d", n-1);
+  if (NA_TYPE(rblapack_dl) != NA_SFLOAT)
+    rblapack_dl = na_change_type(rblapack_dl, NA_SFLOAT);
+  dl = NA_PTR_TYPE(rblapack_dl, real*);
   if (!NA_IsNArray(rblapack_dlf))
     rb_raise(rb_eArgError, "dlf (5th argument) must be NArray");
   if (NA_RANK(rblapack_dlf) != 1)
@@ -130,15 +129,6 @@ rblapack_sgtrfs(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_dlf) != NA_SFLOAT)
     rblapack_dlf = na_change_type(rblapack_dlf, NA_SFLOAT);
   dlf = NA_PTR_TYPE(rblapack_dlf, real*);
-  if (!NA_IsNArray(rblapack_duf))
-    rb_raise(rb_eArgError, "duf (7th argument) must be NArray");
-  if (NA_RANK(rblapack_duf) != 1)
-    rb_raise(rb_eArgError, "rank of duf (7th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_duf) != (n-1))
-    rb_raise(rb_eRuntimeError, "shape 0 of duf must be %d", n-1);
-  if (NA_TYPE(rblapack_duf) != NA_SFLOAT)
-    rblapack_duf = na_change_type(rblapack_duf, NA_SFLOAT);
-  duf = NA_PTR_TYPE(rblapack_duf, real*);
   if (!NA_IsNArray(rblapack_du2))
     rb_raise(rb_eArgError, "du2 (8th argument) must be NArray");
   if (NA_RANK(rblapack_du2) != 1)
@@ -148,15 +138,6 @@ rblapack_sgtrfs(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_du2) != NA_SFLOAT)
     rblapack_du2 = na_change_type(rblapack_du2, NA_SFLOAT);
   du2 = NA_PTR_TYPE(rblapack_du2, real*);
-  if (!NA_IsNArray(rblapack_dl))
-    rb_raise(rb_eArgError, "dl (2th argument) must be NArray");
-  if (NA_RANK(rblapack_dl) != 1)
-    rb_raise(rb_eArgError, "rank of dl (2th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_dl) != (n-1))
-    rb_raise(rb_eRuntimeError, "shape 0 of dl must be %d", n-1);
-  if (NA_TYPE(rblapack_dl) != NA_SFLOAT)
-    rblapack_dl = na_change_type(rblapack_dl, NA_SFLOAT);
-  dl = NA_PTR_TYPE(rblapack_dl, real*);
   if (!NA_IsNArray(rblapack_du))
     rb_raise(rb_eArgError, "du (4th argument) must be NArray");
   if (NA_RANK(rblapack_du) != 1)
@@ -166,6 +147,25 @@ rblapack_sgtrfs(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_du) != NA_SFLOAT)
     rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
   du = NA_PTR_TYPE(rblapack_du, real*);
+  if (!NA_IsNArray(rblapack_b))
+    rb_raise(rb_eArgError, "b (10th argument) must be NArray");
+  if (NA_RANK(rblapack_b) != 2)
+    rb_raise(rb_eArgError, "rank of b (10th argument) must be %d", 2);
+  ldb = NA_SHAPE0(rblapack_b);
+  if (NA_SHAPE1(rblapack_b) != nrhs)
+    rb_raise(rb_eRuntimeError, "shape 1 of b must be the same as shape 1 of x");
+  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
+    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
+  b = NA_PTR_TYPE(rblapack_b, real*);
+  if (!NA_IsNArray(rblapack_duf))
+    rb_raise(rb_eArgError, "duf (7th argument) must be NArray");
+  if (NA_RANK(rblapack_duf) != 1)
+    rb_raise(rb_eArgError, "rank of duf (7th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_duf) != (n-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of duf must be %d", n-1);
+  if (NA_TYPE(rblapack_duf) != NA_SFLOAT)
+    rblapack_duf = na_change_type(rblapack_duf, NA_SFLOAT);
+  duf = NA_PTR_TYPE(rblapack_duf, real*);
   {
     int shape[1];
     shape[0] = nrhs;

data/ext/sgtsv.c

@@ -53,6 +53,14 @@ rblapack_sgtsv(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  if (!NA_IsNArray(rblapack_d))
+    rb_raise(rb_eArgError, "d (2th argument) must be NArray");
+  if (NA_RANK(rblapack_d) != 1)
+    rb_raise(rb_eArgError, "rank of d (2th argument) must be %d", 1);
+  n = NA_SHAPE0(rblapack_d);
+  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
+    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
+  d = NA_PTR_TYPE(rblapack_d, real*);
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (4th argument) must be NArray");
   if (NA_RANK(rblapack_b) != 2)
@@ -62,14 +70,6 @@ rblapack_sgtsv(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_b) != NA_SFLOAT)
     rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
   b = NA_PTR_TYPE(rblapack_b, real*);
-  if (!NA_IsNArray(rblapack_d))
-    rb_raise(rb_eArgError, "d (2th argument) must be NArray");
-  if (NA_RANK(rblapack_d) != 1)
-    rb_raise(rb_eArgError, "rank of d (2th argument) must be %d", 1);
-  n = NA_SHAPE0(rblapack_d);
-  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
-    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
-  d = NA_PTR_TYPE(rblapack_d, real*);
   if (!NA_IsNArray(rblapack_dl))
     rb_raise(rb_eArgError, "dl (1th argument) must be NArray");
   if (NA_RANK(rblapack_dl) != 1)

data/ext/sgtsvx.c

@@ -87,62 +87,35 @@ rblapack_sgtsvx(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_ipiv))
-    rb_raise(rb_eArgError, "ipiv (10th argument) must be NArray");
-  if (NA_RANK(rblapack_ipiv) != 1)
-    rb_raise(rb_eArgError, "rank of ipiv (10th argument) must be %d", 1);
-  n = NA_SHAPE0(rblapack_ipiv);
-  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
-    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
-  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
-  if (!NA_IsNArray(rblapack_b))
-    rb_raise(rb_eArgError, "b (11th argument) must be NArray");
-  if (NA_RANK(rblapack_b) != 2)
-    rb_raise(rb_eArgError, "rank of b (11th argument) must be %d", 2);
-  ldb = NA_SHAPE0(rblapack_b);
-  nrhs = NA_SHAPE1(rblapack_b);
-  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
-    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
-  b = NA_PTR_TYPE(rblapack_b, real*);
   fact = StringValueCStr(rblapack_fact)[0];
-  trans = StringValueCStr(rblapack_trans)[0];
+  if (!NA_IsNArray(rblapack_d))
+    rb_raise(rb_eArgError, "d (4th argument) must be NArray");
+  if (NA_RANK(rblapack_d) != 1)
+    rb_raise(rb_eArgError, "rank of d (4th argument) must be %d", 1);
+  n = NA_SHAPE0(rblapack_d);
+  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
+    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
+  d = NA_PTR_TYPE(rblapack_d, real*);
   if (!NA_IsNArray(rblapack_df))
     rb_raise(rb_eArgError, "df (7th argument) must be NArray");
   if (NA_RANK(rblapack_df) != 1)
     rb_raise(rb_eArgError, "rank of df (7th argument) must be %d", 1);
   if (NA_SHAPE0(rblapack_df) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of df must be the same as shape 0 of ipiv");
+    rb_raise(rb_eRuntimeError, "shape 0 of df must be the same as shape 0 of d");
   if (NA_TYPE(rblapack_df) != NA_SFLOAT)
     rblapack_df = na_change_type(rblapack_df, NA_SFLOAT);
   df = NA_PTR_TYPE(rblapack_df, real*);
-  if (!NA_IsNArray(rblapack_d))
-    rb_raise(rb_eArgError, "d (4th argument) must be NArray");
-  if (NA_RANK(rblapack_d) != 1)
-    rb_raise(rb_eArgError, "rank of d (4th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_d) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of d must be the same as shape 0 of ipiv");
-  if (NA_TYPE(rblapack_d) != NA_SFLOAT)
-    rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
-  d = NA_PTR_TYPE(rblapack_d, real*);
-  if (!NA_IsNArray(rblapack_dlf))
-    rb_raise(rb_eArgError, "dlf (6th argument) must be NArray");
-  if (NA_RANK(rblapack_dlf) != 1)
-    rb_raise(rb_eArgError, "rank of dlf (6th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_dlf) != (n-1))
-    rb_raise(rb_eRuntimeError, "shape 0 of dlf must be %d", n-1);
-  if (NA_TYPE(rblapack_dlf) != NA_SFLOAT)
-    rblapack_dlf = na_change_type(rblapack_dlf, NA_SFLOAT);
-  dlf = NA_PTR_TYPE(rblapack_dlf, real*);
-  if (!NA_IsNArray(rblapack_duf))
-    rb_raise(rb_eArgError, "duf (8th argument) must be NArray");
-  if (NA_RANK(rblapack_duf) != 1)
-    rb_raise(rb_eArgError, "rank of duf (8th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_duf) != (n-1))
-    rb_raise(rb_eRuntimeError, "shape 0 of duf must be %d", n-1);
-  if (NA_TYPE(rblapack_duf) != NA_SFLOAT)
-    rblapack_duf = na_change_type(rblapack_duf, NA_SFLOAT);
-  duf = NA_PTR_TYPE(rblapack_duf, real*);
+  if (!NA_IsNArray(rblapack_ipiv))
+    rb_raise(rb_eArgError, "ipiv (10th argument) must be NArray");
+  if (NA_RANK(rblapack_ipiv) != 1)
+    rb_raise(rb_eArgError, "rank of ipiv (10th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_ipiv) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 0 of d");
+  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
+    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
+  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
   ldx = MAX(1,n);
+  trans = StringValueCStr(rblapack_trans)[0];
   if (!NA_IsNArray(rblapack_du))
     rb_raise(rb_eArgError, "du (5th argument) must be NArray");
   if (NA_RANK(rblapack_du) != 1)
@@ -152,15 +125,24 @@ rblapack_sgtsvx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_du) != NA_SFLOAT)
     rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
   du = NA_PTR_TYPE(rblapack_du, real*);
-  if (!NA_IsNArray(rblapack_du2))
-    rb_raise(rb_eArgError, "du2 (9th argument) must be NArray");
-  if (NA_RANK(rblapack_du2) != 1)
-    rb_raise(rb_eArgError, "rank of du2 (9th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_du2) != (n-2))
-    rb_raise(rb_eRuntimeError, "shape 0 of du2 must be %d", n-2);
-  if (NA_TYPE(rblapack_du2) != NA_SFLOAT)
-    rblapack_du2 = na_change_type(rblapack_du2, NA_SFLOAT);
-  du2 = NA_PTR_TYPE(rblapack_du2, real*);
+  if (!NA_IsNArray(rblapack_duf))
+    rb_raise(rb_eArgError, "duf (8th argument) must be NArray");
+  if (NA_RANK(rblapack_duf) != 1)
+    rb_raise(rb_eArgError, "rank of duf (8th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_duf) != (n-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of duf must be %d", n-1);
+  if (NA_TYPE(rblapack_duf) != NA_SFLOAT)
+    rblapack_duf = na_change_type(rblapack_duf, NA_SFLOAT);
+  duf = NA_PTR_TYPE(rblapack_duf, real*);
+  if (!NA_IsNArray(rblapack_b))
+    rb_raise(rb_eArgError, "b (11th argument) must be NArray");
+  if (NA_RANK(rblapack_b) != 2)
+    rb_raise(rb_eArgError, "rank of b (11th argument) must be %d", 2);
+  ldb = NA_SHAPE0(rblapack_b);
+  nrhs = NA_SHAPE1(rblapack_b);
+  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
+    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
+  b = NA_PTR_TYPE(rblapack_b, real*);
   if (!NA_IsNArray(rblapack_dl))
     rb_raise(rb_eArgError, "dl (3th argument) must be NArray");
   if (NA_RANK(rblapack_dl) != 1)
@@ -170,6 +152,24 @@ rblapack_sgtsvx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_dl) != NA_SFLOAT)
     rblapack_dl = na_change_type(rblapack_dl, NA_SFLOAT);
   dl = NA_PTR_TYPE(rblapack_dl, real*);
+  if (!NA_IsNArray(rblapack_du2))
+    rb_raise(rb_eArgError, "du2 (9th argument) must be NArray");
+  if (NA_RANK(rblapack_du2) != 1)
+    rb_raise(rb_eArgError, "rank of du2 (9th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_du2) != (n-2))
+    rb_raise(rb_eRuntimeError, "shape 0 of du2 must be %d", n-2);
+  if (NA_TYPE(rblapack_du2) != NA_SFLOAT)
+    rblapack_du2 = na_change_type(rblapack_du2, NA_SFLOAT);
+  du2 = NA_PTR_TYPE(rblapack_du2, real*);
+  if (!NA_IsNArray(rblapack_dlf))
+    rb_raise(rb_eArgError, "dlf (6th argument) must be NArray");
+  if (NA_RANK(rblapack_dlf) != 1)
+    rb_raise(rb_eArgError, "rank of dlf (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_dlf) != (n-1))
+    rb_raise(rb_eRuntimeError, "shape 0 of dlf must be %d", n-1);
+  if (NA_TYPE(rblapack_dlf) != NA_SFLOAT)
+    rblapack_dlf = na_change_type(rblapack_dlf, NA_SFLOAT);
+  dlf = NA_PTR_TYPE(rblapack_dlf, real*);
   {
     int shape[2];
     shape[0] = ldx;

data/ext/sgttrs.c

@@ -56,33 +56,24 @@ rblapack_sgttrs(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_ipiv))
-    rb_raise(rb_eArgError, "ipiv (6th argument) must be NArray");
-  if (NA_RANK(rblapack_ipiv) != 1)
-    rb_raise(rb_eArgError, "rank of ipiv (6th argument) must be %d", 1);
-  n = NA_SHAPE0(rblapack_ipiv);
-  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
-    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
-  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
-  if (!NA_IsNArray(rblapack_b))
-    rb_raise(rb_eArgError, "b (7th argument) must be NArray");
-  if (NA_RANK(rblapack_b) != 2)
-    rb_raise(rb_eArgError, "rank of b (7th argument) must be %d", 2);
-  ldb = NA_SHAPE0(rblapack_b);
-  nrhs = NA_SHAPE1(rblapack_b);
-  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
-    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
-  b = NA_PTR_TYPE(rblapack_b, real*);
   trans = StringValueCStr(rblapack_trans)[0];
   if (!NA_IsNArray(rblapack_d))
     rb_raise(rb_eArgError, "d (3th argument) must be NArray");
   if (NA_RANK(rblapack_d) != 1)
     rb_raise(rb_eArgError, "rank of d (3th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_d) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of d must be the same as shape 0 of ipiv");
+  n = NA_SHAPE0(rblapack_d);
   if (NA_TYPE(rblapack_d) != NA_SFLOAT)
     rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
   d = NA_PTR_TYPE(rblapack_d, real*);
+  if (!NA_IsNArray(rblapack_ipiv))
+    rb_raise(rb_eArgError, "ipiv (6th argument) must be NArray");
+  if (NA_RANK(rblapack_ipiv) != 1)
+    rb_raise(rb_eArgError, "rank of ipiv (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_ipiv) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 0 of d");
+  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
+    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
+  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
   if (!NA_IsNArray(rblapack_dl))
     rb_raise(rb_eArgError, "dl (2th argument) must be NArray");
   if (NA_RANK(rblapack_dl) != 1)
@@ -110,6 +101,15 @@ rblapack_sgttrs(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_du) != NA_SFLOAT)
     rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
   du = NA_PTR_TYPE(rblapack_du, real*);
+  if (!NA_IsNArray(rblapack_b))
+    rb_raise(rb_eArgError, "b (7th argument) must be NArray");
+  if (NA_RANK(rblapack_b) != 2)
+    rb_raise(rb_eArgError, "rank of b (7th argument) must be %d", 2);
+  ldb = NA_SHAPE0(rblapack_b);
+  nrhs = NA_SHAPE1(rblapack_b);
+  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
+    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
+  b = NA_PTR_TYPE(rblapack_b, real*);
   {
     int shape[2];
     shape[0] = ldb;

data/ext/sgtts2.c

@@ -54,32 +54,24 @@ rblapack_sgtts2(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_ipiv))
-    rb_raise(rb_eArgError, "ipiv (6th argument) must be NArray");
-  if (NA_RANK(rblapack_ipiv) != 1)
-    rb_raise(rb_eArgError, "rank of ipiv (6th argument) must be %d", 1);
-  n = NA_SHAPE0(rblapack_ipiv);
-  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
-    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
-  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
-  if (!NA_IsNArray(rblapack_b))
-    rb_raise(rb_eArgError, "b (7th argument) must be NArray");
-  if (NA_RANK(rblapack_b) != 2)
-    rb_raise(rb_eArgError, "rank of b (7th argument) must be %d", 2);
-  ldb = NA_SHAPE0(rblapack_b);
-  nrhs = NA_SHAPE1(rblapack_b);
-  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
-    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
-  b = NA_PTR_TYPE(rblapack_b, real*);
+  itrans = NUM2INT(rblapack_itrans);
   if (!NA_IsNArray(rblapack_d))
     rb_raise(rb_eArgError, "d (3th argument) must be NArray");
   if (NA_RANK(rblapack_d) != 1)
     rb_raise(rb_eArgError, "rank of d (3th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_d) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of d must be the same as shape 0 of ipiv");
+  n = NA_SHAPE0(rblapack_d);
   if (NA_TYPE(rblapack_d) != NA_SFLOAT)
     rblapack_d = na_change_type(rblapack_d, NA_SFLOAT);
   d = NA_PTR_TYPE(rblapack_d, real*);
+  if (!NA_IsNArray(rblapack_ipiv))
+    rb_raise(rb_eArgError, "ipiv (6th argument) must be NArray");
+  if (NA_RANK(rblapack_ipiv) != 1)
+    rb_raise(rb_eArgError, "rank of ipiv (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_ipiv) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of ipiv must be the same as shape 0 of d");
+  if (NA_TYPE(rblapack_ipiv) != NA_LINT)
+    rblapack_ipiv = na_change_type(rblapack_ipiv, NA_LINT);
+  ipiv = NA_PTR_TYPE(rblapack_ipiv, integer*);
   if (!NA_IsNArray(rblapack_dl))
     rb_raise(rb_eArgError, "dl (2th argument) must be NArray");
   if (NA_RANK(rblapack_dl) != 1)
@@ -107,7 +99,15 @@ rblapack_sgtts2(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_du) != NA_SFLOAT)
     rblapack_du = na_change_type(rblapack_du, NA_SFLOAT);
   du = NA_PTR_TYPE(rblapack_du, real*);
-  itrans = NUM2INT(rblapack_itrans);
+  if (!NA_IsNArray(rblapack_b))
+    rb_raise(rb_eArgError, "b (7th argument) must be NArray");
+  if (NA_RANK(rblapack_b) != 2)
+    rb_raise(rb_eArgError, "rank of b (7th argument) must be %d", 2);
+  ldb = NA_SHAPE0(rblapack_b);
+  nrhs = NA_SHAPE1(rblapack_b);
+  if (NA_TYPE(rblapack_b) != NA_SFLOAT)
+    rblapack_b = na_change_type(rblapack_b, NA_SFLOAT);
+  b = NA_PTR_TYPE(rblapack_b, real*);
   {
     int shape[2];
     shape[0] = ldb;