when_exe 0.3.4 → 0.3.5

data/lib/when_exe/coordinates.rb

@@ -22,18 +22,18 @@ module When::Coordinates
                     When::HOUR=>'PT1H', When::MINUTE=>'PT1M', When::SECOND=>'PT1S'}
   MATCH     = {'NS'=>/(N|S|北緯|南緯)/, 'EW'=>/(E|W|東経|西経)/}
 
-  # 60進->10進変換
+  # 60進->10進変換(1/225度単位)
   #
   # @param [String]   src 60進法で表した方向付きの数値
   # @param [String]   dir 方向 ('NS' または 'EW')
   #
-  # @return [Numeric] 10進変換した数値 (src が nil なら0.0を、Numeric ならそのままsrcを返す)
+  # @return [Numeric] 10進変換した数値 (src が nil なら0.0を、Numeric なら 225*src を返す)
   #
-  def self.to_deg(src, dir)
+  def self.to_deg_225(src, dir)
     case src
     when String
       src = src.gsub(/_+/,'').gsub(/@/, '.')
-      return src.to_f unless (src =~ MATCH[dir])
+      return src.to_r * Spatial::DEGREE if (src =~ /E[-+]/ || src !~ MATCH[dir])
       sign  = ($1 == dir[1..1]) ? -1 : +1
       value = src.gsub(MATCH[dir], '').strip
       if ((value + "00000") =~ /^(\d+)\.(\d{2})(\d{2})(\d+)$/)
@@ -42,29 +42,53 @@ module When::Coordinates
       else
         deg, min, sec = value.split(/[^\d.]+/)
       end
-      return sign * (deg.to_i + (min||0).to_f/60 +  (sec||0).to_f/3600)
+      return sign * (deg.to_i *  Spatial::DEGREE +
+                (min||0).to_f * (Spatial::DEGREE/60.0) +
+                (sec||0).to_f * (Spatial::DEGREE/3600.0))
     when NilClass
       0.0
     when Numeric
-      src
+      src * Spatial::DEGREE
     else
       raise TypeError, "Invalid Location Type"
     end
   end
 
+  # 60進->10進変換(度単位)
+  #
+  # @param [String]   src 60進法で表した方向付きの数値
+  # @param [String]   dir 方向 ('NS' または 'EW')
+  #
+  # @return [Numeric] 10進変換した数値 (src が nil なら0.0を、Numeric ならそのままsrcを返す)
+  #
+  def self.to_deg(src, dir)
+    to_deg_225(src, dir) / Spatial::DEGREE
+  end
+
   # 10進->60進変換
   #
   # @param [Numeric] src 数値
   # @param [String]  dir 方向 ('NS' または 'EW')
+  # @param [Integer] round 秒の小数点以下最大桁数
   #
   # @return [String] 60進変換した数値
   #
-  def self.to_dms(src, dir)
+  def self.to_dms(src, dir, round=6)
     dir      = (src >= 0) ? dir[0..0] : dir[1..1]
-    deg, min  =  src.abs.divmod(1)
-    min, sec  = (60*min).divmod(1)
-    sec       = (60*sec).floor
-    (['N','S'].include?(dir) ? "%02d.%02d%02d%s" : "%03d.%02d%02d%s") % [deg, min, sec, dir]
+    deg, min  =     src.abs.divmod(1)
+    min, sec  =    (60*min).divmod(1)
+    sec       = (60*10**round*sec).round
+    fig = round + 2
+    round.times do
+      div, mod = sec.divmod(10)
+      if mod == 0
+        fig -= 1
+        sec  = div
+      else
+        break
+      end
+    end
+    (['N','S'].include?(dir) ? "%02d.%02d%0#{fig}d%s" : "%03d.%02d%0#{fig}d%s") % [deg, min, sec, dir]
   end
 
   #
@@ -503,17 +527,21 @@ module When::Coordinates
       #
       def succ
         value = @current
