zakuro 0.2.0 → 0.3.0

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

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

@@ -0,0 +1,44 @@
+# frozen_string_literal: true
+
+require_relative '../../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Lunar
+      #
+      # Localization 入暦特定
+      #
+      module Localization
+        # @return [Integer] 1近点月
+        ANOMALISTIC_MONTH = Const::Number::Cycle::ANOMALISTIC_MONTH
+        # @return [Integer] 積年
+        TOTAL_YEAR = Const::Number::Stack::TOTAL_YEAR
+        # @return [Integer] 開始年
+        BEGIN_YEAR = Const::Number::Stack::BEGIN_YEAR
+        # @return [Integer] 年
+        YEAR = Const::Number::Cycle::YEAR
+
+        #
+        # 対象年の最初の入暦を求める
+        #
+        # @param [Remainder] lunar_age 天正閏余
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Cycle::LunarRemainder] 入暦
+        #
+        def self.first_remainder(lunar_age:, western_year:)
+          # 積年の開始から対象年までの年数
+          total_year = TOTAL_YEAR + western_year - BEGIN_YEAR
+
+          # 通積分 - 天正閏余
+          total_day = total_year * YEAR - lunar_age.to_minute
+
+          Cycle::LunarRemainder.new(total: (total_day % ANOMALISTIC_MONTH))
+        end
+      end
+    end
+  end
+end

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

@@ -0,0 +1,80 @@
+# frozen_string_literal: true
+
+require_relative '../../../../calculation/stella/lunar/abstract_location'
+
+require_relative '../../const/remainder'
+
+require_relative './localization'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Lunar
+      #
+      # Location 入暦
+      #
+      class Location < Calculation::Lunar::AbstractLocation
+        # @return [Cycle::LunarRemainder] 1近点月
+        ANOMALISTIC_MONTH = Const::Remainder::Lunar::ANOMALISTIC_MONTH
+        # @return [Cycle::LunarRemainder] 弦
+        QUARTER = Const::Remainder::Lunar::QUARTER
+
+        #
+        # 初期化
+        #
+        # @param [Cycle::LunarRemainder] lunar_age 天正閏余（大余小余）
+        # @param [Integer] western_year 西暦年
+        #
+        def initialize(lunar_age:, western_year:)
+          super(lunar_age: lunar_age, western_year: western_year)
+        end
+
+        #
+        # 入暦を計算する
+        #
+        def run
+          if calculated
+            decrease(limit: ANOMALISTIC_MONTH)
+            return
+          end
+
+          first
+        end
+
+        #
+        # 弦の分だけ月地点を進める
+        #
+        def add_quarter
+          remainder.add!(QUARTER)
+        end
+
+        private
+
+        #
+        # 初回計算
+        #
+        def first
+          @remainder = Localization.first_remainder(
+            lunar_age: remainder, western_year: western_year
+          )
+          decrease(limit: ANOMALISTIC_MONTH)
+
+          @calculated = true
+        end
+
+        #
+        # 大余小余に合わせて減算する（折り返す）
+        #
+        # @param [Cycle::LunarRemainder] limit 上限
+        #
+        def decrease(limit:)
+          return if remainder < limit
+
+          remainder.sub!(limit)
+        end
+      end
+    end
+  end
+end