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