zakuro 0.1.5 → 0.5.0

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

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

@@ -0,0 +1,75 @@
+# frozen_string_literal: true
+
+require_relative '../../../../calculation/stella/lunar/choukei_value'
+
+require_relative '../../cycle/remainder'
+
+require_relative './adjustment'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Lunar
+      #
+      # Value 月補正値
+      #
+      module Value
+        # @return [Output::Logger] ロガー
+        LOGGER = Output::Logger.new(location: 'lunar_value')
+
+        # :reek:TooManyStatements { max_statements: 9 }
+
+        #
+        # 月の運行による補正値を算出する
+        #
+        # @param [Cycle::LunarRemainder] remainder 月の大余小余
+        #
+        # @return [Integer] 補正値
+        #
+        def self.get(remainder:)
+          valid?(remainder: remainder)
+
+          day, minute = Calculation::Lunar::ChoukeiValue.remainder_without_second(
+            remainder: remainder
+          )
+
+          # 引き当て
+          row = Adjustment.specify(day: day, minute: minute)
+
+          value = row.value
+
+          # LOGGER.debug("value.per: #{value.per}")
+          # LOGGER.debug("row.denominator: #{row.denominator}")
+          # LOGGER.debug("value.stack: #{value.stack}")
+
+          minus_minute = Adjustment.minus_minute(day: day, minute: minute)
+
+          day = Calculation::Lunar::ChoukeiValue.rounded_day(
+            per: value.per, denominator: row.denominator, minute: minus_minute
+          )
+
+          # LOGGER.debug("day: #{day}")
+
+          value.stack + day
+        end
+
+        #
+        # 大余小余を検証する
+        #
+        # @param [Cycle::LunarRemainder] remainder 大余小余
+        #
+        # @return [True] 正しい（月の位相計算に使う大余小余）
+        # @return [True] 正しくない
+        #
+        def self.valid?(remainder:)
+          return if remainder.is_a?(Cycle::LunarRemainder)
+
+          raise ArgumentError, "unmatch parameter type: #{remainder.class}"
+        end
+        private_class_method :valid?
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/origin/average_november.rb

@@ -0,0 +1,34 @@
+# frozen_string_literal: true
+
+require_relative './lunar_age'
+require_relative './winter_solstice'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Origin
+      #
+      # AverageNovember 11月経
+      #
+      module AverageNovember
+        #
+        # 11月経朔（冬至が含まれる月の1日）を求める
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Remainder] 11月経朔
+        #
+        def self.get(western_year:)
+          # 冬至
+          winter_solstice = WinterSolstice.get(western_year: western_year)
+          # 天正閏余
+          lunar_age = LunarAge.get(western_year: western_year)
+
+          winter_solstice.sub(lunar_age)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/origin/lunar_age.rb

@@ -0,0 +1,62 @@
+# frozen_string_literal: true
+
+require_relative '../../../../output/logger'
+
+require_relative '../../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Origin
+      #
+      # LunarAge 天正閏余
+      #
+      module LunarAge
+        # @return [Integer] 朔望月
+        SYNODIC_MONTH = Const::Number::Cycle::SYNODIC_MONTH
+        # @return [Integer] 一年
+        YEAR = Const::Number::Cycle::YEAR
+        # @return [Integer] 積年
+        TOTAL_YEAR = Const::Number::Stack::TOTAL_YEAR
+        # @return [Integer] 暦の開始年
+        BEGIN_YEAR = Const::Number::Stack::BEGIN_YEAR
+
+        # @return [Output::Logger] ロガー
+        LOGGER = Output::Logger.new(location: 'lunar_age')
+
+        # :reek:TooManyStatements { max_statements: 7 }
+
+        #
+        # 対象年の天正閏余（冬至より前にある11月経朔との差 = 月齢）を算出する
+        # 太陽と月の運動による補正値を算出し、その補正結果を返す
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Remainder] 天正閏余
+        #
+        def self.get(western_year:)
+          # 積年の開始から対象年までの年数
+          total = TOTAL_YEAR + western_year - BEGIN_YEAR
+
+          # 12朔望月に対する1年の余り（単位:分）
+          remainder_minute = YEAR - (SYNODIC_MONTH * 12)
+
+          # 朔望月に含まれなかった余り（単位:年）
+          remainder_year = total % SYNODIC_MONTH
+
+          LOGGER.debug("[01]: #{remainder_year}")
+
+          # 天正閏余
+          lunar_age = remainder_minute * remainder_year % SYNODIC_MONTH
+
+          LOGGER.debug("[02]: #{lunar_age}")
+
+          # 大余・小余に変換する
+          Cycle::Remainder.new(total: lunar_age)
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/origin/winter_solstice.rb

