Class: When::Coordinates::Residue::BestRationalApproximations

Inherits:

Parts::Enumerator

• Object
• Enumerator
• Parts::Enumerator
• When::Coordinates::Residue::BestRationalApproximations

Defined in:

lib/when_exe/coordinates.rb

Overview

最良近似分数の系列を生成する Enumerator

Instance Attribute Summary

Attributes inherited from Parts::Enumerator

#count, #count_limit, #current, #direction, #exdate, #first, #index, #last, #object, #options, #parent, #processed

Instance Method Summary

• #_rewind ⇒ void

  Enumerator の巻き戻し.

• #initialize(parent, options = {}) ⇒ BestRationalApproximations constructor

  オブジェクトの生成.

• #succ ⇒ Array<Numeric>

  最良近似分数を生成する.

Methods inherited from Parts::Enumerator

_options, #each, #has_next?, #next, #with_index, #with_object

Constructor Details

#initialize(parent, options = {}) ⇒ BestRationalApproximations

オブジェクトの生成

Parameters:

• parent (When::Coordinates::Residue) —

  生成元の剰余類オブジェクト

• options (Hash) (defaults to: {}) —

  下記の通り

Options Hash (options):

• :error (Numeric) —

  収束とみなす誤差(デフォルト 1E-5)

• :count_limit (Integer) —

  最大繰り返し回数

# File 'lib/when_exe/coordinates.rb', line 530

def initialize(*args)
  @y = args[0].divisor
  @x = args[0].remainder + @y * args[0].carry
  super
  @error       = @options[:error] || 1e-15
  @count_limit = @options[:count_limit]
end

Instance Method Details

#_rewind ⇒ void

This method returns an undefined value.

Enumerator の巻き戻し

# File 'lib/when_exe/coordinates.rb', line 489

def _rewind
  @z = @x/@y
  @k = @z.floor
  @p = [1,@k]
  @q = [0, 1]
  super
end

#succ ⇒ Array<Numeric>

最良近似分数を生成する

Returns:

• (Array<Numeric>) —

  ( remainder, divisor, error )

  remainder (Integer) 分子

  divisor (Integer) 分母

  error (Float) 誤差

# File 'lib/when_exe/coordinates.rb', line 504

def succ
  value = @current
  if (@z==@k) ||
     (@count_limit.kind_of?(Numeric) && @count >= @count_limit) ||
     (@error.kind_of?(Numeric) && @e && @error >= @e.abs)
    @current = nil
  else
    @z = 1.0/(@z-@k)
    @k = @z.floor
    @e = @p[1].to_f/@q[1]-@x.to_f/@y
    @current = [@p[1], @q[1], @e]
    @p = [@p[1], @p[1]*@k + @p[0]]
    @q = [@q[1], @q[1]*@k + @q[0]]
    @count += 1
  end
  return value
end