when_exe 0.3.7 → 0.3.8

data/lib/when_exe/calendartypes.rb

@@ -86,7 +86,18 @@ module When::CalendarTypes
     def universal_time(sdn=nil)
       return super - @time_standard.localtime_difference unless sdn
       time = When::TM::JulianDate._d_to_t(sdn-0.5)
-      @time_standard.to_dynamical_time(time) - When::TimeStandard.to_dynamical_time(time) - @time_standard.localtime_difference
+      @time_standard.to_dynamical_time(time) - When::TimeStandard.to_dynamical_time(time)
+    end
+
+    # この時法の時刻を128秒単位の実数に変換する
+    #
+    # @param [Array<Numeric>] clk_time
+    # @param [Integer] sdn 参照事象の通し番号(ダミー)
+    #
+    # @return [Numeric]
+    #
+    def to_local_time(clk_time, sdn=nil)
+      super - universal_time(sdn)
     end
 
     #
@@ -177,13 +188,13 @@ module When::CalendarTypes
         'rule_table'    => table
       }
       list  = ''
-      while range.include?(date[YEAR])
-        month = date[MONTH] * 1
+      while range.include?(date[When::YEAR])
+        month = date[When::MONTH] * 1
         char  = Pattern[month..month]
-        char  = char.downcase unless date.length(MONTH) == length
+        char  = char.downcase unless date.length(When::MONTH) == length
         list += char
         succ  = date + duration
-        unless date[YEAR] == succ[YEAR]
+        unless date[When::YEAR] == succ[When::YEAR]
           table << list
           list = ''
         end
@@ -845,7 +856,7 @@ module When::CalendarTypes
 
     # オブジェクトの正規化
     #
-    #   YearLengthTableBased+オブジェクトの性質定義を初期設定します。
+    #   PatternTableBasedLuniSolarWithEphemeris オブジェクトの性質定義を初期設定します。
     #
     def _normalize(args=[], options={})
       @engine   = When.Calendar(@engine)
@@ -1050,20 +1061,56 @@ module When::CalendarTypes
   end
 
   #
-  # 年初を太陽黄経で決定する暦
+  # 年初を太陽黄経または別の暦で決定する暦
   #
   class YearLengthTableBased < TableBased
 
+    module CalendarBased
+
+      # 年初の通日(別の暦使用)
+      #
+      # @param [Numeric] year 年
+      #
+      # @return [Integer] 年初の通日
+      #
+      def _new_year_sdn(year)
+        @engine._coordinates_to_number(year + @diff_to_CE - @engine.diff_to_CE, @engine_month, @engine_day)
+      end
+
+      # オブジェクトの正規化
+      def _normalize_engine
+        Rational
+        @engine_month = @engine_month ? @engine_month.to_i - @indices[-2].base : 0
+        @engine_day   = @engine_day   ? @engine_day.to_i   - @indices[-1].base : 0
+        @engine       = When.Calendar(@engine || When::Gregorian)
+      end
+    end
+
     # 天体暦アルゴリズム
     #
     # @return [Array<When::Ephemeris::Formula>]
     #
     attr_reader :formula
 
-    #protected
-
     private
 
+    # 年初の通日(天体暦使用)
+    #
+    # @param [Numeric] year 年
+    #
+    # @return [Integer] 年初の通日
+    #
+    def _new_year_sdn(year)
+      solar_sdn(@formula[0].cn_to_time(year.to_f + @cycle_offset) + @day_offset)
+    end
+
+    # オブジェクトの正規化
+    def _normalize_engine
+      Rational
+      @cycle_offset = (@cycle_offset||0).to_r
+      @day_offset   = (@day_offset||0).to_r
+    end
+
     # 年初の通日
     #
     # @param [Array<Numeric>] date ( y )
@@ -1073,26 +1120,61 @@ module When::CalendarTypes
     # @return [Integer] 年初の通日
     #
     def _sdn_(date)
-      y = +date[0]
-      t = @formula[0].cn_to_time(y.to_f + @cycle_offset)
-      return solar_sdn(t + @day_offset)
+      _new_year_sdn(+date[0])
     end
 
     # オブジェクトの正規化
     #
-    #   YearLengthTableBased+オブジェクトの性質定義を初期設定します。
+    #   YearLengthTableBased オブジェクトの性質定義を初期設定します。
     #
     def _normalize(args=[], options={})
 
-      Rational
-      @cycle_offset = (@cycle_offset||0).to_r
-      @day_offset   = (@day_offset||0).to_r
-      @formula      = 'Formula?formula=1S'
-
+      extend CalendarBased unless @formula || @location || @long || @lat || @alt || @time_basis
+      _normalize_engine
       super
     end
   end
 