@@ -0,0 +1,55 @@
+# frozen_string_literal: true
+
+require_relative '../../../../output/logger'
+
+require_relative '../../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Origin
+      #
+      # WinterSolstice 冬至
+      #
+      module WinterSolstice
+        # @return [Integer] 通余
+        REMAINDER_ALL_YEAR = Const::Number::Derivation::REMAINDER_ALL_YEAR
+        # @return [Integer] 60日
+        SIXTY_DAYS = Const::Number::Derivation::SIXTY_DAYS
+        # @return [Integer] 積年
+        TOTAL_YEAR = Const::Number::Stack::TOTAL_YEAR
+        # @return [Integer] 暦の開始年
+        BEGIN_YEAR = Const::Number::Stack::BEGIN_YEAR
+
+        # @return [Output::Logger] ロガー
+        LOGGER = Output::Logger.new(location: 'winter_solstice')
+
+        # :reek:TooManyStatements { max_statements: 6 }
+
+        #
+        # 対象年の前年の冬至を求める
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Remainder] 前年の冬至
+        #
+        def self.get(western_year:)
+          # 積年の開始から対象年までの年数
+          total = TOTAL_YEAR + western_year - BEGIN_YEAR
+          remainder_year = total % SIXTY_DAYS
+
+          LOGGER.debug("[01]:#{remainder_year}")
+
+          # 通余を使う
+          winter_solstice_minute = (remainder_year * REMAINDER_ALL_YEAR) % SIXTY_DAYS
+
+          LOGGER.debug("[02]:#{winter_solstice_minute}")
+
+          Cycle::Remainder.new(total: winter_solstice_minute)
+        end
+      end
+    end
+  end
+end

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

