zakuro 0.1.4 → 0.4.0

data/lib/zakuro/version/context.rb

@@ -11,7 +11,7 @@ module Zakuro
     attr_reader :version_name
     attr_reader :resolver
 
-    VERSION_NAMES = %w[Genka Gihou Taien Senmyou Joukyou Kansei Tenpou Gregorio].freeze
+    VERSION_NAMES = %w[Genka Gihou Daien Senmyou Joukyou Kansei Tenpou Gregorio].freeze
 
     def initialize(version_name:)
       raise ArgumentError.new, 'invalid version' unless VERSION_NAMES.include?(version_name)

data/lib/zakuro/version/daien/const/number.rb

@@ -0,0 +1,55 @@
+# frozen_string_literal: true
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    #
+    # Const 定数
+    #
+    module Const
+      #
+      # Number 数値
+      #
+      module Number
+        #
+        # Cycle 周期
+        #
+        module Cycle
+          # @return [Integer] 通法（1日=3040分）
+          DAY = 3040
+          #
+          # @return [Float] 転日（1近点月 = 27日1685分79秒（1分=80秒））
+          # @note 小数点以下の値によって大幅に結果が変わる。除算結果を設定した
+          #
+          ANOMALISTIC_MONTH = 83_765 + (79.0 / 80)
+          # @return [Integer] 朔望月
+          SYNODIC_MONTH = 89_773
+          # @return [Integer] 一年
+          YEAR = 1_110_343
+        end
+
+        #
+        # Derivation 導出
+        #
+        module Derivation
+          # @return [Integer] 通余:  (YEAR - DAY * 12 * 30)
+          #   1110343 - 1094400
+          REMAINDER_ALL_YEAR = 15_943
+          # @return [Integer] 旬周（60日） DAY * 60
+          SIXTY_DAYS = 182_400
+        end
+
+        #
+        # 累積
+        #
+        module Stack
+          # @return [Integer] 積年（甲子夜半朔旦冬至〜暦の開始前）
+          TOTAL_YEAR = 96_961_740
+          # @return [Integer] 暦の開始年（開元12年）
+          BEGIN_YEAR = 724
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/const/remainder.rb

@@ -0,0 +1,53 @@
+# frozen_string_literal: true
+
+require_relative '../cycle/remainder'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    #
+    # Const 定数
+    #
+    module Const
+      #
+      # Remainder 大余小余
+      #
+      module Remainder
+        #
+        # Solar 太陽
+        #
+        module Solar
+          # @return [Remainder] 気策（24分の1年）
+          # 1年（1_110_343） / 24 / 1日（3040）= 15 余り 664.2916666666642
+          #   * 0.2916666666642 * 24（1分=24秒） = 6.999999999940799 ≒ 7
+          SOLAR_TERM_AVERAGE = Cycle::Remainder.new(day: 15, minute: 664, second: 7)
+          #
+          # @note 揲法 89773 = 29-1613
+          #   * 89773 / 4 = 22443.25 / 3040 = 7 余り 1163.25
+          #   * 0.25 * 24（1分=24秒） = 6
+          #
+          # @return [Cycle::Remainder] 弦（1分=80秒）
+          QUARTER = Cycle::Remainder.new(day: 7, minute: 1163, second: 6)
+        end
+
+        #
+        # Lunar 月
+        #
+        module Lunar
+          # @return [Cycle::LunarRemainder] 転日（1近点月）
+          ANOMALISTIC_MONTH = \
+            Cycle::LunarRemainder.new(day: 27, minute: 1685, second: 79)
+          #
+          # @note 揲法 89773 = 29-1613
+          #   * 89773 / 4 = 22443.25 / 3040 = 7 余り 1163.25
+          #   * 0.25 * 80（1分=80秒） = 20
+          #
+          # @return [Cycle::LunarRemainder] 弦（1分=80秒）
+          # TODO: 秒が20では通らない。少なくとも0.01247は必要になる
+          QUARTER = Cycle::LunarRemainder.new(day: 7, minute: 1163, second: 20.01247)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/cycle/remainder.rb