+  #
+  # 年初を特定の日の日の出で決定する暦
+  #
+  class TableBasedWithSunrise < YearLengthTableBased
+
+    # 年初の通日(天体暦使用)
+    #
+    # @param [Numeric] year 年
+    #
+    # @return [Integer] 年初の通日
+    #
+    def _new_year_sdn(year)
+      event_time  = @formula[0].cn_to_time(year.to_f + @cycle_offset)
+      event_date  = (event_time + 0.5 + @formula[0].long/360.0).floor
+      sunset_time = @formula[0].sunrise(event_date)
+      event_date -= 1 if sunrise_time > event_time
+      event_date + @day_offset
+    end
+  end
+
+  #
+  # 年初を特定の日の日の入りで決定する暦
+  #
+  class TableBasedWithSunset < YearLengthTableBased
+
+    # 年初の通日(天体暦使用)
+    #
+    # @param [Numeric] year 年
+    #
+    # @return [Integer] 年初の通日
+    #
+    def _new_year_sdn(year)
+      event_time  = @formula[0].cn_to_time(year.to_f + @cycle_offset)
+      event_date  = (event_time + 0.5 + @formula[0].long/360.0).floor
+      sunset_time = @formula[0].sunset(event_date)
+      event_date += 1 if sunset_time <= event_time
+      event_date + @day_offset
+    end
+  end
+
   # 月日の配当が太陽または月の位置によって決定される暦
   #
   #   Calendar based on the ephemeris of the Sun or the Moon
@@ -1154,13 +1236,13 @@ module When::CalendarTypes
     #
     def _length(date)
       y, m = date
-      if (m)
+      if m
         #  指定した月に含まれる日の数を返します。
         m += @months_in_year * +y
-        return _new_month(m+1) - _new_month(m)
+        _new_month(m+1) - _new_month(m)
       else
         #  指定した年に含まれる月の数を返します。
-        return @months_in_year
+        @months_in_year
       end
     end
 
@@ -1175,7 +1257,7 @@ module When::CalendarTypes
     def _sum_(date)
       y, = date
       m = @months_in_year * +y
-      return _new_month(m+@months_in_year) - _new_month(m)
+      _new_month(m+@months_in_year) - _new_month(m)
     end
 
     private
@@ -1206,7 +1288,7 @@ module When::CalendarTypes
     # @return [Integer] 月初の通日
     #
     def _new_month_(m)
-      return solar_sdn(@formula[0].cn_to_time(m + @cycle_offset))
+      solar_sdn(@formula[0].cn_to_time(m + @cycle_offset))
     end
 
     private
@@ -1239,7 +1321,7 @@ module When::CalendarTypes
     # @return [Integer] 月初の通日
     #
     def _new_month_(m)
-      return lunar_sdn(@formula[-1].cn_to_time(m + @cycle_offset))
+      lunar_sdn(@formula[-1].cn_to_time(m + @cycle_offset))
     end
 
     private
@@ -1307,13 +1389,13 @@ module When::CalendarTypes
     #
     def _length(date)
       y, m = date
-      if (m)
+      if m
         #  指定した月に含まれる日の数を返します。
         m += _new_year_month(+y)
-        return _new_month(m+1) - _new_month(m)
+        _new_month(m+1) - _new_month(m)
       else
         #  指定した年に含まれる月の数を返します。
-        return _ids([y]).length
+        _ids([y]).length
       end
     end
 
@@ -1329,7 +1411,7 @@ module When::CalendarTypes
     #
     def _sum_(date)
       y = +date[0]
-      return _new_month(_new_year_month(y+1)) - _new_month(_new_year_month(y))
+      _new_month(_new_year_month(y+1)) - _new_month(_new_year_month(y))
     end
 
     # 太陽月初の通日
@@ -1368,216 +1450,4 @@ module When::CalendarTypes
       super
     end
   end
