BEGIN { $| = 1; print "1..5\n"; }
END {print "not ok 1\n" unless $loaded;}
use Math::MatrixReal;
$loaded = 1;
print "ok 1\n";
my $DEBUG = 0;
do 'funcs.pl';
###############################
## compute a 2x2 inverse
$matrix = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 3.0 7.0 ]
[ 2.0 5.0 ]
MATRIX
$inverse = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 5.0 -7.0 ]
[ -2.0 3.0 ]
MATRIX
ok_matrix(2, $matrix ** -1 , $inverse );
################################
## A*A^-1 should = indentity
$matrix = random_matrix(20);
$one = $matrix->clone();
$one->one();
ok_matrix(3, $matrix * $matrix ** -1, $one );
#################################
## inverse of identity is identity
ok_matrix(4, $one, $one ** -1 );
#############
ok_matrix(5, $matrix->inverse->inverse, $matrix );
### numerically instable
#my $a = random_matrix(10);
#my $b = random_matrix(10);
#ok_matrix(6, ($a*$b)->inverse, ($b->inverse * $a->inverse) );