rand.t

use Test::More tests => 16;
use File::Spec;
use lib File::Spec->catfile("..","lib");
use Math::MatrixReal;
my ($e,$res) = (0,0);
my $eps = 1e-8;

do './funcs.pl';

{
    my $matrix = Math::MatrixReal->new_random( 10,10, { integer => 1 } );
    ok ( ref $matrix eq 'Math::MatrixReal' , 'new_random returns the correct object' );

    my ($rows,$cols) = $matrix->dim;
    ok( $rows == 10 && $cols == 10, 'new_random returns the correct size' );
    for my $r ( 1 .. $rows ){
        for my $c ( 1 .. $cols ) {
            $e = $matrix->element($r,$c);
            $res += abs( $e-int($e) );
        }
    }
    ok( $res < $eps, 'new_random option type integer works' );
}
{
    $matrix = Math::MatrixReal->new_random( 5 );
    ($rows,$cols) = $matrix->dim;
    ok( $rows == 5 && $cols == 5, 'new_random is square if called with one argument' );
}

{
    ($rows,$cols) = (1+int(rand(10)), 1+int(rand(10)) );
    my $matrix = Math::MatrixReal->new_random( $rows,$cols, { bounded_by => [-$rows, $rows] } );
    my $min = $matrix->element(1,1); 
    my $max = $min;
    for my $r ( 1 .. $rows ){
        for my $c ( 2 .. $cols ) {
            $e = $matrix->element($r,$c);
            $e < $min ? $min = $e :  $e > $max ? $max = $e : 0  ;
        }
    }
    ok( $min >= -$rows && $max <= $rows, 'new_random option bounded_by works' );

}
{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random },
                 q{new_random fails with no args} 
    );
}
{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(0, 17.5) },
                 q{new_random fails with invalid args}
    );
}

{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(10,20, { bounded_by => [] } ) }, 
                 q{new_random fails with invalid bounded_by}
    );
}
{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(10,20, { bounded_by => [1,-1] } ) },
                 q{new_random fails with invalid bounded_by range}
    );
}
{ 
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(10,20, { symmetric => 1 } ) },
                 q{new_random fails with rectangular + symmetric}
    );
}

{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(10,20, { tridiag => 1 } ) },
	             q{new_random fails with nonsquare tridiag}
    );
}


{
    assert_dies( sub { my $matrix = Math::MatrixReal->new_random(10,20, { diag => 1 } ) },
                 q{new_random fails with nonsquare diag},
    );
}

{
    ok( Math::MatrixReal->new_random(10, { symmetric => 1 } )->is_symmetric, 
        'new_random can do symmetric');
}
{
    ok( Math::MatrixReal->new_random(5, { tridiag => 1, integer => 1 } )->is_tridiagonal, 
        'new_random with tridiag works');
}
{
    my $a = Math::MatrixReal->new_random(5, { tridiag => 1, symmetric => 1 } );
    ok( $a->is_tridiagonal && $a->is_symmetric,
       'new_random with tridiag+symmetric works');
}

{
    ok( Math::MatrixReal->new_random(5, { diag => 1, integer => 1 } )->is_diagonal,
	    'new_random with diag works');
}