@@ -0,0 +1,59 @@
+# frozen_string_literal: true
+
+require_relative '../../../calculation/cycle/abstract_remainder'
+
+require_relative '../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Cycle
+      #
+      # Remainder 時刻情報（大余小余）
+      #
+      class Remainder < Calculation::Cycle::AbstractRemainder
+        # 『歴代天文律暦等志彙編　七』中華書房 p.2056
+        # 「象統;二十四」「其秒盈象統,従小餘」
+        # 上記は舊唐志の記述だが、新唐志（p.2218）でも同様であることを確認した
+        # @return [Integer] 分（1分=24秒）
+        MINUTE = 24
+
+        #
+        # 初期化
+        #
+        # @param [Integer] day 大余（"日"に相当）
+        # @param [Integer] minute 小余（"分"に相当）
+        # @param [Integer] second 秒
+        # @param [Integer] total 繰り上げなしの小余
+        #
+        def initialize(day: -1, minute: -1, second: -1, total: -1)
+          super(base_day: Const::Number::Cycle::DAY, base_mitune: MINUTE,
+                day: day, minute: minute, second: second, total: total)
+        end
+      end
+
+      #
+      # LunarRemainder 月の位相計算向け時刻情報（大余小余）
+      #
+      class LunarRemainder < Calculation::Cycle::AbstractRemainder
+        # @return [Integer] 分（1分=80秒）
+        MINUTE = 80
+
+        #
+        # 初期化
+        #
+        # @param [Integer] day 大余（"日"に相当）
+        # @param [Integer] minute 小余（"分"に相当）
+        # @param [Integer] second 秒
+        # @param [Integer] total 繰り上げなしの小余
+        #
+        def initialize(day: -1, minute: -1, second: -1, total: -1)
+          super(base_day: Const::Number::Cycle::DAY, base_mitune: MINUTE,
+                day: day, minute: minute, second: second, total: total)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/cycle/solar_term.rb

@@ -0,0 +1,34 @@
+# frozen_string_literal: true
+
+require_relative '../../../calculation/cycle/abstract_solar_term'
+
+require_relative '../const/remainder'
+
+require_relative './remainder'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Cycle
+      #
+      # SolarTerm 二十四節気
+      #
+      class SolarTerm < Calculation::Cycle::AbstractSolarTerm
+        # @return [Remainder] 気策（24分の1年）
+        SOLAR_TERM_AVERAGE = Const::Remainder::Solar::SOLAR_TERM_AVERAGE
+
+        #
+        # 初期化
+        #
+        # @param [Integer] index 連番
+        # @param [Remainder] remainder 時刻情報（大余小余）
+        #
+        def initialize(index: -1, remainder: Remainder.new)
+          super(index: index, remainder: remainder, average: SOLAR_TERM_AVERAGE)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/daien.rb

@@ -0,0 +1,36 @@
+# frozen_string_literal: true
+
+require_relative '../abstract_version'
+
+# :nodoc:
+module Zakuro
+  #
+  # Daien 大衍暦
+  #
+  module Daien
+    #
+    # Gateway アクセサメソッド群
+    #
+    class Gateway < AbstractVersion
+      # @return [True] リリースあり
+      RELEASE = true
+
+      # @return [String] 暦クラス名
+      VERSION_NAME = 'Daien'
+
+      #
+      # 西暦日から和暦日に変換する
+      #
+      # @param [Date] western_date 西暦日
+      #
+      # @return [Result::Single] 和暦日
+      #
+      def self.to_japan_date(western_date:)
+        date = Western::Calendar.create(date: western_date)
+
+        context = Context.new(version_name: VERSION_NAME)
+        Calculation::Summary::Single.get(context: context, date: date)
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/monthly/lunar_phase.rb