-        if (@z==@k) ||
-           (@count_limit.kind_of?(Numeric) && @count >= @count_limit) ||
+        if (@count_limit.kind_of?(Numeric) && @count >= @count_limit) ||
            (@error.kind_of?(Numeric) && @e && @error >= @e.abs)
           @current = nil
         else
-          @z = 1.0/(@z-@k)
-          @k = @z.floor
-          @e = @p[1].to_f/@q[1]-@x.to_f/@y
-          @current = [@p[1], @q[1], @e]
-          @p = [@p[1], @p[1]*@k + @p[0]]
-          @q = [@q[1], @q[1]*@k + @q[0]]
+          if @z==@k
+            @e = 0
+            @current = [@p[1], @q[1], 0]
+          else
+            @z = 1.0/(@z-@k)
+            @k = @z.floor
+            @e = @p[1].to_f/@q[1]-@x.to_f/@y
+            @current = [@p[1], @q[1], @e]
+            @p = [@p[1], @p[1]*@k + @p[0]]
+            @q = [@q[1], @q[1]*@k + @q[0]]
+          end
           @count += 1
         end
         return value
@@ -1244,6 +1272,14 @@ module When::Coordinates
         end
       end
 
+      # 設定情報を取得する
+      #
+      # @return [Hash] 設定情報
+      #
+      def _setup_info
+        {:location => @default_location}
+      end
+
       # デフォルトの空間位置を読みだす
       #
       # @return [When::Coordinates::Spatial]
@@ -1277,6 +1313,9 @@ module When::Coordinates
       end
     end
 
+    # @private
+    HashProperty = [:label, [:alt, 0.0], [:datum, When::Ephemeris::Earth], :ref]
+
     # Degree / Internal Location Unit(16")
     #
     #   (3600 を 2 の因数で割りつくした値を単位とする)
@@ -1336,6 +1375,32 @@ module When::Coordinates
     #
     attr_reader :tz
 
+    # 緯度文字列
+    #
+    # @param [Integer] round 秒の小数点以下最大桁数
+    #
+    # @return [String] 緯度文字列(DD.MMSSsss[NS])
+    #
+    def lat_s
+      When::Coordinates.to_dms(lat / When::Coordinates::Spatial::DEGREE, 'NS', round=6)
+    end
+
+    # 経度文字列
+    #
+    # @param [Integer] round 秒の小数点以下最大桁数
+    #
+    # @return [String] 経度文字列(DDD.MMSSsss[EW])
+    #
+    def long_s
+      When::Coordinates.to_dms(long / When::Coordinates::Spatial::DEGREE, 'EW', round=6)
+    end
+
+    # 高度 / m
+    #
+    # @return [Numeric]
+    # @return [:Center] 天体の中心の場合
+    #
+    attr_reader :alt
     # 観測地の地心距離 / kmを返します。
     #
     # @return [Numeric]
@@ -1395,17 +1460,20 @@ module When::Coordinates
 
     # 要素の正規化
     def _normalize(args=[], options={})
-      @long  = When::Coordinates.to_deg(@long, 'EW') * DEGREE if @long
-      @lat   = When::Coordinates.to_deg(@lat,  'NS') * DEGREE if @lat
-      @datum = When.Resource(@datum || 'Earth', '_ep:')
-      if @tz.kind_of?(String)
-        @label ||= @tz
-        @tz      = When::Parts::Timezone.tz_info[@tz]
-      end
+
+      # 時間帯による指定
+      @tz = When::Parts::Timezone.tz_info[@tz] if @tz.kind_of?(String)
       if @tz
-        @long  ||= @tz.longitude * DEGREE
-        @lat   ||= @tz.latitude  * DEGREE
+        @label ||= @tz.label
+        @long  ||= @tz.longitude
+        @lat   ||= @tz.latitude
       end
+
+      # データの整形
+      @label = When::BasicTypes::M17n.new(@label)        if @label.kind_of?(Hash)
+      @long  = When::Coordinates.to_deg_225(@long, 'EW') if @long
+      @lat   = When::Coordinates.to_deg_225(@lat,  'NS') if @lat
+      @datum = When.Resource(@datum || 'Earth', '_ep:')
       @long ||= 0.0
       @lat  ||= 0.0
       @alt    =
