gps_pvt 0.2.0 → 0.3.0

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(", ")
@@ -3667,6 +3675,15 @@ SWIGINTERN Matrix< double,Array2D_Dense< double > >::self_t &Matrix_Sl_double_Sc
     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 +4097,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(", ")
@@ -4132,6 +4146,15 @@ SWIGINTERN Matrix< Complex< double >,Array2D_Dense< Complex< double > > >::self_
     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 +6794,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 +6922,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 +6979,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 +9076,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 +9132,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 +9185,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 +9833,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 +9961,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 +10018,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 +12115,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 +12171,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 +12224,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);
@@ -12308,6 +12359,8 @@ _wrap_MatrixD_swap_rowsN___(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< double,Array2D_Dense< double > >::self_t *) &(arg1)->swapRows((unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -12370,6 +12423,8 @@ _wrap_MatrixD_swap_columnsN___(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< double,Array2D_Dense< double > >::self_t *) &(arg1)->swapColumns((unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -12430,6 +12485,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 +12531,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 +12585,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 +12796,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 +12852,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 +12897,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);
@@ -13295,6 +13360,81 @@ fail:
 }
 
 
+/*
+  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 +14018,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 +14148,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 +14205,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 +16290,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 +16346,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 +16399,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 +17057,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 +17187,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 +17244,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 +19329,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 +19385,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 +19438,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);
@@ -19405,6 +19573,8 @@ _wrap_MatrixComplexD_swap_rowsN___(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< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &(arg1)->swapRows((unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -19467,6 +19637,8 @@ _wrap_MatrixComplexD_swap_columnsN___(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< Complex< double >,Array2D_Dense< Complex< double > > >::self_t *) &(arg1)->swapColumns((unsigned int const &)*arg2,(unsigned int const &)*arg3);
@@ -19527,6 +19699,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 +19745,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 +19799,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 +20008,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 +20062,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 +20107,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);
@@ -20388,6 +20570,81 @@ fail:
 }
 
 
+/*
+  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 +20675,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 +20693,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 +20711,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 +20729,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,
 };
 
 
@@ -20903,6 +21164,7 @@ SWIGEXPORT void Init_SylphideMath(void) {
   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;
@@ -21036,6 +21298,7 @@ SWIGEXPORT void Init_SylphideMath(void) {
   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;