@@ -0,0 +1,106 @@
+# frozen_string_literal: true
+
+require_relative '../../../calculation/monthly/abstract_lunar_phase'
+
+require_relative '../const/remainder'
+
+require_relative '../stella/solar/location'
+require_relative '../stella/solar/value'
+require_relative '../stella/lunar/location'
+require_relative '../stella/lunar/value'
+
+require_relative '../stella/origin/lunar_age'
+require_relative '../stella/origin/average_november'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Monthly
+      #
+      # LunarPhase 月の位相
+      #
+      class LunarPhase < Calculation::Monthly::AbstractLunarPhase
+        # @return [Cycle::Remainder] 弦
+        QUARTER = Const::Remainder::Solar::QUARTER
+
+        #
+        # 初期化
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        def initialize(western_year:)
+          # 天正閏余
+          lunar_age = Origin::LunarAge.get(western_year: western_year)
+
+          super(
+            quater: QUARTER,
+            average_remainder: Origin::AverageNovember.get(western_year: western_year),
+            solar_location: Solar::Location.new(lunar_age: lunar_age),
+            lunar_location: Lunar::Location.new(
+              western_year: western_year,
+              lunar_age: Cycle::LunarRemainder.new(total: 0).add!(lunar_age)
+            )
+          )
+        end
+
+        private
+
+        # :reek:TooManyStatements { max_statements: 6 }
+
+        #
+        # 現在の定朔を取得する
+        #
+        # @return [Remainder] 定朔
+        #
+        def current_remainder
+          # debug("@average_remainder.format: #{@average_remainder.format(form: '%d-%d-%.5f')}")
+
+          sum = correction_value
+          adjusted = @average_remainder.add(
+            Cycle::Remainder.new(day: 0, minute: sum, second: 0)
+          )
+          # TODO: 大衍暦では進朔するか
+          # adjusted.up_on_new_moon!
+
+          debug("result: #{adjusted.format}")
+
+          adjusted
+        end
+
+        #
+        # 太陽運動の補正値を得る
+        #
+        # @return [Integer] 太陽運動の補正値
+        #
+        def correction_solar_value
+          @solar_location.run
+          # debug("@solar_term.remainder: #{@solar_location.remainder.format(form: '%d-%d-%.5f')}")
+          # debug("@solar_term.index: #{@solar_location.index}")
+
+          Solar::Value.get(solar_location: @solar_location)
+        end
+
+        # :reek:TooManyStatements { max_statements: 6 }
+
+        #
+        # 月運動の補正値を得る
+        #
+        # @return [Integer] 月運動の補正値
+        #
+        def correction_moon_value
+          @lunar_location.run
+
+          remainder = @lunar_location.adjusted_remainder
+
+          # debug("[lunar]remainder.format: #{remainder.format(form: '%d-%d-%.5f')}")
+          # debug("[lunar]remainder.day: #{remainder.day}")
+          # debug("[lunar]remainder.minute: #{remainder.minute}")
+
+          Lunar::Value.get(remainder: remainder)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/range/annual_range.rb

@@ -0,0 +1,39 @@
+# frozen_string_literal: true
+
+require_relative '../../../calculation/range/medieval_annual_range'
+require_relative '../monthly/lunar_phase'
+require_relative '../stella/solar/average'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Range
+      #
+      # AnnualRange 年間範囲
+      #
+      module AnnualRange
+        #
+        # 一覧取得する
+        #
+        #   * 対象年に対して、前年11月-当年11月までを出力する
+        #   * 対象年（西暦）と計算年（元号x年）の紐付けは行わない
+        #
+        # @param [Context] context 暦コンテキスト
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Array<Month>] 1年データ
+        #
+        def self.get(context:, western_year:)
+          lunar_phase = Monthly::LunarPhase.new(western_year: western_year)
+          solar_average = Solar::Average.new(western_year: western_year)
+
+          Calculation::Range::MedievalAnnualRange.get(
+            context: context, lunar_phase: lunar_phase, solar_average: solar_average
+          )
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/lunar/adjustment.rb

