RubyGems - when_exe - Versions diffs - 0.4.2 → 0.5.0 - Mend

when_exe 0.4.2 → 0.5.0

Files changed (144) hide show

checksums.yaml +5 -5
data/LICENSE.ja.txt +8 -21
data/LICENSE.txt +5 -27
data/README.md +25 -79
data/bin/.pryrc +9 -0
data/bin/make_ttl.rb +3 -2
data/bin/make_ttl.rb.config +3 -3
data/lib/when_exe.rb +101 -5
data/lib/when_exe/basictypes.rb +47 -16
data/lib/when_exe/calendarnote.rb +27 -18
data/lib/when_exe/calendartypes.rb +12 -50
data/lib/when_exe/coordinates.rb +37 -391
data/lib/when_exe/ephemeris.rb +53 -14
data/lib/when_exe/ephemeris/eclipse.rb +8 -8
data/lib/when_exe/ephemeris/moon.rb +333 -333
data/lib/when_exe/ephemeris/notes.rb +29 -17
data/lib/when_exe/ephemeris/phase_table.rb +825 -0
data/lib/when_exe/ephemeris/term_table.rb +603 -0
data/lib/when_exe/events.rb +1888 -0
data/lib/when_exe/google_api.rb +26 -31
data/lib/when_exe/icalendar.rb +129 -12
data/lib/when_exe/inspect.rb +5 -5
data/lib/when_exe/linkeddata.rb +43 -23
data/lib/when_exe/locales/akt.rb +63 -0
data/lib/when_exe/locales/encoding_conversion.rb +134 -134
data/lib/when_exe/locales/iast.rb +90 -90
data/lib/when_exe/locales/locale.rb +88 -18
data/lib/when_exe/locales/transliteration_table.rb +62 -62
data/lib/when_exe/mini_application.rb +14 -8
data/lib/when_exe/namespace.rb +42 -0
data/lib/when_exe/parts/enumerator.rb +16 -5
data/lib/when_exe/parts/geometric_complex.rb +52 -2
data/lib/when_exe/parts/method_cash.rb +229 -224
data/lib/when_exe/parts/resource.rb +65 -35
data/lib/when_exe/parts/timezone.rb +25 -34
data/lib/when_exe/region/balinese.rb +5 -3
data/lib/when_exe/region/chinese.rb +39 -11
data/lib/when_exe/region/chinese/calendars.rb +1016 -1016
data/lib/when_exe/region/chinese/epochs.rb +7 -7
data/lib/when_exe/region/chinese/notes.rb +3 -3
data/lib/when_exe/region/chinese/twins.rb +134 -11
data/lib/when_exe/region/dee.rb +4 -4
data/lib/when_exe/region/discordian.rb +1 -1
data/lib/when_exe/region/geologicalage.rb +139 -128
data/lib/when_exe/region/hanke_henry.rb +4 -4
data/lib/when_exe/region/indian.rb +77 -8
data/lib/when_exe/region/international_fixed.rb +3 -3
data/lib/when_exe/region/islamic.rb +1 -1
data/lib/when_exe/region/japanese.rb +12 -8
data/lib/when_exe/region/japanese/calendars.rb +397 -397
data/lib/when_exe/region/japanese/eclipses.rb +1194 -1194
data/lib/when_exe/region/japanese/epochs.rb +8 -4
data/lib/when_exe/region/japanese/location.rb +93 -0
data/lib/when_exe/region/japanese/nihon_shoki.rb +70 -70
data/lib/when_exe/region/japanese/notes.rb +57 -18
data/lib/when_exe/region/japanese/residues.rb +33 -11
data/lib/when_exe/region/japanese/twins.rb +27 -14
data/lib/when_exe/region/japanese/weeks.rb +7 -7
data/lib/when_exe/region/korean.rb +384 -384
data/lib/when_exe/region/location.rb +40 -0
data/lib/when_exe/region/martian.rb +1 -1
data/lib/when_exe/region/mayan.rb +2 -2
data/lib/when_exe/region/roman.rb +2 -2
data/lib/when_exe/region/ryukyu.rb +101 -101
data/lib/when_exe/region/saudi_arabian.rb +57 -0
data/lib/when_exe/region/shire.rb +1 -1
data/lib/when_exe/region/world.rb +3 -3
data/lib/when_exe/region/zoroastrian.rb +1 -1
data/lib/when_exe/spatial.rb +611 -0
data/lib/when_exe/timestandard.rb +9 -7
data/lib/when_exe/tmobjects.rb +34 -2
data/lib/when_exe/tmposition.rb +225 -1322
data/lib/when_exe/tmptypes.rb +1279 -0
data/lib/when_exe/tmreference.rb +42 -7
data/lib/when_exe/version.rb +3 -3
data/test/events/example-datasets +6 -0
data/test/events/history-dataset.csv +22 -0
data/test/events/japanese-holiday-index.csv +28 -0
data/test/events/japanese-holiday.csv +94 -0
data/test/events/japanese-holiday.ttl +948 -0
data/test/events/make_events_ttl.rb +18 -0
data/test/events/mori_wikichoshi.csv +14 -0
data/test/events/ndl_koyomi.csv +220 -0
data/test/events/ndl_koyomi_index.csv +44 -0
data/test/events/primeminister-dataset.csv +19 -0
data/test/events/shogun-dataset.csv +22 -0
data/test/events/test-history-dataset-edge-sparql.csv +26 -0
data/test/events/test-history-dataset-edge.csv +27 -0
data/test/events/test-history-dataset-sparql.csv +22 -0
data/test/events/test-history-dataset.csv +23 -0
data/test/events/test-history-events-edge.ttl +89 -0
data/test/events/test-history-events.csv +6 -0
data/test/examples/JapanHolidays.ics +2 -2
data/test/examples/JapanHolidaysRFC6350.ics +3 -3
data/test/examples/Residue.m17n +2 -2
data/test/examples/Spatial.m17n +3 -3
data/test/examples/Terms.m17n +3 -3
data/test/scripts/8.ext.rb +6 -6
data/test/scripts/8.rb +6 -6
data/test/scripts/geometric_complex.rb +41 -41
data/test/scripts/korea.rb +59 -59
data/test/scripts/thai.rb +36 -36
data/test/test.rb +6 -0
data/test/test/basictypes.rb +431 -431
data/test/test/calendarnote.rb +86 -86
data/test/test/calendartypes.rb +97 -97
data/test/test/coordinates.rb +399 -397
data/test/test/ephemeris.rb +115 -115
data/test/test/ephemeris/moon.rb +14 -14
data/test/test/ephemeris/planets.rb +14 -14
data/test/test/ephemeris/sun.rb +14 -14
data/test/test/events.rb +32 -0
data/test/test/google_api.rb +45 -32
data/test/test/inspect.rb +153 -153
data/test/test/parts.rb +488 -488
data/test/test/region/armenian.rb +20 -20
data/test/test/region/bahai.rb +58 -58
data/test/test/region/balinese.rb +34 -34
data/test/test/region/chinese.rb +229 -229
data/test/test/region/christian.rb +226 -226
data/test/test/region/coptic.rb +27 -27
data/test/test/region/discordian.rb +20 -20
data/test/test/region/french.rb +33 -33
data/test/test/region/geologicalage.rb +17 -17
data/test/test/region/indian.rb +14 -1
data/test/test/region/iran.rb +54 -54
data/test/test/region/islamic.rb +64 -54
data/test/test/region/japanese.rb +281 -261
data/test/test/region/jewish.rb +63 -63
data/test/test/region/m17n.rb +21 -21
data/test/test/region/mayan.rb +17 -16
data/test/test/region/reforms.rb +2 -2
data/test/test/region/residue.rb +7 -7
data/test/test/region/shire.rb +58 -58
data/test/test/region/swedish.rb +45 -45
data/test/test/region/thai.rb +116 -116
data/test/test/region/tibetan.rb +30 -30
data/test/test/region/vietnamese.rb +102 -102
data/test/test/region/zoroastrian.rb +58 -58
data/test/test/timestandard.rb +81 -81
data/test/test/tmobjects.rb +402 -402
data/test/test/tmposition.rb +66 -4
data/test/test/tmreference.rb +157 -157
metadata +36 -93

