ruby-lapack 1.4.1a → 1.5

data/ext/sporfsx.c

@@ -78,6 +78,7 @@ rblapack_sporfsx(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (3th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -87,16 +88,6 @@ rblapack_sporfsx(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_b))
-    rb_raise(rb_eArgError, "b (6th argument) must be NArray");
-  if (NA_RANK(rblapack_b) != 2)
-    rb_raise(rb_eArgError, "rank of b (6th 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*);
-  n_err_bnds = 3;
   if (!NA_IsNArray(rblapack_s))
     rb_raise(rb_eArgError, "s (5th argument) must be NArray");
   if (NA_RANK(rblapack_s) != 1)
@@ -106,25 +97,27 @@ rblapack_sporfsx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_s) != NA_SFLOAT)
     rblapack_s = na_change_type(rblapack_s, NA_SFLOAT);
   s = NA_PTR_TYPE(rblapack_s, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_x))
     rb_raise(rb_eArgError, "x (7th argument) must be NArray");
   if (NA_RANK(rblapack_x) != 2)
     rb_raise(rb_eArgError, "rank of x (7th argument) must be %d", 2);
   ldx = NA_SHAPE0(rblapack_x);
-  if (NA_SHAPE1(rblapack_x) != nrhs)
-    rb_raise(rb_eRuntimeError, "shape 1 of x must be the same as shape 1 of b");
+  nrhs = NA_SHAPE1(rblapack_x);
   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_params))
-    rb_raise(rb_eArgError, "params (8th argument) must be NArray");
-  if (NA_RANK(rblapack_params) != 1)
-    rb_raise(rb_eArgError, "rank of params (8th argument) must be %d", 1);
-  nparams = NA_SHAPE0(rblapack_params);
-  if (NA_TYPE(rblapack_params) != NA_SFLOAT)
-    rblapack_params = na_change_type(rblapack_params, NA_SFLOAT);
-  params = NA_PTR_TYPE(rblapack_params, real*);
+  n_err_bnds = 3;
+  equed = StringValueCStr(rblapack_equed)[0];
+  if (!NA_IsNArray(rblapack_b))
+    rb_raise(rb_eArgError, "b (6th argument) must be NArray");
+  if (NA_RANK(rblapack_b) != 2)
+    rb_raise(rb_eArgError, "rank of b (6th 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_af))
     rb_raise(rb_eArgError, "af (4th argument) must be NArray");
   if (NA_RANK(rblapack_af) != 2)
@@ -135,7 +128,14 @@ rblapack_sporfsx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_af) != NA_SFLOAT)
     rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
   af = NA_PTR_TYPE(rblapack_af, real*);
