when_exe 0.3.4 → 0.3.5

data/lib/when_exe/region/chinese.rb

@@ -12,7 +12,7 @@ module When
     ChineseTerms = [self, [
       "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
       "locale:[=ja:, en=en:, alias]",
-      "names:[ChineseTerms]",
+      "names:[ChineseTerms=]",
       "[中国太陽暦(節月)=ja:%%<二十四節気>#%.<暦月と節月>, *ChineseSolar=en:Solar_term]",
       "[中国太陰太陽暦=ja:%%<中国暦>, *ChineseLuniSolar=en:Chinese_calendar]",
       "[彝暦=ja:%%<イ族>, *Yi=en:Yi_people]",
@@ -69,7 +69,7 @@ module When
     YiTerms = [self, [
       "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
       "locale:[=ja:, en=en:, alias]",
-      "names:[YiTerms]",
+      "names:[YiTerms=]",
 
       [self,
         "names:[月=ja:%%<月_(暦)>, *Month]",
@@ -117,129 +117,559 @@ module When
     #
     class ChineseTrueLunation < MeanLunation
 
-      # 太陽の位置補正表
-      # @return [Array<Array< 入気定日加減数, 朓朒数, 損益率, 損益率増分 >>]
-      attr_reader :s
+      #
+      # 唐代のアルゴリズム
+      #
+      # @private
+      module TangMethods
 
-      # 月の位置補正表
-      # @return [Array<Array< 区間の時間／分, 損益率 >>]
-      attr_reader :m
+        # 太陽の位置補正表
+        # @return [Array<Array< 入気定日加減数, 朓朒数, 損益率, 損益率増分 >>]
+        #attr_reader :s
 
-      # 近点月
-      # @return [Numeric]
-      attr_reader :anomalistic_month_length
+        # 月の位置補正表
+        # @return [Array<Array< 区間の時間／分, 損益率 >>]
+        #attr_reader :m
 
-      # 元期の近点離隔
-      # @return [Numeric]
-      attr_reader :anomalistic_month_shift
+        #
+        # 立成の初期化
+        #
+        def _initialize_rissei
 
-      # 当該日付のthiti の変化範囲(唐代の定朔の暦法用 cn_to_time(1L) を使用する)
-      #
-      # @param [When::TM::TemporalPosition] date 日付
-      #
-      # @return [Range] 当該日付のthiti の変化範囲(朔を含む場合 nil)
-      #
-      def thiti_range(date)
-        date = date.floor
-        p0, p1 = [date, date.succ].map {|d|
-          t  = d.to_f
-          c  = (30.0 * ((t - CYCLE_0M) * @cycle_number_1m + @cycle_number_0m)).floor
-          t0 = t1 = nil
-          loop do
-            t0 = cn_to_time( c   / 30.0)
-            t1 = cn_to_time((c+1)/ 30.0)
-            if t0 > t
-              c -= 1
-            elsif t1 <= t
-              c += 1
+          # 近点月
+          @anomalistic_month_length =  @anomalistic_month_length.to_r
+
+          # 通法
+          @denominator = @year_length.denominator if @year_length.kind_of?(Rational)
+          @denominator = [@denominator||0, @lunation_length.denominator].max if @lunation_length.kind_of?(Rational)
+
+          # 計算式 C) のための繰り返し回数
+          @repeat_count ||= @method.upcase == 'C' ? 2 : 1
+          @repeat_count   = @repeat_count.to_i
+
+          # 太陽の盈縮表の生成
+          #   [先後數, 朓朒数, (立成b, 立成c)]
+          #      ↓
+          #   [区間の始めの冬至からの経過日数, 朓朒数, 立成b, 立成c]
+          @s = @s.map {|item| item.dup}
+          if @rissei # 立成の計算・比較
+            puts "\n  ["+ @denominator.to_s + "]" if @rissei == @rissei.upcase
+            (0...@s.size).each do |i|
+              bc = send('_rissei_' + @rissei.downcase, i)
+              @s[i] += bc if @s[i].size <= 2
+              if @rissei == @rissei.upcase
+                sq = (bc[0]-@s[i][2])**2 + (bc[1]-@s[i][3])**2
+                puts '%9.1f,%5d,%9.4f,%8.4f,%6.2f' %
+                  (@s[i] + [sq == 0 ? -Float::INFINITY : Math.log(sq)/Math.log(100)])
+              end
+            end
+          end
+          (0...@s.size).each do |i|
+            @s[i-1][1,0] = @year_length / @s.size + (@s[i][0]-@s[i-1][0]) / @denominator
+          end
+          (0...@s.size).each do |i|
+            @s[i].shift
+          end
+
+          # 月の遅速表の生成
+          #   [変日差, 損益率]
+          #      ↓
+          #   [変日(区間の終わりの遠/近地点からの経過日数), 区間の終わりの朓朒積, 区間の変日差/日, 区間の損益率]
+          @m = @m.map {|item| item.dup}
+          sum_t = sum_v = 0
+          (0...@m.size).each do |i|
+            sum_t += @m[i][0]
+            sum_v += @m[i][1]
+            @m[i]  = [sum_t / @denominator, sum_v, @m[i][0] / @denominator, @m[i][1]]
+          end
+        end
+
+        # 戊寅暦 立成 １次補間
+        def _rissei_a(k)
+          tv = [0,1].map {|i| _tv_s(k+i)}
+          [(tv[1][1]-tv[0][1]) / (tv[1][0]-tv[0][0]), 0.0]
+        end
+
+        # 儀鳳暦 立成
+        def _rissei_b(k)
+          v =
+            if k % 3 < 2
+              t   = [0,1,2].map {|n| @s[(k+n) % @s.size][1]}
+              t01 = t[1] - t[0]
+              t12 = t[2] - t[1]
+              (t01 + t12) / 2.0 + t01 - t12
             else
-              break
+              t   = [-1,0,1].map {|n| @s[(k+n) % @s.size][1]}
+              t01 = t[1] - t[0]
+              t12 = t[2] - t[1]
+              t02 = t[2] - t[0]
+              t02 / 2.0
             end
+          r = v < 0 ? 12.0/(17*11) : 12.0/(16*11)
+          [v * r, (t12-t01) * r * r]
+        end
+
+        # 宣明暦 立成 ２次補間
+        def _rissei_c(k)
+          i = [k, k+1, k+2, k % 6 == 5 ? -1 : 0]
+          t0, v0 = _tv_s(i[0]+i[3])
+          t1, v1 = _tv_s(i[1]+i[3])
+          t2, v2 = _tv_s(i[2]+i[3])
+          t01    = t1 - t0
+          t02    = t2 - t0
+          t12    = t2 - t1
+          v01    = v1 - v0
+          v02    = v2 - v0
+          c      = (v02*t01 - v01*t02) / (t01*t02*t12) * 2
+          b      = (v01 - t01*(t01-1)/2*c) / t01
+          [b-t01*c*i[3], c]
+        end
+
+        # 日行盈縮
+        #
+        # @param [Integer] i 区間の番号
+        #
+        # @return [区間の始めの冬至からの経過日数, 朓朒数]
+        #
+        def _tv_s(i)
+          [@year_length * i / @s.size + @s[i % @s.size][0] / @denominator, @s[i % @s.size][1]]
+        end
+
+        # 月行遅速
+        #
+        # @param [Numeric] t 直前の遠/近地点からの経過日数
+        #
+        # @return [区間の始めからの経過日数, 区間の始めの朓朒積, 区間の変日差/日, 区間の損益率, 次の区間の損益率]
+        #
+        def _tv_m(t)
+          (0...@m.size).each do |i|
+            next if t > @m[i][0]
+            return [t - (@m[i][0] - @m[i][2]), @m[i][1] - @m[i][3], @m[i][2], @m[i][3], @m[(i+1) % @m.size][3]]
           end
-          (c + (t-t0) / (t1-t0)) % 30.0
-        }
-        p0 >= p1 ? nil : p0...p1
+        end
       end
 
-      private
-
-      # 周期番号 -> 日時
       #
-      # @param [Numeric] cn 周期番号
+      # 日本暦日原典 計算 A)
       #
-      # @return [Numeric] ユリウス日
+      # @private
+      module MethodA
+
+        include TangMethods
+
+        private
+
+        # 周期番号 -> 日時
+        #
+        # @param [Numeric] cn 周期番号
+        #
+        # @return [Numeric] ユリウス日
+        #
+        def cn_to_time_(cn, time0=nil)
+          time = super
+          t = time - @day_epoch
+          time + (delta_s(t.divmod(@year_length)[1]) + delta_m(t.divmod(@anomalistic_month_length)[1])).to_r / @denominator
+        end
+
+        #
+        # 朔の日時の太陽運動の不斉による補正
+        #
+        # @param [Numeric] t 直前の冬至からの日数
+        #
+        # @return [Numeric] 補正量 / @denominator
+        #
+        def delta_s(t)
+
+          t0, a, b, c = nil
+          @s.each do |v|
+            t0, a, b, c = v
+            break if t <= t0
+            t -= t0
+          end
+
+          # ２４気からの日数 (大余, 小余)
+          t, dt = t.divmod(1)
+          dt    = (dt * @denominator).to_i
+
+          # 損益率
+          b0 = (b + c * t).to_i
+
+          # 朓朒数
+          a0 = (a + b * t + c * t * (t-1.0)/2).to_i
+
+          # 補正値
+          a0 + (b0.to_f * dt / @denominator + 0.5).floor
+        end
+
+        #
+        # 朔の日時の月運動の不斉による補正
+        #
+        # @param [Numeric] t 直前の遠/近地点からの日数
+        #
+        # @return [Numeric] 補正量 / @denominator
+        #
+        def delta_m(t)
+
+          dt, a0, t0, b0, b1 = _tv_m(t)
+
+          # 補正値 (A式)
+          a0 + (b0.to_f * dt / t0 + 0.5).floor
+        end
+      end
+
       #
-      def cn_to_time_(cn, time0=nil)
-        time = super
-        return time unless @formula == '1L'
-        t = time - @day_epoch
-        time + (delta_s(t.divmod(@year_length)[1]) + delta_m(t.divmod(@anomalistic_month_length)[1])).to_r / @denominator
+      # 日本暦日原典 計算 C)
+      #
+      # @private
+      module MethodC
+
+        include TangMethods
+
+        private
+
+        # 周期番号 -> 日時
+        #
+        # @param [Numeric] cn 周期番号
+        #
+        # @return [Numeric] ユリウス日
+        #
+        def cn_to_time_(cn, time0=nil)
+          time = super
+          t    = time - @day_epoch
+          t   += delta_s(t.divmod(@year_length)[1]) / @denominator
+          d    = 0
+          @repeat_count.times do
+            d  = delta_m((t+d/@denominator).divmod(@anomalistic_month_length)[1])
+          end
+          t + @day_epoch + d / @denominator
+        end
+
+        #
+        # 朔の日時の太陽運動の不斉による補正
+        #
+        # @param [Numeric] t 直前の冬至からの日数
+        #
+        # @return [Numeric] 補正量 / @denominator
+        #
+        def delta_s(t)
+
+          t0, a, b, c = nil
+          @s.each do |v|
+            t0, a, b, c = v
+            break if t <= t0
+            t -= t0
+          end
+
+          # ２４気からの日数 (大余, 小余)
+          t, dt = t.divmod(1)
+
+          # 損益率
+          b0 = b + c * t
+
+          # 朓朒数
+          a0 = a + b * t + c * t * (t-1)/2
+
+          # 補正値
+          a0 + b0 * dt
+        end
+
+        #
+        # 朔の日時の月運動の不斉による補正
+        #
+        # @param [Numeric] t 直前の遠/近地点からの日数
+        #
+        # @return [Numeric] 補正量 / @denominator
+        #
+        def delta_m(t)
+
+          #Sm, 遅速積, 区間長, Δ1, Δ2
+          dt, a0, t0, b0, b1 = _tv_m(t)
+
+          # 補正値 (B,C式)
+          a0 + dt / t0 * ((b0+b1)/2.0 + 0.5*(1+t0-dt)*(b0-b1))
+        end
       end
 
       #
-      # 朔の日時の太陽運動の不斉による補正
+      # 日本暦日原典 計算 B)
       #
-      # @param [Numeric] t 直前の冬至からの日数
+      # @private
+      MethodB = MethodC
+
       #
-      # @return [Numeric] 補正量 / @denominator
+      # 元明代のアルゴリズム
       #
-      def delta_s(t)
+      # @private
+      module JujiMethods
 
-        t0, a, b, c = nil
-        @s.each do |v|
-          t0, a, b, c = v
-          break if t <= t0
-          t -= t0
+        # 授時暦・大統暦
+        module C
+          def _shift_l(year); @year_span * (year / @year_span) end
+          alias :_shift_s :_shift_l
         end
 
-        # ２４気からの日数 (大余, 小余)
-        t, dt = t.divmod(1)
-        dt    = (dt * @denominator).to_i
+        # 貞享暦・宝暦暦
+        module Y
+          def _shift_l(year); year   end
+          def _shift_s(year); year-1 end
+        end
 
-        # 損益率
-        b0 = (b + c * t).to_i
+        # 寛政暦?
+        module D
+          def _shift_l(year); 5 + year.div(10) * 10 end
+          alias :_shift_s :_shift_l
+        end
 
-        # 朓朒数
-        a0 = (a + b * t + c * t * (t-1.0)/2).to_i
+        # 暦元天正冬至から当該年の天正冬至までの日数
+        def _winter_solstice_(year)
+          year * (@year_length - @year_delta  * _shift_s(year))
+        end
 
-        # 補正値
-        a0 + (b0.to_f * dt / @denominator + 0.5).floor
+        # 暦元天正冬至から当該年の近日点通過までの日数
+        def _perihelion_(year)
+          _winter_solstice(year) + @anomalistic_year_shift
+        end
+
+        # 歳周(当該年の日数)
+        def _year_length_(year)
+          @year_length - 2 * @year_delta * _shift_l(year)
+        end
       end
 
       #
-      # 朔の日時の月運動の不斉による補正
+      # 授時暦の平均太陽の計算
       #
-      # @param [Numeric] t 直前の遠/近地点からの日数
+      # @private
+      module MethodS
+
+        include JujiMethods
+
+        # 日時 -> 周期番号
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric] 周期番号
+        #
+        def time_to_cn(t, cn0=nil)
+          cn0 ||= (t.to_f - @day_epoch) / @year_length + @year_epoch + @longitude_shift
+          root(cn0 * 12, t.to_f) {|cn| cn_to_time(cn) }
+        end
+
+        # 周期番号 -> 日時(平気)
+        #
+        # @param [Numeric] cn 周期番号
+        # @param [Numeric] time0 日時の初期近似値(ダミー)
+        #
+        # @return [Numeric] ユリウス日
+        #
+        def cn_to_time_(cn, time0=nil)
+          year, mean_term = (cn / 12.0 - @longitude_shift - @year_epoch).divmod(1)
+          @day_epoch + _winter_solstice(year) + _year_length_(year) * mean_term
+        end
+      end
+
       #
-      # @return [Numeric] 補正量 / @denominator
+      # 授時暦の定朔の計算
       #
-      def delta_m(t)
+      # @private
+      module MethodJ
+
+        include JujiMethods
+
+        # 周期番号 -> 日時(定朔)
+        #
+        # @param [Numeric] cn 周期番号
+        # @param [Numeric] time0 日時の初期近似値(ダミー)
+        #
+        # @return [Numeric] ユリウス日
+        #
+        def cn_to_time_(cn, time0=nil)
+          # 暦元天正冬至から当該経朔までの日数
+          mean_lunation  = cn * @lunation_length - @lunation_shift
+
+          # 当該経朔を含む近点年の暦元からの年数
+          year  = (mean_lunation / @year_length).floor
+          until (_perihelion(year)..._perihelion(year+1)).include?(mean_lunation)
+            year += _perihelion(year) > mean_lunation ? -1 : +1
+          end
+
+          # 定朔
+          solar_unit  = _year_length(year) / @year_length
+          mean_motion = (@lunar_mean_motion - @solar_weight / solar_unit) * 10000_0000
+          @day_epoch + mean_lunation - send('_anomaly_' + @anomaly_method.downcase, mean_lunation, year, solar_unit, mean_motion)
+        end
+
+        # 経朔 - 定朔 ( A 方式 - 差分)
+        def _anomaly_a(mean_lunation, year, solar_unit, mean_motion)
+
+          # 盈縮差(太陽の中心差) / (日 / 10000_0000)
+          solar_anomalies = (0..@solar_weight).to_a.map {|day|
+            solar_unit * equation_of_centre(((mean_lunation + day - _perihelion(year)) / solar_unit) % @year_length, @s)
+          }
+          solar_anomalies[1] ||= solar_anomalies[0]
+
+          # 遅速差(月の中心差) / (日 / 10000_0000)
+          lunar_anomalies = [0,1].map {|day|
+            equation_of_centre(((mean_lunation + day + @anomalistic_month_shift) % @anomalistic_month_length) / @lunar_unit, @m)
+          }
+
+          # 経朔 - 定朔
+          (lunar_anomalies[0] - solar_anomalies[0]) / ((lunar_anomalies[1] - lunar_anomalies[0]) -
+                                                       (solar_anomalies[1] - solar_anomalies[0]) + mean_motion)
+        end
+
+        # 経朔 - 定朔 ( B 方式 - 微分)
+        def _anomaly_b(mean_lunation, year, solar_unit, mean_motion)
+
+          # 盈縮差(太陽の中心差) / (日 / 10000_0000)
+          solar_anomalies = (0..@solar_weight).to_a.map {|diff|
+            solar_unit * equation_of_centre(((mean_lunation - _perihelion(year)) / solar_unit) % @year_length, @s, diff)
+          }
+          solar_anomalies[1] ||= 0
+
+          # 遅速差(月の中心差) / (日 / 10000_0000)
+          lunar_anomalies = [0,1].map {|diff|
+            equation_of_centre(((mean_lunation + @anomalistic_month_shift) % @anomalistic_month_length) / @lunar_unit, @m, diff)
+          }
+
+          # 経朔 - 定朔
+          (lunar_anomalies[0] - solar_anomalies[0]) / (lunar_anomalies[1] / @lunar_unit -
+                                                       solar_anomalies[1] /  solar_unit + mean_motion)
+        end
+
+        # 経朔 - 定朔 ( C 方式 - 幾何学的補正)
+        def _anomaly_c(mean_lunation, year, solar_unit, mean_motion)
+          diff = 0
+          loop do
+
+            # 盈縮差(太陽の中心差) / (日 / 10000_0000)
+            solar_anomaly = solar_unit * equation_of_centre(((mean_lunation - diff - _perihelion(year)) / solar_unit) % @year_length, @s)
 
-        t *= @denominator
-        a0 = b0 = t0 = 0
-        @m.each do |v|
-          t0, b0 = v
-          break if t <= t0
-          t  -= t0
-          a0 += b0
+            # 遅速差(月の中心差) / (日 / 10000_0000)
+            lunar_anomaly = equation_of_centre(((mean_lunation - diff + @anomalistic_month_shift) % @anomalistic_month_length) / @lunar_unit, @m)
+
+            # 次の差分
+            next_diff = (lunar_anomaly - solar_anomaly) / mean_motion
+            return next_diff if (next_diff - diff).abs < @anomaly_precision
+            diff = next_diff
+          end 
         end
 
-        # 補正値
-        a0 + (b0.to_f * t.to_i / t0.to_i + 0.5).floor
+        # 中心差およびその時間微分
+        def equation_of_centre(mean_anomaly, table, differential=0)
+          table.each do |range, base, *coefficients|
+            if range.include?(mean_anomaly)
+              diff = mean_anomaly - base
+              diff = diff.abs if coefficients[0].size[0] == 0
+              return coefficients[differential].inject(0) {|sum, coefficient| sum * diff + coefficient}
+            end
+          end
+          raise RangeError, 'Mean anomaly out of range: ' + mean_anomaly.to_s
+        end
+
+        # 立成の作成
+        def _initialize_rissei
+          @year_length              =  @year_length.to_f                  # 暦元の冬至年 / 日
+          @year_span                = (@year_span || 1).to_i              # 冬至年の改訂周期 / 年
+          @anomalistic_year_shift   = (@anomalistic_year_shift || 0).to_f # 暦應(冬至から近日点通過までの日数)
+          @lunation_length          =  @lunation_length.to_f              # 朔実(朔望月)
+          @lunation_shift           =  @lunation_shift.to_f               # 閏應(暦元前経朔から暦元天正冬至までの日数)
+          @lunar_mean_motion        =  @lunar_mean_motion.to_f            # 月平行(恒星天に対する月の平均運動 / 日)
+          @anomalistic_month_length =  @anomalistic_month_length.to_f     # 転終(近点月)
+          @anomalistic_month_shift  =  @anomalistic_month_shift.to_f      # 転應(暦元前近/遠地点通過から暦元天正冬至までの日数)
+          @anomaly_method           =  @anomaly_method || 'a'             # (経朔-定朔)の計算方法
+          @anomaly_precision        = (@anomaly_precision || 1.0E-5).to_f # c 方式 での収束判定誤差 / 日
+          @solar_weight             =  @solar_weight   || 0               # (経朔-定朔)の計算で用いる実行差での太陽盈縮の重み(0:非考慮,1:考慮)
+          @lunar_unit               =  @lunar_unit.to_f                   # 太陰遅速計算用招差法定数の時間の単位(限)
+          @m                        =  _rissei_j(@m)                      # 太陰遅速計算用招差法定数
+          @s                        =  _rissei_j(@s)                      # 太陽盈縮計算用招差法定数
+        end
+
+        # 招差法用の表の生成
+        def _rissei_j(table)
+          table.map {|range, base, coefficients|
+            sign = range.last == base ? -1 : +1
+            [range, base, coefficients.reverse,
+             (1...coefficients.size).to_a.reverse.map {|i| sign * i * coefficients[i]}]
+          }
+        end
+      end
+
+      # 近点月
+      # @return [Numeric]
+      #attr_reader :anomalistic_month_length
+
+      # 元期の近点角
+      # @return [Numeric]
+      #attr_reader :anomalistic_month_shift
+
+      # 日時 -> 周期番号(唐代の定朔の暦法用 cn_to_time(1L) を使用する)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [Numeric] 周期番号
+      #
+      def time_to_cn(t)
+        return super unless @cycle_number_1m
+        time = @is_dynamical ? +t : t.to_f
+        cn0  = time * @cycle_number_1m + @cycle_number_0m
+        root(cn0, time) {|cn| cn_to_time(cn)}
       end
 
+      # 当該日付の月の位相の変化範囲(唐代の定朔の暦法用 cn_to_time(1L) を使用する)
+      #
+      # @param [When::TM::TemporalPosition] date 日付
+      #
+      # @return [Array<Numeric>] 当該日付の月の位相の変化範囲
+      #
+      def phase_range(date)
+        date = date.floor
+        [date, date.succ].map {|d|
+          t  = d.to_f
+          c  = (60.0 * ((t - CYCLE_0M) * @cycle_number_1m + @cycle_number_0m)).floor
+          t0 = t1 = nil
+          loop do
+            t0 = cn_to_time( c   / 60.0)
+            t1 = cn_to_time((c+1)/ 60.0)
+            if t0 > t
+              c -= 1
+            elsif t1 <= t
+              c += 1
+            else
+              break
+            end
+          end
+          (c + (t-t0) / (t1-t0)) / 60.0
+        }
+      end
+
+      private
+
       # オブジェクトの正規化
       def _normalize(args=[], options={})
         super
         if @formula == '1L'
-          @anomalistic_month_length =  @anomalistic_month_length.to_r
-          @anomalistic_month_shift  = (@anomalistic_month_shift||0).to_r
-          @s = @s.map {|item| item.dup}
-          (0...@s.size).each do |i|
-            @s[i-1][1,0] = @year_length / @s.size + (@s[i][0]-@s[i-1][0]) / @denominator
-          end
-          (0...@s.size).each do |i|
-            @s[i].shift
+
+          # 月の位相の計算
+          @method ||= @year_span ? 'J' : 'A'
+          instance_eval("class << self; include const_get('Method#{@method.upcase}'); end")
+
+          # 立成の初期化
+          _initialize_rissei
+
+        elsif @year_span
+          # 太陽黄経の計算(消長あり)
+          class << self; include MethodS; end
+          @year_span = @year_span.to_i
+        end
+
+        if self.kind_of?(JujiMethods)
+          case @year_span
+          when 0,1; class << self; include JujiMethods::Y; end
+        # when 10 ; class << self; include JujiMethods::D; end
+          else    ; class << self; include JujiMethods::C; end
           end
         end
       end
@@ -253,6 +683,42 @@ module When
     #
     class ChineseSolar < EphemerisBasedSolar
 
+      class << self
+        #
+        # 太陰太陽暦の定義から太陰太陽暦と太陽暦(節月)の組を作る
+        #
+        # @private
+        def twin(area, definition)
+          definition.inject([]) {|list, cal|
+            if cal.kind_of?(Array) && cal[0] == ChineseLuniSolar
+              solar_name  = cal[1].sub(/=?\]/, '(節月)=]')
+              lunisolar   = cal.dup << "twin:#{area}::" + solar_name.gsub(/(name:\[|=?\])/,'')
+              solar       = cal.dup << "twin:#{area}::" + cal[1].gsub(/(name:\[|=?\])/,'')
+              solar[0..1] = [ChineseSolar, solar_name]
+              list << lunisolar << solar
+            else
+              list << cal
+            end
+          }
+        end
+
+        #
+        # 盈縮差の表の時間の単位を調整する
+        #
+        # @private
+        def change_unit(unit, definition)
+          definition.map {|line|
+            line.map {|item|
+              case item
+              when Range   ; Range.new(item.first*unit, item.last*unit, item.exclude_end?)
+              when Numeric ; item*unit
+              else         ; item
+              end
+            }
+          }
+        end
+      end
+
       #
       # @return [When::CalendarTypes::ChineseLuniSolar] 対で用いる太陰太陽暦の名前
       #
@@ -297,12 +763,20 @@ module When
       #
       # @param [Numeric] y 年
       #
-      # @return [Array<Numeric:月番号>, Array<Numeric:中気のない月の月番号>]
+      # @return [Array<Numeric:月番号>, Hash<Numeric:含む中気の数=>Numeric:月番号>]
       #
       def intercalary_pattern(y)
-        m = _base_month(y)
-        l = _base_ids(y)
-        [l, (0...l.size).to_a.map {|i| _intercalary?(m+i) ? l[i] : nil}.compact]
+        m  = _base_month(y)
+        l  = _base_ids(y)
+        c  = {0=>[], 1=>[], 2=>[]}
+        d0 = Residue.mod(_new_month(m)-1) {|n| _new_epoch(n)}[0]
+        l.size.times do |i|
+          d1 = Residue.mod(_new_month(m+i+1)-1) {|n| _new_epoch(n)}[0]
+          c[d1-d0] << l[i]
+          d0 = d1
+        end
+        c.delete(1)
+        [l, c]
       end
 
       private
@@ -427,480 +901,6 @@ module When
       end
     end
 
-    Chinese = [{}, When::BasicTypes::M17n, [
-      "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
-      "locale:[=ja:, en=en:, alias]",
-      "area:[中国,China]",
-
-      [ChineseLuniSolar,
-        'name:[黄帝暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1228331',
-        'intercalary_month:12'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[顓頊暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1171396&longitude_shift=-1/8',
-        'time_basis:+00,+#{P:12}',
-        'intercalary_month:9'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[夏暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1328411',
-        'intercalary_month:12'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[殷暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1149071',
-        'intercalary_month:12'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[周暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1128251',
-        'intercalary_month:12'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[魯暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1048991',
-        'intercalary_month:12'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[太初暦]',
-        'formula:MeanLunation?year_length=562120/1539&lunation_length=2392/81&day_epoch=1683431'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[四分暦]',
-        'formula:MeanLunation?year_length=1461/4&lunation_length=27759/940&day_epoch=1662611'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[乾象暦]',
-        'formula:MeanLunation?year_length=215130/589&lunation_length=43026/1457&day_epoch=-898129'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[景初暦]',
-        'formula:MeanLunation?year_length=673150/1843&lunation_length=134630/4559&day_epoch=+330191'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[三紀暦]',
-        'formula:MeanLunation?year_length=895220/2451&lunation_length=179044/6063&day_epoch=-28760989'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[玄始暦]',
-        'formula:MeanLunation?year_length=2629759/7200&lunation_length=2629759/89052&day_epoch=-20568349'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[元嘉暦]',
-        'formula:MeanLunation?year_length=111035/304&lunation_length=22207/752&day_epoch=-200089&longitude_shift=-1/12' # 春分の1ヶ月前
-      ],
-
-      [ChineseLuniSolar,
-        'name:[大明暦]',
-        'formula:MeanLunation?year_length=14423804/39491&lunation_length=116321/3939&day_epoch=-17080189'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[正光暦]',
-        'formula:MeanLunation?year_length=2213377/6060&lunation_length=2213377/74952&day_epoch=-59357929'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[興和暦]',
-        'formula:MeanLunation?year_length=6158017/16860&lunation_length=6158017/208530&day_epoch=-105462049'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[天保暦]',
-        'formula:MeanLunation?year_length=8641687/23660&lunation_length=8641687/292635&day_epoch=-38447089'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[天和暦]',
-        'formula:MeanLunation?year_length=8568631/23460&lunation_length=8568631/290160&day_epoch=-317950249'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[大象暦]',
-        'formula:MeanLunation?year_length=4745247/12992&lunation_length=1581749/53563&day_epoch=-13244449'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[開皇暦]',
-        'formula:MeanLunation?year_length=37605463/102960&lunation_length=5372209/181920&day_epoch=-1506155749'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[大業暦]',
-        'formula:MeanLunation?year_length=15573963/42640&lunation_length=33783/1144&day_epoch=-519493909'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[平朔儀鳳暦]',
-        'formula:MeanLunation?year_length=122357/335&lunation_length=39571/1340&day_epoch=-96608689'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[麟徳暦]',
-        'time_basis:+00,+#{P:00}',
-        {'formula'=>['12S', '1L'].map {|f|
-          Ephemeris::ChineseTrueLunation.new({
-            'formula'                  => f,
-            'day_epoch'                => -96608689,
-            'year_length'              => '122357/335',
-            'lunation_length'          =>  '39571/1340',
-            'anomalistic_month_length' =>  '443077/16080', # 27.0 + (743.0+1.0/12)/1340,
-            's'                        => [
-              [    0.0,     0,  +3.9546, -0.0372], # 冬至
-              [ -722.0,   +54,  +3.4091, -0.0372], # 大寒
-              [-1340.0,  +100,  +2.8636, -0.0372], # 小寒
-              [-1854.0,  +138,  +2.3181, +0.0372], # 立春
-              [-2368.0,  +176,  +2.8636, +0.0372], # 雨水
-              [-2986.0,  +222,  +3.4091, +0.0372], # 啓蟄
-              [-3708.0,  +276,  -3.7220, +0.0329], # 春分
-              [-2986.0,  +222,  -3.2086, +0.0329], # 清明
-              [-2368.0,  +176,  -2.6952, +0.0329], # 穀雨
-              [-1854.0,  +138,  -2.1818, -0.0329], # 立夏
-              [-1340.0,  +100,  -2.6952, -0.0329], # 小満
-              [ -722.0,   +54,  -3.2086, -0.0329], # 芒種
-              [    0.0,     0,  -3.7220, +0.0329], # 夏至
-              [ +722.0,   -54,  -3.2086, +0.0329], # 小暑
-              [+1340.0,  -100,  -2.6952, +0.0329], # 大暑
-              [+1854.0,  -138,  -2.1818, -0.0329], # 立秋
-              [+2368.0,  -176,  -2.6952, -0.0329], # 処暑
-              [+2986.0,  -222,  -3.2086, -0.0329], # 白露
-              [+3708.0,  -276,  +3.9546, -0.0372], # 秋分
-              [+2986.0,  -222,  +3.4091, -0.0372], # 寒露
-              [+2368.0,  -176,  +2.8636, -0.0372], # 霜降
-              [+1854.0,  -138,  +2.3181, +0.0372], # 立冬
-              [+1340.0,  -100,  +2.8636, +0.0372], # 小雪
-              [ +722.0,   -54,  +3.4091, +0.0372], # 大雪
-            ],
-            'm'                       => [
-              [1340.0,   -134],                    #   １日
-              [1340.0,   -117],                    #   ２日
-              [1340.0,    -99],                    #   ３日
-              [1340.0,    -78],                    #   ４日
-              [1340.0,    -56],                    #   ５日
-              [1340.0,    -33],                    #   ６日
-              [1191.0,     -9],                    #   ７日
-              [ 149.0,      0],                    #   ７日
-              [1340.0,    +14],                    #   ８日
-              [1340.0,    +38],                    #   ９日
-              [1340.0,    +62],                    #   10日
-              [1340.0,    +85],                    #   11日
-              [1340.0,   +104],                    #   12日
-              [1340.0,   +121],                    #   13日
-              [1042.0,   +102],                    #   14日
-              [ 298.0,    +29],                    #   14日
-              [1340.0,   +128],                    #   15日
-              [1340.0,   +115],                    #   16日
-              [1340.0,    +95],                    #   17日
-              [1340.0,    +74],                    #   18日
-              [1340.0,    +52],                    #   19日
-              [1340.0,    +28],                    #   20日
-              [ 892.0,     +4],                    #   21日
-              [ 448.0,      0],                    #   21日
-              [1340.0,    -20],                    #   22日
-              [1340.0,    -44],                    #   23日
-              [1340.0,    -68],                    #   24日
-              [1340.0,    -89],                    #   25日
-              [1340.0,   -108],                    #   26日
-              [1340.0,   -125],                    #   27日
-              [ 743.0+1.0/12, -71]                 #   28日
-            ]
-         })
-        },
-        'doyo'   => (Rational( 4,15) +  244) / 1340
-       }
-      ],
-
-      [ChineseLuniSolar,
-        'name:[大衍暦]',
-        'time_basis:+00,+#{P:03}',
-        {'formula'=>['12S', '1L'].map {|f|
-          Ephemeris::ChineseTrueLunation.new({
-            'formula'                  => f,
-            'day_epoch'                => -35412747829,
-            'year_length'              => '1110343/3040',
-            'lunation_length'          =>   '89773/3040',
-            'anomalistic_month_length' =>  '6701279/243200', # 27.0 +(1685.0+79.0/80)/3040,
-            'anomalistic_month_shift'  =>  '1/2',
-            's'                        => [
-              [    0.0,     0, +13.4524, -0.1886], # 冬至
-              [-2353.0,  +176, +10.5564, -0.1634], # 大寒
-              [-4198.0,  +314,  +8.0408, -0.1446], # 小寒
-              [-5588.0,  +418,  +5.8160, -0.1318], # 立春
-              [-6564.0,  +491,  +3.7987, -0.1240], # 雨水
-              [-7152.0,  +535,  +1.9265, -0.1240], # 啓蟄
-              [-7366.0,  +551,  -0.2048, -0.1178], # 春分
-              [-7152.0,  +535,  -1.9968, -0.1190], # 清明
-              [-6564.0,  +491,  -3.7956, -0.1240], # 穀雨
-              [-5588.0,  +418,  -5.6626, -0.1324], # 立夏
-              [-4198.0,  +314,  -7.6555, -0.1436], # 小満
-              [-2353.0,  +176,  -9.9405, -0.1436], # 芒種
-              [    0.0,     0, -12.0819, +0.1436], # 夏至
-              [+2353.0,  -176,  -9.7018, +0.1324], # 小暑
-              [+4198.0,  -314,  -7.5450, +0.1240], # 大暑
-              [+5588.0,  -418,  -5.5634, +0.1190], # 立秋
-              [+6564.0,  -491,  -3.7038, +0.1178], # 処暑
-              [+7152.0,  -535,  -1.8954, +0.1178], # 白露
-              [+7366.0,  -551,  +0.1783, +0.1240], # 秋分
-              [+7152.0,  -535,  +2.0042, +0.1318], # 寒露
-              [+6564.0,  -491,  +3.8950, +0.1446], # 霜降
-              [+5588.0,  -418,  +5.9214, +0.1634], # 立冬
-              [+4198.0,  -314,  +8.1610, +0.1886], # 小雪
-              [+2353.0,  -176, +10.9010, +0.1886], # 大雪
-            ],
-            'm'                       => [
-              [3040.0,   +297],                    #   １日
-              [3040.0,   +259],                    #   ２日
-              [3040.0,   +220],                    #   ３日
-              [3040.0,   +180],                    #   ４日
-              [3040.0,   +139],                    #   ５日
-              [3040.0,    +97],                    #   ６日
-              [2701.0,    +48],                    #   ７日
-              [ 339.0,     -6],                    #   ７日
-              [3040.0,    -64],                    #   ８日
-              [3040.0,   -106],                    #   ９日
-              [3040.0,   -148],                    #   10日
-              [3040.0,   -189],                    #   11日
-              [3040.0,   -229],                    #   12日
-              [3040.0,   -267],                    #   13日
-              [2363.0,   -231],                    #   14日
-              [ 677.0,    -66],                    #   14日
-              [3040.0,   -289],                    #   15日
-              [3040.0,   -250],                    #   16日
-              [3040.0,   -211],                    #   17日
-              [3040.0,   -171],                    #   18日
-              [3040.0,   -130],                    #   19日
-              [3040.0,    -87],                    #   20日
-              [2024.0,    -36],                    #   21日
-              [1016.0,    +18],                    #   21日
-              [3040.0,    +73],                    #   22日
-              [3040.0,   +116],                    #   23日
-              [3040.0,   +157],                    #   24日
-              [3040.0,   +198],                    #   25日
-              [3040.0,   +237],                    #   26日
-              [3040.0,   +276],                    #   27日
-              [1686.0,   +165]                     #   28日
-            ]
-         })
-        },
-        'doyo'   => (Rational(13,30) +  531) / 3040
-       }
-      ],
-
-      [ChineseLuniSolar,
-        'name:[五紀暦]',
-        'time_basis:+00,+#{P:06}',
-        {'formula'=>['12S', '1L'].map {|f|
-          Ephemeris::ChineseTrueLunation.new({
-            'formula'                  => f,
-            'day_epoch'                => -96608689,
-            'year_length'              => '122357/335',
-            'lunation_length'          =>  '39571/1340',
-            'anomalistic_month_length' =>  '1366156/49580', # 27.0 + (743.0+5.0/37)/1340,
-            's'                        => [
-              [    0.0,     0,  +5.9668, -0.0843], # 冬至
-              [-1037.0,   +78,  +4.6652, -0.0721], # 大寒
-              [-1850.0,  +139,  +3.5656, -0.0653], # 小寒
-              [-2463.0,  +185,  +2.5583, -0.0590], # 立春
-              [-2893.0,  +217,  +1.6375, -0.0532], # 雨水
-              [-3152.0,  +236,  +0.8384, -0.0532], # 啓蟄
-              [-3246.0,  +243,  -0.0972, -0.0505], # 春分
-              [-3152.0,  +236,  -0.8480, -0.0534], # 清明
-              [-2893.0,  +217,  -1.6517, -0.0561], # 穀雨
-              [-2463.0,  +185,  -2.5057, -0.0584], # 立夏
-              [-1850.0,  +139,  -3.3781, -0.0643], # 小満
-              [-1037.0,   +78,  -4.3954, -0.0643], # 芒種
-              [    0.0,     0,  -5.3592, +0.0643], # 夏至
-              [+1037.0,   -78,  -4.2877, +0.0584], # 小暑
-              [+1850.0,  -139,  -3.3459, +0.0561], # 大暑
-              [+2463.0,  -185,  -2.4475, +0.0534], # 立秋
-              [+2893.0,  -217,  -1.5966, +0.0505], # 処暑
-              [+3152.0,  -236,  -0.8185, +0.0505], # 白露
-              [+3246.0,  -243,  +0.0858, +0.0532], # 秋分
-              [+3152.0,  -236,  +0.8505, +0.0590], # 寒露
-              [+2893.0,  -217,  +1.6943, +0.0653], # 霜降
-              [+2463.0,  -185,  +2.6205, +0.0721], # 立冬
-              [+1850.0,  -139,  +3.6007, +0.0843], # 小雪
-              [+1037.0,   -78,  +4.8330, +0.0843], # 大雪
-            ],
-            'm'                       => [
-              [1340.0,   -135],                    #   １日
-              [1340.0,   -117],                    #   ２日
-              [1340.0,    -99],                    #   ３日
-              [1340.0,    -78],                    #   ４日
-              [1340.0,    -56],                    #   ５日
-              [1340.0,    -33],                    #   ６日
-              [1191.0,     -8],                    #   ７日
-              [ 149.0,     +1],                    #   ７日
-              [1340.0,    +14],                    #   ８日
-              [1340.0,    +38],                    #   ９日
-              [1340.0,    +62],                    #   10日
-              [1340.0,    +85],                    #   11日
-              [1340.0,   +103],                    #   12日
-              [1340.0,   +118],                    #   13日
-              [1042.0,   +105],                    #   14日
-              [ 298.0,    +30],                    #   14日
-              [1340.0,   +128],                    #   15日
-              [1340.0,   +115],                    #   16日
-              [1340.0,    +95],                    #   17日
-              [1340.0,    +74],                    #   18日
-              [1340.0,    +52],                    #   19日
-              [1340.0,    +28],                    #   20日
-              [ 892.0,     +6],                    #   21日
-              [ 448.0,     -3],                    #   21日
-              [1340.0,    -20],                    #   22日
-              [1340.0,    -42],                    #   23日
-              [1340.0,    -65],                    #   24日
-              [1340.0,    -89],                    #   25日
-              [1340.0,   -109],                    #   26日
-              [1340.0,   -125],                    #   27日
-              [ 743.0+5.0/37, -75]                 #   28日
-            ]
-         })
-        },
-        'doyo'   => (Rational( 4,15) +  244) / 1340
-       }
-      ],
-
-      [ChineseSolar,
-        'name:[正元暦(節月)]',
-        'formula:MeanLunation?year_length=399943/1095&lunation_length=32336/1095&day_epoch=-145149709'
-      ],
-
-      [ChineseLuniSolar,
-        'name:[宣明暦]',
-        'time_basis:+00,+#{P:06}',
-        {'formula'=>['12S', '1L'].map {|f|
-          Ephemeris::ChineseTrueLunation.new({
-            'formula'                  => f,
-            'day_epoch'                => -2580308749,
-            'year_length'              => '3068055/8400',
-            'lunation_length'          =>  '248057/8400',
-            'anomalistic_month_length' =>  '23145819/840000', # 27.0 + 4658.19 / 8400,
-            'anomalistic_month_shift'  =>  '1/2',
-            's'                        => [
-              [     0.0,     0, +33.4511, -0.3695], # 冬至
-              [ -6000.0,  +449, +28.0389, -0.3606], # 大寒
-              [-11000.0,  +823, +22.6998, -0.3519], # 小寒
-              [-15000.0, +1122, +17.8923, -0.4068], # 立春
-              [-18000.0, +1346, +11.7966, -0.3998], # 雨水
-              [-19800.0, +1481,  +5.7986, -0.3998], # 啓蟄
-              [-20400.0, +1526,  -0.2433, -0.3779], # 春分
-              [-19800.0, +1481,  -6.1254, -0.3634], # 清明
-              [-18000.0, +1346, -12.2048, -0.2987], # 穀雨
-              [-15000.0, +1122, -16.9060, -0.2919], # 立夏
-              [-11000.0,  +823, -21.5362, -0.2854], # 小満
-              [ -6000.0,  +449, -26.0498, -0.2854], # 芒種
-              [     0.0,     0, -30.3119, +0.2854], # 夏至
-              [ +6000.0,  -449, -25.8126, +0.2919], # 小暑
-              [+11000.0,  -823, -21.2454, +0.2987], # 大暑
-              [+15000.0, -1122, -17.0296, +0.3634], # 立秋
-              [+18000.0, -1346, -11.4744, +0.3779], # 処暑
-              [+19800.0, -1481,  -5.6429, +0.3779], # 白露
-              [+20400.0, -1526,  +0.1432, +0.3998], # 秋分
-              [+19800.0, -1481,  +6.1488, +0.4068], # 寒露
-              [+18000.0, -1346, +12.6336, +0.3519], # 霜降
-              [+15000.0, -1122, +17.8043, +0.3606], # 立冬
-              [+11000.0,  -823, +23.0590, +0.3695], # 小雪
-              [ +6000.0,  -449, +28.4618, +0.3695], # 大雪
-            ],
-            'm'                       => [
-              [ 8400.0,   +830],                    # 進１日
-              [ 8400.0,   +726],                    #   ２日
-              [ 8400.0,   +606],                    #   ３日
-              [ 8400.0,   +471],                    #   ４日
-              [ 8400.0,   +337],                    #   ５日
-              [ 8400.0,   +202],                    #   ６日
-              [ 7465.0,    +53],                    #   ７日
-              [  935.0,     -7],                    #   ７日
-              [ 8400.0,    -82],                    #   ８日
-              [ 8400.0,   -224],                    #   ９日
-              [ 8400.0,   -366],                    #   10日
-              [ 8400.0,   -509],                    #   11日
-              [ 8400.0,   -643],                    #   12日
-              [ 8400.0,   -748],                    #   13日
-              [ 6529.095, -646],                    #   14日
-              [ 8400.0,   -830],                    # 退１日
-              [ 8400.0,   -726],                    #   ２日
-              [ 8400.0,   -598],                    #   ３日
-              [ 8400.0,   -464],                    #   ４日
-              [ 8400.0,   -329],                    #   ５日
-              [ 8400.0,   -195],                    #   ６日
-              [ 7465.0,    -53],                    #   ７日
-              [  935.0,     +7],                    #   ７日
-              [ 8400.0,    +82],                    #   ８日
-              [ 8400.0,   +225],                    #   ９日
-              [ 8400.0,   +366],                    #   10日
-              [ 8400.0,   +501],                    #   11日
-              [ 8400.0,   +628],                    #   12日
-              [ 8400.0,   +740],                    #   13日
-              [ 6529.095, +646]                     #   14日
-            ]
-         })
-        },
-        'doyo'   => (Rational( 1, 2) + 1468) / 8400
-       }
-      ],
-
-      [ChineseSolar,
-        'name:[前貞享暦(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.675000&year_epoch=1684&year_length=365.241696&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[貞享暦(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.689990&year_epoch=1684&year_length=365.241696&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[貞享補暦A(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.903800&year_epoch=1684&year_length=365.241696&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[貞享補暦B(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.622300&year_epoch=1684&year_length=365.241696&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[宝暦暦(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.622100&year_epoch=1684&year_length=365.241696&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[修正宝暦暦(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.762200&year_epoch=1684&year_length=365.241766&year_delta=1'
-      ],
-
-      [ChineseSolar,
-        'name:[寛政暦(節月)]',
-        'formula:VariableYearLengthMethod?day_epoch=2336118.720200&year_epoch=1684&year_length=365.242360&year_delta=0'
-      ]
-    ].inject([]) {|list, cal|
-      if cal.kind_of?(Array) && cal[0] == ChineseLuniSolar
-        solar_name  = cal[1].sub(/\]/, '(節月)]')
-        lunisolar   = cal.dup << 'twin:Chinese::' + solar_name.gsub(/(name:\[|\])/,'')
-        solar       = cal.dup << 'twin:Chinese::' + cal[1].gsub(/(name:\[|\])/,'')
-        solar[0..1] = [ChineseSolar, solar_name]
-        list << lunisolar << solar
-      else
-        list << cal
-      end
-    }]
-
     #
     # 太平天国 2.1.1-3.2.30
     #