Algorithm::Simplex::Types - Types into which we coerce matrix input


Algorithm-Simplex documentation Contained in the Algorithm-Simplex distribution.

Index


Code Index:

Name

Top

Algorithm::Simplex::Types - Types into which we coerce matrix input for PDL and Rational models

Methods

Top

make_fractions

Make each rational entry a Math::Cephes::Fraction object with the help of Math::BigRat

display_fractions

Convert each fraction object entry into a string.

display_piddle

Convert a PDL into an ArrayRef[ArrayRef[Num]]


Algorithm-Simplex documentation Contained in the Algorithm-Simplex distribution.
package Algorithm::Simplex::Types;
use strict;
use warnings;
use Moose::Util::TypeConstraints;
use PDL::Lite;
use Math::BigRat;
use Math::Cephes::Fraction qw(:fract);
use Data::Dumper;

subtype 'Piddle' 
    => as 'PDL' 
    => where { $_->isa('PDL') }
    => message { "This thingy $_ is not a Piddle!" };
coerce 'Piddle' 
    => from 'ArrayRef[ArrayRef[Num]]' 
    => via { PDL->pdl($_) };
    
subtype 'PiddleDisplay'
    => as 'ArrayRef[ArrayRef[Str]]',
    => where { 1 }
    => message { "This thingys $_ is not a PiddleDisplay!"};
coerce 'PiddleDisplay'
    => from 'Piddle'
    => via { &display_piddle($_) };      

subtype 'FractionMatrix' 
    => as 'ArrayRef[ArrayRef[Math::Cephes::Fraction]]' 
    => where { 1 } 
    => message { "This thingy $_ is not a matrix of Fraction objects.  It is a " . Dumper $_ };
coerce 'FractionMatrix' 
    => from 'ArrayRef[ArrayRef[Num]]' 
    => via { &make_fractions($_) };
        
subtype 'FractionDisplay'
    => as 'ArrayRef[ArrayRef[Str]]'
    => where { 1 }
    => message { "This thingy $_ is not a FractDisplay!" };
coerce 'FractionDisplay'
    => from 'FractionMatrix'    
    => via { &display_fractions($_) };

sub make_fractions {
    my $tableau = shift;
    
    for my $i ( 0 ..  scalar @{ $tableau } - 1 ) {
        for my $j ( 0 .. scalar @{ $tableau->[0] } - 1 ) {
            # Using Math::BigRat to make fraction from decimal
            my $x = Math::BigRat->new( $tableau->[$i]->[$j] );
            $tableau->[$i]->[$j] = fract( $x->numerator, $x->denominator);
        }
    }
    return $tableau;
}

sub display_fractions {
    my $fraction_tableau = shift;

    my $display_tableau;
    for my $i ( 0 .. scalar @{ $fraction_tableau } - 1 ) {
        for my $j ( 0 .. scalar @{ $fraction_tableau->[0] } - 1  ) {
            $display_tableau->[$i]->[$j] = $fraction_tableau->[$i]->[$j]->as_string;
        }
    }
    return $display_tableau;

}

sub display_piddle {
    my $piddle_tableau = shift;

    my @display_tableau;
    my ($number_of_columns, $number_of_rows) = ($piddle_tableau->dims);
    my $number_of_zero_based_rows = $number_of_rows - 1;
    my $number_of_zero_based_columns = $number_of_columns - 1;
    for my $i ( 0 .. $number_of_zero_based_rows ) {
        my $row = $piddle_tableau->slice("0:$number_of_zero_based_columns,($i)");
        my @row   = $row->list;
        push @display_tableau, \@row;
    }
    
    return \@display_tableau;
}

1;