when_exe 0.2.100 → 0.3.1

data/lib/when_exe/region/sun.rb

@@ -0,0 +1,210 @@
+# -*- coding: utf-8 -*-
+=begin
+  Copyright (C) 2011-2013 Takashi SUGA
+
+  You may use and/or modify this file according to the license described in the LICENSE.txt file included in this archive.
+=end
+
+module When
+  module Ephemeris
+
+    #------------------------------------------
+    #   太陽・地球 位置計算用データ
+    #------------------------------------------
+
+    # 太陽の地心黄経
+    # @private
+    P3L    =
+      [[LIN   ,  279.0357964 ,     360.007693787, 3.02500E-8],
+       [SIN   ,  353.06      ,     719.981008 ,   0.02009 ,-2.98832E-8],
+       [SIN   ,  248.64      ,     -19.34136  ,  -0.0048  ],
+       [SIN   ,  285.0       ,     329.6446699,   0.0020  ,+8.0E-9    ],
+       [SIN   ,  334.2       ,   -4452.671109 ,   0.0018  ,-0.1634E-6 ],
+
+       [SIN   ,  293.7       ,      -0.2020   ,   0.00178 ],
+       [SIN   ,  242.4       ,     450.3688   ,   0.0015  ,+1.4E-9    ],
+       [SIN   ,  211.1       ,     225.1844   ,   0.0013  ,+0.7E-9    ],
+       [SIN   ,  208.0       ,     659.2893   ,   0.0008  ,+1.6E-8    ],
+       [SIN   ,   53.5       ,      90.3783   ,   0.0007  ],
+
+       [SIN   ,   12.1       ,     -30.35     ,   0.0007  ],
+       [SIN   ,  239.1       ,     337.18     ,   0.0006  ],
+       [SIN   ,   10.1       ,      -1.50     ,   0.0005  ],
+       [SIN   ,   99.1       ,     -22.81     ,   0.0005  ],
+       [SIN   ,  264.8       ,     315.56     ,   0.0004  ],
+
+       [SIN   ,  233.8       ,     299.30     ,   0.0004  ],
+       [SIN   ,  198.1       ,     720.02     ,  -0.0004  ],
+       [SIN   ,  349.6       ,    1079.97     ,   0.0003  ],
+       [SIN   ,  241.2       ,     -44.43     ,   0.0003  ]]
+    #  [SIN   ,   65.0       ,     -44.43     ,   0.0003  ]], 旧版
+    # 太陽の動径の対数
+    # @private
+    P3Q    =
+      [[COS   ,  356.53      ,     359.990504 ,  -0.007261 ,-1.49416E-8],
+       [COST  ,  356.53      ,     359.990504 ,   0.0000002,-1.49416E-8],
+       [COS   ,  353.1       ,     719.981008 ,  -0.000091 ,-2.98832E-8],
+       [COS   ,  205.8       ,    4452.671109 ,   0.000013 ,-0.1634E-6 ],
+       [COS   ,   62.0       ,     450.4      ,   0.000007 ,+1.4E-9    ],
+       [COS   ,  105.0       ,     329.6      ,   0.000007 ,+8.0E-9    ]]
+
+    #
+    # The Sun
+    #
+    class Sun < Datum
+      class << self
+
+        include Ephemeris
+
+        S0=[356.531, 359.990504, -1.49416E-8] # 太陽の中心差の位相
+        S1=[  1.9159, -4.8E-5  , -1.44444E-9] # 太陽の中心差の振幅
+
+        # 視黄経 / CIRCLE
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def true_longitude(t)
+          y = julian_year_from_1975(+t)
+          return (sind(S0[0]+y*(S0[1]+y*S0[2]))*(S1[0]+y*(S1[1]+y*S1[2])) +
+                  trigonometric(y, P3L)) / 360.0 + 1974.0
+        end
+
+        # 距離 / AU
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def pi(t)
+          return 10.0**(0.000030+trigonometric(julian_year_from_1975(+t), P3Q))
+        end
+
+        # 平均黄経 / CIRCLE
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def mean_longitude(t)
+          return trigonometric(julian_year_from_1975(+t), P3L, 0.0, 1) / 360.0 + 1974.0
+        end
+      end
+
+      # 位置 (黄道座標)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [When::Ephemeris::Coords]
+      #
+      def _coords(t)
+        Coords.polar(0, 0, 0, 0)
+      end
+
+      # 視黄経 / CIRCLE
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [0]
+      #
+      def true_longitude(t)
+        0
+      end
+
+      # 平均黄経 / CIRCLE
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [0]
+      #
+      def mean_longitude(t)
+        0
+      end
+
+      # 平均運動 / (DEG/YEAR)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [0]
+      #
+      def mean_motion
+        0
+      end
+
+      # オブジェクトの生成
+      # @private
+      def initialize(*args)
+        options = [args.pop] if args[-1].kind_of?(Hash)
+        surface_radius, aberration, luminosity, *rest = args
+        surface_radius ||= 696000.0
+        aberration     ||=      0.00000
+        luminosity     ||=      4.58
+        args  = [surface_radius, aberration, luminosity] + rest
+        args += options if options
+        super(*args)
+      end
+    end
+
+    #
+    # The Earth
+    #
+    class Earth < Datum
+
+      # 地球の位置 (黄道座標)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [When::Ephemeris::Coords]
+      #
+      def _coords(t)
+        t  = +t
+        radius = Sun.pi(t)
+        dl = 0.5 + @aberration / radius / 360 # 略算式に含まれる光行差をキャンセルする
+        Coords.polar(Sun.true_longitude(t)+dl, 0, radius, Sun.mean_longitude(t)+dl)
+      end
+
+      # 地球の平均運動 / (DEG/YEAR)
+      #
+      # @return [Numeric]
+      #
+      def mean_motion
+        P3L[0][2]
+      end
+
+      # オブジェクトの生成
+      # @private
+      def initialize(*args)
+        options = {
+          'shape' => [0.998327112, +0.001676399, -0.000003519, -11.514/60.0],
+          'sid'   => [6.697375, 2400.0513369, 0.0000259],
+          'zeros' => {'Z'=>0, 'A'=>-0.58555, '0'=>-0.85255, 'T'=>-7.36111},
+          'air'   => [10.0, 0.40, 20.0]
+        }
+        options.update(args.pop) if args[-1].kind_of?(Hash)
+        surface_radius, aberration, luminosity, *rest = args
+        surface_radius ||= 6378.14
+        aberration     ||=    0.0056932
+        luminosity     ||=   -3.50
+        args = [surface_radius, aberration, luminosity] + rest + [options]
+        super(*args)
+      end
+    end
+
+    #
+    # Typical Geometrical Datum
+    #
+    # provisionally same as When::Ephemeris::Earth
+    #
+    class JGD2000 < Earth # TODO
+    end
+  end
+end

