# -*-perl-*-
use strict;
use Test;
BEGIN {
plan tests => 6;
}
use PDL::LiteF;
sub tapprox {
my($a,$b) = @_;
my $d = max( abs($a-$b) );
$d < 1.0e-6;
}
# test rvals
# synonym for centre/center
# squared option
#
my $x0 = pdl( [ 2, 1, 2 ], [ 1, 0, 1 ], [ 2, 1, 2 ] );
my $a1 = rvals(3,3);
print "\na1: $a1\n";
ok( tapprox( $x0->sqrt, $a1 ), 1 ); # 1
my $a2 = rvals(3,3,{squared=>1});
print "\na2: $a2\n";
ok( tapprox( $x0, $a2 ), 1 ); # 2
my $x1 = pdl( [ 8, 5, 4 ], [ 5, 2, 1 ], [ 4, 1, 0 ] );
my $a3 = rvals(3,3,{centre=>[2,2]});
print "\na3: $a3\n";
ok( tapprox( $x1->sqrt, $a3 ), 1 ); # 3
my $a4 = rvals(3,3,{center=>[2,2]});
print "\na4: $a4\n";
ok( tapprox( $x1->sqrt, $a4 ), 1 ); # 4
ok( tapprox( $x1->sqrt, rvals(3,3,{ceNteR=>[2,2]}) ), 1 ); # 5
ok( tapprox( $x1, rvals(3,3,{center=>[2,2],squared=>1}) ), 1 ); # 6