File: physical.t

package info (click to toggle)
pdl 1%3A2.4.7%2Bdfsg-2
  • links: PTS
  • area: main
  • in suites: squeeze
  • size: 10,128 kB
  • ctags: 5,821
  • sloc: perl: 26,328; fortran: 13,113; ansic: 9,378; makefile: 71; sh: 50; sed: 6
file content (44 lines) | stat: -rw-r--r-- 605 bytes parent folder | download | duplicates (7) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
 
# Test ->*physical*(). This is not yet good enough: we need
# nasty test cases,

use PDL::LiteF;

sub ok {
	my $no = shift ;
	my $result = shift ;
	print "not " unless $result ;
	print "ok $no\n" ;
}

sub tapprox {
	my($a,$b) = @_;
	$c = abs($a-$b);
	$d = max($c);
	$d < 0.01;
}

# Cheat
print "1..1\n";
print "ok 1\n";

if(0) {

print "1..6\n";

$a = zeroes(4,4);
ok(1,$a->isphysical());

$b = xvals $a + 0.1 * yvals $a;
$c = $b->slice("1:3:2,:");
ok(2,! $c->isphysical());

$d = $b->physical();
ok(3, $d == $b);

$e = $c->physical();
ok(4, $e != $c);
ok(5, $e->isphysical());

ok(6, tapprox($c,$e));

}