module Zakuro::Calculation::Solar::ChoukeiValue

ChoukeiValue 再考長慶宣明暦算法

Constants

LOGGER

@return [Output::Logger] ロガー

Public Class Methods

get(remainder:, row:) click to toggle source

補正値を返す

@param [Cycle::AbstractRemainder] remainder 大余小余 @param [Adjustment::Row] row 24気損益眺朒(ちょうじく)数

@return [Integer] 補正値

# File lib/zakuro/calculation/stella/solar/choukei_value.rb, line 26
def self.get(remainder:, row:)
  # 損益率/眺朒(ちょうじく)数
  # パラメータ:
  #  a: 眺朒(ちょうじく)数の初日の値
  #  b: 損益率初日の値
  #  c: 損益率の毎日の差
  #  n: 定気の日から数えた日数

  # LOGGER.debug("row.per_day: #{row.per_day}")
  # LOGGER.debug("row.stack: #{row.stack}")

  day_stack = calc_day_stack(remainder: remainder, row: row)

  # LOGGER.debug("day_stack: #{day_stack}")

  month_stack = calc_month_stack(row: row, day: remainder.day)

  # LOGGER.debug("month_stack: #{month_stack}")

  # 冬至であれば眺朒数がプラスになり続けて損益率が「益」で、小雪であればマイナスの眺朒数がプラスされ続けて「損」
  month_stack + day_stack
end

Private Class Methods

calc_day_stack(remainder:, row:) click to toggle source

損益率を求める

@param [Remainder] remainder 入定気 @param [Adjustment::Row] row 24気損益眺朒(ちょうじく)数

@return [Integer] 損益率

# File lib/zakuro/calculation/stella/solar/choukei_value.rb, line 57
def self.calc_day_stack(remainder:, row:)
  ratio = calc_ratio(day: remainder.day, per_term: row.per_term, per_day: row.per_day)

  # LOGGER.debug("ratio.sign: #{ratio.sign}")
  # LOGGER.debug("ratio.abs: #{ratio.abs}")

  calc_day_stack_from_ratio(
    ratio: ratio, minute: remainder.minute, limit: remainder.base_day
  )
end
calc_day_stack_from_ratio(ratio:, minute:, limit:) click to toggle source

小余を含めた損益率を求める

@param [Integer] sign 正負(大余に対応する損益率) @param [Integer] ratio 大余に対応する損益率 @param [Integer] minute 小余

@return [Integer] 小余を含めた損益率

# File lib/zakuro/calculation/stella/solar/choukei_value.rb, line 97
def self.calc_day_stack_from_ratio(ratio:, minute:, limit:)
  minute_stack = ratio.abs * minute
  day_stack = (minute_stack / limit).floor
  # 四捨五入
  # NOTE 資料では「この余りが4200をこえていれば切り上げる」とあり「>=」とした
  # 1612年の7月(慶長17年7月)が境界値4200だが、繰り上げを行なっていたため
  day_stack += 1 if minute_stack % limit >= (limit / 2)
  day_stack *= ratio.sign

  day_stack
end
calc_month_stack(row:, day:) click to toggle source

眺朒(ちょうじく)数を求める

眺朒(ちょうじく)数 = a + (n * b) + (1/2)n(n-1)c

@param [Adjustment::Row] row 24気損益眺朒(ちょうじく)数 @param [Integer] day 大余

@return [Integer] 眺朒(ちょうじく)数

# File lib/zakuro/calculation/stella/solar/choukei_value.rb, line 119
def self.calc_month_stack(row:, day:)
  # row.stack: 眺朒(ちょうじく)積
  # row.per_term: 眺朒(ちょうじく)数
  # row.per_day: 毎日差
  month_stack = Type::OldFloat.new(
    row.stack + day * row.per_term + \
    (1 / 2.0) * (day * (day - 1) * row.per_day)
  )
  # 切り捨て(プラスマイナスに関わらず小数点以下切り捨て)
  month_stack.floor!

  month_stack.get
end
calc_ratio(day:, per_term:, per_day:) click to toggle source

大余に対応する損益率を求める

損益率 = b + n * c

@param [Integer] day 大余 @param [Integer] per_term 眺朒(ちょうじく)数 @param [Integer] per_day 毎日差

@return [Type::OldFloat] 大余に対応する損益率

# File lib/zakuro/calculation/stella/solar/choukei_value.rb, line 79
def self.calc_ratio(day:, per_term:, per_day:)
  ratio = Type::OldFloat.new(per_term + day * per_day)
  # 小数点以下は無視する
  ratio.floor!

  ratio
end