@@ -0,0 +1,93 @@
+# frozen_string_literal: true
+
+require_relative '../../cycle/solar_term'
+
+require_relative '../../const/number'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Solar
+      #
+      # Adjustment 24気損益眺朒（ちょうじく）数
+      #
+      module Adjustment
+        #
+        # Row 行データ
+        #
+        class Row
+          # @return [Integer] 眺朒（ちょうじく）積
+          attr_reader :stack
+          # @return [Integer] 眺朒（ちょうじく）数
+          attr_reader :per_term
+          # @return [Integer] 毎日差
+          attr_reader :per_day
+
+          #
+          # 初期化
+          #
+          # @param [Integer] stack 眺朒（ちょうじく）積
+          # @param [Integer] per_term 眺朒（ちょうじく）数
+          # @param [Integer] per_day 毎日差
+          #
+          def initialize(stack:, per_term:, per_day:)
+            @stack = stack
+            @per_term = per_term
+            @per_day = per_day
+          end
+
+          #
+          # 文字化
+          #
+          # @return [String] 文字
+          #
+          def to_s
+            "stack:#{@stack}, per_term:#{@per_term}, per_day:#{@per_day}"
+          end
+        end
+
+        # @return [Array<Row>] 24気損益眺朒（ちょうじく）数
+        MAP = {
+          touji: Row.new(stack: 0.0, per_term: +13.4524, per_day: -0.1886),        # 冬至（とうじ）
+          shoukan: Row.new(stack: +176.0, per_term: +10.5564, per_day: -0.1634),   # 小寒（しょうかん）
+          daikan: Row.new(stack: +314.0, per_term: +8.0408, per_day: -0.1446),     # 大寒（だいかん）
+          risshun: Row.new(stack: +418.0, per_term: +5.8160, per_day: -0.1318),    # 立春（りっしゅん）
+          usui: Row.new(stack: +491.0, per_term: +3.7987, per_day: -0.1240),       # 雨水（うすい）
+          keichitsu: Row.new(stack: +535.0, per_term: +1.9265, per_day: -0.1240),  # 啓蟄（けいちつ）
+          shunbun: Row.new(stack: +551.0, per_term: -0.2048, per_day: -0.1178),    # 春分（しゅんぶん）
+          seimei: Row.new(stack: +535.0, per_term: -1.9968, per_day: -0.1190),     # 清明（せいめい）
+          kokuu: Row.new(stack: +491.0, per_term: -3.7956, per_day: -0.1240),      # 穀雨（こくう）
+          rikka: Row.new(stack: +418.0, per_term: -5.6626, per_day: -0.1324),      # 立夏（りっか）
+          shouman: Row.new(stack: +314.0, per_term: -7.6555, per_day: -0.1436),    # 小満（しょうまん）
+          boushu: Row.new(stack: +176.0, per_term: -9.9405, per_day: -0.1436),     # 芒種（ぼうしゅ）
+          geshi: Row.new(stack: 0.0, per_term: -12.0819, per_day: +0.1436),        # 夏至（げし）
+          shousho: Row.new(stack: -176.0, per_term: -9.7018, per_day: +0.1324),    # 小暑（しょうしょ）
+          taisho: Row.new(stack: -314.0, per_term: -7.5450, per_day: +0.1240),     # 大暑（たいしょ）
+          risshuu: Row.new(stack: -418.0, per_term: -5.5634, per_day: +0.1190),    # 立秋（りっしゅう）
+          shosho: Row.new(stack: -491.0, per_term: -3.7038, per_day: +0.1178),     # 処暑（しょしょ）
+          hakuro: Row.new(stack: -535.0, per_term: -1.8954, per_day: +0.1178),     # 白露（はくろ）
+          shuubun: Row.new(stack: -551.0, per_term: +0.1783, per_day: +0.1240),    # 秋分（しゅうぶん）
+          kanro: Row.new(stack: -535.0, per_term: +2.0042, per_day: +0.1318),      # 寒露（かんろ）
+          soukou: Row.new(stack: -491.0, per_term: +3.8950, per_day: +0.1446),     # 霜降（そうこう）
+          rittou: Row.new(stack: -418.0, per_term: +5.9214, per_day: +0.1634),     # 立冬（りっとう）
+          shousetsu: Row.new(stack: -314.0, per_term: +8.1610, per_day: +0.1886),  # 小雪（しょうせつ）
+          taisetsu: Row.new(stack: -176.0, per_term: +10.9010, per_day: +0.1886)   # 大雪（たいせつ）
+        }.freeze
+
+        #
+        # 24気損益眺朒（ちょうじく）数の行データを特定する
+        #
+        # @param [Integer] index 連番（二十四節気）
+        #
+        # @return [Row] 行データ
+        #
+        def self.specify(index:)
+          key = Cycle::SolarTerm::ORDER[index]
+          Adjustment::MAP[key].clone
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/solar/average.rb