data/lib/when_exe/ephemeris.rb CHANGED Viewed

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 =begin
-  Copyright (C) 2011-2015 Takashi SUGA
+  Copyright (C) 2011-2021 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
@@ -26,6 +26,9 @@ module When::Ephemeris
   autoload :V50,     'when_exe/ephemeris/v50'
   autoload :Hindu,   'when_exe/region/indian'
+  autoload :SolarFormulaWithTable, 'when_exe/ephemeris/term_table'
+  autoload :LunarFormulaWithTable, 'when_exe/ephemeris/phase_table'
   include Math
   DAY       =   86400.0            # 日 / 秒
@@ -642,7 +645,7 @@ module When::Ephemeris
     # @return [When::Ephemeris::Coords]
     #
     def rh_to_h(t, loc)
-      rotate_y(loc.lat / (360.0*When::Coordinates::Spatial::DEGREE) - 0.25)
+      rotate_y(loc.lat / (360.0*loc.degree) - 0.25)
     end
     # 赤道座標 -> 地平座標
@@ -655,7 +658,7 @@ module When::Ephemeris
     #
     def r_to_h(t, loc)
       rotate_z(-loc.local_sidereal_time(t) / 24).
-      rotate_y(loc.lat / (360.0*When::Coordinates::Spatial::DEGREE) - 0.25)
+      rotate_y(loc.lat / (360.0*loc.degree) - 0.25)
     end
     # 球面の余弦 spherical law of cosines
@@ -1342,7 +1345,7 @@ module When::Ephemeris
     #
     def day_event(t, event, target=When.Resource('_ep:Sun'), height=nil)
       # 時刻の初期値
-      dl  = @location.long / (360.0 * When::Coordinates::Spatial::DEGREE)
+      dl  = @location.long / (360.0 * @location.degree)
       t0  = (+t + dl).round - dl
       dt  = _coords(t0, When::Coordinates::Spatial::EQUATORIAL, target).phi -
             @location.local_sidereal_time(t0) / 24.0 + 0.25 * (event||0)
@@ -1526,8 +1529,8 @@ module When::Ephemeris
       # difference of ecliptic longitude between zenith and target star
       # when the event happens
       pos  = _coords(+t, When::Coordinates::Spatial::EQUATORIAL, target)
-      long = @location.long / When::Coordinates::Spatial::DEGREE
-      lat  = @location.lat  / When::Coordinates::Spatial::DEGREE
+      long = @location.long / @location.degree
+      lat  = @location.lat  / @location.degree
       dp   = (sind(hs) - sind(lat) * sinc(pos.theta)) /
                         (cosd(lat) * cosc(pos.theta))
@@ -1574,7 +1577,7 @@ module When::Ephemeris
       seed[:precision] = nil
       seed[:clock]   ||= When::TM::Clock.local_time
       t = When::TM::JulianDate._d_to_t(d)
