File Coverage

blib/lib/Math/Rational/Approx/ContFrac.pm

Criterion	Covered	Total	%
statement	33	36	91.6
branch	2	2	100.0
condition			n/a
subroutine	11	13	84.6
pod	3	3	100.0
total	49	54	90.7

# --8<--8<--8<--8<--
#
# Copyright (C) 2012 Smithsonian Astrophysical Observatory
#
# This file is part of Math::Rational::Approx::ContFrac
#
# Math::Rational::Approx::ContFrac is free software: you can
# redistribute it and/or modify it under the terms of the GNU General
# Public License as published by the Free Software Foundation, either
# version 3 of the License, or (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program.  If not, see .
#
# -->8-->8-->8-->8--

package Math::Rational::Approx::ContFrac;

use strict;
use warnings;
use Carp;

our $VERSION = '0.01';

use Math::BigFloat;

use Moo;
use MooX::Types::MooseLike::Numeric ':all';

use Params::Validate qw[ validate_pos ARRAYREF ];

use Math::Rational::Approx qw[ contfrac contfrac_nd ];

has x => (
  is  => 'ro',
  isa => sub { die( "must be a positive number\n" )
                 unless is_PositiveNum($_[0]) },
  required => 1,
);

has n => (
  is  => 'rwp',
  isa => PositiveInt,
  required => 1,
);


has _terms => (
  is  => 'rwp',
  init_arg => undef,
  default => sub { [] },
);

sub terms { [ @{$_[0]->_terms} ] }

has _resid => (
  is => 'rwp',
  init_arg => undef,
  lazy => 1,
  builder => '_build_resid',
);

sub resid { $_[0]->x->copy }


sub _build_resid  { Math::BigFloat->new( $_[0]->x )  }

sub approx {

	my $self = shift;

	 my ( $n ) = validate_pos( @_,
	                         { optional => 1,
	                           callbacks => {
	                              'positive integer' => sub { is_PositiveInt($_[0]) },
	                                        },
	                           });

	$self->_set_n( $self->n + $n )
	  if defined $n;

	my ( undef, $x ) = contfrac( $self->_resid, $self->n - @{$self->_terms}, $self->_terms );
	$self->_set__resid( $x );

	return contfrac_nd( $self->_terms );
}

1;


__END__