@@ -1416,6 +1484,21 @@ module When::Coordinates
         end
     end
 
+    #
+    # 空間位置オブジェクトの内容を Hash 化
+    #
+    # @param [Object] options When::Parts::Resource#to_h を参照
+    # @option options [Symbol] :method :to_h を label の内容の Hash 書き出しのために追加
+    #
+    # @return [Hash] Hash 化した空間位置オブジェクト
+    #
+    def _to_h(options={})
+      hash = super
+      hash[:long] = long_s
+      hash[:lat]  = lat_s
+      hash
+    end
+
     # 観測地の惑星中心を原点とする三次元座標
     #
     # @param [Numeric] t ユリウス日(Terrestrial Time)
@@ -1502,7 +1585,7 @@ module When::Coordinates
         # Formula
         instance_eval('class << self; attr_reader :formula; end') if @location && @border
         if respond_to?(:formula)
-          instance_eval('class << self; include When::Ephemeris::Formula::Methods; end')
+          instance_eval('class << self; include When::Ephemeris::Formula::ForwardedFormula; end')
           @formula ||= When::Ephemeris::Formula.new({:location=>@location})
           @formula   = When.Resource(Array(@formula), '_ep:')
         end
@@ -1587,10 +1670,10 @@ module When::Coordinates
 
     # 代表暦注
     #
-    # @return [When::CalendarTypes::CalendarNote]
+    # @return [When::CalendarNote]
     # @return [Array<Array<klass, Array<klass, method, block>>>] 最外側のArray要素は年・月・日に対応
     #
-    #   klass  [String, When::CalendarTypes::CalendarNote, When::Coordinates::Residue]
+    #   klass  [String, When::CalendarNote, When::Coordinates::Residue]
     #
     #   method [String, Symbol] (デフォルト 'day', 'month' or 'year' (対応する桁による))
     #
@@ -1644,6 +1727,13 @@ module When::Coordinates
       @_diff_to_CE ||= @epoch_in_CE ? @epoch_in_CE - @origin_of_MSC : 0
     end
 
+    # @private
+    #
+    # 対応する ::Date の start 属性
+    def _default_start
+      ::Date::GREGORIAN
+    end
+
     # protected
 
     #
@@ -1946,7 +2036,13 @@ module When::Coordinates
 
     # その他のメソッド
     #   When::Coordinates::Temporal で定義されていないメソッドは
-    #   処理を @note or @formula[0] (When::Ephemeris::Formula) に委譲する
+    #   処理を下記に移譲する(番号は優先順位)
+    #     When::CalendarNote
+    #       (1) @note
+    #       (2) SolarTerms
+    #       (3) LunarPhases
+    #     When::Ephemeris::Formula
+    #       (4)@formula[0]
     #
     def method_missing(name, *args, &block)
       unless When::Parts::MethodCash::Escape.key?(name)
@@ -1957,9 +2053,20 @@ module When::Coordinates
             end
           } unless When::Parts::MethodCash.escape(name)
           return @note.send(name, *(args + [self]), &block)
-        elsif When::Ephemeris::Formula.method_defined?(name)
+        end
+        ['SolarTerms', 'LunarPhases'].each do |note|
+          if When.CalendarNote(note).respond_to?(name)
+            instance_eval %Q{
+              def #{name}(*args, &block)
+                When.CalendarNote("#{note}").send("#{name}", *args, &block)
+              end
+            } unless When::Parts::MethodCash.escape(name)
+            return When.CalendarNote(note).send(name, *args, &block)
+          end
+        end
+        if When::Ephemeris::Formula.method_defined?(name)
           unless respond_to?(:forwarded_formula, true)
-            instance_eval('class << self; include When::Ephemeris::Formula::Methods; end')
+            instance_eval('class << self; include When::Ephemeris::Formula::ForwardedFormula; end')
             @formula ||= When::Ephemeris::Formula.new({:location=>@location})
             @formula   = When.Resource(Array(@formula), '_ep:')
           end

