File Coverage

blib/lib/Math/Approx/Symbolic.pm

Criterion	Covered	Total	%
statement	34	34	100.0
branch	4	6	66.6
condition			n/a
subroutine	7	7	100.0
pod	1	1	100.0
total	46	48	95.8

=head1 NAME

Math::Approx::Symbolic - Symbolic representation of interpolated polynomials

=head1 SYNOPSIS

  use Math::Approx::Symbolic;
  # ... use as you would use Math::Approx ...

  my $symbolic = $approximation->symbolic();

  # ... $symbolic is now a Math::Symbolic object.

=head1 DESCRIPTION

This module is a thin wrapper around the Math::Approx module. It subclasses
Math::Approx and adds the "symbolic" subroutine that returns a
Math::Symbolic object representing the calculated approximation.

=head2 EXPORT

None. Ever.

=cut

package Math::Approx::Symbolic;

use 5.006;
use strict;
use warnings;

use Math::Approx;
use Math::Symbolic;

use base 'Math::Approx';

our $VERSION = '0.100';

=head2 symbolic() method

This is the only method added to the ones from Math::Approx.
It takes an optional argument indicating the variable name to use
for the symbolic representation of the approximation polynomial.

It returns a Math::Symbolic object representing the approximation
polynomial.

=cut

sub symbolic {
    my $self = shift;
    my $var  = shift;
    $var = 'x' unless defined $var;
    $var = Math::Symbolic::Variable->new($var)
      unless ref($var) =~ /^Math::Symbolic::Variable/;
    my $degree = $self->{N};
    my @coeff  = @{ $self->{A} };

    my $constant = shift @coeff;

    my $cur_degree = $degree;
    my @exps;
    foreach ( reverse @coeff ) {
        push @exps,
          Math::Symbolic::Operator->new(
            '*',
            Math::Symbolic::Constant->new($_),
            (
                  $cur_degree == 1
                ? $var
                : Math::Symbolic::Operator->new(
                    '^', $var, Math::Symbolic::Constant->new($cur_degree)
                )
            )
          );
        $cur_degree--;
    }
    my $symbolic = shift @exps;
    $symbolic += $_ foreach @exps;
    return $symbolic + Math::Symbolic::Constant->new($constant);
}

1;
__END__