@@ -4,11 +4,7 @@ use PDL;
4
4
use PDL::FFT;
5
5
use Test::More;
6
6
use Test::Exception;
7
-
8
- sub tapprox {
9
- my ($pa ,$pb ) = @_ ;
10
- all approx $pa , $pb , 0.01;
11
- }
7
+ use Test::PDL -atol => 0.01, -require_equal_types => 0;
12
8
13
9
foreach my $type (double,float,cdouble,cfloat){
14
10
my $pa = pdl($type ,1,-1,1,-1);
@@ -29,35 +25,35 @@ fft($pb,$pc);
29
25
ifft($pb ,$pc );
30
26
fft($pd );
31
27
ifft($pd );
32
- ok (tapprox( $pc ,0) , " fft zeroes" ) ;
33
- ok (tapprox( $pd -> im,0) , " fft zeroes using complex ndarrays" ) ;
34
- ok (tapprox( $pa ,$pb ) , " original image recovered" ) ;
28
+ is_pdl $pc , $pb -> zeroes , " fft zeroes" ;
29
+ is_pdl $pd -> im, $pb -> zeroes , " fft zeroes using complex ndarrays" ;
30
+ is_pdl $pa , $pb , " original image recovered" ;
35
31
}
36
32
37
33
{
38
34
my $pb = $pa -> copy;
39
35
my $pc = $pb -> zeroes;
40
36
my $pd =czip($pb , $pc );
41
37
fftnd($pb ,$pc ); ifftnd($pb ,$pc );
42
- ok ( tapprox( $pc ,0) , " fftnd zeroes" ) ;
43
- ok ( tapprox( $pa ,$pb ) , " fftnd real image" ) ;
38
+ is_pdl $pc , $pb -> zeroes , " fftnd zeroes" ;
39
+ is_pdl $pa , $pb , " fftnd real image" ;
44
40
fftnd($pd ); ifftnd($pd );
45
- ok ( tapprox( $pd ,$pb ) , " fftnd native complex image with imag zeroes" ) ;
41
+ is_pdl $pd , $pb , " fftnd native complex image with imag zeroes" ;
46
42
}
47
43
48
44
{
49
45
my $pb = $pa -> slice(" 1:35,1:69" );
50
46
my $pc = $pb -> copy; fftnd($pb ,$pc ); ifftnd($pb ,$pc );
51
- ok ( tapprox( $pc ,$pb ) , " fftnd real and imaginary" ) ;
52
- ok ( tapprox( $pa -> slice(" 1:35,1:69" ),$pb ) , " fftnd original restored" ) ;
47
+ is_pdl $pc , $pb , " fftnd real and imaginary" ;
48
+ is_pdl $pa -> slice(" 1:35,1:69" ), $pb , " fftnd original restored" ;
53
49
}
54
50
55
51
{
56
52
my $pb = $pa -> copy;
57
53
# Test real ffts
58
54
realfft($pb );
59
55
realifft($pb );
60
- ok( tapprox( $pa ,$pb ) , " realfft" ) ;
56
+ is_pdl $pa , $pb , " realfft" ;
61
57
}
62
58
63
59
# Test that errors are properly caught
0 commit comments