@@ -0,0 +1,240 @@
+# frozen_string_literal: true
+
+require_relative '../../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Lunar
+      #
+      # Adjustment 補正値情報
+      #
+      module Adjustment
+        # 『歴代天文律暦等志彙編　七』中華書房 p.2230
+        #
+        # 遠/近の地点での中間
+        #
+        # @return [Hash<Integer>] 遠/近の地点での中間
+        DAY_LIMIT = {
+          7 => 2701,
+          14 => 2363,
+          21 => 2024,
+          28 => 1686
+        }.freeze
+
+        #
+        # Row 行情報
+        #
+        class Row
+          # @return [Integer] 入暦（1-14）
+          attr_reader :day
+          # @return [Range] 小余範囲
+          attr_reader :range
+          # @return [Value] 補正値
+          attr_reader :value
+
+          #
+          # 初期化
+          #
+          # @param [Integer] day 入暦（1-14）
+          # @param [Range] range 小余範囲
+          # @param [Value] value 補正値
+          #
+          def initialize(day:, range:, value:)
+            @day = day
+            @range = range
+            @value = value
+          end
+
+          # :reek:ControlParameter
+
+          #
+          # 一致するか
+          #
+          # @param [Integer] day 入暦（1-14）
+          # @param [Integer] minute 小余
+          #
+          # @return [True] 一致
+          # @return [False] 不一致
+          #
+          def match?(day:, minute:)
+            return false unless @day == day
+
+            return false unless @range.include?(minute: minute)
+
+            true
+          end
+
+          #
+          # 分母を返す
+          #
+          # @return [Integer] 分母
+          #
+          def denominator
+            @range.denominator
+          end
+        end
+
+        #
+        # Range 小余範囲
+        #
+        class Range
+          # @return [Integer] 下限
+          MIN = 0
+          # @return [Integer] 上限
+          MAX = Const::Number::Cycle::DAY
+
+          # @return [Integer] 下限
+          attr_reader :min
+          # @return [Integer] 上限
+          attr_reader :max
+
+          #
+          # 初期化
+          #
+          # @param [Integer] min 下限
+          # @param [Integer] max 上限
+          #
+          def initialize(min: MIN, max: MAX)
+            @min = min
+            @max = max
+          end
+
+          #
+          # 含まれるか
+          #
+          # @param [Integer] minute 小余
+          #
+          # @return [True] 含まれる
+          # @return [False] 含まれない
+          #
+          def include?(minute:)
+            minute >= @min && minute <= @max
+          end
+
+          #
+          # 分母を返す
+          #
+          # @return [Integer] 分母
+          #
+          def denominator
+            @max - @min
+          end
+        end
+
+        #
+        # Value 補正値
+        #
+        class Value
+          # @return [Integer] 増減率
+          attr_reader :per
+          # @return [Integer] 遅速積
+          attr_reader :stack
+
+          #
+          # 初期化
+          #
+          # @param [Integer] per 増減率
+          # @param [Integer] stack 遅速積
+          #
+          def initialize(per:, stack:)
+            @per = per
+            @stack = stack
+          end
+
+          #
+          # 文字化
+          #
+          # @return [String] 文字
+          #
+          def to_s
+            "per:#{@per}, stack:#{@stack}"
+          end
+        end
+
+        # rubocop:disable Layout/LineLength
+
+        # 『歴代天文律暦等志彙編　七』中華書房 p.2228-2230
+        #
+        # @note 7日、14日、21日、28日の小余は DAY_LIMIT を参照のこと
+        #
+        # @return [Array<Row>] 月の補正値情報
+        #
+        LIST = [
+          Row.new(day: 1, range: Range.new, value: Value.new(per: +297, stack: 0)),
+          Row.new(day: 2, range: Range.new, value: Value.new(per: +259, stack: +297)),
+          Row.new(day: 3, range: Range.new, value: Value.new(per: +220, stack: +556)),
+          Row.new(day: 4, range: Range.new, value: Value.new(per: +180, stack: +776)),
+          Row.new(day: 5, range: Range.new, value: Value.new(per: +139, stack: +956)),
+          Row.new(day: 6, range: Range.new, value: Value.new(per: +97, stack: +1095)),
+          Row.new(day: 7, range: Range.new(max: DAY_LIMIT[7]), value: Value.new(per: +48, stack: +1192)),
+          Row.new(day: 7, range: Range.new(min: DAY_LIMIT[7]), value: Value.new(per: -6, stack: +1240)), # stack: +1192 + 48
+          Row.new(day: 8, range: Range.new, value: Value.new(per: -64, stack: +1234)),
+          Row.new(day: 9, range: Range.new, value: Value.new(per: -106, stack: +1170)),
+          Row.new(day: 10, range: Range.new, value: Value.new(per: -148, stack: +1064)),
+          Row.new(day: 11, range: Range.new, value: Value.new(per: -189, stack: +916)),
+          Row.new(day: 12, range: Range.new, value: Value.new(per: -229, stack: +727)),
+          Row.new(day: 13, range: Range.new, value: Value.new(per: -267, stack: +498)),
+          Row.new(day: 14, range: Range.new(max: DAY_LIMIT[14]), value: Value.new(per: -231, stack: +231)),
+          Row.new(day: 14, range: Range.new(min: DAY_LIMIT[14]), value: Value.new(per: -66, stack: 0)), # stack: +232 - 231
+          Row.new(day: 15, range: Range.new, value: Value.new(per: -289, stack: -66)),
+          Row.new(day: 16, range: Range.new, value: Value.new(per: -250, stack: -355)),
+          Row.new(day: 17, range: Range.new, value: Value.new(per: -211, stack: -605)),
+          Row.new(day: 18, range: Range.new, value: Value.new(per: -171, stack: -816)),
+          Row.new(day: 19, range: Range.new, value: Value.new(per: -130, stack: -987)),
+          Row.new(day: 20, range: Range.new, value: Value.new(per: -87, stack: -1117)),
+          Row.new(day: 21, range: Range.new(max: DAY_LIMIT[21]), value: Value.new(per: -36, stack: -1204)),
+          Row.new(day: 21, range: Range.new(min: DAY_LIMIT[21]), value: Value.new(per: +18, stack: -1240)), # stack: -1204 - 36
+          Row.new(day: 22, range: Range.new, value: Value.new(per: +73, stack: -1222)),
+          Row.new(day: 23, range: Range.new, value: Value.new(per: +116, stack: -1149)),
+          Row.new(day: 24, range: Range.new, value: Value.new(per: +157, stack: -1033)),
+          Row.new(day: 25, range: Range.new, value: Value.new(per: +198, stack: -876)),
+          Row.new(day: 26, range: Range.new, value: Value.new(per: +237, stack: -678)),
+          Row.new(day: 27, range: Range.new, value: Value.new(per: +276, stack: -441)),
+          Row.new(day: 28, range: Range.new(max: DAY_LIMIT[28]), value: Value.new(per: +165, stack: -165))
+        ].freeze
+        # rubocop:enable Layout/LineLength
+
+        #
+        # 月軌道の補正に必要な基本値を引き当てる
+        #
+        # @param [True, False] forward 進（遠地点より数える）/退（近地点より数える）
+        # @param [Integer] day 大余
+        # @param [Integer] minute 小余
+        #
+        # @return [Row] 補正値
+        #
+        def self.specify(day:, minute:)
+          LIST.each do |row|
+            # NOTE: 範囲が重複している場合、最初に引き当てたほうを優先する
+            return row if row.match?(day: day, minute: minute)
+          end
+
+          raise ArgumentError.new, "invalid parameter: #{day}/#{minute}"
+        end
+
+        # :reek:ControlParameter
+
+        #
+        # 小余の下げ幅を求める
+        #
+        # @param [Integer] day 大余
+        # @param [Integer] minute 小余
+        #
+        # @return [Integer] 小余の下げ幅
+        #
+        def self.minus_minute(day:, minute:)
+          limit = DAY_LIMIT.fetch(day, -1)
+          # 該当なし
+          return minute if limit == -1
+
+          return minute unless minute > limit
+
+          minute - limit
+        end
+      end
+    end
+  end
+end