data/lib/when_exe/core/compatibility.rb

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 =begin
-  Copyright (C) 2013 Takashi SUGA
+  Copyright (C) 2013-2014 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
@@ -33,7 +33,7 @@ class String
     def to_r
       case self
       when /\.|E/i
-        to_f
+        to_f.to_r
       when /\//
         Rational
         Rational(*split(/\//).map {|v| v.to_i})
@@ -54,6 +54,26 @@ class String
   end
 end
 
+#
+# 浮動小数点数
+#
+class Float
+  unless const_defined?(:INFINITY)
+    # @private
+    INFINITY = MAX
+  end
+
+  unless method_defined?(:to_r)
+    Rational
+    # @private
+    def to_r
+      it  = When::Coordinates::Residue.new(self, 1).enum_for
+      res = it.succ while it.has_next?
+      Rational(*res[0..1])
+    end
+  end
+end
+
 #
 # 例外クラス
 #

data/lib/when_exe/core/extension.rb

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 =begin
-  Copyright (C) 2013 Takashi SUGA
+  Copyright (C) 2013-2014 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
@@ -19,7 +19,7 @@ class Time
   #
   # 対応する When::TM::JulianDate を生成
   #
-  # @param [Hash] options
+  # @param [Hash] options 以下の通り
   # @option options [When::TM::Clock] :clock
   # @option options [When::Parts::Timezone] :tz
   #
@@ -55,7 +55,27 @@ end
 #
 class Date
 
-  include When::TM::TemporalPosition::Conversion if When::TM.const_defined?(:TemporalPosition)
+  if When::TM.const_defined?(:TemporalPosition)
+
+    include When::TM::TemporalPosition::Conversion
+
+    #
+    # 対応する When::TM::CalDate or DateAndTime を生成
+    #
+    # @param [Hash] options 暦法や時法などの指定
+    #   see {When::TM::TemporalPosition._instance}
+    #
+    # @return [When::TM::CalDate, When::TM::DateAndTime]
+    #
+    # @note 暦法の指定がない場合、start メソッドの値によって
+    #       Julian / Gregorian / Civil 暦法を選択する
+    #
+    def tm_position(options={})
+      options[:frame] ||= When::CalendarTypes::Christian._default_start(self)
+      super(options)
+    end
+    alias :to_tm_position :tm_position
+  end
 
   alias :__method_missing :method_missing
 
@@ -130,6 +150,18 @@ class Numeric
   end
   alias :pair :to_pair
 
+  #
+  # self を秒数とみなして When::Parts::Timezone::Base を取得
+  #
+  # @return [When::Parts::Timezone::Base]
+  #
+  # @note core/extension
+  #
+  def clock
+    When.Clock(self)
+  end
+  alias :to_clock :clock
+
   # メソッド名に相当する単位で表した When::TM::IntervalLength
   # @method unit_interval_length
   # @return [When::TM::IntervalLength]
@@ -233,6 +265,18 @@ class String
   end
   alias :to_calendar :calendar
 
+  #
+  # self を時間帯文字列とみなして When::Parts::Timezone::Base を取得
+  #
+  # @return [When::Parts::Timezone::Base]
+  #
+  # @note core/extension
+  #
+  def clock
+    When.Clock(self)
+  end
+  alias :to_clock :clock
+
   #
   # self をプレフィクス '_c:' を省略した IRI とみなして When::TM::CalendarNote を取得
   #

data/lib/when_exe/ephemeris.rb

@@ -1147,35 +1147,7 @@ module When::Ephemeris
     #
     # 天体暦機能を When::TM::Calendar クラスに提供する
     #
-    module Methods
-
-      # 月の位相が指定の周期番号パターンになる最も近い過去の日時
-      #
-      # @param [Numeric] date ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] date
-      # @param [Numeric] n 相対周期番号(n=0 なら date または date の直前が基準)
-      # @param [Numeric] d 単位周期数
-      #
-      # @return [Numeric, When::TM::TemporalPosition] 周期番号が d で割って n になる日時
-      #
-      def nearest_past_from_lunar_phase(date, n=0, d=1)
-        @formula ||= When.Resource(['_ep:Formula?formula=12S', '_ep:Formula?formula=1L'])
-        @formula[-1].nearest_past(date, n, d)
-      end
-
-      # 二十四節気が指定の周期番号パターンになる最も近い過去の日時
-      #
-      # @param [Numeric] date ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] date
-      # @param [Numeric] n 相対周期番号(n=0 なら date または date の直前が基準)
-      # @param [Numeric] d 単位周期数
-      #
-      # @return [Numeric, When::TM::TemporalPosition] 周期番号が d で割って n になる日時
-      #
-      def nearest_past_from_solar_term(date, n=0, d=1)
-        @formula ||= When.Resource(['_ep:Formula?formula=12S', '_ep:Formula?formula=1L'])
-        @formula[0].nearest_past(date, n, d)
-      end
+    module ForwardedFormula
 
       private
 
@@ -1262,18 +1234,17 @@ module When::Ephemeris
       root(time0, cn) {|t| time_to_cn(t)}
     end
 
-    # 当該日付のthiti の変化範囲
+    # 当該日付の月の位相の変化範囲
     #
     # @param [When::TM::TemporalPosition] date 日付
     #
-    # @return [Range] 当該日付のthiti の変化範囲(朔を含む場合 nil)
+    # @return [Array<Numeric>] 当該日付の月の位相の変化範囲
     #
-    def thiti_range(date)
-      date   = date.floor
-      p0, p1 = [date, date.succ].map {|d|
-        (30.0 * time_to_cn(d)) % 30.0
+    def phase_range(date)
+      date = date.floor
+      [date, date.succ].map {|d|
+        time_to_cn(d)
       }
-      p0 >= p1 ? nil : p0...p1
     end
 
     # 指定の周期番号パターンになる最も近い過去の日時
@@ -1774,6 +1745,52 @@ module When::Ephemeris
   #
   class MeanLunation < Formula
 
+    #
+    # Lunar Calendar Formula
+    #
+    module LunarMethod
+
+      private
+
+      # 周期番号 -> 日時
+      #
+      # @param [Numeric] cn 周期番号
+      # @param [Numeric] time0 日時の初期近似値
+      #
+      # @return [Numeric] ユリウス日
+      #
+      # @note 半ティティの日時の丸め誤差に配慮
+      #
+      def cn_to_time_(cn, time0=nil)
+        time0 ||= (cn - @cycle_number_0m) / @cycle_number_1m
+        return time0 if (cn * 60 - (cn * 60).round).abs > @cycle_precision
+        ((time0 + 1.0/256 - @day_epoch) / @half_tithi).floor  * @half_tithi  + @day_epoch
+      end
+    end
+
+    #
+    # Solar Calendar Formula for Fixed Year Length Method
+    #
+    module SolarMethod
+
+      private
+
+      # 周期番号 -> 日時
+      #
+      # @param [Numeric] cn 周期番号
+      # @param [Numeric] time0 日時の初期近似値
+      #
+      # @return [Numeric] ユリウス日
+      #
+      # @note 太陽黄経が整数になる日時の丸め誤差に配慮
+      #
+      def cn_to_time_(cn, time0=nil)
+        time0 ||= (cn - @cycle_number_0m) / @cycle_number_1m
+        return time0 if (cn * 360 - (cn * 360).round).abs > @cycle_precision
+        ((time0 + 1.0/256 - @day_epoch) / @solar_degree).floor * @solar_degree + @day_epoch
+      end
+    end
+
     # 計算の基準経度 / 度
     # @return [Numeric]
     attr_reader :long
@@ -1851,88 +1868,44 @@ module When::Ephemeris
 
     private
 
-    # 周期番号 -> 日時
-    #
-    # @param [Numeric] cn 周期番号
-    # @param [Numeric] time0 日時の初期近似値
-    #
-    # @return [Numeric] ユリウス日
-    #
-    def cn_to_time_(cn, time0=nil)
-      time0 ||= (cn - @cycle_number_0m) / @cycle_number_1m
-      case @formula
-      when /S/ ; ((time0 + 1.0/256 - @day_epoch) / @solar_terms).floor * @solar_terms + @day_epoch
-      when /L/ ; ((time0 + 1.0/256 - @day_epoch) / @month_tithi).floor * @month_tithi + @day_epoch
-      end
-    end
-
     # オブジェクトの正規化
     def _normalize(args=[], options={})
       Rational
       @time_standard   ||= 'universal'
       @epoch_shift     ||= 1721139        # 西暦 0 年  春分
       @day_shift       ||= Rational(-1,2) # 夜半 -1/2, 日出 -1/4
-      @day_shift          = @day_shift.to_r
-      @longitude_shift ||= Rational(-1,4) # 冬至 -1/4, 立春 -1/8
-      @longitude_shift   = @longitude_shift.to_r
-      @day_epoch         = @day_epoch.to_i + @day_shift
-      @year_length       = @year_length.to_r
-      @lunation_length   = @lunation_length.to_r
-      @month_length      = 1 / (1.to_r/@year_length + 1.to_r/@lunation_length)
-      @denominator       = [@year_length.denominator, @lunation_length.denominator].max
-      @solar_terms       = @year_length / 24
-      @month_tithi       = @lunation_length / 30
-      @year_epoch        = 0
-      @year_epoch        = @longitude_shift -_mean_sun_(@epoch_shift).to_i
-      @month_epoch       = 0
-      @month_epoch       = @longitude_shift -_mean_moon_(@epoch_shift).to_i
-      super
-    end
-  end
-
-  #
-  # Solar Calendar Formula for Variable Year Length Method
-  #
-  class VariableYearLengthMethod < Formula
-
-    # 日時 -> 周期番号
-    #
-    # @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)
-      t, n = (cn / 12.0 - @longitude_shift - @year_epoch).divmod(1)
-      @day_epoch + @year_length * t - @year_delta * t * (t-1) + (@year_length - 2 * @year_delta * t) * n
-    end
-
-    private
-
-    # オブジェクトの正規化
-    def _normalize(args=[], options={})
-      Rational
-      @time_standard   ||= 'universal'
-      @day_shift       ||= Rational(-1,2) # 夜半 -1/2, 日出 -1/4
       @day_shift         = @day_shift.to_r
       @longitude_shift ||= Rational(-1,4) # 冬至 -1/4, 立春 -1/8
       @longitude_shift   = @longitude_shift.to_r
-      @day_epoch         = @day_epoch.to_f + @day_shift
-      @year_epoch        = @year_epoch.to_f
-      @year_length       = @year_length.to_f
-      @year_delta        = @year_delta.to_f * 1.0E-6
+      @day_epoch         = (@day_epoch.to_f == @day_epoch.to_i ? @day_epoch.to_i : @day_epoch.to_f) + @day_shift
+      @year_length       = @year_length.to_r
+      @year_delta        = @year_delta.to_f * 1.0E-6 if @year_delta
+      if @year_epoch
+        @year_epoch      = @year_epoch.to_f
+      else
+        @year_epoch      = 0
+        @year_epoch      = @longitude_shift -_mean_sun_(@epoch_shift).to_i
+      end
+      @cycle_precision ||= 1.0E-8
+      @cycle_precision   = @cycle_precision.to_f
+
+      if @lunation_length && /S/i !~ @formula
+        # 月の位相の計算
+        @lunation_length = @lunation_length.to_r
+        @month_length    = 1 / (1.to_r/@year_length + 1.to_r/@lunation_length)
+        @half_tithi      = @lunation_length / 60
+        if @month_epoch
+          @month_epoch   = @month_epoch.to_f
+        else
+          @month_epoch   = 0
+          @month_epoch   = @longitude_shift -_mean_moon_(@epoch_shift).to_i
+        end
+        class << self; include LunarMethod; end
+      else
+        # 太陽黄経の計算
+        @solar_degree    = @year_length / 360
+        class << self; include SolarMethod;    end
+      end
       super
     end
   end