@@ -0,0 +1,98 @@
+# frozen_string_literal: true
+
+require_relative '../../../../calculation/stella/solar/abstract_average'
+
+require_relative '../../const/remainder'
+
+require_relative '../../cycle/solar_term'
+
+require_relative '../origin/winter_solstice'
+
+require_relative './location'
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Solar
+      #
+      # Average 平気（太陽軌道平均）
+      #
+      class Average < Calculation::Solar::AbstractAverage
+        #
+        # 初期化
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        def initialize(western_year:)
+          solar_term = Average.first_solar_term(western_year: western_year)
+          super(solar_term: solar_term)
+        end
+
+        #
+        # 冬至から数えた1年データの月ごとに二十四節気を割り当てる
+        #
+        # @param [Array<Month>] annual_range 1年データ
+        #
+        # @return [Array<Month>] 1年データ
+        #
+        def set(annual_range:)
+          super(annual_range: annual_range)
+        end
+
+        #
+        # 計算開始する二十四節気を求める
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [SolarTerm] 二十四節気
+        #
+        def self.first_solar_term(western_year:)
+          # 天正冬至
+          winter_solstice = Origin::WinterSolstice.get(western_year: western_year)
+
+          # 二十四節気（冬至）
+          solar_term = Cycle::SolarTerm.new(index: 0, remainder: winter_solstice)
+
+          first_solar_term_index = Average.calc_first_solar_term_index(western_year: western_year)
+
+          # 対象の二十四節気まで戻す
+          solar_term.prev_by_index(first_solar_term_index)
+
+          solar_term
+        end
+
+        # :reek:TooManyStatements { max_statements: 6 }
+
+        #
+        # 計算開始する二十四節気番号を求める
+        #
+        # * 前提として入定気は冬至の手前にある
+        # * 例えば、定気が大雪であれば入定気は大雪の範囲内にある
+        # * 入定気は、定気の開始位置に重複しない限り、常に定気より後にある
+        # * 基本的に定気の一つ前から起算すれば、当時から求めた11月(閏10/閏11月)に二十四節気を割り当てられる
+        #
+        # @param [Integer] western_year 西暦年
+        #
+        # @return [Integer] 二十四節気番号
+        #
+        def self.calc_first_solar_term_index(western_year:)
+          # 天正閏余
+          lunar_age = Origin::LunarAge.get(western_year: western_year)
+
+          solar_location = Solar::Location.new(lunar_age: lunar_age)
+          solar_location.run
+
+          solar_term_index = solar_location.index
+
+          # 入定気の一つ後の二十四節気まで戻す（ただし11月経朔が二十四節気上にある場合は戻さない）
+          # NOTE: 789年で大余0のパターンあり
+          solar_term_index += 1 unless solar_location.remainder.day.zero?
+
+          solar_term_index
+        end
+      end
+    end
+  end
+end

data/lib/zakuro/version/daien/stella/solar/interval.rb

