gps_pvt 0.2.1 → 0.3.3

data/ext/gps_pvt/SylphideMath/SylphideMath_wrap.cxx

@@ -1862,8 +1862,9 @@ int SWIG_Ruby_arity( VALUE proc, int minimal )
 #define SWIGTYPE_p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void swig_types[12]
 #define SWIGTYPE_p_self_t swig_types[13]
 #define SWIGTYPE_p_swig__GC_VALUE swig_types[14]
-static swig_type_info *swig_types[16];
-static swig_module_info swig_module = {swig_types, 15, 0, 0, 0, 0};
+#define SWIGTYPE_p_unsigned_int swig_types[15]
+static swig_type_info *swig_types[17];
+static swig_module_info swig_module = {swig_types, 16, 0, 0, 0, 0};
 #define SWIG_TypeQuery(name) SWIG_TypeQueryModule(&swig_module, &swig_module, name)
 #define SWIG_MangledTypeQuery(name) SWIG_MangledTypeQueryModule(&swig_module, &swig_module, name)
 
@@ -3202,6 +3203,16 @@ SWIG_AsVal_unsigned_SS_int (VALUE obj, unsigned int *val)
   return res;
 }
 
+
+  inline bool is_lt_zero_after_asval(const unsigned int &i){
+    return ((sizeof(unsigned int) == sizeof(unsigned long)) && ((UINT_MAX >> 1) <= i));
+  }
+  void raise_if_lt_zero_after_asval(const unsigned int &i){
+    if(is_lt_zero_after_asval(i)){
+      SWIG_exception(SWIG_ValueError, "Expected positive value.");
+    }
+  } 
+
 SWIGINTERN double Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg____getitem__(Matrix_Frozen< double,Array2D_Dense< double >,MatViewBase > const *self,unsigned int const &row,unsigned int const &column){
     return (self)->operator()(row, column);
   }
@@ -3615,13 +3626,10 @@ SWIGINTERN Matrix< double,Array2D_Dense< double > > *new_Matrix_Sl_double_Sc_Arr
       MatrixUtil::replace(res, replacer);
       return new Matrix<double, Array2D_Dense< double >, MatrixViewBase< >>(res);
     }else if(value && rb_respond_to(*value, id_r) && rb_respond_to(*value, id_c)){
-      /* "unsigned" is remove because SWIG_AsVal(unsigned int)
-       * can not detect less than zero in Windows Ruby devkit.
-       */
-      int r, c; 
+      unsigned int r, c; 
       VALUE v_r(rb_funcall(*value, id_r, 0, 0)), v_c(rb_funcall(*value, id_c, 0, 0));
-      if(!SWIG_IsOK(SWIG_AsVal_int (v_r, &r)) || (r < 0)
-          || !SWIG_IsOK(SWIG_AsVal_int (v_c, &c)) || (c < 0)){
+      if(!SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (v_r, &r)) || is_lt_zero_after_asval(r)
+          || !SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (v_c, &c)) || is_lt_zero_after_asval(c)){
         throw std::runtime_error(
             std::string("Unexpected length [")
               .append(inspect_str(v_r)).append(", ")
@@ -3645,27 +3653,39 @@ SWIGINTERN Matrix< double,Array2D_Dense< double > > Matrix_Sl_double_Sc_Array2D_
     return Matrix<double, Array2D_Dense<double> >(
         Matrix_Frozen<double, Array2D_Dense< double >, MatrixViewBase< >>::getI(size));
   }
-SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(Matrix< double,Array2D_Dense< double > > *self,Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &matrix){
-    return self->replace(matrix);
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__swap_rows(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,unsigned int const &r1,unsigned int const &r2){
+    self->swapRows(r1, r2);
+  }
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__swap_columns(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,unsigned int const &c1,unsigned int const &c2){
+    self->swapColumns(c1, c2);
   }
-SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(Matrix< double,Array2D_Dense< double > > *self,Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &matrix){
-    return self->replace(matrix);
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &matrix){
+    self->replace(matrix);
   }
-SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(Matrix< double,Array2D_Dense< double > > *self,void const *replacer=NULL){
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &matrix){
+    self->replace(matrix);
+  }
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,void const *replacer=NULL){
     if(!MatrixUtil::replace(*self, replacer)){
       throw std::runtime_error("Unsupported replacement");
     }
-    return *self;
   }
-SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_5(Matrix< double,Array2D_Dense< double > > *self,double const *serialized){
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_5(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,double const *serialized){
     if(!MatrixUtil::replace(*self, serialized)){
       throw std::runtime_error("Unsupported replacement");
     }
-    return *self;
   }
-SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(Matrix< double,Array2D_Dense< double > > *self,void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &),MatrixUtil::each_which_t const &each_which=MatrixUtil::EACH_ALL){
+SWIGINTERN void Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(Matrix< double,Array2D_Dense< double > > *self,Matrix< double,Array2D_Dense< double > >::self_t *self_p,void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &),MatrixUtil::each_which_t const &each_which=MatrixUtil::EACH_ALL){
     MatrixUtil::each(*self, each_func, each_which, self);
-    return *self;
+  }
+SWIGINTERN Matrix< double,Array2D_Dense< double > > &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__resize(Matrix< double,Array2D_Dense< double > > *self,unsigned int const *r_p,unsigned int const *c_p){
+    unsigned int r(r_p ? *r_p : self->rows()), c(c_p ? *c_p : self->columns());
+    Matrix<double, Array2D_Dense<double > > mat_new(r, c);
+    unsigned int r_min(r), c_min(c);
+    if(r_min > self->rows()){r_min = self->rows();}
+    if(c_min > self->columns()){c_min = self->columns();}
+    mat_new.partial(r_min, c_min).replace(self->partial(r_min, c_min), false);
+    return (*(self) = mat_new);
   }
 SWIGINTERN Complex< double > Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg____getitem__(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatViewBase > const *self,unsigned int const &row,unsigned int const &column){
     return (self)->operator()(row, column);
@@ -4080,13 +4100,10 @@ SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *new_M
       MatrixUtil::replace(res, replacer);
       return new Matrix<Complex< double >, Array2D_Dense< Complex< double > >, MatrixViewBase< >>(res);
     }else if(value && rb_respond_to(*value, id_r) && rb_respond_to(*value, id_c)){
-      /* "unsigned" is remove because SWIG_AsVal(unsigned int)
-       * can not detect less than zero in Windows Ruby devkit.
-       */
-      int r, c; 
+      unsigned int r, c; 
       VALUE v_r(rb_funcall(*value, id_r, 0, 0)), v_c(rb_funcall(*value, id_c, 0, 0));
-      if(!SWIG_IsOK(SWIG_AsVal_int (v_r, &r)) || (r < 0)
-          || !SWIG_IsOK(SWIG_AsVal_int (v_c, &c)) || (c < 0)){
+      if(!SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (v_r, &r)) || is_lt_zero_after_asval(r)
+          || !SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (v_c, &c)) || is_lt_zero_after_asval(c)){
         throw std::runtime_error(
             std::string("Unexpected length [")
               .append(inspect_str(v_r)).append(", ")
@@ -4110,27 +4127,39 @@ SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > > Matrix
     return Matrix<Complex< double >, Array2D_Dense<Complex< double >> >(
         Matrix_Frozen<Complex< double >, Array2D_Dense< Complex< double > >, MatrixViewBase< >>::getI(size));
   }
-SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &matrix){
-    return self->replace(matrix);
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__swap_rows(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,unsigned int const &r1,unsigned int const &r2){
+    self->swapRows(r1, r2);
+  }
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__swap_columns(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,unsigned int const &c1,unsigned int const &c2){
+    self->swapColumns(c1, c2);
   }
-SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &matrix){
-    return self->replace(matrix);
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &matrix){
+    self->replace(matrix);
   }
-SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,void const *replacer=NULL){
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &matrix){
+    self->replace(matrix);
+  }
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,void const *replacer=NULL){
     if(!MatrixUtil::replace(*self, replacer)){
       throw std::runtime_error("Unsupported replacement");
     }
-    return *self;
   }
-SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_5(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Complex< double > const *serialized){
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_5(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,Complex< double > const *serialized){
     if(!MatrixUtil::replace(*self, serialized)){
       throw std::runtime_error("Unsupported replacement");
     }
-    return *self;
   }
-SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &),MatrixUtil::each_which_t const &each_which=MatrixUtil::EACH_ALL){
+SWIGINTERN void Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p,void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &),MatrixUtil::each_which_t const &each_which=MatrixUtil::EACH_ALL){
     MatrixUtil::each(*self, each_func, each_which, self);
-    return *self;
+  }
+SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > > &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__resize(Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *self,unsigned int const *r_p,unsigned int const *c_p){
+    unsigned int r(r_p ? *r_p : self->rows()), c(c_p ? *c_p : self->columns());
+    Matrix<Complex<double>, Array2D_Dense<Complex<double> > > mat_new(r, c);
+    unsigned int r_min(r), c_min(c);
+    if(r_min > self->rows()){r_min = self->rows();}
+    if(c_min > self->columns()){c_min = self->columns();}
+    mat_new.partial(r_min, c_min).replace(self->partial(r_min, c_min), false);
+    return (*(self) = mat_new);
   }
 static swig_class SwigClassGC_VALUE;
 
@@ -6771,6 +6800,8 @@ _wrap_Matrix_FrozenD___getitem__(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = (double)Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg____getitem__((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -6897,6 +6928,10 @@ _wrap_Matrix_FrozenD_circular__SWIG_0(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg__circular__SWIG_0((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -6950,6 +6985,8 @@ _wrap_Matrix_FrozenD_circular__SWIG_1(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg__circular__SWIG_1((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -9045,6 +9082,10 @@ _wrap_Matrix_FrozenD_partial(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg__partial((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -9097,6 +9138,7 @@ _wrap_Matrix_FrozenD_row_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg__row_vector((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2);
@@ -9149,6 +9191,7 @@ _wrap_Matrix_FrozenD_column_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatViewBase_Sg__column_vector((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2);
@@ -9796,6 +9839,8 @@ _wrap_Matrix_FrozenD_f___getitem__(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = (double)Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg____getitem__((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -9922,6 +9967,10 @@ _wrap_Matrix_FrozenD_f_circular__SWIG_0(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg__circular__SWIG_0((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -9975,6 +10024,8 @@ _wrap_Matrix_FrozenD_f_circular__SWIG_1(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg__circular__SWIG_1((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -12070,6 +12121,10 @@ _wrap_Matrix_FrozenD_f_partial(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg__partial((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -12122,6 +12177,7 @@ _wrap_Matrix_FrozenD_f_row_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg__row_vector((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2);
@@ -12174,6 +12230,7 @@ _wrap_Matrix_FrozenD_f_column_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg__column_vector((Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2);
@@ -12263,130 +12320,6 @@ free_Matrix_Frozen_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sc_MatView_f_Sg_(voi
 */
 static swig_class SwigClassMatrixD;
 
-/*
-  Document-method: GPS_PVT::SylphideMath::MatrixD.swap_rows
-
-  call-seq:
-    swap_rows(unsigned int const & row1, unsigned int const & row2) -> MatrixD
-
-An instance method.
-
-*/
-SWIGINTERN VALUE
-_wrap_MatrixD_swap_rowsN___(int argc, VALUE *argv, VALUE self) {
-  Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  unsigned int *arg2 = 0 ;
-  unsigned int *arg3 = 0 ;
-  void *argp1 = 0 ;
-  int res1 = 0 ;
-  unsigned int temp2 ;
-  unsigned int val2 ;
-  int ecode2 = 0 ;
-  unsigned int temp3 ;
-  unsigned int val3 ;
-  int ecode3 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
-  VALUE vresult = Qnil;
-  
-  if ((argc < 2) || (argc > 2)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
-  }
-  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
-  if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","swapRows", 1, self )); 
-  }
-  arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  ecode2 = SWIG_AsVal_unsigned_SS_int(argv[0], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), Ruby_Format_TypeError( "", "unsigned int","swapRows", 2, argv[0] ));
-  } 
-  temp2 = static_cast< unsigned int >(val2);
-  arg2 = &temp2;
-  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[1], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swapRows", 3, argv[1] ));
-  } 
-  temp3 = static_cast< unsigned int >(val3);
-  arg3 = &temp3;
-  {
-    try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &(arg1)->swapRows((unsigned int const &)*arg2,(unsigned int const &)*arg3);
-    } catch (const native_exception &e) {
-      e.regenerate();
-      SWIG_fail;
-    } catch (const std::exception& e) {
-      SWIG_exception_fail(SWIG_RuntimeError, e.what());
-    }
-  }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
-  return vresult;
-fail:
-  return Qnil;
-}
-
-
-/*
-  Document-method: GPS_PVT::SylphideMath::MatrixD.swap_columns
-
-  call-seq:
-    swap_columns(unsigned int const & column1, unsigned int const & column2) -> MatrixD
-
-An instance method.
-
-*/
-SWIGINTERN VALUE
-_wrap_MatrixD_swap_columnsN___(int argc, VALUE *argv, VALUE self) {
-  Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  unsigned int *arg2 = 0 ;
-  unsigned int *arg3 = 0 ;
-  void *argp1 = 0 ;
-  int res1 = 0 ;
-  unsigned int temp2 ;
-  unsigned int val2 ;
-  int ecode2 = 0 ;
-  unsigned int temp3 ;
-  unsigned int val3 ;
-  int ecode3 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
-  VALUE vresult = Qnil;
-  
-  if ((argc < 2) || (argc > 2)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
-  }
-  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
-  if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","swapColumns", 1, self )); 
-  }
-  arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  ecode2 = SWIG_AsVal_unsigned_SS_int(argv[0], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), Ruby_Format_TypeError( "", "unsigned int","swapColumns", 2, argv[0] ));
-  } 
-  temp2 = static_cast< unsigned int >(val2);
-  arg2 = &temp2;
-  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[1], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swapColumns", 3, argv[1] ));
-  } 
-  temp3 = static_cast< unsigned int >(val3);
-  arg3 = &temp3;
-  {
-    try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &(arg1)->swapColumns((unsigned int const &)*arg2,(unsigned int const &)*arg3);
-    } catch (const native_exception &e) {
-      e.regenerate();
-      SWIG_fail;
-    } catch (const std::exception& e) {
-      SWIG_exception_fail(SWIG_RuntimeError, e.what());
-    }
-  }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
-  return vresult;
-fail:
-  return Qnil;
-}
-
-
 /*
   Document-method: GPS_PVT::SylphideMath::MatrixD.new
 
@@ -12430,6 +12363,8 @@ _wrap_new_MatrixD__SWIG_0(int argc, VALUE *argv, VALUE self) {
   {
     arg3 = &argv[2];
   }
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< double,Array2D_Dense< double > > *)new_Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg___SWIG_0((unsigned int const &)*arg1,(unsigned int const &)*arg2,(void const *)arg3);
@@ -12474,6 +12409,8 @@ _wrap_new_MatrixD__SWIG_1(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< double,Array2D_Dense< double > > *)new_Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg___SWIG_0((unsigned int const &)*arg1,(unsigned int const &)*arg2);
@@ -12526,6 +12463,8 @@ _wrap_new_MatrixD__SWIG_2(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "double const *","Matrix<(double,Array2D_Dense<(double)>)>", 3, argv[2] )); 
   }
   arg3 = reinterpret_cast< double * >(argp3);
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< double,Array2D_Dense< double > > *)new_Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg___SWIG_2((unsigned int const &)*arg1,(unsigned int const &)*arg2,(double const *)arg3);
@@ -12735,6 +12674,8 @@ _wrap_MatrixD___setitem__(int argc, VALUE *argv, VALUE self) {
   } 
   temp4 = static_cast< double >(val4);
   arg4 = &temp4;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = (double *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg____setitem__(arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(double const &)*arg4);
@@ -12789,6 +12730,7 @@ _wrap_MatrixD_scalar(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< double >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg1);
   {
     try {
       result = Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__getScalar((unsigned int const &)*arg1,(double const &)*arg2);
@@ -12833,6 +12775,7 @@ _wrap_MatrixD_I(int argc, VALUE *argv, VALUE self) {
   } 
   temp1 = static_cast< unsigned int >(val1);
   arg1 = &temp1;
+  raise_if_lt_zero_after_asval(*arg1);
   {
     try {
       result = Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__getI((unsigned int const &)*arg1);
@@ -12851,47 +12794,56 @@ fail:
 
 
 /*
-  Document-method: GPS_PVT::SylphideMath::MatrixD.replace!
+  Document-method: GPS_PVT::SylphideMath::MatrixD.swap_rows
 
   call-seq:
-    replace!(matrix) -> MatrixD
-    replace!(matrix) -> MatrixD
-    replace!(replacer=nil) -> MatrixD
-    replace! -> MatrixD
-    replace!(serialized) -> MatrixD
+    swap_rows(unsigned int const & r1, unsigned int const & r2)
+
+An instance method.
 
-Replace all or a portion of the MatrixD.
 */
 SWIGINTERN VALUE
-_wrap_MatrixD_replaceN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
+_wrap_MatrixD_swap_rowsN___(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > *arg2 = 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  unsigned int *arg3 = 0 ;
+  unsigned int *arg4 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 ;
-  int res2 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
+  unsigned int temp3 ;
+  unsigned int val3 ;
+  int ecode3 = 0 ;
+  unsigned int temp4 ;
+  unsigned int val4 ;
+  int ecode4 = 0 ;
   VALUE vresult = Qnil;
   
-  if ((argc < 1) || (argc > 1)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
   }
   res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
   if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 1, self )); 
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","swap_rows", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2, SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t,  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 2, argv[0] )); 
-  }
-  if (!argp2) {
-    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 2, argv[0])); 
-  }
-  arg2 = reinterpret_cast< Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > * >(argp2);
+  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[0], &val3);
+  if (!SWIG_IsOK(ecode3)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swap_rows", 3, argv[0] ));
+  } 
+  temp3 = static_cast< unsigned int >(val3);
+  arg3 = &temp3;
+  ecode4 = SWIG_AsVal_unsigned_SS_int(argv[1], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), Ruby_Format_TypeError( "", "unsigned int","swap_rows", 4, argv[1] ));
+  } 
+  temp4 = static_cast< unsigned int >(val4);
+  arg4 = &temp4;
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(arg1,(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &)*arg2);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__swap_rows(arg1,arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -12906,36 +12858,57 @@ fail:
 }
 
 
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixD.swap_columns
+
+  call-seq:
+    swap_columns(unsigned int const & c1, unsigned int const & c2)
+
+An instance method.
+
+*/
 SWIGINTERN VALUE
-_wrap_MatrixD_replaceN_____SWIG_2(int argc, VALUE *argv, VALUE self) {
+_wrap_MatrixD_swap_columnsN___(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > *arg2 = 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  unsigned int *arg3 = 0 ;
+  unsigned int *arg4 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 ;
-  int res2 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
+  unsigned int temp3 ;
+  unsigned int val3 ;
+  int ecode3 = 0 ;
+  unsigned int temp4 ;
+  unsigned int val4 ;
+  int ecode4 = 0 ;
   VALUE vresult = Qnil;
   
-  if ((argc < 1) || (argc > 1)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
   }
   res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
   if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 1, self )); 
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","swap_columns", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2, SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t,  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 2, argv[0] )); 
-  }
-  if (!argp2) {
-    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 2, argv[0])); 
-  }
-  arg2 = reinterpret_cast< Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > * >(argp2);
+  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[0], &val3);
+  if (!SWIG_IsOK(ecode3)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swap_columns", 3, argv[0] ));
+  } 
+  temp3 = static_cast< unsigned int >(val3);
+  arg3 = &temp3;
+  ecode4 = SWIG_AsVal_unsigned_SS_int(argv[1], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), Ruby_Format_TypeError( "", "unsigned int","swap_columns", 4, argv[1] ));
+  } 
+  temp4 = static_cast< unsigned int >(val4);
+  arg4 = &temp4;
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(arg1,(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &)*arg2);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__swap_columns(arg1,arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -12950,29 +12923,49 @@ fail:
 }
 
 
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixD.replace!
+
+  call-seq:
+    replace!(matrix)
+    replace!(matrix)
+    replace!(replacer=nil)
+    replace!
+    replace!(serialized)
+
+Replace all or a portion of the MatrixD.
+*/
 SWIGINTERN VALUE
-_wrap_MatrixD_replaceN_____SWIG_3(int argc, VALUE *argv, VALUE self) {
+_wrap_MatrixD_replaceN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  void *arg2 = (void *) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > *arg3 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
+  void *argp3 ;
+  int res3 = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
   res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
   if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace", 1, self )); 
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  {
-    arg2 = &argv[0];
+  res3 = SWIG_ConvertPtr(argv[0], &argp3, SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t,  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 3, argv[0] )); 
+  }
+  if (!argp3) {
+    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 3, argv[0])); 
   }
+  arg3 = reinterpret_cast< Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > * >(argp3);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(arg1,(void const *)arg2);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(arg1,arg2,(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &)*arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -12988,24 +12981,36 @@ fail:
 
 
 SWIGINTERN VALUE
-_wrap_MatrixD_replaceN_____SWIG_4(int argc, VALUE *argv, VALUE self) {
+_wrap_MatrixD_replaceN_____SWIG_2(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > *arg3 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
+  void *argp3 ;
+  int res3 = 0 ;
   VALUE vresult = Qnil;
   
-  if ((argc < 0) || (argc > 0)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
+  
+  if ((argc < 1) || (argc > 1)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
   res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
   if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace", 1, self )); 
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
+  res3 = SWIG_ConvertPtr(argv[0], &argp3, SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t,  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 3, argv[0] )); 
+  }
+  if (!argp3) {
+    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 3, argv[0])); 
+  }
+  arg3 = reinterpret_cast< Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > * >(argp3);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(arg1);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(arg1,arg2,(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &)*arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -13021,16 +13026,15 @@ fail:
 
 
 SWIGINTERN VALUE
-_wrap_MatrixD_replaceN_____SWIG_5(int argc, VALUE *argv, VALUE self) {
+_wrap_MatrixD_replaceN_____SWIG_3(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  double *arg2 = (double *) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  void *arg3 = (void *) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 = 0 ;
-  int res2 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
@@ -13039,14 +13043,12 @@ _wrap_MatrixD_replaceN_____SWIG_5(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2,SWIGTYPE_p_double, 0 |  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "double const *","replace", 2, argv[0] )); 
+  {
+    arg3 = &argv[0];
   }
-  arg2 = reinterpret_cast< double * >(argp2);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_5(arg1,(double const *)arg2);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(arg1,arg2,(void const *)arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -13061,13 +13063,89 @@ fail:
 }
 
 
-SWIGINTERN VALUE _wrap_MatrixD_replaceN___(int nargs, VALUE *args, VALUE self) {
-  int argc;
-  VALUE argv[3];
-  int ii;
+SWIGINTERN VALUE
+_wrap_MatrixD_replaceN_____SWIG_4(int argc, VALUE *argv, VALUE self) {
+  Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  VALUE vresult = Qnil;
   
-  argc = nargs + 1;
-  argv[0] = self;
+  
+  if ((argc < 0) || (argc > 0)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
+  {
+    try {
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_3(arg1,arg2);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = self;
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
+SWIGINTERN VALUE
+_wrap_MatrixD_replaceN_____SWIG_5(int argc, VALUE *argv, VALUE self) {
+  Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  double *arg3 = (double *) 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  void *argp3 = 0 ;
+  int res3 = 0 ;
+  VALUE vresult = Qnil;
+  
+  
+  if ((argc < 1) || (argc > 1)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","replace", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
+  res3 = SWIG_ConvertPtr(argv[0], &argp3,SWIGTYPE_p_double, 0 |  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "double const *","replace", 3, argv[0] )); 
+  }
+  arg3 = reinterpret_cast< double * >(argp3);
+  {
+    try {
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__replace__SWIG_5(arg1,arg2,(double const *)arg3);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = self;
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
+SWIGINTERN VALUE _wrap_MatrixD_replaceN___(int nargs, VALUE *args, VALUE self) {
+  int argc;
+  VALUE argv[3];
+  int ii;
+  
+  argc = nargs + 1;
+  argv[0] = self;
   if (argc > 3) SWIG_fail;
   for (ii = 1; (ii < argc); ++ii) {
     argv[ii] = args[ii-1];
@@ -13143,11 +13221,11 @@ SWIGINTERN VALUE _wrap_MatrixD_replaceN___(int nargs, VALUE *args, VALUE self) {
   
 fail:
   Ruby_Format_OverloadedError( argc, 3, "replace!", 
-    "    Matrix< double,Array2D_Dense< double > >::self_t replace!(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &matrix)\n"
-    "    Matrix< double,Array2D_Dense< double > >::self_t replace!(Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &matrix)\n"
-    "    Matrix< double,Array2D_Dense< double > >::self_t replace!(void const *replacer)\n"
-    "    Matrix< double,Array2D_Dense< double > >::self_t replace!()\n"
-    "    Matrix< double,Array2D_Dense< double > >::self_t & replace!(double const *serialized)\n");
+    "    void replace!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewBase< > > const &matrix)\n"
+    "    void replace!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, Matrix_Frozen< double,Array2D_Dense< double >,MatrixViewFilter< MatrixViewBase< > > > const &matrix)\n"
+    "    void replace!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, void const *replacer)\n"
+    "    void replace!(Matrix< double,Array2D_Dense< double > >::self_t *self_p)\n"
+    "    void replace!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, double const *serialized)\n");
   
   return Qnil;
 }
@@ -13157,8 +13235,8 @@ fail:
   Document-method: GPS_PVT::SylphideMath::MatrixD.map!
 
   call-seq:
-    map!(MatrixUtil::each_which_t const & each_which=MatrixUtil::EACH_ALL) -> MatrixD
-    map! -> MatrixD
+    map!(MatrixUtil::each_which_t const & each_which=MatrixUtil::EACH_ALL)
+    map!
 
 An instance method.
 
@@ -13166,18 +13244,19 @@ An instance method.
 SWIGINTERN VALUE
 _wrap_MatrixD_mapN_____SWIG_0(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  void (*arg2)(double const &,double *,unsigned int const &,unsigned int const &) = (void (*)(double const &,double *,unsigned int const &,unsigned int const &)) 0 ;
-  MatrixUtil::each_which_t *arg3 = 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  void (*arg3)(double const &,double *,unsigned int const &,unsigned int const &) = (void (*)(double const &,double *,unsigned int const &,unsigned int const &)) 0 ;
+  MatrixUtil::each_which_t *arg4 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   {
     if(!rb_block_given_p()){
       return rb_enumeratorize(self, ID2SYM(rb_frame_callee()), argc, argv);
     }
-    arg2 = matrix_each((const double *)0);
+    arg3 = matrix_each((const double *)0);
   }
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -13189,14 +13268,14 @@ _wrap_MatrixD_mapN_____SWIG_0(int argc, VALUE *argv, VALUE self) {
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
   {
     try{
-      arg3 = &const_cast<typename MatrixUtil::each_which_t &>(MatrixUtil::sym2each_which(argv[0]));
+      arg4 = &const_cast<typename MatrixUtil::each_which_t &>(MatrixUtil::sym2each_which(argv[0]));
     }catch(std::runtime_error &e){
       SWIG_exception(SWIG_TypeError, e.what());
     }
   }
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(arg1,arg2,(MatrixUtil::each_which_t const &)*arg3);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(arg1,arg2,arg3,(MatrixUtil::each_which_t const &)*arg4);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -13214,17 +13293,18 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixD_mapN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
-  void (*arg2)(double const &,double *,unsigned int const &,unsigned int const &) = (void (*)(double const &,double *,unsigned int const &,unsigned int const &)) 0 ;
+  Matrix< double,Array2D_Dense< double > >::self_t *arg2 = (Matrix< double,Array2D_Dense< double > >::self_t *) 0 ;
+  void (*arg3)(double const &,double *,unsigned int const &,unsigned int const &) = (void (*)(double const &,double *,unsigned int const &,unsigned int const &)) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< double,Array2D_Dense< double > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   {
     if(!rb_block_given_p()){
       return rb_enumeratorize(self, ID2SYM(rb_frame_callee()), argc, argv);
     }
-    arg2 = matrix_each((const double *)0);
+    arg3 = matrix_each((const double *)0);
   }
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
@@ -13236,7 +13316,7 @@ _wrap_MatrixD_mapN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
   {
     try {
-      result = (Matrix< double,Array2D_Dense< double > >::self_t *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(arg1,arg2);
+      Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__map_bang__SWIG_0(arg1,arg2,arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -13288,13 +13368,88 @@ SWIGINTERN VALUE _wrap_MatrixD_mapN___(int nargs, VALUE *args, VALUE self) {
   
 fail:
   Ruby_Format_OverloadedError( argc, 3, "map!", 
-    "    Matrix< double,Array2D_Dense< double > >::self_t map!(void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &), MatrixUtil::each_which_t const &each_which)\n"
-    "    Matrix< double,Array2D_Dense< double > >::self_t & map!(void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &))\n");
+    "    void map!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &), MatrixUtil::each_which_t const &each_which)\n"
+    "    void map!(Matrix< double,Array2D_Dense< double > >::self_t *self_p, void (*each_func)(double const &,double *,unsigned int const &,unsigned int const &))\n");
   
   return Qnil;
 }
 
 
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixD.resize
+
+  call-seq:
+    resize(r_p, c_p) -> MatrixD
+
+Resize the size of the MatrixD.
+*/
+SWIGINTERN VALUE
+_wrap_MatrixD_resizeN___(int argc, VALUE *argv, VALUE self) {
+  Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *) 0 ;
+  unsigned int *arg2 = (unsigned int *) 0 ;
+  unsigned int *arg3 = (unsigned int *) 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  unsigned int temp2 ;
+  unsigned int temp3 ;
+  Matrix< double,Array2D_Dense< double > > *result = 0 ;
+  VALUE vresult = Qnil;
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< double,Array2D_Dense< double > > *","resize", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< double,Array2D_Dense< double > > * >(argp1);
+  {
+    if(SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (argv[0], &temp2))){
+      raise_if_lt_zero_after_asval(temp2);
+      
+      arg2 = &temp2;
+    }
+    
+    else if(NIL_P(argv[0])){
+      arg2 = NULL;
+    }
+    
+    else{
+      SWIG_exception(SWIG_TypeError, "unsigned int is expected");
+    }
+  }
+  {
+    if(SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (argv[1], &temp3))){
+      raise_if_lt_zero_after_asval(temp3);
+      
+      arg3 = &temp3;
+    }
+    
+    else if(NIL_P(argv[1])){
+      arg3 = NULL;
+    }
+    
+    else{
+      SWIG_exception(SWIG_TypeError, "unsigned int is expected");
+    }
+  }
+  {
+    try {
+      result = (Matrix< double,Array2D_Dense< double > > *) &Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg__resize(arg1,(unsigned int const *)arg2,(unsigned int const *)arg3);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
 SWIGINTERN void
 free_Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg_(void *self) {
     Matrix< double,Array2D_Dense< double > > *arg1 = (Matrix< double,Array2D_Dense< double > > *)self;
@@ -13878,6 +14033,8 @@ _wrap_Matrix_FrozenComplexD___getitem__(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg____getitem__((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -14006,6 +14163,10 @@ _wrap_Matrix_FrozenComplexD_circular__SWIG_0(int argc, VALUE *argv, VALUE self)
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg__circular__SWIG_0((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -14059,6 +14220,8 @@ _wrap_Matrix_FrozenComplexD_circular__SWIG_1(int argc, VALUE *argv, VALUE self)
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg__circular__SWIG_1((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -16142,6 +16305,10 @@ _wrap_Matrix_FrozenComplexD_partial(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg__partial((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -16194,6 +16361,7 @@ _wrap_Matrix_FrozenComplexD_row_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg__row_vector((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2);
@@ -16246,6 +16414,7 @@ _wrap_Matrix_FrozenComplexD_column_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatViewBase_Sg__column_vector((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const *)arg1,(unsigned int const &)*arg2);
@@ -16903,6 +17072,8 @@ _wrap_Matrix_FrozenComplexD_f___getitem__(int argc, VALUE *argv, VALUE self) {
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg____getitem__((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -17031,6 +17202,10 @@ _wrap_Matrix_FrozenComplexD_f_circular__SWIG_0(int argc, VALUE *argv, VALUE self
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg__circular__SWIG_0((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -17084,6 +17259,8 @@ _wrap_Matrix_FrozenComplexD_f_circular__SWIG_1(int argc, VALUE *argv, VALUE self
   } 
   temp3 = static_cast< unsigned int >(val3);
   arg3 = &temp3;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg__circular__SWIG_1((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -19167,6 +19344,10 @@ _wrap_Matrix_FrozenComplexD_f_partial(int argc, VALUE *argv, VALUE self) {
   } 
   temp5 = static_cast< unsigned int >(val5);
   arg5 = &temp5;
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  raise_if_lt_zero_after_asval(*arg5);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg__partial((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4,(unsigned int const &)*arg5);
@@ -19219,6 +19400,7 @@ _wrap_Matrix_FrozenComplexD_f_row_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg__row_vector((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2);
@@ -19271,6 +19453,7 @@ _wrap_Matrix_FrozenComplexD_f_column_vector(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sc_MatView_f_Sg__column_vector((Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const *)arg1,(unsigned int const &)*arg2);
@@ -19360,130 +19543,6 @@ free_Matrix_Frozen_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_doubl
 */
 static swig_class SwigClassMatrixComplexD;
 
-/*
-  Document-method: GPS_PVT::SylphideMath::MatrixComplexD.swap_rows
-
-  call-seq:
-    swap_rows(unsigned int const & row1, unsigned int const & row2) -> MatrixComplexD
-
-An instance method.
-
-*/
-SWIGINTERN VALUE
-_wrap_MatrixComplexD_swap_rowsN___(int argc, VALUE *argv, VALUE self) {
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  unsigned int *arg2 = 0 ;
-  unsigned int *arg3 = 0 ;
-  void *argp1 = 0 ;
-  int res1 = 0 ;
-  unsigned int temp2 ;
-  unsigned int val2 ;
-  int ecode2 = 0 ;
-  unsigned int temp3 ;
-  unsigned int val3 ;
-  int ecode3 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
-  VALUE vresult = Qnil;
-  
-  if ((argc < 2) || (argc > 2)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
-  }
-  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
-  if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","swapRows", 1, self )); 
-  }
-  arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
-  ecode2 = SWIG_AsVal_unsigned_SS_int(argv[0], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), Ruby_Format_TypeError( "", "unsigned int","swapRows", 2, argv[0] ));
-  } 
-  temp2 = static_cast< unsigned int >(val2);
-  arg2 = &temp2;
-  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[1], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swapRows", 3, argv[1] ));
-  } 
-  temp3 = static_cast< unsigned int >(val3);
-  arg3 = &temp3;
-  {
-    try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &(arg1)->swapRows((unsigned int const &)*arg2,(unsigned int const &)*arg3);
-    } catch (const native_exception &e) {
-      e.regenerate();
-      SWIG_fail;
-    } catch (const std::exception& e) {
-      SWIG_exception_fail(SWIG_RuntimeError, e.what());
-    }
-  }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
-  return vresult;
-fail:
-  return Qnil;
-}
-
-
-/*
-  Document-method: GPS_PVT::SylphideMath::MatrixComplexD.swap_columns
-
-  call-seq:
-    swap_columns(unsigned int const & column1, unsigned int const & column2) -> MatrixComplexD
-
-An instance method.
-
-*/
-SWIGINTERN VALUE
-_wrap_MatrixComplexD_swap_columnsN___(int argc, VALUE *argv, VALUE self) {
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  unsigned int *arg2 = 0 ;
-  unsigned int *arg3 = 0 ;
-  void *argp1 = 0 ;
-  int res1 = 0 ;
-  unsigned int temp2 ;
-  unsigned int val2 ;
-  int ecode2 = 0 ;
-  unsigned int temp3 ;
-  unsigned int val3 ;
-  int ecode3 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
-  VALUE vresult = Qnil;
-  
-  if ((argc < 2) || (argc > 2)) {
-    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
-  }
-  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
-  if (!SWIG_IsOK(res1)) {
-    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","swapColumns", 1, self )); 
-  }
-  arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
-  ecode2 = SWIG_AsVal_unsigned_SS_int(argv[0], &val2);
-  if (!SWIG_IsOK(ecode2)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode2), Ruby_Format_TypeError( "", "unsigned int","swapColumns", 2, argv[0] ));
-  } 
-  temp2 = static_cast< unsigned int >(val2);
-  arg2 = &temp2;
-  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[1], &val3);
-  if (!SWIG_IsOK(ecode3)) {
-    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swapColumns", 3, argv[1] ));
-  } 
-  temp3 = static_cast< unsigned int >(val3);
-  arg3 = &temp3;
-  {
-    try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &(arg1)->swapColumns((unsigned int const &)*arg2,(unsigned int const &)*arg3);
-    } catch (const native_exception &e) {
-      e.regenerate();
-      SWIG_fail;
-    } catch (const std::exception& e) {
-      SWIG_exception_fail(SWIG_RuntimeError, e.what());
-    }
-  }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
-  return vresult;
-fail:
-  return Qnil;
-}
-
-
 /*
   Document-method: GPS_PVT::SylphideMath::MatrixComplexD.new
 
@@ -19527,6 +19586,8 @@ _wrap_new_MatrixComplexD__SWIG_0(int argc, VALUE *argv, VALUE self) {
   {
     arg3 = &argv[2];
   }
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *)new_Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg___SWIG_0((unsigned int const &)*arg1,(unsigned int const &)*arg2,(void const *)arg3);
@@ -19571,6 +19632,8 @@ _wrap_new_MatrixComplexD__SWIG_1(int argc, VALUE *argv, VALUE self) {
   } 
   temp2 = static_cast< unsigned int >(val2);
   arg2 = &temp2;
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *)new_Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg___SWIG_0((unsigned int const &)*arg1,(unsigned int const &)*arg2);
@@ -19623,6 +19686,8 @@ _wrap_new_MatrixComplexD__SWIG_2(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Complex< double > const *","Matrix<(Complex<(double)>,Array2D_Dense<(Complex<(double)>)>)>", 3, argv[2] )); 
   }
   arg3 = reinterpret_cast< Complex< double > * >(argp3);
+  raise_if_lt_zero_after_asval(*arg1);
+  raise_if_lt_zero_after_asval(*arg2);
   {
     try {
       result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *)new_Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg___SWIG_2((unsigned int const &)*arg1,(unsigned int const &)*arg2,(Complex< double > const *)arg3);
@@ -19830,6 +19895,8 @@ _wrap_MatrixComplexD___setitem__(int argc, VALUE *argv, VALUE self) {
       SWIG_exception(SWIG_TypeError, "in method '__setitem__', expecting type Complex< double >");
     }
   }
+  raise_if_lt_zero_after_asval(*arg2);
+  raise_if_lt_zero_after_asval(*arg3);
   {
     try {
       result = (Complex< double > *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg____setitem__(arg1,(unsigned int const &)*arg2,(unsigned int const &)*arg3,(Complex< double > const &)*arg4);
@@ -19882,6 +19949,7 @@ _wrap_MatrixComplexD_scalar(int argc, VALUE *argv, VALUE self) {
       SWIG_exception(SWIG_TypeError, "in method 'Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__getScalar', expecting type Complex< double >");
     }
   }
+  raise_if_lt_zero_after_asval(*arg1);
   {
     try {
       result = Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__getScalar((unsigned int const &)*arg1,(Complex< double > const &)*arg2);
@@ -19926,6 +19994,7 @@ _wrap_MatrixComplexD_I(int argc, VALUE *argv, VALUE self) {
   } 
   temp1 = static_cast< unsigned int >(val1);
   arg1 = &temp1;
+  raise_if_lt_zero_after_asval(*arg1);
   {
     try {
       result = Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__getI((unsigned int const &)*arg1);
@@ -19943,29 +20012,160 @@ fail:
 }
 
 
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixComplexD.swap_rows
+
+  call-seq:
+    swap_rows(unsigned int const & r1, unsigned int const & r2)
+
+An instance method.
+
+*/
+SWIGINTERN VALUE
+_wrap_MatrixComplexD_swap_rowsN___(int argc, VALUE *argv, VALUE self) {
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  unsigned int *arg3 = 0 ;
+  unsigned int *arg4 = 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  unsigned int temp3 ;
+  unsigned int val3 ;
+  int ecode3 = 0 ;
+  unsigned int temp4 ;
+  unsigned int val4 ;
+  int ecode4 = 0 ;
+  VALUE vresult = Qnil;
+  
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","swap_rows", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
+  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[0], &val3);
+  if (!SWIG_IsOK(ecode3)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swap_rows", 3, argv[0] ));
+  } 
+  temp3 = static_cast< unsigned int >(val3);
+  arg3 = &temp3;
+  ecode4 = SWIG_AsVal_unsigned_SS_int(argv[1], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), Ruby_Format_TypeError( "", "unsigned int","swap_rows", 4, argv[1] ));
+  } 
+  temp4 = static_cast< unsigned int >(val4);
+  arg4 = &temp4;
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  {
+    try {
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__swap_rows(arg1,arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = self;
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixComplexD.swap_columns
+
+  call-seq:
+    swap_columns(unsigned int const & c1, unsigned int const & c2)
+
+An instance method.
+
+*/
+SWIGINTERN VALUE
+_wrap_MatrixComplexD_swap_columnsN___(int argc, VALUE *argv, VALUE self) {
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  unsigned int *arg3 = 0 ;
+  unsigned int *arg4 = 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  unsigned int temp3 ;
+  unsigned int val3 ;
+  int ecode3 = 0 ;
+  unsigned int temp4 ;
+  unsigned int val4 ;
+  int ecode4 = 0 ;
+  VALUE vresult = Qnil;
+  
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","swap_columns", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
+  ecode3 = SWIG_AsVal_unsigned_SS_int(argv[0], &val3);
+  if (!SWIG_IsOK(ecode3)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode3), Ruby_Format_TypeError( "", "unsigned int","swap_columns", 3, argv[0] ));
+  } 
+  temp3 = static_cast< unsigned int >(val3);
+  arg3 = &temp3;
+  ecode4 = SWIG_AsVal_unsigned_SS_int(argv[1], &val4);
+  if (!SWIG_IsOK(ecode4)) {
+    SWIG_exception_fail(SWIG_ArgError(ecode4), Ruby_Format_TypeError( "", "unsigned int","swap_columns", 4, argv[1] ));
+  } 
+  temp4 = static_cast< unsigned int >(val4);
+  arg4 = &temp4;
+  raise_if_lt_zero_after_asval(*arg3);
+  raise_if_lt_zero_after_asval(*arg4);
+  {
+    try {
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__swap_columns(arg1,arg2,(unsigned int const &)*arg3,(unsigned int const &)*arg4);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = self;
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
 /*
   Document-method: GPS_PVT::SylphideMath::MatrixComplexD.replace!
 
   call-seq:
-    replace!(matrix) -> MatrixComplexD
-    replace!(matrix) -> MatrixComplexD
-    replace!(replacer=nil) -> MatrixComplexD
-    replace! -> MatrixComplexD
-    replace!(serialized) -> MatrixComplexD
+    replace!(matrix)
+    replace!(matrix)
+    replace!(replacer=nil)
+    replace!
+    replace!(serialized)
 
 Replace all or a portion of the MatrixComplexD.
 */
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_replaceN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > *arg2 = 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > *arg3 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 ;
-  int res2 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
+  void *argp3 ;
+  int res3 = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
@@ -19974,17 +20174,17 @@ _wrap_MatrixComplexD_replaceN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2, SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t,  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 2, argv[0] )); 
+  res3 = SWIG_ConvertPtr(argv[0], &argp3, SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t,  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 3, argv[0] )); 
   }
-  if (!argp2) {
-    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 2, argv[0])); 
+  if (!argp3) {
+    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &","replace<(T,Array2D_Dense<(T)>,MatViewBase)>", 3, argv[0])); 
   }
-  arg2 = reinterpret_cast< Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > * >(argp2);
+  arg3 = reinterpret_cast< Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > * >(argp3);
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(arg1,(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &)*arg2);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatViewBase_Sg___SWIG_1(arg1,arg2,(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &)*arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20002,14 +20202,15 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_replaceN_____SWIG_2(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > *arg2 = 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > *arg3 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 ;
-  int res2 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
+  void *argp3 ;
+  int res3 = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
@@ -20018,17 +20219,17 @@ _wrap_MatrixComplexD_replaceN_____SWIG_2(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2, SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t,  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 2, argv[0] )); 
+  res3 = SWIG_ConvertPtr(argv[0], &argp3, SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t,  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 3, argv[0] )); 
   }
-  if (!argp2) {
-    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 2, argv[0])); 
+  if (!argp3) {
+    SWIG_exception_fail(SWIG_ValueError, Ruby_Format_TypeError("invalid null reference ", "Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &","replace<(T,Array2D_Dense<(T)>,MatView_f)>", 3, argv[0])); 
   }
-  arg2 = reinterpret_cast< Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > * >(argp2);
+  arg3 = reinterpret_cast< Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > * >(argp3);
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(arg1,(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &)*arg2);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace_Sl_T_Sc_Array2D_Dense_Sl_T_Sg__Sc_MatView_f_Sg___SWIG_2(arg1,arg2,(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &)*arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20046,12 +20247,13 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_replaceN_____SWIG_3(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  void *arg2 = (void *) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  void *arg3 = (void *) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
@@ -20061,11 +20263,11 @@ _wrap_MatrixComplexD_replaceN_____SWIG_3(int argc, VALUE *argv, VALUE self) {
   }
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
   {
-    arg2 = &argv[0];
+    arg3 = &argv[0];
   }
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(arg1,(void const *)arg2);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(arg1,arg2,(void const *)arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20083,11 +20285,12 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_replaceN_____SWIG_4(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
   }
@@ -20098,7 +20301,7 @@ _wrap_MatrixComplexD_replaceN_____SWIG_4(int argc, VALUE *argv, VALUE self) {
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(arg1);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_3(arg1,arg2);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20116,14 +20319,15 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_replaceN_____SWIG_5(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  Complex< double > *arg2 = (Complex< double > *) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  Complex< double > *arg3 = (Complex< double > *) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  void *argp2 = 0 ;
-  int res2 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
+  void *argp3 = 0 ;
+  int res3 = 0 ;
   VALUE vresult = Qnil;
   
+  
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
   }
@@ -20132,14 +20336,14 @@ _wrap_MatrixComplexD_replaceN_____SWIG_5(int argc, VALUE *argv, VALUE self) {
     SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","replace", 1, self )); 
   }
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
-  res2 = SWIG_ConvertPtr(argv[0], &argp2,SWIGTYPE_p_ComplexT_double_t, 0 |  0 );
-  if (!SWIG_IsOK(res2)) {
-    SWIG_exception_fail(SWIG_ArgError(res2), Ruby_Format_TypeError( "", "Complex< double > const *","replace", 2, argv[0] )); 
+  res3 = SWIG_ConvertPtr(argv[0], &argp3,SWIGTYPE_p_ComplexT_double_t, 0 |  0 );
+  if (!SWIG_IsOK(res3)) {
+    SWIG_exception_fail(SWIG_ArgError(res3), Ruby_Format_TypeError( "", "Complex< double > const *","replace", 3, argv[0] )); 
   }
-  arg2 = reinterpret_cast< Complex< double > * >(argp2);
+  arg3 = reinterpret_cast< Complex< double > * >(argp3);
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_5(arg1,(Complex< double > const *)arg2);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__replace__SWIG_5(arg1,arg2,(Complex< double > const *)arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20236,11 +20440,11 @@ SWIGINTERN VALUE _wrap_MatrixComplexD_replaceN___(int nargs, VALUE *args, VALUE
   
 fail:
   Ruby_Format_OverloadedError( argc, 3, "replace!", 
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t replace!(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &matrix)\n"
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t replace!(Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &matrix)\n"
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t replace!(void const *replacer)\n"
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t replace!()\n"
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t & replace!(Complex< double > const *serialized)\n");
+    "    void replace!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewBase< > > const &matrix)\n"
+    "    void replace!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, Matrix_Frozen< Complex< double >,Array2D_Dense< Complex< double > >,MatrixViewFilter< MatrixViewBase< > > > const &matrix)\n"
+    "    void replace!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, void const *replacer)\n"
+    "    void replace!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p)\n"
+    "    void replace!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, Complex< double > const *serialized)\n");
   
   return Qnil;
 }
@@ -20250,8 +20454,8 @@ fail:
   Document-method: GPS_PVT::SylphideMath::MatrixComplexD.map!
 
   call-seq:
-    map!(MatrixUtil::each_which_t const & each_which=MatrixUtil::EACH_ALL) -> MatrixComplexD
-    map! -> MatrixComplexD
+    map!(MatrixUtil::each_which_t const & each_which=MatrixUtil::EACH_ALL)
+    map!
 
 An instance method.
 
@@ -20259,18 +20463,19 @@ An instance method.
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_mapN_____SWIG_0(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  void (*arg2)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &) = (void (*)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &)) 0 ;
-  MatrixUtil::each_which_t *arg3 = 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  void (*arg3)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &) = (void (*)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &)) 0 ;
+  MatrixUtil::each_which_t *arg4 = 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   {
     if(!rb_block_given_p()){
       return rb_enumeratorize(self, ID2SYM(rb_frame_callee()), argc, argv);
     }
-    arg2 = matrix_each((const Complex< double > *)0);
+    arg3 = matrix_each((const Complex< double > *)0);
   }
   if ((argc < 1) || (argc > 1)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 1)",argc); SWIG_fail;
@@ -20282,14 +20487,14 @@ _wrap_MatrixComplexD_mapN_____SWIG_0(int argc, VALUE *argv, VALUE self) {
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
   {
     try{
-      arg3 = &const_cast<typename MatrixUtil::each_which_t &>(MatrixUtil::sym2each_which(argv[0]));
+      arg4 = &const_cast<typename MatrixUtil::each_which_t &>(MatrixUtil::sym2each_which(argv[0]));
     }catch(std::runtime_error &e){
       SWIG_exception(SWIG_TypeError, e.what());
     }
   }
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(arg1,arg2,(MatrixUtil::each_which_t const &)*arg3);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(arg1,arg2,arg3,(MatrixUtil::each_which_t const &)*arg4);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20307,17 +20512,18 @@ fail:
 SWIGINTERN VALUE
 _wrap_MatrixComplexD_mapN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
-  void (*arg2)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &) = (void (*)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &)) 0 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *arg2 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) 0 ;
+  void (*arg3)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &) = (void (*)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &)) 0 ;
   void *argp1 = 0 ;
   int res1 = 0 ;
-  Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *result = 0 ;
   VALUE vresult = Qnil;
   
+  
   {
     if(!rb_block_given_p()){
       return rb_enumeratorize(self, ID2SYM(rb_frame_callee()), argc, argv);
     }
-    arg2 = matrix_each((const Complex< double > *)0);
+    arg3 = matrix_each((const Complex< double > *)0);
   }
   if ((argc < 0) || (argc > 0)) {
     rb_raise(rb_eArgError, "wrong # of arguments(%d for 0)",argc); SWIG_fail;
@@ -20329,7 +20535,7 @@ _wrap_MatrixComplexD_mapN_____SWIG_1(int argc, VALUE *argv, VALUE self) {
   arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
   {
     try {
-      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(arg1,arg2);
+      Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__map_bang__SWIG_0(arg1,arg2,arg3);
     } catch (const native_exception &e) {
       e.regenerate();
       SWIG_fail;
@@ -20381,13 +20587,88 @@ SWIGINTERN VALUE _wrap_MatrixComplexD_mapN___(int nargs, VALUE *args, VALUE self
   
 fail:
   Ruby_Format_OverloadedError( argc, 3, "map!", 
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t map!(void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &), MatrixUtil::each_which_t const &each_which)\n"
-    "    Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t & map!(void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &))\n");
+    "    void map!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &), MatrixUtil::each_which_t const &each_which)\n"
+    "    void map!(Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *self_p, void (*each_func)(Complex< double > const &,Complex< double > *,unsigned int const &,unsigned int const &))\n");
   
   return Qnil;
 }
 
 
+/*
+  Document-method: GPS_PVT::SylphideMath::MatrixComplexD.resize
+
+  call-seq:
+    resize(r_p, c_p) -> MatrixComplexD
+
+Resize the size of the MatrixComplexD.
+*/
+SWIGINTERN VALUE
+_wrap_MatrixComplexD_resizeN___(int argc, VALUE *argv, VALUE self) {
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) 0 ;
+  unsigned int *arg2 = (unsigned int *) 0 ;
+  unsigned int *arg3 = (unsigned int *) 0 ;
+  void *argp1 = 0 ;
+  int res1 = 0 ;
+  unsigned int temp2 ;
+  unsigned int temp3 ;
+  Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *result = 0 ;
+  VALUE vresult = Qnil;
+  
+  if ((argc < 2) || (argc > 2)) {
+    rb_raise(rb_eArgError, "wrong # of arguments(%d for 2)",argc); SWIG_fail;
+  }
+  res1 = SWIG_ConvertPtr(self, &argp1,SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  if (!SWIG_IsOK(res1)) {
+    SWIG_exception_fail(SWIG_ArgError(res1), Ruby_Format_TypeError( "", "Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *","resize", 1, self )); 
+  }
+  arg1 = reinterpret_cast< Matrix< Complex< double >,Array2D_Dense< Complex< double > > > * >(argp1);
+  {
+    if(SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (argv[0], &temp2))){
+      raise_if_lt_zero_after_asval(temp2);
+      
+      arg2 = &temp2;
+    }
+    
+    else if(NIL_P(argv[0])){
+      arg2 = NULL;
+    }
+    
+    else{
+      SWIG_exception(SWIG_TypeError, "unsigned int is expected");
+    }
+  }
+  {
+    if(SWIG_IsOK(SWIG_AsVal_unsigned_SS_int (argv[1], &temp3))){
+      raise_if_lt_zero_after_asval(temp3);
+      
+      arg3 = &temp3;
+    }
+    
+    else if(NIL_P(argv[1])){
+      arg3 = NULL;
+    }
+    
+    else{
+      SWIG_exception(SWIG_TypeError, "unsigned int is expected");
+    }
+  }
+  {
+    try {
+      result = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *) &Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg__resize(arg1,(unsigned int const *)arg2,(unsigned int const *)arg3);
+    } catch (const native_exception &e) {
+      e.regenerate();
+      SWIG_fail;
+    } catch (const std::exception& e) {
+      SWIG_exception_fail(SWIG_RuntimeError, e.what());
+    }
+  }
+  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  return vresult;
+fail:
+  return Qnil;
+}
+
+
 SWIGINTERN void
 free_Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg_(void *self) {
     Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *arg1 = (Matrix< Complex< double >,Array2D_Dense< Complex< double > > > *)self;
@@ -20418,6 +20699,7 @@ static swig_type_info _swigt__p_f_r_q_const__Complex__double___p_Complex__double
 static swig_type_info _swigt__p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void = {"_p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void", "void (*)(double const &,double *,unsigned int const &,unsigned int const &)", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_self_t = {"_p_self_t", "self_t *", 0, 0, (void*)0, 0};
 static swig_type_info _swigt__p_swig__GC_VALUE = {"_p_swig__GC_VALUE", "swig::GC_VALUE *", 0, 0, (void*)0, 0};
+static swig_type_info _swigt__p_unsigned_int = {"_p_unsigned_int", "unsigned int *", 0, 0, (void*)0, 0};
 
 static swig_type_info *swig_type_initial[] = {
   &_swigt__p_ComplexT_double_t,
@@ -20435,6 +20717,7 @@ static swig_type_info *swig_type_initial[] = {
   &_swigt__p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void,
   &_swigt__p_self_t,
   &_swigt__p_swig__GC_VALUE,
+  &_swigt__p_unsigned_int,
 };
 
 static swig_cast_info _swigc__p_ComplexT_double_t[] = {  {&_swigt__p_ComplexT_double_t, 0, 0, 0},{0, 0, 0, 0}};
@@ -20452,6 +20735,7 @@ static swig_cast_info _swigc__p_f_r_q_const__Complex__double___p_Complex__double
 static swig_cast_info _swigc__p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void[] = {  {&_swigt__p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_self_t[] = {  {&_swigt__p_self_t, 0, 0, 0},{0, 0, 0, 0}};
 static swig_cast_info _swigc__p_swig__GC_VALUE[] = {  {&_swigt__p_swig__GC_VALUE, 0, 0, 0},{0, 0, 0, 0}};
+static swig_cast_info _swigc__p_unsigned_int[] = {  {&_swigt__p_unsigned_int, 0, 0, 0},{0, 0, 0, 0}};
 
 static swig_cast_info *swig_cast_initial[] = {
   _swigc__p_ComplexT_double_t,
@@ -20469,6 +20753,7 @@ static swig_cast_info *swig_cast_initial[] = {
   _swigc__p_f_r_q_const__double_p_double_r_q_const__unsigned_int_r_q_const__unsigned_int__void,
   _swigc__p_self_t,
   _swigc__p_swig__GC_VALUE,
+  _swigc__p_unsigned_int,
 };
 
 
@@ -20893,16 +21178,17 @@ SWIGEXPORT void Init_SylphideMath(void) {
   SWIG_TypeClientData(SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, (void *) &SwigClassMatrixD);
   rb_define_alloc_func(SwigClassMatrixD.klass, _wrap_MatrixD_allocate);
   rb_define_method(SwigClassMatrixD.klass, "initialize", VALUEFUNC(_wrap_new_MatrixD), -1);
-  rb_define_method(SwigClassMatrixD.klass, "swap_rows!", VALUEFUNC(_wrap_MatrixD_swap_rowsN___), -1);
-  rb_define_method(SwigClassMatrixD.klass, "swap_columns!", VALUEFUNC(_wrap_MatrixD_swap_columnsN___), -1);
   rb_define_method(SwigClassMatrixD.klass, "[]=", VALUEFUNC(_wrap_MatrixD___setitem__), -1);
   rb_define_singleton_method(SwigClassMatrixD.klass, "scalar", VALUEFUNC(_wrap_MatrixD_scalar), -1);
   rb_define_singleton_method(SwigClassMatrixD.klass, "I", VALUEFUNC(_wrap_MatrixD_I), -1);
+  rb_define_method(SwigClassMatrixD.klass, "swap_rows!", VALUEFUNC(_wrap_MatrixD_swap_rowsN___), -1);
+  rb_define_method(SwigClassMatrixD.klass, "swap_columns!", VALUEFUNC(_wrap_MatrixD_swap_columnsN___), -1);
   rb_define_method(SwigClassMatrixD.klass, "replace!", VALUEFUNC(_wrap_MatrixD_replaceN___), -1);
   rb_define_method(SwigClassMatrixD.klass, "map!", VALUEFUNC(_wrap_MatrixD_mapN___), -1);
   rb_define_alias(SwigClassMatrixD.klass, "collect!", "map!");
   rb_define_alias(SwigClassMatrixD.klass, "map_with_index!", "map!");
   rb_define_alias(SwigClassMatrixD.klass, "collect_with_index!", "map!");
+  rb_define_method(SwigClassMatrixD.klass, "resize!", VALUEFUNC(_wrap_MatrixD_resizeN___), -1);
   SwigClassMatrixD.mark = 0;
   SwigClassMatrixD.destroy = (void (*)(void *)) free_Matrix_Sl_double_Sc_Array2D_Dense_Sl_double_Sg__Sg_;
   SwigClassMatrixD.trackObjects = 0;
@@ -21026,16 +21312,17 @@ SWIGEXPORT void Init_SylphideMath(void) {
   SWIG_TypeClientData(SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, (void *) &SwigClassMatrixComplexD);
   rb_define_alloc_func(SwigClassMatrixComplexD.klass, _wrap_MatrixComplexD_allocate);
   rb_define_method(SwigClassMatrixComplexD.klass, "initialize", VALUEFUNC(_wrap_new_MatrixComplexD), -1);
-  rb_define_method(SwigClassMatrixComplexD.klass, "swap_rows!", VALUEFUNC(_wrap_MatrixComplexD_swap_rowsN___), -1);
-  rb_define_method(SwigClassMatrixComplexD.klass, "swap_columns!", VALUEFUNC(_wrap_MatrixComplexD_swap_columnsN___), -1);
   rb_define_method(SwigClassMatrixComplexD.klass, "[]=", VALUEFUNC(_wrap_MatrixComplexD___setitem__), -1);
   rb_define_singleton_method(SwigClassMatrixComplexD.klass, "scalar", VALUEFUNC(_wrap_MatrixComplexD_scalar), -1);
   rb_define_singleton_method(SwigClassMatrixComplexD.klass, "I", VALUEFUNC(_wrap_MatrixComplexD_I), -1);
+  rb_define_method(SwigClassMatrixComplexD.klass, "swap_rows!", VALUEFUNC(_wrap_MatrixComplexD_swap_rowsN___), -1);
+  rb_define_method(SwigClassMatrixComplexD.klass, "swap_columns!", VALUEFUNC(_wrap_MatrixComplexD_swap_columnsN___), -1);
   rb_define_method(SwigClassMatrixComplexD.klass, "replace!", VALUEFUNC(_wrap_MatrixComplexD_replaceN___), -1);
   rb_define_method(SwigClassMatrixComplexD.klass, "map!", VALUEFUNC(_wrap_MatrixComplexD_mapN___), -1);
   rb_define_alias(SwigClassMatrixComplexD.klass, "collect!", "map!");
   rb_define_alias(SwigClassMatrixComplexD.klass, "map_with_index!", "map!");
   rb_define_alias(SwigClassMatrixComplexD.klass, "collect_with_index!", "map!");
+  rb_define_method(SwigClassMatrixComplexD.klass, "resize!", VALUEFUNC(_wrap_MatrixComplexD_resizeN___), -1);
   SwigClassMatrixComplexD.mark = 0;
   SwigClassMatrixComplexD.destroy = (void (*)(void *)) free_Matrix_Sl_Complex_Sl_double_Sg__Sc_Array2D_Dense_Sl_Complex_Sl_double_Sg__Sg__Sg_;
   SwigClassMatrixComplexD.trackObjects = 0;