ruby-mpfi 0.0.6 → 0.0.7

data/ext/mpfi_complex/mpfi/extconf.rb

@@ -20,5 +20,16 @@ dir_config("mpfr")
 dir_config("mpfi")
 dir_config("gmp")
 if have_header('mpfr.h') && have_library('mpfr') && have_header('gmp.h') && have_library('gmp')
+  ruby_mpfr_header_dir = nil
+  begin
+    require "rubygems"
+    spec = Gem::Specification.find_by_name("ruby-mpfr")
+    ruby_mpfr_header_dir = File.join(spec.full_gem_path, 'ext/mpfr')
+  rescue LoadError
+  end
+  unless find_header("ruby_mpfr.h", ruby_mpfr_header_dir)
+    header_not_found("ruby_mpfr.h")
+  end
+  find_header("ruby_mpfi.h", File.join(File.dirname(__FILE__), "../../mpfi"))
   create_makefile("mpfi/complex")
 end

data/ext/mpfi_complex/mpfi/ruby_mpfi_complex.c

@@ -66,9 +66,9 @@ static VALUE r_mpfi_complex_initialize_copy (VALUE self, VALUE other){
 /* Return real part. */
 static VALUE r_mpfi_complex_real (VALUE self){
   MPFIComplex *ptr_self;
-  r_mpfi_get_complex_struct(ptr_self, self);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_make_struct_init(ret, ptr_ret);
   mpfi_set(ptr_ret, ptr_self->re);
   return ret;
@@ -77,9 +77,9 @@ static VALUE r_mpfi_complex_real (VALUE self){
 /* Return imaginary part. */
 static VALUE r_mpfi_complex_imaginary (VALUE self){
   MPFIComplex *ptr_self;
-  r_mpfi_get_complex_struct(ptr_self, self);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_make_struct_init(ret, ptr_ret);
   mpfi_set(ptr_ret, ptr_self->im);
   return ret;
@@ -88,9 +88,9 @@ static VALUE r_mpfi_complex_imaginary (VALUE self){
 /* If _p1_ is 0, return real part. If _p1_ is 1, return imaginary part. */
 static VALUE r_mpfi_complex_element (VALUE self, VALUE arg){
   MPFIComplex *ptr_self;
-  r_mpfi_get_complex_struct(ptr_self, self);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_make_struct_init(ret, ptr_ret);
   switch(NUM2INT(arg)){
   case 0:
@@ -109,8 +109,8 @@ static VALUE r_mpfi_complex_element (VALUE self, VALUE arg){
 /* Return conjugate complex number. */
 static VALUE r_mpfi_complex_conjugate (VALUE self){
   MPFIComplex *ptr_self, *ptr_ret;
-  r_mpfi_get_complex_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_make_complex_struct_init(ret, ptr_ret);
   mpfi_complex_conjugate(ptr_ret, ptr_self);
   return ret;
@@ -130,15 +130,16 @@ static VALUE r_mpfi_complex_conjugate (VALUE self){
 /* String for inspect. */
 static VALUE r_mpfi_complex_inspect(VALUE self){
   MPFIComplex *ptr_s;
-  r_mpfi_get_complex_struct(ptr_s, self);
   char *ret_str;
+  VALUE ret_val;
+  r_mpfi_get_complex_struct(ptr_s, self);
   if (!mpfr_asprintf(&ret_str, "#<MPFI:%lx,['%.Re %.Re', '%.Re %.Re'], [%d, %d]>",
 		     NUM2LONG(rb_funcall(self, object_id, 0)), r_mpfi_left_ptr(ptr_s->re),
 		     r_mpfi_right_ptr(ptr_s->re), r_mpfi_left_ptr(ptr_s->im), r_mpfi_right_ptr(ptr_s->im),
 		     mpfi_get_prec(ptr_s->re), mpfi_get_prec(ptr_s->im))) {
     rb_raise(rb_eFatal, "Can not allocate a string by mpfr_asprintf.");
   }
-  VALUE ret_val = rb_str_new2(ret_str);
+  ret_val = rb_str_new2(ret_str);
   mpfr_free_str(ret_str);
   return ret_val;
 }
@@ -146,9 +147,10 @@ static VALUE r_mpfi_complex_inspect(VALUE self){
 /* Return array including strings which the elements is converted to by _p1_. */
 static VALUE r_mpfi_complex_str_ary(VALUE self, VALUE format_str){
   MPFIComplex *ptr_self;
-  r_mpfi_get_complex_struct(ptr_self, self);
-  char *format = StringValuePtr(format_str), *tmp_str;
+  char *format, *tmp_str;
   VALUE ret_val[2];
+  r_mpfi_get_complex_struct(ptr_self, self);
+  format = StringValuePtr(format_str);
   gmp_asprintf(&tmp_str, format, ptr_self->re);
   ret_val[0] = rb_str_new2(tmp_str);
   free(tmp_str);
@@ -161,9 +163,9 @@ static VALUE r_mpfi_complex_str_ary(VALUE self, VALUE format_str){
 /* Return _self_ + _p1_.*/
 static VALUE r_mpfi_complex_add (VALUE self, VALUE other){
   MPFIComplex *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_get_complex_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_complex_struct_init(ret, ptr_ret);
   mpfi_complex_add(ptr_ret, ptr_self, ptr_other);
   return ret;
@@ -172,9 +174,9 @@ static VALUE r_mpfi_complex_add (VALUE self, VALUE other){
 /* Return _self_ - _p1_.*/
 static VALUE r_mpfi_complex_sub (VALUE self, VALUE other){
   MPFIComplex *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_get_complex_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_complex_struct_init(ret, ptr_ret);
   mpfi_complex_sub(ptr_ret, ptr_self, ptr_other);
   return ret;
@@ -183,9 +185,9 @@ static VALUE r_mpfi_complex_sub (VALUE self, VALUE other){
 /* Return _self_ * _p1_.*/
 static VALUE r_mpfi_complex_mul (VALUE self, VALUE other){
   MPFIComplex *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_get_complex_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_complex_struct_init(ret, ptr_ret);
   mpfi_complex_mul(ptr_ret, ptr_self, ptr_other);
   return ret;
@@ -194,9 +196,9 @@ static VALUE r_mpfi_complex_mul (VALUE self, VALUE other){
 /* Return _self_ / _p1_.*/
 static VALUE r_mpfi_complex_div (VALUE self, VALUE other){
   MPFIComplex *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_complex_struct(ptr_self, self);
   r_mpfi_get_complex_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_complex_struct_init(ret, ptr_ret);
   mpfi_complex_div(ptr_ret, ptr_self, ptr_other);
   return ret;

data/ext/mpfi_matrix/mpfi/extconf.rb

@@ -19,5 +19,19 @@ dir_config('mpfi')
 dir_config('gmp')
 
 if have_header('mpfr.h') && have_library('mpfr') && have_header('mpfi.h') && have_library('mpfi') && have_header('gmp.h') && have_library('gmp')
+  ruby_mpfr_dir = nil
+  begin
+    require "rubygems"
+    spec = Gem::Specification.find_by_name("ruby-mpfr")
+    ruby_mpfr_dir = File.join(spec.full_gem_path, 'ext')
+  rescue LoadError
+  end
+  unless find_header("ruby_mpfr.h", File.join(ruby_mpfr_dir, 'mpfr'))
+    header_not_found("ruby_mpfr.h")
+  end
+  unless find_header("ruby_mpfr_matrix.h", File.join(ruby_mpfr_dir, 'mpfr_matrix/mpfr'))
+    header_not_found("ruby_mpfr_matrix.h")
+  end
+  find_header("ruby_mpfi.h", File.join(File.dirname(__FILE__), "../../mpfi"))
   create_makefile("mpfi/matrix")
 end

data/ext/mpfi_matrix/mpfi/func_mpfi_matrix.c

@@ -1,12 +1,12 @@
 #include "ruby_mpfi_matrix.h"
 
 void mpfi_matrix_init(MPFIMatrix *mat, int row, int column){
+  int i;
   mat->row = row;
   mat->column = column;
   mat->size = row * column;
   /* mat->data = (MPFI *)malloc(sizeof(MPFI) * mat->size); */
   mat->data = ALLOC_N(MPFI, mat->size);
-  int i;
   for(i = 0; i < mat->size; i++){
     mpfi_init(mat->data + i);
   }
@@ -305,8 +305,8 @@ int mpfi_matrix_union(MPFIMatrix *z, MPFIMatrix *x, MPFIMatrix *y){
 
 void mpfi_matrix_inner_product(MPFI *pr, MPFIMatrix *x, MPFIMatrix *y){
   MPFI *tmp;
-  r_mpfi_temp_alloc_init(tmp);
   int i;
+  r_mpfi_temp_alloc_init(tmp);
   mpfi_set_si(pr, 0);
   for(i = 0; i < x->size; i++){
     mpfi_mul(tmp, x->data + i, y->data + i);
@@ -317,8 +317,8 @@ void mpfi_matrix_inner_product(MPFI *pr, MPFIMatrix *x, MPFIMatrix *y){
 
 void mpfi_matrix_vector_distance(MPFI *distance, MPFIMatrix *x, MPFIMatrix *y){
   MPFI *tmp;
-  r_mpfi_temp_alloc_init(tmp);
   int i;
+  r_mpfi_temp_alloc_init(tmp);
   mpfi_set_si(distance, 0);
   for(i = 0; i < x->size; i++){
     mpfi_sub(tmp, x->data + i, y->data + i);
@@ -344,8 +344,8 @@ void mpfi_matrix_vector_distance_center_pts(MPFR *distance, MPFIMatrix *x, MPFIM
 
 void mpfi_matrix_vector_norm(MPFI *norm, MPFIMatrix *x){
   MPFI *tmp;
-  r_mpfi_temp_alloc_init(tmp);
   int i;
+  r_mpfi_temp_alloc_init(tmp);
   mpfi_set_si(norm, 0);
   for(i = 0; i < x->size; i++){
     mpfi_mul(tmp, x->data + i, x->data + i);
@@ -357,9 +357,9 @@ void mpfi_matrix_vector_norm(MPFI *norm, MPFIMatrix *x){
 
 void mpfi_matrix_max_norm(MPFI *norm, MPFIMatrix *x){
   MPFI *tmp, *abs;
+  int i;
   r_mpfi_temp_alloc_init(tmp);
   r_mpfi_temp_alloc_init(abs);
-  int i;
   mpfi_set_si(norm, 0);
   for(i = 0; i < x->size; i++){
     mpfi_abs(abs, x->data + i);
@@ -395,24 +395,24 @@ void mpfi_matrix_max_diam_abs(MPFR *diam, MPFIMatrix *x){
 /* ------------------- vector --------------------- */
 
 void mpfi_col_vector_init(MPFIMatrix *mat, int row){
+  int i;
   mat->row = row;
   mat->column = 1;
   mat->size = row;
   /* mat->data = (MPFI *)malloc(sizeof(MPF) * mat->size); */
   mat->data = ALLOC_N(MPFI, mat->size);
-  int i;
   for(i = 0; i < mat->size; i++){
     mpfi_init(mat->data + i);
   }
 }
 
 void mpfi_row_vector_init(MPFIMatrix *mat, int column){
+  int i;
   mat->row = 1;
   mat->column = column;
   mat->size = column;
   /* mat->data = (MPFI *)malloc(sizeof(MPF) * mat->size); */
   mat->data = ALLOC_N(MPFI, mat->size);
-  int i;
   for(i = 0; i < mat->size; i++){
     mpfi_init(mat->data + i);
   }
@@ -421,12 +421,12 @@ void mpfi_row_vector_init(MPFIMatrix *mat, int column){
 /* If length of MPFIMatrix *x is zero, return 1. Otherwise return 0. */
 int mpfi_vector_normalize(MPFIMatrix *new, MPFIMatrix *x){
   MPFRMatrix *fr_mat;
+  MPFR *norm;
+  int i, j, index, ret = 0;
   r_mpfr_matrix_temp_alloc_init(fr_mat, x->row, x->column);
   mpfi_matrix_mid(fr_mat, x);
-  MPFR *norm;
   r_mpfr_temp_alloc_init(norm);
   mpfr_matrix_vector_norm(norm, fr_mat);
-  int i, j, index, ret = 0;
   if(mpfr_cmp_ui(norm, 0) > 0){
     for(j = 0; j < x->column; j++){
       index = j * x->row;
@@ -454,11 +454,11 @@ void mpfi_vector_midpoint(MPFIMatrix *new, MPFIMatrix *x, MPFIMatrix *y){
 int mpfi_vector_set_length(MPFIMatrix *new, MPFIMatrix *x, MPFR *length){
   MPFI *norm_i;
   MPFR *factor_r;
+  int i, j, index, ret = 0;
   r_mpfi_temp_alloc_init(norm_i);
   r_mpfr_temp_alloc_init(factor_r);
   mpfi_matrix_vector_norm(norm_i, x);
   mpfi_mid(factor_r, norm_i);
-  int i, j, index, ret = 0;
   if(mpfr_cmp_ui(factor_r, 0) > 0){
     mpfr_ui_div(factor_r, 1, factor_r, GMP_RNDN);
     mpfr_mul(factor_r, factor_r, length, GMP_RNDN);
@@ -629,8 +629,9 @@ void mpfi_square_matrix_determinant(MPFI *det, MPFIMatrix *x){
     mpfi_3d_square_matrix_determinant(det, x);
   }else{
     MPFIMatrix *ptr_lu;
+    int *indx, i;
     r_mpfi_matrix_temp_alloc_init(ptr_lu, x->row, x->column);
-    int indx[x->row], i;
+    indx = (int *) malloc(sizeof(int) * x->row);
     if((i = mpfi_square_matrix_lu_decomp (ptr_lu, indx, x)) >= 0){
       if (i == 0) {
 	mpfi_set_si(det, 1);
@@ -644,22 +645,22 @@ void mpfi_square_matrix_determinant(MPFI *det, MPFIMatrix *x){
       mpfi_set_ui(det, 0);
     }
     r_mpfi_matrix_temp_free(ptr_lu);
+    free(indx);
   }
-
 }
 
 void mpfi_square_matrix_qr_decomp(MPFIMatrix *q, MPFIMatrix *r, MPFIMatrix *x){
   MPFIMatrix *q_mat, *r_mat;
+  int size, i, j, k, ind1, ind2, ind3;
+  MPFIMatrix *ary;
+  MPFI *tmp;
   r_mpfi_matrix_temp_alloc_init(q_mat, q->row, q->column);
   r_mpfi_matrix_temp_alloc_init(r_mat, r->row, r->column);
 
-  int size = x->row;
-  MPFIMatrix *ary;
+  size = x->row;
   r_mpfi_matrix_temp_alloc_init(ary, size, size);
   mpfi_matrix_set(ary, x);
-  MPFI *tmp;
   r_mpfi_temp_alloc_init(tmp);
-  int i, j, k, ind1, ind2, ind3;
   for (i = 0; i < size; i++) {
     ind1 = i * size;
     ind2 = i + ind1;
@@ -707,5 +708,4 @@ void mpfi_square_matrix_identity(MPFIMatrix *id){
       }
     }
   }
-  
 }

data/ext/mpfi_matrix/mpfi/ruby_mpfi_matrix.c

@@ -273,10 +273,10 @@ static VALUE r_mpfi_row_vector_initialize (VALUE self, VALUE arg){
 static VALUE r_mpfi_matrix_marshal_dump(VALUE self)
 {
   MPFIMatrix *ptr;
-  r_mpfi_get_matrix_struct(ptr, self);
   int i;
   char *tmp_str;
   VALUE ret_ary;
+  r_mpfi_get_matrix_struct(ptr, self);
   ret_ary = rb_ary_new();
   rb_ary_push(ret_ary, INT2FIX(ptr->row));
   rb_ary_push(ret_ary, INT2FIX(ptr->column));
@@ -293,15 +293,15 @@ static VALUE r_mpfi_matrix_marshal_dump(VALUE self)
 static VALUE r_mpfi_matrix_marshal_load(VALUE self, VALUE dump_ary)
 {
   MPFIMatrix *ptr;
+  int i;
+  char *dump;
+  VALUE dump_element;
   r_mpfi_get_matrix_struct(ptr, self);
 
   ptr->row = NUM2INT(rb_ary_entry(dump_ary, 0));
   ptr->column = NUM2INT(rb_ary_entry(dump_ary, 1));
   ptr->size = ptr->row * ptr->column;
   ptr->data = ALLOC_N(MPFI, ptr->size);
-  int i;
-  char *dump;
-  VALUE dump_element;
 
   for(i = 0; i < ptr->size; i++){
     dump_element = rb_ary_entry(dump_ary, i + 2);
@@ -336,8 +336,9 @@ static VALUE r_mpfi_matrix_row_size (VALUE self){
 /* Return element at _arg_.*/
 static VALUE r_mpfi_matrix_at (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self;
+  int i;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  int i = NUM2INT(arg);
+  i = NUM2INT(arg);
   if (i < ptr_self->size) {
     VALUE ret;
     MPFI *ptr_ret;
@@ -352,8 +353,10 @@ static VALUE r_mpfi_matrix_at (VALUE self, VALUE arg){
 /* Return element with _row_ and _column_.*/
 static VALUE r_mpfi_matrix_element (VALUE self, VALUE row, VALUE column){
   MPFIMatrix *ptr_self;
+  int i, j;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  int i = NUM2INT(row), j = NUM2INT(column);
+  i = NUM2INT(row);
+  j = NUM2INT(column);
   if (i < ptr_self->row && j < ptr_self->column) {
     VALUE ret;
     MPFI*ptr_ret;
@@ -368,8 +371,10 @@ static VALUE r_mpfi_matrix_element (VALUE self, VALUE row, VALUE column){
 /* Set _robj_ to element of which row is _row_ and column is _column_. */
 static VALUE r_mpfi_matrix_set_element (VALUE self, VALUE row, VALUE column, VALUE robj){
   MPFIMatrix *ptr_self;
+  int i, j;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  int i = NUM2INT(row), j = NUM2INT(column);
+  i = NUM2INT(row);
+  j = NUM2INT(column);
   if (i < ptr_self->row && j < ptr_self->column) {
     r_mpfi_set_robj(ptr_self->data + i + j * ptr_self->row, robj);
   }
@@ -379,9 +384,10 @@ static VALUE r_mpfi_matrix_set_element (VALUE self, VALUE row, VALUE column, VAL
 /* Evaluate block with each element of matrix. */
 static VALUE r_mpfi_matrix_each_element (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil;
+  VALUE ret;
   int i, j;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  ret = Qnil;
   for (i = 0; i < ptr_self->row; i++) {
     for (j = 0; j < ptr_self->column; j++) {
       volatile VALUE arg = r_mpfi_make_new_fi_obj(ptr_self->data + i + j * ptr_self->row);
@@ -394,9 +400,10 @@ static VALUE r_mpfi_matrix_each_element (VALUE self){
 /* Evaluate block with each element and its index. */
 static VALUE r_mpfi_matrix_each_element_with_index (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil, tmp_i;
+  VALUE ret, tmp_i;
   int i, j;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  ret = Qnil;
   for (i = 0; i < ptr_self->row; i++) {
     tmp_i = INT2NUM(i);
     for (j = 0; j < ptr_self->column; j++) {
@@ -410,8 +417,8 @@ static VALUE r_mpfi_matrix_each_element_with_index (VALUE self){
 /* Return one dimensinal array which has strings converted elements to. */
 static VALUE r_mpfi_matrix_str_ary_for_inspect(VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   char *tmp_str;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret_val[ptr_self->size];
   int i;
   for (i = 0; i < ptr_self->size; i++) {
@@ -428,10 +435,10 @@ static VALUE r_mpfi_matrix_str_ary_for_inspect(VALUE self){
 /* Return two dimensinal array which has strings converted elements to. */
 static VALUE r_mpfi_matrix_str_ary_for_inspect2(VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   char *tmp_str;
-  VALUE ary[ptr_self->row];
   int i, j;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  VALUE ary[ptr_self->row];
   for(i = 0; i < ptr_self->row; i++){
     ary[i] = rb_ary_new();
   }
@@ -450,10 +457,10 @@ static VALUE r_mpfi_matrix_str_ary_for_inspect2(VALUE self){
 
 static VALUE r_mpfi_matrix_to_str_ary(VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   char *tmp_str1, *tmp_str2;
-  VALUE ary[ptr_self->row];
   int i, j;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  VALUE ary[ptr_self->row];
   for(i = 0; i < ptr_self->row; i++){
     ary[i] = rb_ary_new();
   }
@@ -475,11 +482,11 @@ static VALUE r_mpfi_matrix_to_str_ary(VALUE self){
 
 static VALUE r_mpfi_matrix_to_strf_ary(VALUE self, VALUE format_str){
   MPFIMatrix *ptr_self;
+  char *tmp_str1, *tmp_str2, *format;
+  int i, j;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  char *tmp_str1, *tmp_str2;
-  char *format = StringValuePtr(format_str);
+  format = StringValuePtr(format_str);
   VALUE ary[ptr_self->row];
-  int i, j;
   for(i = 0; i < ptr_self->row; i++){
     ary[i] = rb_ary_new();
   }
@@ -501,9 +508,9 @@ static VALUE r_mpfi_matrix_to_strf_ary(VALUE self, VALUE format_str){
 
 static VALUE r_mpfi_matrix_to_array(VALUE self){
   MPFIMatrix *ptr_self;
+  int i;
   r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret_val[ptr_self->size];
-  int i;
   for (i = 0; i < ptr_self->size; i++) {
     ret_val[i] = r_mpfi_make_new_fi_obj(ptr_self->data + i);
   }
@@ -512,9 +519,9 @@ static VALUE r_mpfi_matrix_to_array(VALUE self){
 
 static VALUE r_mpfi_matrix_to_array2(VALUE self){
   MPFIMatrix *ptr_self;
+  int i, j;
   r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ary[ptr_self->row];
-  int i, j;
   for(i = 0; i < ptr_self->row; i++){
     ary[i] = rb_ary_new();
   }
@@ -528,9 +535,10 @@ static VALUE r_mpfi_matrix_to_array2(VALUE self){
 
 static VALUE r_mpfi_matrix_row (VALUE self, VALUE arg) {
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  int num = NUM2INT(arg);
+  int num;
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  num = NUM2INT(arg);
   if(num < ptr_self->row){
     r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, 1, ptr_self->column);
     mpfi_matrix_row(ptr_ret, ptr_self, num);
@@ -542,9 +550,10 @@ static VALUE r_mpfi_matrix_row (VALUE self, VALUE arg) {
 
 static VALUE r_mpfi_matrix_column (VALUE self, VALUE arg) {
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  int num = NUM2INT(arg);
+  int num;
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  num = NUM2INT(arg);
   if(num < ptr_self->column){
     r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, 1);
     mpfi_matrix_column(ptr_ret, ptr_self, num);
@@ -557,8 +566,8 @@ static VALUE r_mpfi_matrix_column (VALUE self, VALUE arg) {
 /* Return transposed matrix which is new object. */
 static VALUE r_mpfi_matrix_transpose (VALUE self){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->column, ptr_self->row);
   mpfi_matrix_transpose(ptr_ret, ptr_self);
   return ret;
@@ -567,6 +576,7 @@ static VALUE r_mpfi_matrix_transpose (VALUE self){
 /* Transpose self. This method is destroying method. */
 static VALUE r_mpfi_matrix_transpose2 (VALUE self){
   MPFIMatrix *ptr_self, tmp;
+  int i;
   r_mpfi_get_matrix_struct(ptr_self, self);
   if (ptr_self->column > 1 && ptr_self->row > 1) {
     mpfi_matrix_init(&tmp, ptr_self->column, ptr_self->row);
@@ -574,7 +584,6 @@ static VALUE r_mpfi_matrix_transpose2 (VALUE self){
     mpfi_matrix_set(ptr_self, &tmp);
     mpfi_matrix_clear(&tmp);
   }
-  int i;
   i = ptr_self->column;
   ptr_self->column = ptr_self->row;
   ptr_self->row = i;
@@ -584,8 +593,8 @@ static VALUE r_mpfi_matrix_transpose2 (VALUE self){
 /* Return transposed matrix which is new object. */
 static VALUE r_mpfi_matrix_neg (VALUE self){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   mpfi_matrix_neg(ptr_ret, ptr_self);
   return ret;
@@ -597,19 +606,18 @@ static VALUE r_mpfi_matrix_equal_p (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE ret = Qtrue;
   if (mpfi_matrix_equal_p(ptr_self, ptr_other) != 0) {
-    ret = Qnil;
+    return Qnil;
   }
-  return ret;
+  return Qtrue;
 }
 
 /* Return _self_ + _other_. */
 static VALUE r_mpfi_matrix_add (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE ret;
   if (ptr_self->column == ptr_other->column && ptr_self->row == ptr_other->row) {
     r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
     mpfi_matrix_add(ptr_ret, ptr_self, ptr_other);
@@ -622,8 +630,8 @@ static VALUE r_mpfi_matrix_add (VALUE self, VALUE other){
 /* Return _self_ + _other_ which is MPFI::Matrix. */
 static VALUE r_mpfi_matrix_add2 (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_make_matrix_struct(ret, ptr_ret);
 
   if(RTEST(rb_funcall(__mpfi_matrix_class__, eqq, 1, other))){
@@ -653,9 +661,9 @@ static VALUE r_mpfi_matrix_add2 (VALUE self, VALUE other){
 /* Return _self_ - _other_. */
 static VALUE r_mpfi_matrix_sub (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_matrix_struct(ret, ptr_ret);
   if (ptr_self->column == ptr_other->column && ptr_self->row == ptr_other->row) {
     r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
@@ -669,8 +677,8 @@ static VALUE r_mpfi_matrix_sub (VALUE self, VALUE other){
 /* Return _self_ - _other_ which is MPFI::Matrix. */
 static VALUE r_mpfi_matrix_sub2 (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_make_matrix_struct(ret, ptr_ret);
 
   if(RTEST(rb_funcall(__mpfi_matrix_class__, eqq, 1, other))){
@@ -699,10 +707,10 @@ static VALUE r_mpfi_matrix_sub2 (VALUE self, VALUE other){
 
 static VALUE r_mpfi_matrix_vector_inner_product (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_arg;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  r_mpfi_get_matrix_struct(ptr_arg, arg);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_get_matrix_struct(ptr_arg, arg);
   r_mpfi_make_struct_init(ret, ptr_ret);
   if(ptr_self->size == ptr_arg->size){
     mpfi_matrix_inner_product(ptr_ret, ptr_self, ptr_arg);
@@ -715,9 +723,9 @@ static VALUE r_mpfi_matrix_vector_inner_product (VALUE self, VALUE arg){
 /* Return _self_ * _other_. */
 static VALUE r_mpfi_matrix_mul_matrix (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE ret;
   r_mpfi_make_matrix_struct(ret, ptr_ret);
   if (ptr_self->column == ptr_other->row) {
     if((ptr_other->column == 1) && (ptr_self->row == 1)){
@@ -735,8 +743,8 @@ static VALUE r_mpfi_matrix_mul_matrix (VALUE self, VALUE other){
 /* Return _self_ * _other_ which is MPFI::Matrix. */
 static VALUE r_mpfi_matrix_mul_matrix2 (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_make_matrix_struct(ret, ptr_ret);
 
   if(RTEST(rb_funcall(__mpfi_matrix_class__, eqq, 1, other))){
@@ -767,8 +775,8 @@ static VALUE r_mpfi_matrix_mul_matrix2 (VALUE self, VALUE other){
 static VALUE r_mpfi_matrix_mul_scalar (VALUE self, VALUE scalar){
   MPFIMatrix *ptr_self, *ptr_ret;
   MPFI *ptr_scalar;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(RTEST(rb_funcall(r_mpfi_class, eqq, 1, scalar))){
     r_mpfi_get_struct(ptr_scalar, scalar);
@@ -785,8 +793,8 @@ static VALUE r_mpfi_matrix_mul_scalar (VALUE self, VALUE scalar){
 static VALUE r_mpfi_matrix_div_scalar (VALUE self, VALUE scalar){
   MPFIMatrix *ptr_self, *ptr_ret;
   MPFI *ptr_scalar;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(RTEST(rb_funcall(r_mpfi_class, eqq, 1, scalar))){
     r_mpfi_get_struct(ptr_scalar, scalar);
@@ -803,8 +811,9 @@ static VALUE r_mpfi_matrix_div_scalar (VALUE self, VALUE scalar){
 
 static VALUE r_mpfi_matrix_include_p (VALUE self, VALUE other){
   MPFIMatrix *ptr_self;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil;
+  ret = Qnil;
   if(RTEST(rb_funcall(r_mpfi_matrix, eqq, 1, other))){
     MPFIMatrix *ptr_other;
     r_mpfi_get_matrix_struct(ptr_other, other);
@@ -825,8 +834,9 @@ static VALUE r_mpfi_matrix_include_p (VALUE self, VALUE other){
 
 static VALUE r_mpfi_matrix_strictly_include_p (VALUE self, VALUE other){
   MPFIMatrix *ptr_self;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil;
+  ret = Qnil;
   if(RTEST(rb_funcall(r_mpfi_matrix, eqq, 1, other))){
     MPFIMatrix *ptr_other;
     r_mpfi_get_matrix_struct(ptr_other, other);
@@ -851,9 +861,9 @@ static VALUE r_mpfi_matrix_bounded_p (VALUE self) {
 
 static VALUE r_mpfi_matrix_mid (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE val_ret;
   MPFRMatrix *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfr_matrix_suitable_matrix_init(&val_ret, &ptr_ret, ptr_self->row, ptr_self->column);
   mpfi_matrix_mid(ptr_ret, ptr_self);
   return val_ret;
@@ -861,8 +871,8 @@ static VALUE r_mpfi_matrix_mid (VALUE self){
 
 static VALUE r_mpfi_matrix_mid_interval (VALUE self){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE val_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init(&val_ret, &ptr_ret, ptr_self->row, ptr_self->column);
   mpfi_matrix_mid_interval(ptr_ret, ptr_self);
   return val_ret;
@@ -870,9 +880,9 @@ static VALUE r_mpfi_matrix_mid_interval (VALUE self){
 
 static VALUE r_mpfi_matrix_intersect (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other, *ptr_ret;
+  VALUE val_ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE val_ret;
   r_mpfi_matrix_suitable_matrix_init(&val_ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(mpfi_matrix_intersect(ptr_ret, ptr_self, ptr_other) == 0){
     return val_ret;
@@ -883,9 +893,9 @@ static VALUE r_mpfi_matrix_intersect (VALUE self, VALUE other){
 
 static VALUE r_mpfi_matrix_union (VALUE self, VALUE other){
   MPFIMatrix *ptr_self, *ptr_other, *ptr_ret;
+  VALUE val_ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_other, other);
-  VALUE val_ret;
   r_mpfi_matrix_suitable_matrix_init(&val_ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(mpfi_matrix_union(ptr_ret, ptr_self, ptr_other) == 0){
     return val_ret;
@@ -897,9 +907,9 @@ static VALUE r_mpfi_matrix_union (VALUE self, VALUE other){
 
 static VALUE r_mpfi_matrix_max_diam_abs (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   MPFR *ptr_ret;
   VALUE ret_val;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfr_make_struct_init(ret_val, ptr_ret);
   mpfi_matrix_max_diam_abs(ptr_ret, ptr_self);
   return ret_val;
@@ -916,13 +926,13 @@ static VALUE r_mpfi_square_matrix_dim (VALUE self){
 static VALUE r_mpfi_square_matrix_lu_decomp (VALUE self){
   MPFIMatrix *ptr_self, *ptr_ret_l, *ptr_ret_u;
   VALUE ret_l, ret_u;
+  int i, j;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret_l, &ptr_ret_l, ptr_self->row, ptr_self->column);
   r_mpfi_matrix_suitable_matrix_init (&ret_u, &ptr_ret_u, ptr_self->row, ptr_self->column);
   VALUE ret_indx[ptr_self->row];
   int indx[ptr_self->row];
   if(mpfi_square_matrix_lu_decomp (ptr_ret_u, indx, ptr_self) >= 0){
-    int i, j;
     for(i = 1; i < ptr_ret_u->row; i++){
       for(j = 0; j < i; j++){
 	mpfi_set(mpfi_matrix_get_element(ptr_ret_l, i, j), mpfi_matrix_get_element(ptr_ret_u, i, j));
@@ -948,9 +958,9 @@ static VALUE r_mpfi_square_matrix_lu_decomp (VALUE self){
 
 static VALUE r_mpfi_square_matrix_determinant (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   MPFI *ptr_ret;
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_make_struct_init(ret, ptr_ret);
   mpfi_square_matrix_determinant(ptr_ret, ptr_self);
   return ret;
@@ -958,8 +968,8 @@ static VALUE r_mpfi_square_matrix_determinant (VALUE self){
 
 static VALUE r_mpfi_square_matrix_qr_decomp (VALUE self){
   MPFIMatrix *ptr_self, *ptr_q, *ptr_r;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE q, r;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&q, &ptr_q, ptr_self->row, ptr_self->column);
   r_mpfi_matrix_suitable_matrix_init (&r, &ptr_r, ptr_self->column, ptr_self->column);
   mpfi_square_matrix_qr_decomp(ptr_q, ptr_r, ptr_self);
@@ -986,9 +996,9 @@ static VALUE r_mpfi_vector_dim (VALUE self){
 /* Take matrix for vector and return length of the vector. */
 static VALUE r_mpfi_vector_abs (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_make_struct_init(ret, ptr_ret);
   mpfi_matrix_vector_norm(ptr_ret, ptr_self);
   return ret;
@@ -997,8 +1007,9 @@ static VALUE r_mpfi_vector_abs (VALUE self){
 /* Return element at _arg_.*/
 static VALUE r_mpfi_vector_element (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self;
+  int i;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  int i = NUM2INT(arg);
+  i = NUM2INT(arg);
   if (i < ptr_self->size) {
     VALUE ret;
     MPFI *ptr_ret;
@@ -1013,8 +1024,9 @@ static VALUE r_mpfi_vector_element (VALUE self, VALUE arg){
 /* Set _robj_ to _arg_ th element. */
 static VALUE r_mpfi_vector_set_element (VALUE self, VALUE arg, VALUE robj){
   MPFIMatrix *ptr_self;
+  int i;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  int i = NUM2INT(arg);
+  i = NUM2INT(arg);
   if (i < ptr_self->size) {
     r_mpfi_set_robj(ptr_self->data + i, robj);
   }
@@ -1024,11 +1036,12 @@ static VALUE r_mpfi_vector_set_element (VALUE self, VALUE arg, VALUE robj){
 /* Evaluate block with each element of vector. */
 static VALUE r_mpfi_vector_each_element (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil, tmp;
+  VALUE ret, tmp;
   MPFI *tmp_ptr;
-  r_mpfi_make_struct_init(tmp, tmp_ptr);
   int i;
+  ret = Qnil;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_make_struct_init(tmp, tmp_ptr);
   for (i = 0; i < ptr_self->size; i++) {
     mpfi_set(tmp_ptr, ptr_self->data + i);
     ret = rb_yield(tmp);
@@ -1039,11 +1052,12 @@ static VALUE r_mpfi_vector_each_element (VALUE self){
 /* Evaluate block with each element and its index. */
 static VALUE r_mpfi_vector_each_element_with_index (VALUE self){
   MPFIMatrix *ptr_self;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = Qnil, tmp;
+  VALUE ret, tmp;
   MPFI *tmp_ptr;
-  r_mpfi_make_struct_init(tmp, tmp_ptr);
   int i;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_make_struct_init(tmp, tmp_ptr);
+  ret = Qnil;
   for (i = 0; i < ptr_self->size; i++) {
     mpfi_set(tmp_ptr, ptr_self->data + i);
     ret = rb_yield(rb_ary_new3(2, tmp, INT2NUM(i)));
@@ -1053,10 +1067,10 @@ static VALUE r_mpfi_vector_each_element_with_index (VALUE self){
 
 static VALUE r_mpfi_vector_inner_product (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_arg;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  r_mpfi_get_matrix_struct(ptr_arg, arg);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_get_matrix_struct(ptr_arg, arg);
   r_mpfi_make_struct_init(ret, ptr_ret);
   if(ptr_self->size == ptr_arg->size){
     mpfi_matrix_inner_product(ptr_ret, ptr_self, ptr_arg);
@@ -1068,10 +1082,10 @@ static VALUE r_mpfi_vector_inner_product (VALUE self, VALUE arg){
 
 static VALUE r_mpfi_vector_distance_from (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_arg;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  r_mpfi_get_matrix_struct(ptr_arg, arg);
   VALUE ret;
   MPFI *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_get_matrix_struct(ptr_arg, arg);
   r_mpfi_make_struct_init(ret, ptr_ret);
   if(ptr_self->size == ptr_arg->size){
     mpfi_matrix_vector_distance(ptr_ret, ptr_self, ptr_arg);
@@ -1083,10 +1097,10 @@ static VALUE r_mpfi_vector_distance_from (VALUE self, VALUE arg){
 
 static VALUE r_mpfi_vector_distance_center_pts (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_arg;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  r_mpfi_get_matrix_struct(ptr_arg, arg);
   VALUE ret;
   MPFR *ptr_ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
+  r_mpfi_get_matrix_struct(ptr_arg, arg);
   r_mpfr_make_struct_init(ret, ptr_ret);
   if(ptr_self->size == ptr_arg->size){
     mpfi_matrix_vector_distance_center_pts(ptr_ret, ptr_self, ptr_arg);
@@ -1098,9 +1112,9 @@ static VALUE r_mpfi_vector_distance_center_pts (VALUE self, VALUE arg){
 
 static VALUE r_mpfi_vector_midpoint (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_arg, *ptr_ret;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_get_matrix_struct(ptr_arg, arg);
-  VALUE ret;
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   mpfi_vector_midpoint(ptr_ret, ptr_self, ptr_arg);
   return ret;
@@ -1108,8 +1122,8 @@ static VALUE r_mpfi_vector_midpoint (VALUE self, VALUE arg){
 
 static VALUE r_mpfi_vector_normalize (VALUE self){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
   VALUE ret;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(mpfi_vector_normalize(ptr_ret, ptr_self) == 0){
     return ret;
@@ -1121,8 +1135,9 @@ static VALUE r_mpfi_vector_normalize (VALUE self){
 /* Return normalized vector of _self_. This method is destroy method. */
 static VALUE r_mpfi_vector_normalize2 (VALUE self){
   MPFIMatrix *ptr_self, tmp;
+  VALUE ret;
   r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret = self;
+  ret = self;
   mpfi_matrix_init(&tmp, ptr_self->column, ptr_self->row);
   if(mpfi_vector_normalize(&tmp, ptr_self) == 0){
     mpfi_matrix_set(ptr_self, &tmp);
@@ -1135,9 +1150,10 @@ static VALUE r_mpfi_vector_normalize2 (VALUE self){
 
 static VALUE r_mpfi_vector_set_length (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret, returned_value = Qnil;
+  VALUE ret, returned_value;
   MPFR *ptr_l;
+  returned_value = Qnil;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(RTEST(rb_funcall(r_mpfr_class, eqq, 1, arg))){
     r_mpfr_get_struct(ptr_l, arg);
@@ -1157,9 +1173,10 @@ static VALUE r_mpfi_vector_set_length (VALUE self, VALUE arg){
 
 static VALUE r_mpfi_vector_set_length2 (VALUE self, VALUE arg){
   MPFIMatrix *ptr_self, *ptr_ret;
-  r_mpfi_get_matrix_struct(ptr_self, self);
-  VALUE ret, returned_value = Qnil;
+  VALUE ret, returned_value;
   MPFR *ptr_l;
+  returned_value = Qnil;
+  r_mpfi_get_matrix_struct(ptr_self, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr_self->row, ptr_self->column);
   if(RTEST(rb_funcall(r_mpfi_class, eqq, 1, arg))){
     r_mpfr_get_struct(ptr_l, arg);
@@ -1180,19 +1197,17 @@ static VALUE r_mpfi_vector_set_length2 (VALUE self, VALUE arg){
 }
 
 
-
 static VALUE r_mpfr_matrix_to_fi_matrix(VALUE self){
   MPFRMatrix *ptr;
-  r_mpfr_get_matrix_struct(ptr, self);
   VALUE ret;
   MPFIMatrix *ptr_ret;
+  r_mpfr_get_matrix_struct(ptr, self);
   r_mpfi_matrix_suitable_matrix_init (&ret, &ptr_ret, ptr->row, ptr->column);
   mpfi_matrix_from_mpfr_matrix(ptr_ret, ptr);
   return ret;
 }
 
 void Init_matrix(){
-
   /* Initialization of MPFI::Matrix, MPFI::SquareMatrix, MPFI::ColumnVector and MPFI::RowVector */
   VALUE tmp_r_mpfi_class = rb_define_class("MPFI", rb_cNumeric);
   r_mpfi_matrix = rb_define_class_under(tmp_r_mpfi_class, "Matrix", rb_cObject);
@@ -1316,6 +1331,4 @@ void Init_matrix(){
   eqq = rb_intern("===");
   __mpfr_matrix_class__ = rb_eval_string("MPFR::Matrix");
   __mpfi_matrix_class__ = rb_eval_string("MPFI::Matrix");
-
 }
-