gps_pvt 0.10.1 → 0.10.3

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: dc12b7e4590d91b7afacd099d36852acdde54d11e8418ff0c99386ac17e68ca4
-  data.tar.gz: f3286d74e35a744a2c451768f8785c259f36f5d6a94d04daaad71857f472429a
+  metadata.gz: 0d6b65fbf88b765e83cf8d2341d42e6919c8b79be5e98d2c493ca23ec5445e8e
+  data.tar.gz: 11bd177fd1cd8ba4ecf5134b64977f6e9a3b0e08cff3a290d1bca6897958f23b
 SHA512:
-  metadata.gz: 0611b095f8e2ccb80baa1dc9ba98291b6a5091dcbb14479b64058541260e5ccfd8ee1589aa40fc73393837b638f23633ac3dcfaa8fe23b396c1edfb960d47d75
-  data.tar.gz: 90731b787b1a62c718a49da481adc87bdedd26d57c207f29d69b35dffb0859bdaf185a6fe739060be68370cf66587feb26ffd18f01a2ed6037900a6c9527e35b
+  metadata.gz: e20cfb1dd15fe7ff0ff7320e6acc32d8966c1ad6d67f3203e27da757fb095e512085bdea4e7a9ac90282df8a4ea92235cb94a3b6210f6d1b8938c17f505f9ee3
+  data.tar.gz: 6897b4cdedcfd23edb23bdaf57f6d008f9c2562107dc7ccbf1788e7249637b5a1a7fc79b0af09d68a967aed5d86811947b4f29b2465b4db8835727adf69af4a0

data/README.md