data/lib/when_exe/region/thai.rb

@@ -0,0 +1,227 @@
+# -*- coding: utf-8 -*-
+=begin
+  Copyright (C) 2011-2012 Takashi SUGA
+
+  You may use and/or modify this file according to the license described in the LICENSE.txt file included in this archive.
+=end
+
+module When
+  class BasicTypes::M17n
+
+    ThaiTerms = [self, [
+      "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
+      "locale:[=en:, ja=ja:, alias]",
+      "names:[ThaiTerms]",
+      "[ThaiLuniSolar=en:Thai_lunar_calendar, タイ太陰太陽暦=ja:%E3%83%81%E3%83%A3%E3%83%B3%E3%83%88%E3%83%A9%E3%82%AB%E3%83%86%E3%82%A3]"
+    ]]
+  end
+
+  module CalendarTypes
+
+    #
+    # The Calendar of Thai people in China
+    #
+    class ThaiC < TableBased
+
+      NormalIDs    =  '1,1<,2,2<,3,3<,4,4<,5,5<,6,6<,7,7<,8,8<,9,9<,10,10<,11,11<,12,12<'
+      LeapIDs      =  '1,1<,2,2<,3,3<,4,4<,5,5<,6,6<,7,7<,8,8<,9&,9*,9,9<,10,10<,11,11<,12,12<'
+
+      # ソンクラーン - 太陽の白羊宮入り
+      #
+      # @param [Integer] y 年
+      #
+      # @return [Integer] ソンクラーンのユリウス日
+      #
+      def songkran_(y)
+        e = _eph(y)
+        return @origin_of_LSC - 1 + e['H'][0]
+      end
+
+      private
+
+      # ThaiC+オブジェクトの性質定義を初期設定します。
+      #
+      # オブジェクトの正規化
+      # @private
+      def _normalize(args=[], options={})
+        @label ||= When.Resource('_m:ThaiTerms::ThaiLuniSolar')
+
+        Rational
+        @mean_month    ||=       29 + Rational(  373,    703)
+        @mean_year     ||=      365 + Rational(25876, 100000)
+        @epoch_in_CE   ||=      638
+        @origin_of_MSC ||=        0 # 638 + 543
+        @origin_of_LSC ||=  1954168
+        @origin_of_MSC   =  @origin_of_MSC.to_i
+        @origin_of_LSC   =  @origin_of_LSC.to_i
+
+        intercalary_month = When.Resource('_m:IndianTerms::IntercalaryMonth::*')
+        @indices ||= [
+          Coordinates::Index.new({:branch=>{0   =>  intercalary_month[0],
+                                            0.5 =>  intercalary_month[1],
+                                           -1.5 =>  intercalary_month[2],
+                                           -1   =>  intercalary_month[3]},
+                                  :trunk =>When.Resource('_m:IndianTerms::LunarMonth::*')}),
+          Coordinates::Index.new
+        ]
+        @rule_table  ||=  {
+          354 => {'Length'=>[15,15,15,14]*6,                            'IDs' => NormalIDs},
+          355 => {'Length'=>[15,15,15,14]*3 + [15]*4 + [15,15,15,14]*2, 'IDs' => NormalIDs},
+          356 => {'Length'=>[15,15,15,14]*3 + [15]*7 + [14] + [15]*4,   'IDs' => NormalIDs},
+          383 => {'Length'=>[15,14,15,15]*6 + [15,14],                  'IDs' => LeapIDs},
+          384 => {'Length'=>[15,15,15,14]*4 + [15]*2 + [15,15,15,14]*2, 'IDs' => LeapIDs}
+        }
+
+        super
+      end
+
+      # 年初の通日
+      #
+      # @param [Array<Numeric>] date ( 年 )
+      #
+      # @return [Integer] 年初の通日
+      #
+      def _sdn_(date)
+        y  = +date[0]
+        e0 = _thai(y,  0)[0]
+        e1 = _thai(y+1,0)[0]
+        e0['T'] -= 1 if (e1['T']-e0['T'] == 353)
+        return @origin_of_LSC - 1 + e0['T'] - 30*3 - 29*2
+      end
+
+      # y で指定した年の暦定数を返します。
+      #
+      # @param [Integer] y 年
+      #
+      # @return [Hash]
+      #
+      def _eph(y) # C
+        h  = (y+4).divmod(9)[0]                   # y
+        h  = (y-h).divmod(3)[0]                   # z
+        h  = (y+1-h).divmod(2)[0]                 # r (2 => h?)
+        h  = (36525876*y+149049-h).divmod(100000) # s
+        a  = (h[0]*11+633-((h[0]+7368).divmod(8878))[0]).divmod(692)
+        m  = (h[0]+a[0]+0).divmod(30)
+        return {'H'=>h, 'A'=>a, 'M'=>m}
+      end
+
+      #  y で指定した年の平閏、8月の大小を判定します。
+      #
+      # @param [Integer] y   年
+      # @param [Integer] dir 前後の閏状態を表すコード
+      #
+      # @return [Array<Hash, String>]
+      #
+      def _thai(y,dir)
+          # 諸元の計算 
+        e = _eph(y)
+        k = 800*125 - e['H'][1]
+        t = e['M'][1]
+
+          # Jyaistha 大の判定 
+        b = ((e['A'][1]<126) && (k<=207*125)) || ((e['A'][1]<137) && (k> 207*125))
+
+          # 閏年の判定
+        c = ((t<6) || (t>24))
+        c = true  if ((t==24) && (_thai(y+1,+1)[1]< 'C'))
+        c = false if ((t==25) && (_thai(y+1,+1)[1]>='C'))
+
+          # Caitra 月 0 日 
+        e['T'] = e['H'][0] - t
+        e['T'] -= 1  if (t==0)
+        e['T'] -= 29 if (c && (t<=6))
+        e['b'] = b
+        e['c'] = c
+        e['t'] = t
+        return e,' '
+      end
+    end
+
+    #
+    # The Calendar of central Thailand
+    #
+    class ThaiB < ThaiC
+
+      YEAR_TYPE = {'A'=>354, 'B'=>355, 'C'=>384, 'D'=>384, 'E'=>384}
+
+      NormalIDs    =  '1,1<,2,2<,3,3<,4,4<,5,5<,6,6<,7,7<,8,8<,9,9<,10,10<,11,11<,12,12<'
+      LeapIDs      =  '1,1<,2,2<,3,3<,4,4<,5,5<,6,6<,7,7<,8&,8*,8,8<,9,9<,10,10<,11,11<,12,12<'
+
+      private
+
+      # オブジェクトの正規化
+      #
+      def _normalize(args=[], options={})
+        @rule_table ||= {
+          354 => {'Length'=>[15,14,15,15]*6,    'IDs' => NormalIDs},
+          355 => {'Length'=>[15,14,15,15]*3 + [15]*4 + [15,14,15,15]*2, 'IDs' => NormalIDs},
+          384 => {'Length'=>[15,14,15,15]*4 + [15]*2 + [15,14,15,15]*2, 'IDs' => LeapIDs  }
+        }
+
+        super
+      end
+
+      # 年初の通日
+      #
+      # @param [Array<Numeric>] date ( 年 )
+      #
+      # @return [Integer] 年初の通日
+      #
+      def _sdn_(date)
+        y  = +date[0]
+        e0,dir0  = _thai(y,  0)
+        e1,dir1  = _thai(y+1,0)
+        e0['T'] -= 1 if (e1['T'] - e0['T'] < YEAR_TYPE[dir0])
+        return @origin_of_LSC + e0['T'] - 29*2 - 30*2
+      end
+
+      # y で指定した年の暦定数を返します。
+      #
+      # @param [Integer] y 年
+      #
+      # @return [Hash]
+      #
+      def _eph(y) # C
+        h = (y*292207+373+800*1).divmod(800)
+      # u = (h[0]+2611).divmod(3232)
+        a = (h[0]*11+650).divmod(692)
+        m = (h[0]+a[0]+0).divmod(30)
+        h[1] *= 125
+        return {'H'=>h, 'A'=>a, 'M'=>m}
+      end
+
+      #  y で指定した年の平閏、8月の大小を判定します。
+      #
+      # @param [Integer] y   年
+      # @param [Integer] dir 前後の閏状態を表すコード
+      #
+      # @return [Array<Hash, String>]
+      #
+      def _thai(y,dir)
+        # 諸元の計算
+        e,w = super(y,dir)
+
+        # カレンダーのタイプ 
+        if (e['c'])
+          if (!e['b'])
+            return  e,'C'   # 普通の閏年   
+          elsif ((0<e['t']) && (e['t']<6))
+            return  e,'D'   # 前年大の閏年 
+          else
+            return  e,'E'   # 翌年大の閏年 
+          end
+        else
+          if (e['b'])
+            return  e,'B'   # 大の平年 
+          elsif ((dir != -1) && (_thai(y+1,+1)[1]=='D'))
+            return  e,'B'   # 大の平年 
+          elsif ((dir != +1) && (_thai(y-1,-1)[1]=='E'))
+            return  e,'B'   # 大の平年 
+          else
+            return  e,'A'   # 小の平年 
+          end
+        end
+      end
+    end
+  end
+end