@@ -0,0 +1,108 @@
+# frozen_string_literal: true
+
+# :nodoc:
+module Zakuro
+  # :nodoc:
+  module Daien
+    # :nodoc:
+    module Solar
+      #
+      # Interval 入気定日加減数（二十四節気の間隔）
+      #
+      module Interval
+        #
+        # 気策の 15-664.7 を加減する
+        #
+        # @note 『歴代天文律暦等志彙編　七』中華書房 p.2224 盈縮分の行で加減した
+        #
+        # @return [Hash<Symbol, Cycle::Remainder>] 一覧
+        MAP = {
+          # 冬至（とうじ）・大雪（たいせつ）
+          touji: Cycle::Remainder.new(day: 14, minute: 1351, second: 7), # -2353
+          taisetsu: Cycle::Remainder.new(day: 14, minute: 1351, second: 7),
+          # 小寒（しょうかん）・小雪（しょうせつ）
+          shoukan: Cycle::Remainder.new(day: 14, minute: 1859, second: 7), # -1845
+          shousetsu: Cycle::Remainder.new(day: 14, minute: 1859, second: 7),
+          # 大寒（だいかん）・立冬（りっとう）
+          daikan: Cycle::Remainder.new(day: 14, minute: 2314, second: 7), # -1390
+          rittou: Cycle::Remainder.new(day: 14, minute: 2314, second: 7),
+          # 立春（りっしゅん）・霜降（そうこう）
+          risshun: Cycle::Remainder.new(day: 14, minute: 2728, second: 7), # -976
+          soukou: Cycle::Remainder.new(day: 14, minute: 2728, second: 7),
+          # 雨水（うすい）・寒露（かんろ）
+          usui: Cycle::Remainder.new(day: 15, minute: 76, second: 7), # -588
+          kanro: Cycle::Remainder.new(day: 15, minute: 76, second: 7),
+          # 啓蟄（けいちつ）・秋分（しゅうぶん）
+          keichitsu: Cycle::Remainder.new(day: 15, minute: 450, second: 7), # -214
+          shuubun: Cycle::Remainder.new(day: 15, minute: 450, second: 7),
+          # 春分（しゅんぶん）・白露（はくろ）
+          shunbun: Cycle::Remainder.new(day: 15, minute: 878, second: 7), # +214
+          hakuro: Cycle::Remainder.new(day: 15, minute: 878, second: 7),
+          # 清明（せいめい）・処暑（しょしょ）
+          seimei: Cycle::Remainder.new(day: 15, minute: 1252, second: 7), # +588
+          shosho: Cycle::Remainder.new(day: 15, minute: 1252, second: 7),
+          # 穀雨（こくう）・立秋（りっしゅう）
+          kokuu: Cycle::Remainder.new(day: 15, minute: 1640, second: 7), # +976
+          risshuu: Cycle::Remainder.new(day: 15, minute: 1640, second: 7),
+          # 立夏（りっか）・大暑（たいしょ）
+          rikka: Cycle::Remainder.new(day: 15, minute: 2054, second: 7), # +1390
+          taisho: Cycle::Remainder.new(day: 15, minute: 2054, second: 7),
+          # 小満（しょうまん）・小暑（しょうしょ）
+          shouman: Cycle::Remainder.new(day: 15, minute: 2509, second: 7), # +1845
+          shousho: Cycle::Remainder.new(day: 15, minute: 2509, second: 7),
+          # 芒種（ぼうしゅ）・夏至（げし）
+          boushu: Cycle::Remainder.new(day: 15, minute: 3017, second: 7), # +2353
+          geshi: Cycle::Remainder.new(day: 15, minute: 3017, second: 7)
+        }.freeze
+
+        # @return [Array<Remainder>] 索引
+        LIST = [
+          MAP[:touji],      #  0
+          MAP[:shoukan],    #  1
+          MAP[:daikan],     #  2
+          MAP[:risshun],    #  3
+          MAP[:usui],       #  4
+          MAP[:keichitsu],  #  5
+          MAP[:shunbun],    #  6
+          MAP[:seimei],     #  7
+          MAP[:kokuu],      #  8
+          MAP[:rikka],      #  9
+          MAP[:shouman],    # 10
+          MAP[:boushu],     # 11
+          MAP[:geshi],      # 12
+          MAP[:shousho],    # 13
+          MAP[:taisho],     # 14
+          MAP[:risshuu],    # 15
+          MAP[:shosho],     # 16
+          MAP[:hakuro],     # 17
+          MAP[:shuubun],    # 18
+          MAP[:kanro],      # 19
+          MAP[:soukou],     # 20
+          MAP[:rittou],     # 21
+          MAP[:shousetsu],  # 22
+          MAP[:taisetsu]    # 23
+        ].freeze
+
+        #
+        # 入気定日加減数を引き当てる
+        #
+        # @param [Integer] index 二十四節気番号
+        #
+        # @return [Cycle::Remainder] 入気定日加減数
+        #
+        def self.index_of(index)
+          LIST[index]
+        end
+
+        #
+        # 入気定日加減数の総数を返す
+        #
+        # @return [Integer] 入気定日加減数の総数
+        #
+        def self.size
+          LIST.size
+        end
+      end
+    end
+  end
+end