use Test::More tests => 1;
use File::Spec;
use lib File::Spec->catfile("..","lib");
use Math::MatrixReal;
do 'funcs.pl';
$eps = 1e-6;
my $A = Math::MatrixReal->new_from_string(<<"MATRIX");
[ 1 2 3 ]
[ 5 7 11 ]
[ 23 19 13 ]
MATRIX
$b = Math::MatrixReal->new_from_cols([[0, 1, 29 ]] );
$x0 = Math::MatrixReal->new_from_cols([[1, 1, -1.1 ]] );
$sol = Math::MatrixReal->new_from_cols([[1, 1, -1 ]] );
SKIP : {
skip 'solve_GSM ? ', 1;
if ( $xn = $A->solve_GSM($x0,$b,$eps) ) {
print $xn;
ok( ($xn - $sol) < $eps, 'solve_GSM seems to work');
} else {
ok( 0, 'solve_GSM' );
}
}