data/lib/when_exe/region/tibetan.rb

@@ -0,0 +1,233 @@
+# -*- coding: utf-8 -*-
+=begin
+  Copyright (C) 2011-2012 Takashi SUGA
+
+  You may use and/or modify this file according to the license described in the LICENSE.txt file included in this archive.
+=end
+
+module When
+  class BasicTypes::M17n
+
+    TibetanTerms = [self, [
+      "namespace:[en=http://en.wikipedia.org/wiki/, ja=http://ja.wikipedia.org/wiki/]",
+      "locale:[=en:, ja=ja:, alias]",
+      "names:[TibetanTerms]",
+      "[Tibetan=en:Tibetan_calendar, 時輪暦=ja:%E3%83%81%E3%83%99%E3%83%83%E3%83%88%E6%9A%A6]",
+
+      [self,
+        "names:[IntercalaryMonth=en:Intercalation, 閏月]",
+        "[%s=,       %s宿月=  ]",
+        "[Adika %s=, 閏%s宿月=]"
+      ],
+
+      [self,
+        "names:[IntercalaryDay=en:Intercalation, 閏日=ja:%E9%96%8F]",
+        "[Double %s=,      欠=]",
+        "[Intercalary %s=, 重=]"
+      ],
+
+      [self,
+        "names:[Month, 月=ja:%E6%9C%88_(%E6%9A%A6)]",
+        "[Margasirsa, 觜=]",
+        "[Pausha,     鬼=]",
+        "[Magha,      星=]",
+        "[Phalguna,   翼=]",
+        "[Chaitra,    角=]",
+        "[Vaisakha,   氐=]",
+        "[Jyeshta,    心=]",
+        "[Ashada,     箕=]",
+        "[Sravana,    牛=]",
+        "[Bhadrapada, 室=]",
+        "[Asvina,     婁=]",
+        "[Kartika,    昴=]"
+      ]
+    ]]
+  end
+
+  module CalendarTypes
+
+    #
+    # Tibetan Calendar
+    #
+    class Tibetan < EphemerisBasedLuniSolar
+
+      #
+      # 定数
+      #
+      Y0 = 1827           # 元期の西暦年
+      D0 = 2388440        # 元期のユリウス日
+      V0 = (6.0+45/60.0)  # 白羊宮入りのオフセット
+      Am = (2.0+1.0/126)  # 月の近地点
+
+      Rem  = 1            # 余りのインデックス
+      Quot = 0            # 商のインデックス
+      Ddd  = 0            # 差分のインデックス
+      Sum  = 1            # 積算値のインデックス
+      Jm = [              # 月の中心差
+        [ 5, 0], [ 5, 5], [ 5,10], [ 5,15], [ 4,19], [ 3,22], [ 2,24],
+        [ 1,25], [-1,24], [-2,22], [-3,19], [-4,15], [-5,10], [-5, 5],
+        [-5, 0], [-5,-5]
+      ]
+      Js = [   # 太陽の中心差
+        [ 6, 0], [ 6, 6], [ 4,10], [ 1,11], [-1,10], [-4, 6], [-6, 0],
+        [-6,-6]
+      ]
+
+      private
+
+      # オブジェクトの正規化
+      #
+      #    @type      = 体系派なら 1, 作用派なら 2
+      #    @parameter = 太陽と月の位相計算用のテーブル
+      #
+      def _normalize(args=[], options={})
+        intercalary_month = When.Resource('_m:TibetanTerms::IntercalaryMonth::*')
+        intercalary_day   = When.Resource('_m:TibetanTerms::IntercalaryDay::*')
+        @indices ||= [
+          Coordinates::Index.new({:branch=>{0=>intercalary_month[0], 1=>intercalary_month[1]},
+                                  :trunk=>When.Resource('_m:TibetanTerms::Month::*'), :shift=>2}),
+          Coordinates::Index.new({:branch=>{-2=>intercalary_day[0], -1=>intercalary_day[1]}})
+        ]
+        @label       ||= When.Resource('_m:TibetanTerms::Tibetan')
+        @epoch_in_CE ||= 0
+        @type        ||= 1
+        @parameter = case @type.to_i
+          when 1
+            {'M0'=>(60+15),
+             'Wo'=>Z(29,31,50,0,480),
+             'W0'=>Z( 3,37,43,2,140),
+             'Wd'=>Z( 0,59, 3,4, 16),
+             'So'=>K( 2,10,58,1, 17),
+             'S0'=>K(24,59, 6,1, 41),
+             'Sd'=>K( 0, 4,21,5, 43),
+             'A0'=>22.0
+            }
+          when 2
+            {'M0'=>(64+15),
+             'Wo'=>Z(29,31,50,0,  0),
+             'W0'=>Z( 3,21,20,0,  0),
+             'Wd'=>Z( 0,59, 3,4,  0),
+             'So'=>L( 2,10,58,2, 10),
+             'S0'=>L(25,42,12,1, 11),
+             'Sd'=>L( 0, 4,21,5,  9),
+             'A0'=>(22.0+28.0/126)
+            }
+          end
+        super
+      end
+
+      # @private
+      def Z(z,q,d,m,s)
+        ((((s/707.0+m)/6.0+d)/60.0+q)/60.0+z)
+      end
+
+      # @private
+      def K(k,q,d,m,s)
+        ((((s/67.0+m)/6.0+d)/60.0+q)/60.0+k)
+      end
+
+      # @private
+      def L(k,q,d,m,s)
+        ((((s/13.0+m)/6.0+d)/60.0+q)/60.0+k)
+      end
+
+      #
+      # 太陽または月の位置の中心差による差分
+      # @private
+      def _sn(t,m,p)
+        t0 = t.floor
+        r  = t0.divmod(m)
+        s  = (1.0/60.0)
+        s  = -s unless r[Quot] % 2 == 0
+        return s * (p[r[Rem]][Sum] + (t-t0) * p[r[Rem]+1][Ddd])
+      end
+
+      # 積月計算
+      #
+      # @param  [Integer] y 年
+      # @param  [Integer] m 年(デフォルト 0)
+      #
+      # @return [Integer] 積月
+      #
+      #   m に 0 以外を指定することで、積月の跳び(閏)を検出できる
+      #
+      def _new_year_month_(y, m=0)
+        mm = (y-Y0) * 12 + (m - 2)
+        mm += (mm * 2 + @parameter['M0']).divmod(65)[Quot]
+        return mm
+      end
+
+      # 月初の通日
+      #
+      # @param  [Integer] m 通月
+      #
+      # @return [Integer] 月初の通日
+      #
+      def _new_month_(m)
+        sdn = [0,-1].map {|d| _new_month_day(m, d)}
+        sdn[0]-sdn[1] == 1 ? sdn[0] : sdn[0]-1
+      end
+
+      #
+      # 積日計算
+      #
+      # @param  [Integer] m 通月
+      # @param  [Integer] d 日
+      #
+      # @return [Integer] 積日
+      #
+      #   積日の跳び(重日)、ダブり(欠日)を検出できる
+      #
+      def _new_month_day_(m, d)
+        d   += 1
+        w    = m * @parameter['Wo'] + @parameter['W0'] + d * @parameter['Wd']  # 中曜
+        s    = m * @parameter['So'] + @parameter['S0'] + d * @parameter['Sd']  # 中日
+        sdn  = D0 + w
+        sdn += _sn(m*Am+@parameter['A0']+d, 14, Jm)
+        sdn -= _sn((s-V0)*(60.0/135.0), 6, Js)
+        return sdn.floor
+      end
+
+      #
+      # 日時要素の翻訳表の取得
+      #
+      # @overload _ids_(date)
+      #   @param [Array<Numeric>] date ( 年 )
+      #   @return [Array<When::Coordinates::Pair>] 1年の月の配置の翻訳表
+      #
+      # @overload _ids_(date)
+      #   @param [Array<Numeric>] date ( 年 月 )
+      #   @return [Array<When::Coordinates::Pair>] 1月の日の配置の翻訳表
+      #   @note 月は 0 始まりの通番
+      #
+      def _ids_(date)
+        y, m = date
+        y    = +y
+        mm   = _new_year_month(y)
+        return (_table(mm, 12, 0, 1) {|i| _new_year_month(y,i)}) unless m
+        mm  += m
+        table = _table(_new_month_day(mm,0), 30, 1, 0) {|i| _new_month_day(mm,i)}
+        table[1,0] =  Pair._force_pair(1,1) unless _new_month(mm) == _new_month_day(mm,0)
+        return table
+      end
+
+      # 日時要素の翻訳表の作成
+      # @private
+      def _table(b0, n, k0, k1)
+        table = [Pair._force_pair(1,0)]
+        (1..n).each do |i|
+          b1 = yield(i)
+          case b1-b0
+          when 0 ; table[-1] = Pair._force_pair(i, -2)
+          when 1 ; table << Pair._force_pair(i+1, 0 )
+          when 2 ; table << Pair._force_pair(i+k0,k1) << Pair._force_pair(i+1, k0)
+          else   ; raise ArgumentError, "Irregal date span: #{b1-b0}"
+          end
+          b0 =  b1
+        end
+        table.pop while table[-1].trunk > n
+        return table
+      end
+    end
+  end
+end