@@ -43,7 +43,7 @@ The format of file is automatically determined with its extension, such as .ubx
 | --antex=file_or_URI | [Antenna Exchange Format](https://igs.org/wg/antenna#files) (supported gps_pvt version >= 0.6.0) |
 | --rinex_clk=file_or_URI | [RINEX clock](https://files.igs.org/pub/data/format/rinex_clock304.txt) file (supported gps_pvt version >= 0.7.0) |
 | <a name=opt_rtcm3>--rtcm3=file_or_URI</a> | [RTCM 10403.x](https://rtcm.myshopify.com/collections/differential-global-navigation-satellite-dgnss-standards). (supported gps_pvt version >= 0.9.0) The latest version uses message type Observation(GPS: 1001..1004; GLONASS: 1009..1012), Epehemris(GPS: 1019; GLOANSS: 1020; SBAS: 1043; QZSS: 1044), MSM(GPS: 1071..1077; GLONASS: 1081..1087; SBAS: 1101..1107; QZSS: 1111..1117) |
-| <a name=opt_supl>--supl=URI</a> | [SUPL, secure user plane location](https://www.openmobilealliance.org/release/SUPL/). (supported gps_pvt version >= 0.10.0) Both [LPP](https://portal.3gpp.org/desktopmodules/Specifications/SpecificationDetails.aspx?specificationId=3710)(default) and [RRLP](https://portal.3gpp.org/desktopmodules/Specifications/SpecificationDetails.aspx?specificationId=2688) are internally used, which can be manually selected by adding ```?protocol=lpp_or_rrlp``` URI query string. |
+| <a name=opt_supl>--supl=URI</a> | [SUPL, secure user plane location](https://www.openmobilealliance.org/release/SUPL/). (supported gps_pvt version >= 0.10.0) Whether [LPP(default)](https://portal.3gpp.org/desktopmodules/Specifications/SpecificationDetails.aspx?specificationId=3710) or [RRLP](https://portal.3gpp.org/desktopmodules/Specifications/SpecificationDetails.aspx?specificationId=2688) are internally used, which can be manually selected by adding ```?protocol=lpp_or_rrlp``` URI query string. |
 
 Since version 0.2.0, SBAS and QZSS are supported in addition to GPS. Since version 0.4.0, GLONASS is also available. QZSS ranging is activated in default, however, SBAS is just utilized for ionospheric correction. GLONASS is also turned off by default. If you want to activate SBAS or GLONASS ranging, "--with=(system or PRN)" options are used with gps_pvt executable like
 
@@ -172,7 +172,7 @@ Utility to get and dump GPS files. After installation of gps_pvt, to type
 
     $ gps_get file_or_URI(s) (options) > output_file
 
-saves data into output_file by using redirection. http(s), ftp, ntrip, and supl can be used as scheme of URI. Serial port is also supported. Note that compressed data is automatically decompressed before output. The following options are available.
+saves data into output_file by using redirection. http(s), ftp, ntrip, and supl can be used as scheme of URI. Serial port is also supported. Note that compressed data is automatically decompressed before output. It is also noted that supl received data are filtered to its LPP/RRLP content, and printed in JSON format. The following options are available.
 
 | key | value | comment | since |
 ----|----|----|----

data/Rakefile

@@ -11,6 +11,14 @@ Rake::ExtensionTask.new("gps_pvt") do |ext|
   ext.lib_dir = "lib/gps_pvt"
 end
 
+require "yard"
+
+YARD::Rake::YardocTask.new do |t|
+  #t.files = ['lib/**/*.rb', 'ext/**/*_wrap.cxx']
+  t.options = ['--embed-mixins']
+  #t.stats_options = ['--list-undoc']
+end
+
 namespace :git do
   task :version do
     @git_version ||= proc{
@@ -83,25 +91,28 @@ end
 
 desc "Generate SWIG wrapper codes"
 task :swig do
-  swig_dir = File::join(File::dirname(__FILE__), 'ext', 'ninja-scan-light', 'tool', 'swig')
   out_base_dir = File::join(File::dirname(__FILE__), 'ext', 'gps_pvt')
-  Dir::chdir(swig_dir){
-    Dir::glob("*.i"){|src|
-      mod_name = File::basename(src, '.*')
-      out_dir = File::join(out_base_dir, mod_name)
-      sh "mkdir -p #{out_dir}"
-      wrapper = File::join(out_dir, "#{mod_name}_wrap.cxx")
-      sh [:make, :clean, wrapper,
-          "BUILD_DIR=#{out_dir}",
-          "SWIGFLAGS='-c++ -ruby -prefix \"GPS_PVT::\"#{" -D__MINGW__" if ENV["MSYSTEM"]}'"].join(' ')
-      open(wrapper, 'r+'){|io|
-        lines = io.read.lines.collect{|line|
-          line.sub(/rb_require\(\"([^\"]+)\"\)/){ # from camel to underscore downcase style
-            "rb_require(\"#{$1.sub('GPS_PVT', 'gps_pvt')}\")"
+  [
+    File::join(File::dirname(__FILE__), 'ext', 'ninja-scan-light', 'tool', 'swig'),
+  ].each{|swig_dir|
+    Dir::chdir(swig_dir){
+      Dir::glob("*.i"){|src|
+        mod_name = File::basename(src, '.*')
+        out_dir = File::join(out_base_dir, mod_name)
+        sh "mkdir -p #{out_dir}"
+        wrapper = File::join(out_dir, "#{mod_name}_wrap.cxx")
+        sh [:make, :clean, wrapper,
+            "BUILD_DIR=#{out_dir}",
+            "SWIGFLAGS='-c++ -ruby -prefix \"GPS_PVT::\"#{" -D__MINGW__" if ENV["MSYSTEM"]}'"].join(' ')
+        open(wrapper, 'r+'){|io|
+          lines = io.read.lines.collect{|line|
+            line.sub(/rb_require\(\"([^\"]+)\"\)/){ # from camel to underscore downcase style
+              "rb_require(\"#{$1.sub('GPS_PVT', 'gps_pvt')}\")"
+            }
           }
+          io.rewind
+          io.write(lines.join)
         }
-        io.rewind
-        io.write(lines.join)
       }
     }
   }

data/ext/gps_pvt/SylphideMath/SylphideMath_wrap.cxx

@@ -8495,7 +8495,7 @@ _wrap_Matrix_FrozenD_each__SWIG_0(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -8531,7 +8531,7 @@ _wrap_Matrix_FrozenD_each__SWIG_1(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -11940,7 +11940,7 @@ _wrap_Matrix_FrozenD_f_each__SWIG_0(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -11976,7 +11976,7 @@ _wrap_Matrix_FrozenD_f_each__SWIG_1(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_double_Array2D_DenseT_double_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -13956,7 +13956,7 @@ _wrap_MatrixD_resizeN___(int argc, VALUE *argv, VALUE self) {
     }
   }
   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);
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_double_Array2D_DenseT_double_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -16448,7 +16448,7 @@ _wrap_Matrix_FrozenComplexD_each__SWIG_0(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -16484,7 +16484,7 @@ _wrap_Matrix_FrozenComplexD_each__SWIG_1(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -19891,7 +19891,7 @@ _wrap_Matrix_FrozenComplexD_f_each__SWIG_0(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -19927,7 +19927,7 @@ _wrap_Matrix_FrozenComplexD_f_each__SWIG_1(int argc, VALUE *argv, VALUE self) {
     (&_e)->regenerate();
     SWIG_fail;
   }
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_Matrix_FrozenT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewFilterT_MatrixViewBaseT_t_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;
@@ -21903,7 +21903,7 @@ _wrap_MatrixComplexD_resizeN___(int argc, VALUE *argv, VALUE self) {
     }
   }
   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);
-  vresult = SWIG_NewPointerObj(SWIG_as_voidptr(result), SWIGTYPE_p_MatrixT_ComplexT_double_t_Array2D_DenseT_ComplexT_double_t_t_MatrixViewBaseT_t_t, 0 |  0 );
+  vresult = self;
   return vresult;
 fail:
   return Qnil;

data/ext/ninja-scan-light/tool/navigation/GPS.h

@@ -95,6 +95,11 @@ class GPS_Signal {
           PRN::content <<= 1;
           PRN::content[0] = tmp;
         }
+        void previous(){
+          bool tmp(PRN::content[3] ^ PRN::content[0]);
+          PRN::content >>= 1;
+          PRN::content[9] = tmp;
+        }
     };
     
     class G2 : public PRN {
@@ -114,6 +119,16 @@ class GPS_Signal {
           PRN::content <<= 1;
           PRN::content[0] = tmp;
         }
+        void previous(){
+          bool tmp(PRN::content[2]
+                     ^ PRN::content[3]
+                     ^ PRN::content[6]
+                     ^ PRN::content[8]
+                     ^ PRN::content[9]
+                     ^ PRN::content[0]);
+          PRN::content >>= 1;
+          PRN::content[9] = tmp;
+        }
         static G2 get_G2(const int &prn){
           switch(prn){
             case  1: return G2(1, 5);
@@ -160,9 +175,9 @@ class GPS_Signal {
     class CA_Code {
       public:
         typedef FloatT float_t;
-        static const float_t FREQENCY;
+        static const float_t FREQUENCY;
         static const float_t length_1chip() {
-          static const float_t res(1. / FREQENCY);
+          static const float_t res(1. / FREQUENCY);
           return res;
         }
       protected:
@@ -177,11 +192,15 @@ class GPS_Signal {
           g1.next();
           g2.next();
         }
+        void previous(){
+          g1.previous();
+          g2.previous();
+        }
     };
 };
 
 template <class FloatT>
-const typename GPS_Signal<FloatT>::float_t GPS_Signal<FloatT>::CA_Code::FREQENCY = 1.023E6;
+const typename GPS_Signal<FloatT>::float_t GPS_Signal<FloatT>::CA_Code::FREQUENCY = 1.023E6;
 
 template <class FloatT = double>
 struct GPS_Time {

data/ext/ninja-scan-light/tool/swig/Coordinate.i

@@ -73,6 +73,8 @@ struct native_exception : public std::exception {
       values[i] = (*self)[i];
     }
   }
+  %typemap(in) FloatT values[3];
+  %typemap(argout) FloatT values[3];
 #endif
   %typemap(in,numinputs=0) const void *check_block {
 #ifdef SWIGRUBY
@@ -98,6 +100,7 @@ struct native_exception : public std::exception {
 #endif
     }
   }
+  %clear const void *check_block;
 };
 #ifdef SWIGRUBY
 %mixin System_3D "Enumerable";
@@ -170,6 +173,9 @@ struct native_exception : public std::exception {
 %template(XYZ) System_XYZ<type, WGS84>;
 %template(LLH) System_LLH<type, WGS84>;
 %template(ENU) System_ENU<type, WGS84>;
+%typemap(out) FloatT &;
+%typemap(in) (type (&res)[3][3]);
+%typemap(argout) (type (&res)[3][3]);
 %enddef
 
 CONCRETIZE(double);

data/ext/ninja-scan-light/tool/swig/SylphideMath.i

@@ -960,6 +960,7 @@ struct MatrixUtil {
           std::string("Unknown enumerate direction: ").append(inspect_str($1)).c_str());
     }
   }
+  %typemap(out) const Matrix_Frozen<T, Array2D_Type, ViewType> & "$result = self;"
   %catches(native_exception) each;
   const Matrix_Frozen<T, Array2D_Type, ViewType> &each(
       void (*each_func)(
@@ -970,6 +971,7 @@ struct MatrixUtil {
     return *$self;
   }
   %alias each "each_with_index";
+  %typemap(out) const Matrix_Frozen<T, Array2D_Type, ViewType> &;
   
   %catches(native_exception, std::invalid_argument) map;
   Matrix<T, Array2D_Dense<T> > map(
@@ -1005,6 +1007,7 @@ struct MatrixUtil {
         (RTEST(idx_selector) ? 1 : 0), &idx_selector,
         (rb_block_call_func_t)rb_equal, value);
   }
+  %typemap(check) VALUE idx_selector;
   %alias index "find_index";
   
   SWIG_Object to_a() const {
@@ -1020,6 +1023,20 @@ struct MatrixUtil {
     return res;
   }
 #endif
+  %typemap(in)
+      Matrix<T, Array2D_Dense<T> > &output_L,
+      Matrix<T, Array2D_Dense<T> > &output_U,
+      Matrix<T, Array2D_Dense<T> > &output_P,
+      Matrix<T, Array2D_Dense<T> > &output_D,
+      Matrix<T, Array2D_Dense<T> > &output_Q,
+      Matrix<T, Array2D_Dense<T> > &output_R;
+  %typemap(argout)
+      Matrix<T, Array2D_Dense<T> > &output_L,
+      Matrix<T, Array2D_Dense<T> > &output_U,
+      Matrix<T, Array2D_Dense<T> > &output_P,
+      Matrix<T, Array2D_Dense<T> > &output_D,
+      Matrix<T, Array2D_Dense<T> > &output_Q,
+      Matrix<T, Array2D_Dense<T> > &output_R;
 };
 
 #if defined(SWIGRUBY)
@@ -1171,6 +1188,15 @@ MAKE_TO_S(Matrix_Frozen)
   %rename("map!") map_bang;
   %alias map_bang "collect!,map_with_index!,collect_with_index!";
 #endif
+
+  %typemap(in)
+      void (*each_func)(const T &src, T *dst, const unsigned int &i, const unsigned int &j);
+  %typemap(typecheck) const typename MatrixUtil::each_which_t each_which;
+  %typemap(in) const typename MatrixUtil::each_which_t each_which;
+  %typemap(typecheck) const void *replacer;
+  %typemap(in) const void *replacer;
+  %typemap(in) self_t *self_p;
+  %typemap(argout) self_t *self_p;
 };
 
 %define INSTANTIATE_MATRIX_TRANSPOSE(type, storage, view_from, view_to)
@@ -1251,6 +1277,12 @@ MAKE_TO_S(Matrix_Frozen)
     }
     output_D = D;
   }
+  %typemap(in)
+      Matrix<ctype, Array2D_Dense<ctype > > &output_D,
+      Matrix<ctype, Array2D_Dense<ctype > > &output_V;
+  %typemap(argout)
+      Matrix<ctype, Array2D_Dense<ctype > > &output_D,
+      Matrix<ctype, Array2D_Dense<ctype > > &output_V;
 };
 %enddef
 %define INSTANTIATE_MATRIX_EIGEN(type, ctype)
@@ -1366,6 +1398,7 @@ INSTANTIATE_MATRIX_PARTIAL(type, Array2D_Dense<type >, MatView_pt, MatView_pt);
 #endif
     else{SWIG_exception(SWIG_TypeError, "$*1_ltype is expected");}
   }
+  %typemap(out) Matrix<type, Array2D_Dense<type > > & "$result = self;"
   Matrix<type, Array2D_Dense<type > > &resize(
       const unsigned int *r_p, const unsigned int *c_p){
     unsigned int r(r_p ? *r_p : $self->rows()), c(c_p ? *c_p : $self->columns());
@@ -1376,6 +1409,10 @@ INSTANTIATE_MATRIX_PARTIAL(type, Array2D_Dense<type >, MatView_pt, MatView_pt);
     mat_new.partial(r_min, c_min).replace($self->partial(r_min, c_min), false);
     return (*($self) = mat_new);
   }
+  %clear Matrix<type, Array2D_Dense<type > > &;
+
+  //%typemap(in) unsigned int *r_p, unsigned int *c_p; // NG; remove custom typemap before Matrix::resize! generation
+  %clear unsigned int *r_p, unsigned int *c_p; // OK, work around version
 };
 
 %template(Matrix ## suffix) Matrix<type, Array2D_Dense<type > >;

data/ext/ninja-scan-light/tool/swig/spec/SylphideMath_spec.rb

@@ -286,7 +286,7 @@ shared_examples 'Matrix' do
       ].each{|arg|
         mat_orig = mat[0].copy
         mat_replaced = arg.kind_of?(Proc) ? mat_orig.send(:replace!, &arg) : mat_orig.send(:replace!, arg)
-        expect(mat_replaced).to equal(mat_orig)
+        expect(mat_replaced).to be(mat_orig)
         expect(mat_replaced).not_to equal(mat[1])
         expect(mat_replaced.to_a).to eq(mat[1].to_a)
       }
@@ -306,7 +306,7 @@ shared_examples 'Matrix' do
             idxs[a], idxs[b] = [b, a]
             mat_builtin = Matrix::columns(mat_builtin.column_vectors.values_at(*idxs))
           end
-          expect(mat_mod).to equal(mat[0])
+          expect(mat_mod).to be(mat[0])
           expect(mat[0].to_a).to eq(mat_builtin.to_a)
         }
       }
@@ -391,11 +391,12 @@ shared_examples 'Matrix' do
           candidates = (func.to_s =~ /with_index$/) \
               ? indices.collect{|i, j| [compare_with[i][j], i, j]} \
               : indices.collect{|i, j| [compare_with[i][j]]}
-          mat.send(*[func, k].compact){|*v|
+          mat2 = mat.send(*[func, k].compact){|*v|
             i = candidates.find_index(v)
             expect(i).not_to be(nil)
             candidates.delete_at(i)
           }
+          expect(mat2).to be(mat)
           expect(candidates.empty?).to be(true)
         }
       }
@@ -449,7 +450,7 @@ shared_examples 'Matrix' do
             v[0] * 2
           }
           expect(candidates.empty?).to be(true)
-          expect(mat2.to_a).to eq(mat.to_a)
+          expect(mat2).to be(mat)
           expect(mat2.to_a).to eq(compare_with.collect.with_index{|values, i|
                 values.collect.with_index{|v, j|
                   indices.include?([i, j]) ? (v * 2) : v
@@ -541,6 +542,7 @@ shared_examples 'Matrix' do
       }
       mat_orig = mat[0].to_a
       r, c = [:rows, :columns].collect{|f| mat[0].send(f)}
+      expect(mat[0].resize!(r, c)).to be(mat[0])
       expect(mat[0].resize!(r, c).to_a).to eq(mat_orig)
       expect(mat[0].resize!(r, nil).to_a).to eq(mat_orig)
       expect(mat[0].resize!(nil, c).to_a).to eq(mat_orig)

data/gps_pvt.gemspec

@@ -64,6 +64,7 @@ Gem::Specification.new do |spec|
   spec.add_development_dependency "rspec", "~> 3.0"
   spec.add_development_dependency "matrix" if GPS_PVT::version_compare(RUBY_VERSION, "3.1") >= 0
   spec.add_development_dependency "racc"
+  spec.add_development_dependency "yard"
   spec.add_development_dependency "github_changelog_generator" unless ((Gem::Platform.local.os =~ /mingw/) && (GPS_PVT::version_compare(RUBY_VERSION, "3.1") >= 0))
 
   # For more information and examples about making a new gem, checkout our

data/lib/gps_pvt/asn1/asn1.rb

@@ -1,3 +1,5 @@
+require 'set'
+
 module GPS_PVT
 end
 
@@ -261,10 +263,10 @@ resolve_tree = proc{|root|
     res
   }
   
-  # This proc is useless because assumed automatic tagging 
-  # does not require to re-order based on manual tags or class numbers
+  # This proc is used to re-order for non-automatic tagging transformation
   # @see https://stackoverflow.com/a/31450137
   get_universal_class_number = proc{|type|
+    # TODO where is CHOICE class number definition?
     next get_universal_class_number.call(type[1][:root][0][:type]) if type[0] == :CHOICE
     { # @see Table.1 unless any comment
       :BOOLEAN => 1,
@@ -274,6 +276,9 @@ resolve_tree = proc{|root|
       :NULL => 5,
       :ENUMERATED => 10,
       :SEQUENCE => 16,
+      :SEQUENCE_OF => 16,
+      :SET => 17,
+      :SET_OF => 17,
       :NumericString => 18, # @see Table.6
       :PrintableString => 19, # @see Table.6
       :IA5String => 22, # @see Table.6
@@ -281,6 +286,23 @@ resolve_tree = proc{|root|
       :UTCTime => 23, # @see 43.3
     }[type[0]]
   }
+
+  reorder_children = proc{|type_opts|
+    # reordering when automatic tagging transformation is not selected
+    tags = Set[]
+    [type_opts[:root], type_opts[:extension]].compact.each{|items|
+      items.collect!{|v|
+        tag = v[:type][1][:tag]
+        tag_class_num = case tag
+        when Array; [{"APPLICATION" => 1, "PRIVATE" => 3}[tag[0]], tag[1]]
+        when Integer; [2, tag] # (Context-specific) 
+        else; [0, get_universal_class_number.call(v[:type])] # UNIVERSAL
+        end
+        raise unless tags.add?(tag_class_num) # Identical tag is already existed!
+        [v, tag_class_num]
+      }.sort!{|a, b| a[1] <=> b[1]}.collect!{|v| v[0]}
+    }
+  }
   
   prepare_coding = proc{|tree|
     next tree.each{|k, v|
@@ -313,7 +335,8 @@ resolve_tree = proc{|root|
       }
     when :BIT_STRING, :OCTET_STRING
       opts[:size_range] = find_range.call(opts, :size)
-    when :SEQUENCE
+    when :SEQUENCE, :SET
+      reorder_children.call(opts) if (opts[:automatic_tagging] == false) # only for SET
       (opts[:root] + (opts[:extension] || [])).each.with_index{|v, i|
         v[:name] = v[:name] ? v[:name].to_sym : i
         v[:type][1][:typename] ||= v[:name] if v[:name].kind_of?(Symbol) && v[:type] # for debugger
@@ -325,11 +348,11 @@ resolve_tree = proc{|root|
           v2[:name]
         } if v[:group]
       }
-    when :SEQUENCE_OF
+    when :SEQUENCE_OF, :SET_OF
       opts[:size_range] = find_range.call(opts, :size)
       prepare_coding.call(opts)
     when :CHOICE
-      # Skip reordering based on automatic tagging assumption
+      reorder_children.call(opts) if (opts[:automatic_tagging] == false)
       opts[:extension] = opts[:extension].collect{|v|
         v[:group] || [v] # 22. Note says "Version brackets have no effect"
       }.flatten(1) if opts[:extension]
@@ -372,7 +395,7 @@ generate_skeleton = proc{|tree, data|
     when :BIT_STRING, :OCTET_STRING
       data || ({:BIT_STRING => [0], :OCTET_STRING => [0xFF]}[type] \
           * (opts[:size_range][:root].first rescue 0))
-    when :SEQUENCE
+    when :SEQUENCE, :SET
       data ||= {}
       Hash[*((opts[:root] + (opts[:extension] || [])).collect{|v|
         if v[:group] then
@@ -386,7 +409,7 @@ generate_skeleton = proc{|tree, data|
           end
         end
       }.compact.flatten(2))]
-    when :SEQUENCE_OF
+    when :SEQUENCE_OF, :SET_OF
       next data.collect{|v| generate_skeleton.call(opts, v)} if data
       v = Marshal::dump(generate_skeleton.call(opts))
       (opts[:size_range][:root].first rescue 0).times.collect{
@@ -488,7 +511,7 @@ encode = proc{|tree, data|
         end
       end
       res += data.collect{|v| "%0#{bits}b"%[v]}.join
-    when :SEQUENCE
+    when :SEQUENCE, :SET
       opt_def_flags, root_encoded = opts[:root].collect{|v| # 18.2
         has_elm = data.include?(v[:name])
         elm = data[v[:name]]
@@ -532,7 +555,7 @@ encode = proc{|tree, data|
       end
 
       "#{ext_bit}#{opt_def_flags.join}#{root_encoded.join}#{ext_encoded}"
-    when :SEQUENCE_OF
+    when :SEQUENCE_OF, :SET_OF
       ext_bit, len_enc = case (cat = opts[:size_range].belong_to(data.size))
       when :additional
         # 19.4 -> 10.9.4.2(semi_constrained_whole_number)
@@ -691,7 +714,7 @@ decode = proc{|tree, str|
         end
       end
       str.slice!(0, bits * len).scan(/.{#{bits}}/).collect{|chunk| chunk.to_i(2)}
-    when :SEQUENCE
+    when :SEQUENCE, :SET
       has_extension = (opts[:extension] && (str.slice!(0) == '1'))
       data = Hash[*(
         opts[:root].collect{|v| [v[:name], v[:default]] if v[:default]}.compact.flatten(1)
@@ -709,7 +732,7 @@ decode = proc{|tree, str|
             v[:group] ? decoded.to_a : [[v[:name], decoded]]
           }.compact.flatten(2))]) if has_extension
       data
-    when :SEQUENCE_OF
+    when :SEQUENCE_OF, :SET_OF
       len_dec = if opts[:size_range][:additional] && (str.slice!(0) == '1') then
         # 19.4 -> 10.9.4.2(semi_constrained_whole_number)
         :length_otherwise