-      seed[:frame].jul_trans(@is_dynamical ? When::TM::JulianDate.dynamical_time(t) :
+      seed[:frame].jul_trans(@is_dynamical ? When::TM::JulianDate.dynamical_time(t, {:time_standard=>seed[:time_standard]}) :
                                              When::TM::JulianDate.universal_time(t), seed)
     end
@@ -1595,8 +1598,8 @@ module When::Ephemeris
     #   @sun_1m          =  t = CYCLE_1M日 までの太陽年番号の比例定数
     def _normalize(args=[], options={})
       @location         = When.Resource(@location, '_l:') if @location.kind_of?(String)
-      @long, @lat, @alt = [@location.long / When::Coordinates::Spatial::DEGREE,
-                           @location.lat  / When::Coordinates::Spatial::DEGREE,
+      @long, @lat, @alt = [@location.long / @location.degree,
+                           @location.lat  / @location.degree,
                            @location.alt] if @location
       @formula        ||= '1L'
       @time_standard  ||= 'dynamical'
@@ -1775,6 +1778,43 @@ module When::Ephemeris
     end
   end
+  #
+  #  一部表引き
+  #
+  class FormulaWithTable < Formula
+    private
+    # 周期番号 -> 日時
+    #
+    # @param [Numeric] cn 周期番号
+    # @param [Numeric] time0 日時の初期近似値
+    #
+    # @return [Numeric] ユリウス日
+    #
+    def cn_to_time_(cn, time0=nil)
+      count = ((cn - self.class::CycleBase) * self.class::CycleUnit / 360.0).floor
+      return super unless count >= 0
+      degree = (cn * self.class::CycleUnit) % 360
+      degree = degree.round if (degree - degree.round).abs < 1.5E-8
+      if degree == 360
+        degree = 0
+        count += 1
+      end
+      epoch = self.class::Epoch[degree]
+      return super unless epoch
+      delta = self.class::Delta[degree][count]
+      return super unless delta
+      epoch + count * self.class::CyclePeriod + delta * 1E-9
+    end
+    # オブジェクトの正規化
+    def _normalize(args=[], options={})
+      @formula = self.class::FormulaType
+      super
+    end
+  end
   #
   # Luni-Solar Calendar Formula for Mean Lunation Type
   #
@@ -1919,7 +1959,7 @@ module When::Ephemeris
         @year_epoch      = @year_epoch.to_f
       else
         @year_epoch      = 0
-        @year_epoch      = @longitude_shift -_mean_sun_(@epoch_shift).to_i
+        @year_epoch      = @longitude_shift - _mean_sun_(@epoch_shift).to_i
       end
       @cycle_precision ||= 1.0E-8
       @cycle_precision   = @cycle_precision.to_f
@@ -1933,7 +1973,7 @@ module When::Ephemeris
           @month_epoch   = @month_epoch.to_f
         else
           @month_epoch   = 0
-          @month_epoch   = @longitude_shift -_mean_moon_(@epoch_shift).to_i
+          @month_epoch   = @longitude_shift - _mean_moon_(@epoch_shift).to_i
         end
         extend LunarMethod
       else
@@ -1943,9 +1983,8 @@ module When::Ephemeris
       end
       # 通法
-      @denominator = @year_length.denominator if @year_length.kind_of?(Rational)
-      @denominator = [@denominator||0, @lunation_length.denominator].max if @lunation_length.kind_of?(Rational)
+      @denominator   = @denominator.to_i if @denominator
+      @denominator ||= [@year_length, @lunation_length].map {|length| length.kind_of?(Rational) ? length.denominator : 0}.max
       super
     end
   end

data/lib/when_exe/ephemeris/eclipse.rb CHANGED Viewed

@@ -1,6 +1,6 @@
 # -*- coding: utf-8 -*-
 =begin
-    Copyright (C) 2014-2015 Takashi SUGA
+    Copyright (C) 2014-2016 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
@@ -16,14 +16,14 @@ class When::Coordinates::Spatial
   # 日食の情報
   #
   # @param [When::TM::TemporalPosition] date
-  # @param [Range<When::TM::TemporalPosition>] date
+  # @param [::Range<When::TM::TemporalPosition>] date
   # @param [Block] block
   #
-  # @return [Array<String, Numeric, Array<Array<Numeric or When::TM::TemporalPosition, String>>>] 食の情報(のArray(dateがRangeの場合))
+  # @return [Array<String, Numeric, Array<Array<Numeric or When::TM::TemporalPosition, String>>>] 食の情報(のArray(dateが::Rangeの場合))
   #   @see When::Coordinates::Spatial#eclipse_info
   #
   def solar_eclipse(date, &block)
-    if date.kind_of?(Range)
+    if date.kind_of?(::Range)
       last  = date.last.to_i
       last -= 1 if date.exclude_end?
       first = date.first.to_i
@@ -51,16 +51,16 @@ class When::Coordinates::Spatial
   # 月食の情報
   #
   # @param [When::TM::TemporalPosition] date
-  # @param [Range<When::TM::TemporalPosition>] date
-  #   @note Rangeの場合夜半より翌日に向けmargin経過時点より前は前日扱い
+  # @param [::Range<When::TM::TemporalPosition>] date
+  #   @note ::Rangeの場合夜半より翌日に向けmargin経過時点より前は前日扱い
   # @param [When::TM::Duration] margin
   # @param [Block] block
   #
-  # @return [Array<String, Numeric, Array<Array<Numeric or When::TM::TemporalPosition, String>>>] 食の情報(のArray(dateがRangeの場合))
+  # @return [Array<String, Numeric, Array<Array<Numeric or When::TM::TemporalPosition, String>>>] 食の情報(のArray(dateが::Rangeの場合))
   #   @see When::Coordinates::Spatial#eclipse_info
   #
   def lunar_eclipse(date, margin=When::PT6H, &block)
-    if date.kind_of?(Range)
+    if date.kind_of?(::Range)
       last  = date.last.to_i
       last -= 1 if date.exclude_end?
       first = date.first.to_i

data/lib/when_exe/ephemeris/moon.rb CHANGED Viewed

@@ -1,333 +1,333 @@
-# -*- coding: utf-8 -*-
-=begin
-  Copyright (C) 2011-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
-require 'when_exe/ephemeris/sun'
-module When
-  module Ephemeris
-    #------------------------------------------
-    #     月と地球の本影 計算用データ
-    #------------------------------------------
-    # 月の地心黄経
-    # @private
-    P0dL =
-      [[SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
-       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
-       [SIN   ,   66.0       ,       0.2      ,  0.0006   ],
-       [SIN   ,  249.0       ,     -19.3      ,  0.0006   ]]
-    # @private
-    P0L =
-      [[LIN   ,  124.8757417 ,  4812.6788201  , -1.330000E-7],
-       [SINL  ,  338.9158263 ,  4771.9886313  ,  6.2887   ,+8.99400E-7], # 中心差
-       [SIN   ,  107.236832  ,   -4133.3536154,  1.2740   ,+1.22612E-6], # 出差
-       [SIN   ,   51.678963  ,    8905.3422467,  0.6583   ,-0.32672E-6], # 二均差
-       [SIN   ,  317.831653  ,    9543.9772627,  0.2136   ,+1.7988E-6 ], # 中心差 x 2
-       [SIN   ,  176.528826  ,     359.990504 ,  0.1856   ,-1.561E-8  ], # 年差
-       [SIN   ,  292.463     ,    9664.0403505,  0.1143   ,-0.6812E-6 ],
-       [SIN   ,   86.16      ,     638.63475182, 0.0588   ,+2.1256E-6 ],
-       [SIN   ,  103.78      ,   -3773.36305278, 0.0572   ,-1.22769E-6],
-       [SIN   ,   30.58      ,   13677.331    ,  0.0533   ],
-       [SIN   ,  124.86      ,   -8545.352    ,  0.0459   ],
-       [SIN   ,  342.38      ,    4411.998    ,  0.0410   ],
-       [SIN   ,   25.83      ,    4452.671    ,  0.0348   ], # 月角差( 二均差 / 2)
-       [SIN   ,  155.45      ,    5131.979    ,  0.0305   ],
-       [SIN   ,  240.79      ,     758.698    ,  0.0153   ],
-       [SIN   ,  271.38      ,   14436.029    ,  0.0125   ],
-       [SIN   ,  226.45      ,   -4892.052    ,  0.0110   ],
-       [SIN   ,   55.58      ,  -13038.696    ,  0.0107   ],
-       [SIN   ,  296.75      ,   14315.966    ,  0.0100   ],
-       [SIN   ,   34.5       ,   -8266.71     ,  0.0085   ],
-       [SIN   ,  290.7       ,   -4493.34     ,  0.0079   ],
-       [SIN   ,  228.2       ,    9265.33     ,  0.0068   ],
-       [SIN   ,  133.1       ,     319.32     ,  0.0052   ],
-       [SIN   ,  202.4       ,    4812.66     ,  0.0050   ],
-       [SIN   ,   68.6       ,     -19.34     ,  0.0048   ],
-       [SIN   ,   34.1       ,   13317.34     ,  0.0040   ],
-       [SIN   ,    9.5       ,   18449.32     ,  0.0040   ],
-       [SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
-       [SIN   ,  103.3       ,   17810.68     ,  0.0039   ],
-       [SIN   ,   65.1       ,    5410.62     ,  0.0037   ],
-       [SIN   ,  321.3       ,    9183.99     ,  0.0027   ],
-       [SIN   ,  174.8       ,  -13797.39     ,  0.0026   ],
-       [SIN   ,   82.7       ,     998.63     ,  0.0024   ],
-       [SIN   ,    4.7       ,    9224.66     ,  0.0024   ],
-       [SIN   ,  121.4       ,   -8185.36     ,  0.0022   ],
-       [SIN   ,  134.4       ,    9903.97     ,  0.0021   ],
-       [SIN   ,  173.1       ,     719.98     ,  0.0021   ],
-       [SIN   ,  100.3       ,   -3413.37     ,  0.0021   ],
-       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
-       [SIN   ,   98.1       ,    4013.29     ,  0.0018   ],
-       [SIN   ,  344.1       ,   18569.38     ,  0.0016   ],
-       [SIN   ,   52.1       ,  -12678.71     ,  0.0012   ],
-       [SIN   ,  250.3       ,   19208.02     ,  0.0011   ],
-       [SIN   ,   81.0       ,   -8586.0      ,  0.0009   ],
-       [SIN   ,  207.0       ,   14037.3      ,  0.0008   ],
-       [SIN   ,   31.0       ,   -7906.7      ,  0.0008   ],
-       [SIN   ,  346.0       ,    4052.0      ,  0.0007   ],
-       [SIN   ,  294.0       ,   -4853.3      ,  0.0007   ],
-       [SIN   ,   90.0       ,     278.6      ,  0.0007   ],
-       [SIN   ,  237.0       ,    1118.7      ,  0.0006   ],
-       [SIN   ,   82.0       ,   22582.7      ,  0.0005   ],
-       [SIN   ,  276.0       ,   19088.0      ,  0.0005   ],
-       [SIN   ,   73.0       ,  -17450.7      ,  0.0005   ],
-       [SIN   ,  112.0       ,    5091.3      ,  0.0005   ],
-       [SIN   ,  116.0       ,    -398.7      ,  0.0004   ],
-       [SIN   ,   25.0       ,    -120.1      ,  0.0004   ],
-       [SIN   ,  181.0       ,    9584.7      ,  0.0004   ],
-       [SIN   ,   18.0       ,     720.0      ,  0.0004   ],
-       [SIN   ,   60.0       ,   -3814.0      ,  0.0003   ],
-       [SIN   ,   13.0       ,   -3494.7      ,  0.0003   ],
-       [SIN   ,   13.0       ,   18089.3      ,  0.0003   ],
-       [SIN   ,  152.0       ,    5492.0      ,  0.0003   ],
-       [SIN   ,  317.0       ,     -40.7      ,  0.0003   ],
-       [SIN   ,  348.0       ,   23221.3      ,  0.0003   ]]
-    # 月の地心黄緯
-    # @private
-    P0dB =
-      [[SIN   ,   68.64      ,     -19.341    ,  0.0267   ],
-       [SIN   ,  342.0       ,     -19.36     ,  0.0043   ],
-       [SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
-       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
-       [SIN   ,  358.0       ,     -19.4      ,  0.0005   ]]
-    # @private
-    P0B =
-      [[SINL  ,  236.231     ,    4832.0201248,  5.1281218,-0.3406E-6 ],
-       [SIN   ,  215.147     ,    9604.0088   ,  0.2806   ],
-       [SIN   ,   77.316     ,      60.0316   ,  0.2777   ],
-       [SIN   ,    4.563     ,   -4073.3220   ,  0.1732   ],
-       [SIN   ,  308.98      ,    8965.374    ,  0.0554   ],
-       [SIN   ,  343.48      ,     698.667    ,  0.0463   ],
-       [SIN   ,  287.90      ,   13737.362    ,  0.0326   ],
-       [SIN   ,  194.06      ,   14375.997    ,  0.0172   ],
-       [SIN   ,   25.6       ,   -8845.31     ,  0.0093   ],
-       [SIN   ,   98.4       ,   -4711.96     ,  0.0088   ],
-       [SIN   ,    1.1       ,   -3713.33     ,  0.0082   ],
-       [SIN   ,  322.4       ,    5470.66     ,  0.0043   ],
-       [SIN   ,  266.8       ,   18509.35     ,  0.0042   ],
-       [SIN   ,  188.0       ,   -4433.31     ,  0.0034   ],
-       [SIN   ,  312.5       ,    8605.38     ,  0.0025   ],
-       [SIN   ,  291.4       ,   13377.37     ,  0.0022   ],
-       [SIN   ,  340.0       ,    1058.66     ,  0.0021   ],
-       [SIN   ,  218.6       ,    9244.02     ,  0.0019   ],
-       [SIN   ,  291.8       ,   -8206.68     ,  0.0018   ],
-       [SIN   ,   52.8       ,    5192.01     ,  0.0018   ],
-       [SIN   ,  168.7       ,   14496.06     ,  0.0017   ],
-       [SIN   ,   73.8       ,     420.02     ,  0.0016   ],
-       [SIN   ,  262.1       ,    9284.69     ,  0.0015   ],
-       [SIN   ,   31.7       ,    9964.00     ,  0.0015   ],
-       [SIN   ,  260.8       ,    -299.96     ,  0.0014   ],
-       [SIN   ,  239.7       ,    4472.03     ,  0.0013   ],
-       [SIN   ,   30.4       ,     379.35     ,  0.0013   ],
-       [SIN   ,  304.9       ,    4812.68     ,  0.0012   ],
-       [SIN   ,   12.4       ,   -4851.36     ,  0.0012   ],
-       [SIN   ,  173.0       ,   19147.99     ,  0.0011   ],
-       [SIN   ,  312.9       ,  -12978.66     ,  0.0010   ],
-       [SIN   ,    1.0       ,   17870.7      ,  0.0008   ],
-       [SIN   ,  190.0       ,    9724.1      ,  0.0008   ],
-       [SIN   ,   22.0       ,   13098.7      ,  0.0007   ],
-       [SIN   ,  117.0       ,    5590.7      ,  0.0006   ],
-       [SIN   ,   47.0       ,  -13617.3      ,  0.0006   ],
-       [SIN   ,   22.0       ,   -8485.3      ,  0.0005   ],
-       [SIN   ,  150.0       ,    4193.4      ,  0.0005   ],
-       [SIN   ,  119.0       ,   -9483.9      ,  0.0004   ],
-       [SIN   ,  246.0       ,   23281.3      ,  0.0004   ],
-       [SIN   ,  301.0       ,   10242.6      ,  0.0004   ],
-       [SIN   ,  126.0       ,    9325.4      ,  0.0004   ],
-       [SIN   ,  104.0       ,   14097.4      ,  0.0004   ],
-       [SIN   ,  340.0       ,   22642.7      ,  0.0003   ],
-       [SIN   ,  270.0       ,   18149.4      ,  0.0003   ],
-       [SIN   ,  358.0       ,   -3353.3      ,  0.0003   ],
-       [SIN   ,  148.0       ,   19268.0      ,  0.0003   ]]
-    # 月の視差
-    # @private
-    P0P =
-      [[COS   ,  338.92      ,    4771.98849108, 0.0518   ,+9.19178E-7],
-       [COS   ,  287.2       ,   -4133.35355678, 0.0095   ,+1.2262E-6 ],
-       [COS   ,   51.7       ,    8905.34223034, 0.0078   ,-0.3268E-6 ],
-       [COS   ,  317.8       ,    9543.97698216, 0.0028   ,+1.7988E-6 ],
-       [COS   ,   31.0       ,   13677.3      ,  0.0009   ],
-       [COS   ,  305.0       ,   -8545.4      ,  0.0005   ],
-       [COS   ,  284.0       ,   -3773.4      ,  0.0004   ],
-       [COS   ,  342.0       ,    4412.0      ,  0.0003   ]]
-    #
-    # The Moon
-    #
-    class Moon < Datum
-      class << self
-        include Ephemeris
-        # 月の真黄経 / CIRCLE
-        #
-        # @param [Numeric] t ユリウス日(Terrestrial Time)
-        # @param [When::TM::TemporalPosition] t
-        #
-        # @return [Numeric]
-        #
-        def true_longitude(t)
-          y = julian_year_from_1975(+t)
-          return trigonometric(y, P0L, sind(trigonometric(y, P0dL))) / 360.0 + 38770.0
-        end
-        # 月の真黄緯/ CIRCLE
-        #
-        # @param [Numeric] t ユリウス日(Terrestrial Time)
-        # @param [When::TM::TemporalPosition] t
-        #
-        # @return [Numeric]
-        #
-        def latitude(t)
-          y = julian_year_from_1975(+t)
-          return trigonometric(y, P0B, trigonometric(y, P0dB)) / 360.0
-        end
-        # 月の距離 / km
-        #
-        # @param [Numeric] t ユリウス日(Terrestrial Time)
-        # @param [When::TM::TemporalPosition] t
-        #
-        # @return [Numeric]
-        #
-        def pi(t)
-          return 6378.14 / ((trigonometric(julian_year_from_1975(+t), P0P) + 0.9507)*DEG)
-        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), P0L, 0.0, 1) / 360.0 + 38770.0
-        end
-      end
-      # 月の位置 (黄道座標)
-      #
-      # @param [Numeric] t ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] t
-      #
-      # @return [When::Ephemeris::Coords]
-      #
-      def _coords(t)
-        t  = +t
-        Coords.polar(Moon.true_longitude(t),
-                     Moon.latitude(t),
-                     Moon.pi(t)/AU,
-                     Moon.mean_longitude(t)) +
-        When.Resource('_ep:Earth')._coords(t)
-      end
-      # 真黄経 / CIRCLE
-      #
-      # @param [Numeric] t ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] t
-      #
-      # @return [Numeric]
-      #
-      def true_longitude(t)
-        Moon.true_longitude(t)
-      end
-      # 平均黄経 / CIRCLE
-      #
-      # @param [Numeric] t ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] t
-      #
-      # @return [Numeric]
-      #
-      def mean_longitude(t)
-        Moon.mean_longitude(t)
-      end
-      # 平均運動 / (DEG/YEAR)
-      #
-      # @return [Numeric]
-      #
-      def mean_motion
-        P0L[0][2]
-      end
-      #
-      # オブジェクトの生成
-      # @private
-      def initialize(*args)
-        options = [args.pop] if args[-1].kind_of?(Hash)
-        surface_radius, aberration, luminosity, *rest = args
-        surface_radius ||= 1738.1
-        aberration     ||=    0.00020
-        luminosity     ||=    0.40
-        args  = [surface_radius, aberration, luminosity] + rest
-        args += options if options
-        super(*args)
-      end
-    end
-    #
-    # The Shadow of the Earth
-    #
-    class Shadow < Earth
-      # 地球の影の位置 (黄道座標)
-      #
-      # @param [Numeric] t ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] t
-      #
-      # @return [When::Ephemeris::Coords]
-      #
-      def _coords(t)
-        t  = +t
-        radius = Sun.pi(t) + Moon.pi(t)/AU
-        dl = 0.5 + @aberration / radius / 360 # 略算式に含まれる光行差をキャンセルする
-        Coords.polar(Sun.true_longitude(t)+dl, 0, radius, Sun.mean_longitude(t)+dl)
-      end
-      # 地球の影の視半径 / CIRCLE
-      #
-      # @param [Numeric] t ユリウス日(Terrestrial Time)
-      # @param [When::TM::TemporalPosition] t
-      # @param [When::Coordinates::Spatial] base 観測地
-      #
-      # @return [Numeric]
-      #
-      def apparent_radius(t, base=nil)
-        t  = +t
-        radius = Sun.pi(t) * AU
-        (asin(surface_radius/Moon.pi(t)) -
-         asin(Ephemeris::Sun::Radius/radius) +
-         asin(surface_radius/radius)) * 1.02 / CIRCLE
-      end
-    end
-  end
-end
+# -*- coding: utf-8 -*-
+=begin
+  Copyright (C) 2011-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
+require 'when_exe/ephemeris/sun'
+module When
+  module Ephemeris
+    #------------------------------------------
+    #     月と地球の本影 計算用データ
+    #------------------------------------------
+    # 月の地心黄経
+    # @private
+    P0dL =
+      [[SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
+       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
+       [SIN   ,   66.0       ,       0.2      ,  0.0006   ],
+       [SIN   ,  249.0       ,     -19.3      ,  0.0006   ]]
+    # @private
+    P0L =
+      [[LIN   ,  124.8757417 ,  4812.6788201  , -1.330000E-7],
+       [SINL  ,  338.9158263 ,  4771.9886313  ,  6.2887   ,+8.99400E-7], # 中心差
+       [SIN   ,  107.236832  ,   -4133.3536154,  1.2740   ,+1.22612E-6], # 出差
+       [SIN   ,   51.678963  ,    8905.3422467,  0.6583   ,-0.32672E-6], # 二均差
+       [SIN   ,  317.831653  ,    9543.9772627,  0.2136   ,+1.7988E-6 ], # 中心差 x 2
+       [SIN   ,  176.528826  ,     359.990504 ,  0.1856   ,-1.561E-8  ], # 年差
+       [SIN   ,  292.463     ,    9664.0403505,  0.1143   ,-0.6812E-6 ],
+       [SIN   ,   86.16      ,     638.63475182, 0.0588   ,+2.1256E-6 ],
+       [SIN   ,  103.78      ,   -3773.36305278, 0.0572   ,-1.22769E-6],
+       [SIN   ,   30.58      ,   13677.331    ,  0.0533   ],
+       [SIN   ,  124.86      ,   -8545.352    ,  0.0459   ],
+       [SIN   ,  342.38      ,    4411.998    ,  0.0410   ],
+       [SIN   ,   25.83      ,    4452.671    ,  0.0348   ], # 月角差( 二均差 / 2)
+       [SIN   ,  155.45      ,    5131.979    ,  0.0305   ],
+       [SIN   ,  240.79      ,     758.698    ,  0.0153   ],
+       [SIN   ,  271.38      ,   14436.029    ,  0.0125   ],
+       [SIN   ,  226.45      ,   -4892.052    ,  0.0110   ],
+       [SIN   ,   55.58      ,  -13038.696    ,  0.0107   ],
+       [SIN   ,  296.75      ,   14315.966    ,  0.0100   ],
+       [SIN   ,   34.5       ,   -8266.71     ,  0.0085   ],
+       [SIN   ,  290.7       ,   -4493.34     ,  0.0079   ],
+       [SIN   ,  228.2       ,    9265.33     ,  0.0068   ],
+       [SIN   ,  133.1       ,     319.32     ,  0.0052   ],
+       [SIN   ,  202.4       ,    4812.66     ,  0.0050   ],
+       [SIN   ,   68.6       ,     -19.34     ,  0.0048   ],
+       [SIN   ,   34.1       ,   13317.34     ,  0.0040   ],
+       [SIN   ,    9.5       ,   18449.32     ,  0.0040   ],
+       [SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
+       [SIN   ,  103.3       ,   17810.68     ,  0.0039   ],
+       [SIN   ,   65.1       ,    5410.62     ,  0.0037   ],
+       [SIN   ,  321.3       ,    9183.99     ,  0.0027   ],
+       [SIN   ,  174.8       ,  -13797.39     ,  0.0026   ],
+       [SIN   ,   82.7       ,     998.63     ,  0.0024   ],
+       [SIN   ,    4.7       ,    9224.66     ,  0.0024   ],
+       [SIN   ,  121.4       ,   -8185.36     ,  0.0022   ],
+       [SIN   ,  134.4       ,    9903.97     ,  0.0021   ],
+       [SIN   ,  173.1       ,     719.98     ,  0.0021   ],
+       [SIN   ,  100.3       ,   -3413.37     ,  0.0021   ],
+       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
+       [SIN   ,   98.1       ,    4013.29     ,  0.0018   ],
+       [SIN   ,  344.1       ,   18569.38     ,  0.0016   ],
+       [SIN   ,   52.1       ,  -12678.71     ,  0.0012   ],
+       [SIN   ,  250.3       ,   19208.02     ,  0.0011   ],
+       [SIN   ,   81.0       ,   -8586.0      ,  0.0009   ],
+       [SIN   ,  207.0       ,   14037.3      ,  0.0008   ],
+       [SIN   ,   31.0       ,   -7906.7      ,  0.0008   ],
+       [SIN   ,  346.0       ,    4052.0      ,  0.0007   ],
+       [SIN   ,  294.0       ,   -4853.3      ,  0.0007   ],
+       [SIN   ,   90.0       ,     278.6      ,  0.0007   ],
+       [SIN   ,  237.0       ,    1118.7      ,  0.0006   ],
+       [SIN   ,   82.0       ,   22582.7      ,  0.0005   ],
+       [SIN   ,  276.0       ,   19088.0      ,  0.0005   ],
+       [SIN   ,   73.0       ,  -17450.7      ,  0.0005   ],
+       [SIN   ,  112.0       ,    5091.3      ,  0.0005   ],
+       [SIN   ,  116.0       ,    -398.7      ,  0.0004   ],
+       [SIN   ,   25.0       ,    -120.1      ,  0.0004   ],
+       [SIN   ,  181.0       ,    9584.7      ,  0.0004   ],
+       [SIN   ,   18.0       ,     720.0      ,  0.0004   ],
+       [SIN   ,   60.0       ,   -3814.0      ,  0.0003   ],
+       [SIN   ,   13.0       ,   -3494.7      ,  0.0003   ],
+       [SIN   ,   13.0       ,   18089.3      ,  0.0003   ],
+       [SIN   ,  152.0       ,    5492.0      ,  0.0003   ],
+       [SIN   ,  317.0       ,     -40.7      ,  0.0003   ],
+       [SIN   ,  348.0       ,   23221.3      ,  0.0003   ]]
+    # 月の地心黄緯
+    # @private
+    P0dB =
+      [[SIN   ,   68.64      ,     -19.341    ,  0.0267   ],
+       [SIN   ,  342.0       ,     -19.36     ,  0.0043   ],
+       [SIN   ,   93.8       ,      -1.33     ,  0.0040   ],
+       [SIN   ,  248.6       ,     -19.34     ,  0.0020   ],
+       [SIN   ,  358.0       ,     -19.4      ,  0.0005   ]]
+    # @private
+    P0B =
+      [[SINL  ,  236.231     ,    4832.0201248,  5.1281218,-0.3406E-6 ],
+       [SIN   ,  215.147     ,    9604.0088   ,  0.2806   ],
+       [SIN   ,   77.316     ,      60.0316   ,  0.2777   ],
+       [SIN   ,    4.563     ,   -4073.3220   ,  0.1732   ],
+       [SIN   ,  308.98      ,    8965.374    ,  0.0554   ],
+       [SIN   ,  343.48      ,     698.667    ,  0.0463   ],
+       [SIN   ,  287.90      ,   13737.362    ,  0.0326   ],
+       [SIN   ,  194.06      ,   14375.997    ,  0.0172   ],
+       [SIN   ,   25.6       ,   -8845.31     ,  0.0093   ],
+       [SIN   ,   98.4       ,   -4711.96     ,  0.0088   ],
+       [SIN   ,    1.1       ,   -3713.33     ,  0.0082   ],
+       [SIN   ,  322.4       ,    5470.66     ,  0.0043   ],
+       [SIN   ,  266.8       ,   18509.35     ,  0.0042   ],
+       [SIN   ,  188.0       ,   -4433.31     ,  0.0034   ],
+       [SIN   ,  312.5       ,    8605.38     ,  0.0025   ],
+       [SIN   ,  291.4       ,   13377.37     ,  0.0022   ],
+       [SIN   ,  340.0       ,    1058.66     ,  0.0021   ],
+       [SIN   ,  218.6       ,    9244.02     ,  0.0019   ],
+       [SIN   ,  291.8       ,   -8206.68     ,  0.0018   ],
+       [SIN   ,   52.8       ,    5192.01     ,  0.0018   ],
+       [SIN   ,  168.7       ,   14496.06     ,  0.0017   ],
+       [SIN   ,   73.8       ,     420.02     ,  0.0016   ],
+       [SIN   ,  262.1       ,    9284.69     ,  0.0015   ],
+       [SIN   ,   31.7       ,    9964.00     ,  0.0015   ],
+       [SIN   ,  260.8       ,    -299.96     ,  0.0014   ],
+       [SIN   ,  239.7       ,    4472.03     ,  0.0013   ],
+       [SIN   ,   30.4       ,     379.35     ,  0.0013   ],
+       [SIN   ,  304.9       ,    4812.68     ,  0.0012   ],
+       [SIN   ,   12.4       ,   -4851.36     ,  0.0012   ],
+       [SIN   ,  173.0       ,   19147.99     ,  0.0011   ],
+       [SIN   ,  312.9       ,  -12978.66     ,  0.0010   ],
+       [SIN   ,    1.0       ,   17870.7      ,  0.0008   ],
+       [SIN   ,  190.0       ,    9724.1      ,  0.0008   ],
+       [SIN   ,   22.0       ,   13098.7      ,  0.0007   ],
+       [SIN   ,  117.0       ,    5590.7      ,  0.0006   ],
+       [SIN   ,   47.0       ,  -13617.3      ,  0.0006   ],
+       [SIN   ,   22.0       ,   -8485.3      ,  0.0005   ],
+       [SIN   ,  150.0       ,    4193.4      ,  0.0005   ],
+       [SIN   ,  119.0       ,   -9483.9      ,  0.0004   ],
+       [SIN   ,  246.0       ,   23281.3      ,  0.0004   ],
+       [SIN   ,  301.0       ,   10242.6      ,  0.0004   ],
+       [SIN   ,  126.0       ,    9325.4      ,  0.0004   ],
+       [SIN   ,  104.0       ,   14097.4      ,  0.0004   ],
+       [SIN   ,  340.0       ,   22642.7      ,  0.0003   ],
+       [SIN   ,  270.0       ,   18149.4      ,  0.0003   ],
+       [SIN   ,  358.0       ,   -3353.3      ,  0.0003   ],
+       [SIN   ,  148.0       ,   19268.0      ,  0.0003   ]]
+    # 月の視差
+    # @private
+    P0P =
+      [[COS   ,  338.92      ,    4771.98849108, 0.0518   ,+9.19178E-7],
+       [COS   ,  287.2       ,   -4133.35355678, 0.0095   ,+1.2262E-6 ],
+       [COS   ,   51.7       ,    8905.34223034, 0.0078   ,-0.3268E-6 ],
+       [COS   ,  317.8       ,    9543.97698216, 0.0028   ,+1.7988E-6 ],
+       [COS   ,   31.0       ,   13677.3      ,  0.0009   ],
+       [COS   ,  305.0       ,   -8545.4      ,  0.0005   ],
+       [COS   ,  284.0       ,   -3773.4      ,  0.0004   ],
+       [COS   ,  342.0       ,    4412.0      ,  0.0003   ]]
+    #
+    # The Moon
+    #
+    class Moon < Datum
+      class << self
+        include Ephemeris
+        # 月の真黄経 / CIRCLE
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def true_longitude(t)
+          y = julian_year_from_1975(+t)
+          return trigonometric(y, P0L, sind(trigonometric(y, P0dL))) / 360.0 + 38770.0
+        end
+        # 月の真黄緯/ CIRCLE
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def latitude(t)
+          y = julian_year_from_1975(+t)
+          return trigonometric(y, P0B, trigonometric(y, P0dB)) / 360.0
+        end
+        # 月の距離 / km
+        #
+        # @param [Numeric] t ユリウス日(Terrestrial Time)
+        # @param [When::TM::TemporalPosition] t
+        #
+        # @return [Numeric]
+        #
+        def pi(t)
+          return 6378.14 / ((trigonometric(julian_year_from_1975(+t), P0P) + 0.9507)*DEG)
+        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), P0L, 0.0, 1) / 360.0 + 38770.0
+        end
+      end
+      # 月の位置 (黄道座標)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [When::Ephemeris::Coords]
+      #
+      def _coords(t)
+        t  = +t
+        Coords.polar(Moon.true_longitude(t),
+                     Moon.latitude(t),
+                     Moon.pi(t)/AU,
+                     Moon.mean_longitude(t)) +
+        When.Resource('_ep:Earth')._coords(t)
+      end
+      # 真黄経 / CIRCLE
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [Numeric]
+      #
+      def true_longitude(t)
+        Moon.true_longitude(t)
+      end
+      # 平均黄経 / CIRCLE
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [Numeric]
+      #
+      def mean_longitude(t)
+        Moon.mean_longitude(t)
+      end
+      # 平均運動 / (DEG/YEAR)
+      #
+      # @return [Numeric]
+      #
+      def mean_motion
+        P0L[0][2]
+      end
+      #
+      # オブジェクトの生成
+      # @private
+      def initialize(*args)
+        options = [args.pop] if args[-1].kind_of?(Hash)
+        surface_radius, aberration, luminosity, *rest = args
+        surface_radius ||= 1738.1
+        aberration     ||=    0.00020
+        luminosity     ||=    0.40
+        args  = [surface_radius, aberration, luminosity] + rest
+        args += options if options
+        super(*args)
+      end
+    end
+    #
+    # The Shadow of the Earth
+    #
+    class Shadow < Earth
+      # 地球の影の位置 (黄道座標)
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      #
+      # @return [When::Ephemeris::Coords]
+      #
+      def _coords(t)
+        t  = +t
+        radius = Sun.pi(t) + Moon.pi(t)/AU
+        dl = 0.5 + @aberration / radius / 360 # 略算式に含まれる光行差をキャンセルする
+        Coords.polar(Sun.true_longitude(t)+dl, 0, radius, Sun.mean_longitude(t)+dl)
+      end
+      # 地球の影の視半径 / CIRCLE
+      #
+      # @param [Numeric] t ユリウス日(Terrestrial Time)
+      # @param [When::TM::TemporalPosition] t
+      # @param [When::Coordinates::Spatial] base 観測地
+      #
+      # @return [Numeric]
+      #
+      def apparent_radius(t, base=nil)
+        t  = +t
+        radius = Sun.pi(t) * AU
+        (asin(surface_radius/Moon.pi(t)) -
+         asin(Ephemeris::Sun::Radius/radius) +
+         asin(surface_radius/radius)) * 1.02 / CIRCLE
+      end
+    end
+  end
+end