-  equed = StringValueCStr(rblapack_equed)[0];
+  if (!NA_IsNArray(rblapack_params))
+    rb_raise(rb_eArgError, "params (8th argument) must be NArray");
+  if (NA_RANK(rblapack_params) != 1)
+    rb_raise(rb_eArgError, "rank of params (8th argument) must be %d", 1);
+  nparams = NA_SHAPE0(rblapack_params);
+  if (NA_TYPE(rblapack_params) != NA_SFLOAT)
+    rblapack_params = na_change_type(rblapack_params, NA_SFLOAT);
+  params = NA_PTR_TYPE(rblapack_params, real*);
   {
     int shape[1];
     shape[0] = nrhs;

data/ext/sposv.c

@@ -47,15 +47,6 @@ rblapack_sposv(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_a))
-    rb_raise(rb_eArgError, "a (2th argument) must be NArray");
-  if (NA_RANK(rblapack_a) != 2)
-    rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
-  lda = NA_SHAPE0(rblapack_a);
-  n = NA_SHAPE1(rblapack_a);
-  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];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (3th argument) must be NArray");
@@ -66,6 +57,15 @@ rblapack_sposv(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_a))
+    rb_raise(rb_eArgError, "a (2th argument) must be NArray");
+  if (NA_RANK(rblapack_a) != 2)
+    rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
+  lda = NA_SHAPE0(rblapack_a);
+  n = NA_SHAPE1(rblapack_a);
+  if (NA_TYPE(rblapack_a) != NA_SFLOAT)
+    rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
+  a = NA_PTR_TYPE(rblapack_a, real*);
   {
     int shape[2];
     shape[0] = lda;

data/ext/sposvx.c

@@ -75,6 +75,7 @@ rblapack_sposvx(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  fact = StringValueCStr(rblapack_fact)[0];
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (3th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -84,17 +85,7 @@ rblapack_sposvx(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*);
-  fact = StringValueCStr(rblapack_fact)[0];
-  if (!NA_IsNArray(rblapack_s))
-    rb_raise(rb_eArgError, "s (6th argument) must be NArray");
-  if (NA_RANK(rblapack_s) != 1)
-    rb_raise(rb_eArgError, "rank of s (6th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_s) != n)
-    rb_raise(rb_eRuntimeError, "shape 0 of s must be the same as shape 1 of a");
-  if (NA_TYPE(rblapack_s) != NA_SFLOAT)
-    rblapack_s = na_change_type(rblapack_s, NA_SFLOAT);
-  s = NA_PTR_TYPE(rblapack_s, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
+  equed = StringValueCStr(rblapack_equed)[0];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (7th argument) must be NArray");
   if (NA_RANK(rblapack_b) != 2)
@@ -104,7 +95,16 @@ rblapack_sposvx(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*);
-  equed = StringValueCStr(rblapack_equed)[0];
+  uplo = StringValueCStr(rblapack_uplo)[0];
+  if (!NA_IsNArray(rblapack_s))
+    rb_raise(rb_eArgError, "s (6th argument) must be NArray");
+  if (NA_RANK(rblapack_s) != 1)
+    rb_raise(rb_eArgError, "rank of s (6th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_s) != n)
+    rb_raise(rb_eRuntimeError, "shape 0 of s must be the same as shape 1 of a");
+  if (NA_TYPE(rblapack_s) != NA_SFLOAT)
+    rblapack_s = na_change_type(rblapack_s, NA_SFLOAT);
+  s = NA_PTR_TYPE(rblapack_s, real*);
   if (!NA_IsNArray(rblapack_af))
     rb_raise(rb_eArgError, "af (4th argument) must be NArray");
   if (NA_RANK(rblapack_af) != 2)

data/ext/sposvxx.c

@@ -86,6 +86,7 @@ rblapack_sposvxx(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  fact = StringValueCStr(rblapack_fact)[0];
   if (!NA_IsNArray(rblapack_a))
     rb_raise(rb_eArgError, "a (3th argument) must be NArray");
   if (NA_RANK(rblapack_a) != 2)
@@ -95,8 +96,18 @@ rblapack_sposvxx(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*);
-  fact = StringValueCStr(rblapack_fact)[0];
+  equed = StringValueCStr(rblapack_equed)[0];
+  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*);
   n_err_bnds = 3;
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_s))
     rb_raise(rb_eArgError, "s (6th argument) must be NArray");
   if (NA_RANK(rblapack_s) != 1)
@@ -106,24 +117,6 @@ rblapack_sposvxx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_s) != NA_SFLOAT)
     rblapack_s = na_change_type(rblapack_s, NA_SFLOAT);
   s = NA_PTR_TYPE(rblapack_s, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
-  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*);
-  if (!NA_IsNArray(rblapack_params))
-    rb_raise(rb_eArgError, "params (8th argument) must be NArray");
-  if (NA_RANK(rblapack_params) != 1)
-    rb_raise(rb_eArgError, "rank of params (8th argument) must be %d", 1);
-  nparams = NA_SHAPE0(rblapack_params);
-  if (NA_TYPE(rblapack_params) != NA_SFLOAT)
-    rblapack_params = na_change_type(rblapack_params, NA_SFLOAT);
-  params = NA_PTR_TYPE(rblapack_params, real*);
   if (!NA_IsNArray(rblapack_af))
     rb_raise(rb_eArgError, "af (4th argument) must be NArray");
   if (NA_RANK(rblapack_af) != 2)
@@ -135,7 +128,14 @@ rblapack_sposvxx(int argc, VALUE *argv, VALUE self){
     rblapack_af = na_change_type(rblapack_af, NA_SFLOAT);
   af = NA_PTR_TYPE(rblapack_af, real*);
   ldx = MAX(1,n);
-  equed = StringValueCStr(rblapack_equed)[0];
+  if (!NA_IsNArray(rblapack_params))
+    rb_raise(rb_eArgError, "params (8th argument) must be NArray");
+  if (NA_RANK(rblapack_params) != 1)
+    rb_raise(rb_eArgError, "rank of params (8th argument) must be %d", 1);
+  nparams = NA_SHAPE0(rblapack_params);
+  if (NA_TYPE(rblapack_params) != NA_SFLOAT)
+    rblapack_params = na_change_type(rblapack_params, NA_SFLOAT);
+  params = NA_PTR_TYPE(rblapack_params, real*);
   {
     int shape[2];
     shape[0] = ldx;

data/ext/spotf2.c

@@ -40,6 +40,7 @@ rblapack_spotf2(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_spotf2(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/spotrf.c

@@ -40,6 +40,7 @@ rblapack_spotrf(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_spotrf(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/spotri.c

@@ -40,6 +40,7 @@ rblapack_spotri(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_spotri(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/spotrs.c

@@ -45,15 +45,6 @@ rblapack_spotrs(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_a))
-    rb_raise(rb_eArgError, "a (2th argument) must be NArray");
-  if (NA_RANK(rblapack_a) != 2)
-    rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
-  lda = NA_SHAPE0(rblapack_a);
-  n = NA_SHAPE1(rblapack_a);
-  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];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (3th argument) must be NArray");
@@ -64,6 +55,15 @@ rblapack_spotrs(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_a))
+    rb_raise(rb_eArgError, "a (2th argument) must be NArray");
+  if (NA_RANK(rblapack_a) != 2)
+    rb_raise(rb_eArgError, "rank of a (2th argument) must be %d", 2);
+  lda = NA_SHAPE0(rblapack_a);
+  n = NA_SHAPE1(rblapack_a);
+  if (NA_TYPE(rblapack_a) != NA_SFLOAT)
+    rblapack_a = na_change_type(rblapack_a, NA_SFLOAT);
+  a = NA_PTR_TYPE(rblapack_a, real*);
   {
     int shape[2];
     shape[0] = ldb;

data/ext/sppcon.c

@@ -45,6 +45,7 @@ rblapack_sppcon(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   anorm = (real)NUM2DBL(rblapack_anorm);
   if (!NA_IsNArray(rblapack_ap))
     rb_raise(rb_eArgError, "ap (2th argument) must be NArray");
@@ -55,7 +56,6 @@ rblapack_sppcon(int argc, VALUE *argv, VALUE self){
     rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
   ap = NA_PTR_TYPE(rblapack_ap, real*);
   n = ((int)sqrtf(ldap*8+1.0f)-1)/2;
-  uplo = StringValueCStr(rblapack_uplo)[0];
   work = ALLOC_N(real, (3*n));
   iwork = ALLOC_N(integer, (n));
 

data/ext/sppequ.c

@@ -44,6 +44,7 @@ rblapack_sppequ(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_ap))
     rb_raise(rb_eArgError, "ap (2th argument) must be NArray");
   if (NA_RANK(rblapack_ap) != 1)
@@ -53,7 +54,6 @@ rblapack_sppequ(int argc, VALUE *argv, VALUE self){
     rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
   ap = NA_PTR_TYPE(rblapack_ap, real*);
   n = ((int)sqrtf(ldap*8+1.0f)-1)/2;
-  uplo = StringValueCStr(rblapack_uplo)[0];
   {
     int shape[1];
     shape[0] = n;

data/ext/spprfs.c

@@ -57,36 +57,17 @@ rblapack_spprfs(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
-  if (!NA_IsNArray(rblapack_x))
-    rb_raise(rb_eArgError, "x (5th argument) must be NArray");
-  if (NA_RANK(rblapack_x) != 2)
-    rb_raise(rb_eArgError, "rank of x (5th argument) must be %d", 2);
-  ldx = NA_SHAPE0(rblapack_x);
-  nrhs = NA_SHAPE1(rblapack_x);
-  if (NA_TYPE(rblapack_x) != NA_SFLOAT)
-    rblapack_x = na_change_type(rblapack_x, NA_SFLOAT);
-  x = NA_PTR_TYPE(rblapack_x, real*);
   uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (4th argument) must be NArray");
   if (NA_RANK(rblapack_b) != 2)
     rb_raise(rb_eArgError, "rank of b (4th 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");
+  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*);
   n = ldb;
-  if (!NA_IsNArray(rblapack_afp))
-    rb_raise(rb_eArgError, "afp (3th argument) must be NArray");
-  if (NA_RANK(rblapack_afp) != 1)
-    rb_raise(rb_eArgError, "rank of afp (3th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_afp) != (n*(n+1)/2))
-    rb_raise(rb_eRuntimeError, "shape 0 of afp must be %d", n*(n+1)/2);
-  if (NA_TYPE(rblapack_afp) != NA_SFLOAT)
-    rblapack_afp = na_change_type(rblapack_afp, NA_SFLOAT);
-  afp = NA_PTR_TYPE(rblapack_afp, real*);
   if (!NA_IsNArray(rblapack_ap))
     rb_raise(rb_eArgError, "ap (2th argument) must be NArray");
   if (NA_RANK(rblapack_ap) != 1)
@@ -96,6 +77,25 @@ rblapack_spprfs(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*);
+  if (!NA_IsNArray(rblapack_x))
+    rb_raise(rb_eArgError, "x (5th argument) must be NArray");
+  if (NA_RANK(rblapack_x) != 2)
+    rb_raise(rb_eArgError, "rank of x (5th argument) must be %d", 2);
+  ldx = NA_SHAPE0(rblapack_x);
+  if (NA_SHAPE1(rblapack_x) != nrhs)
+    rb_raise(rb_eRuntimeError, "shape 1 of x must be the same as shape 1 of b");
+  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_afp))
+    rb_raise(rb_eArgError, "afp (3th argument) must be NArray");
+  if (NA_RANK(rblapack_afp) != 1)
+    rb_raise(rb_eArgError, "rank of afp (3th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_afp) != (n*(n+1)/2))
+    rb_raise(rb_eRuntimeError, "shape 0 of afp must be %d", n*(n+1)/2);
+  if (NA_TYPE(rblapack_afp) != NA_SFLOAT)
+    rblapack_afp = na_change_type(rblapack_afp, NA_SFLOAT);
+  afp = NA_PTR_TYPE(rblapack_afp, real*);
   {
     int shape[1];
     shape[0] = nrhs;

data/ext/sppsv.c

@@ -48,6 +48,7 @@ rblapack_sppsv(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (4th argument) must be NArray");
   if (NA_RANK(rblapack_b) != 2)
@@ -57,7 +58,6 @@ rblapack_sppsv(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*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
   n = NUM2INT(rblapack_n);
   if (!NA_IsNArray(rblapack_ap))
     rb_raise(rb_eArgError, "ap (3th argument) must be NArray");

data/ext/sppsvx.c

@@ -73,6 +73,8 @@ rblapack_sppsvx(int argc, VALUE *argv, VALUE self){
   } else {
   }
 
+  fact = StringValueCStr(rblapack_fact)[0];
+  equed = StringValueCStr(rblapack_equed)[0];
   if (!NA_IsNArray(rblapack_b))
     rb_raise(rb_eArgError, "b (7th argument) must be NArray");
   if (NA_RANK(rblapack_b) != 2)
@@ -82,6 +84,7 @@ rblapack_sppsvx(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*);
+  uplo = StringValueCStr(rblapack_uplo)[0];
   if (!NA_IsNArray(rblapack_s))
     rb_raise(rb_eArgError, "s (6th argument) must be NArray");
   if (NA_RANK(rblapack_s) != 1)
@@ -90,10 +93,16 @@ rblapack_sppsvx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_s) != NA_SFLOAT)
     rblapack_s = na_change_type(rblapack_s, NA_SFLOAT);
   s = NA_PTR_TYPE(rblapack_s, real*);
-  uplo = StringValueCStr(rblapack_uplo)[0];
+  if (!NA_IsNArray(rblapack_ap))
+    rb_raise(rb_eArgError, "ap (3th argument) must be NArray");
+  if (NA_RANK(rblapack_ap) != 1)
+    rb_raise(rb_eArgError, "rank of ap (3th argument) must be %d", 1);
+  if (NA_SHAPE0(rblapack_ap) != (n*(n+1)/2))
+    rb_raise(rb_eRuntimeError, "shape 0 of ap must be %d", n*(n+1)/2);
+  if (NA_TYPE(rblapack_ap) != NA_SFLOAT)
+    rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
+  ap = NA_PTR_TYPE(rblapack_ap, real*);
   ldx = MAX(1,n);
-  fact = StringValueCStr(rblapack_fact)[0];
-  equed = StringValueCStr(rblapack_equed)[0];
   if (!NA_IsNArray(rblapack_afp))
     rb_raise(rb_eArgError, "afp (4th argument) must be NArray");
   if (NA_RANK(rblapack_afp) != 1)
@@ -103,15 +112,6 @@ rblapack_sppsvx(int argc, VALUE *argv, VALUE self){
   if (NA_TYPE(rblapack_afp) != NA_SFLOAT)
     rblapack_afp = na_change_type(rblapack_afp, NA_SFLOAT);
   afp = NA_PTR_TYPE(rblapack_afp, real*);
-  if (!NA_IsNArray(rblapack_ap))
-    rb_raise(rb_eArgError, "ap (3th argument) must be NArray");
-  if (NA_RANK(rblapack_ap) != 1)
-    rb_raise(rb_eArgError, "rank of ap (3th argument) must be %d", 1);
-  if (NA_SHAPE0(rblapack_ap) != (n*(n+1)/2))
-    rb_raise(rb_eRuntimeError, "shape 0 of ap must be %d", n*(n+1)/2);
-  if (NA_TYPE(rblapack_ap) != NA_SFLOAT)
-    rblapack_ap = na_change_type(rblapack_ap, NA_SFLOAT);
-  ap = NA_PTR_TYPE(rblapack_ap, real*);
   {
     int shape[2];
     shape[0] = ldx;