-
-  #
-  # 日時要素の境界 - Border
-  #
-  class Border < When::TM::ReferenceSystem
-    #
-    # 境界の振舞
-    #
-    # @return [Numeric]
-    #
-    #  Pair(-1,+1) - 暦年/暦日が進む(境界が前年/日にあり、境界後が当年/日の扱いになる)
-    #
-    #  Pair( 0, 0) - 暦年/暦日が戻る(境界が当年/日にあり、境界前が前年/日の扱いになる)
-    #
-    def behavior
-      @border[0]
-    end
-
-    # 境界の取得
-    #
-    # @param [Array<Numeric>] date 境界を計算する年/日
-    # @param [When::TM::ReferenceSystem] frame 使用する暦法/時法
-    #
-    # @return [Array<Numeric>] その年/日の境界
-    #
-    def border(date=[], frame=nil)
-      last = date.length-1
-      return @border if last<0
-      bDate  = date[0..last] + @border[(last+1)..-1]
-      branch = @border[last] * 0
-      bDate[last] = When::Coordinates::Pair.new(date[last] * 1, branch)
-      return bDate
-    end
-
-    # 境界の正規化
-    #
-    # @param [Array<Numeric>] date 境界を計算する年/日
-    # @param [When::TM::ReferenceSystem] frame 使用する暦法/時法
-    #
-    # @return [Array<Numeric>] その年/日の境界
-    #
-    def _adjust_epoch(date, frame=nil)
-      s_date  = date.dup
-      e_date  = border([+date[0]], frame)
-      branch  = behavior * 0
-      branch += 1 if (s_date[1..-1] <=> e_date[1..-1]) < 0
-      s_date[0] = When::Coordinates::Pair.new(+s_date[0]-branch, branch)
-      return s_date
-    end
-
-    # 日付の補正
-    # @private
-    def _date_adjust(source)
-      source
-    end
-
-    private
-
-    # 要素の正規化
-    def _normalize(args=[], options={})
-      @border = When::Coordinates::Pair._en_pair_date_time(@border) if (@border.kind_of?(String))
-    end
-  end
-
-  #
-  # 日時要素の境界 - 年/日によって、異なる境界を使用する場合
-  #
-  class MultiBorder < Border
-
-    #
-    # 境界の配列
-    # @return [Array<When::CalendarTypes::Border>]
-    attr_reader :borders
-
-    #
-    # 境界の振舞
-    #
-    # @return [Numeric]
-    #
-    #   Pair(-1,+1) - 暦年/暦日が進む(境界が前年/日にあり、境界後が当年/日の扱いになる)
-    #
-    #   Pair( 0, 0) - 暦年/暦日が戻る(境界が当年/日にあり、境界前が前年/日の扱いになる)
-    #
-    def behavior
-      @borders[0][:border].behavior
-    end
-
-    # 境界の取得
-    #
-    # @param [Array<Numeric>] date 境界を計算する年/日
-    # @param [When::TM::ReferenceSystem] frame 使用する暦法/時法
-    #
-    # @return [Array<Numeric>] その年/日の境界
-    #
-    def border(date=[], frame=nil)
-      last = date.length-1
-      return @borders[0][:boder] if (last<0)
-      @borders.each do |border|
-        return border[:border].border(date, frame) if date[0] >= border[:key]
-      end
-      date[0..last]
-    end
-
-    # 境界の正規化
-    #
-    # @param [Array<Numeric>] date 境界を計算する年/日
-    # @param [When::TM::ReferenceSystem] frame 使用する暦法/時法
-    #
-    # @return [Array<Numeric>] その年/日の境界
-    #
-    def _adjust_epoch(date, frame=nil)
-      @borders.each do |border|
-        next unless date[0] >= border[:key]
-        s_date  = date.dup
-        e_date  = border[:border].border(date[0..0], frame)
-        branch  = border[:border].behavior * 0
-        branch += 1 if ((s_date[1..-1] <=> e_date[1..-1]) < 0)
-        s_date[0] = When::Coordinates::Pair.new(+s_date[0]-branch, branch)
-        return s_date
-      end
-      date
-    end
-
-    private
-
-    # 要素の正規化
-    def _normalize(args=[], options={})
-      if @borders.kind_of?(String)
-        list = @borders.split(/(\(.+?\))/)
-        list.shift if list[0]==''
-        list.unshift(-Float::MAX) unless list[0] =~ /\(/
-        list.push('0-1-1') if list[-1] =~ /\(/
-        @borders = []
-        loop do
-          key, border, *list = list
-          break unless key
-          key = $1.to_i if key.kind_of?(String) && /\((.+?)\)/ =~ key
-          border = "_c:Border?border=#{border}" unless border =~ /\A[A-Z_]/i
-          border = When.Calendar(border)
-          @borders << {:key=>key, :border=>border}
-        end
-      end
-      @borders = @borders.sort_by {|border| -border[:key]}
-    end
-  end
-
-  #
-  # 日時要素の境界 - 日の出,日の入り
-  #
-  class DayBorder < Border
-
-    # 境界の取得
-    #
-    # @param [Array<Numeric>] date 境界を計算する日
-    # @param [When::TM::ReferenceSystem] clock 使用する時法
-    #
-    # @return [Array<Numeric>] その日の境界
-    #
-    # @note 属性 @event によって境界を計算する (see {When::Ephemeris::Formula#day_event})
-    #
-    def border(date=[], clock=When::UTC)
-      return @border unless date[0] && clock.formula
-
-      time =
-        clock._number_to_coordinates(clock.second *
-          clock.time_standard.from_dynamical_time(
-            When::TM::JulianDate._d_to_t(
-              clock.formula.first.day_event(
-                clock.time_standard.to_dynamical_date(date[0] + @border[0]*0), @event, When.Resource('_ep:Sun'), @height
-              ))))
-
-      time[0] += When::TM::JulianDate::JD19700101
-      time[0]  = When::Coordinates::Pair.new(time[0]-@border[0]*0, @border[0]*0) unless @border[0]*0 == 0
-      clock._encode(time, false)
-    end
-
-    # 日付の補正
-    # @private
-    def _date_adjust(source)
-      source * 1 + @border[0] * 0
-    end
-  end
-
-  #
-  # 日時要素の境界 - 日の出
-  #
-  class Sunrise < DayBorder
-
-    private
-
-    # 要素の正規化
-    def _normalize(args=[], options={})
-      @border   = [0,0,0,0]
-      @event    = -1
-      @height ||= '0'
-    end
-  end
-
-  #
-  # 日時要素の境界 - 日の入り
-  #
-  class Sunset < DayBorder
-
-    private
-
-    # 要素の正規化
-    def _normalize(args=[], options={})
-      @border   = [When::Coordinates::Pair.new(+1,-1),0,0,0]
-      @event    = +1
-      @height ||= '0'
-    end
-  end
 end