-
Notifications
You must be signed in to change notification settings - Fork 135
/
unpolar.t
82 lines (67 loc) · 2.63 KB
/
unpolar.t
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
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
use v6;
use Test;
plan 199;
# L<S32::Numeric/Numeric/=item cis>
my $pi = 312689/99532;
{
is-approx(cis(0), 1 + 0i, "cis(0) == 1");
is-approx(cis($pi), -1 + 0i, "cis(pi) == -1");
is-approx(cis($pi / 2), 1i, "cis(pi/2) == i");
is-approx(cis(3*$pi / 2),-1i, "cis(3pi/2) == -i");
}
# Test that 1.unpolar == cis
# L<S32::Numeric/Numeric/=item cis>
# L<S32::Numeric/Numeric/=item unpolar>
{
for 1...20 -> $i {
my $angle = 2 * $pi * $i / 20;
is-approx(cis($i), 1.unpolar($i), "cis(x) == 1.unpolar(x) No. $i");
is-approx($i.cis, 1.unpolar($i), "x.cis == 1.unpolar(x) No. $i");
is-approx($i.Rat.cis, 1.unpolar($i), "x.Rat.cis == 1.unpolar(x) No. $i");
is-approx($i.Num.cis, 1.unpolar($i), "x.Num.cis == 1.unpolar(x) No. $i");
}
}
# L<S32::Numeric/Numeric/=item abs>
# L<S32::Numeric/Numeric/=item unpolar>
#
# Test that unpolar() doesn't change the absolute value
{
my $counter = 1;
for 1...10 -> $abs {
for 1...10 -> $a {
my $angle = 2 * $pi * $a / 10;
is-approx($abs.unpolar($angle).abs, $abs,
"unpolar doesn't change the absolute value (No. $counter)");
$counter++;
}
}
}
# L<S32::Numeric/Numeric/=item unpolar>
{
# Basic tests for unpolar()
my $s = 2 * sqrt(2);
is-approx(4.unpolar(0), 4, "4.unpolar(0) == 4");
is-approx(4.unpolar($pi/4), $s + ($s)i ,"4.unpolar(pi/4) == 2+2i");
is-approx(4.unpolar($pi/2), 4i, "4.unpolar(pi/2) == 4i");
is-approx(4.unpolar(3.Num*$pi/4), -$s + ($s)i,"4.unpolar(3*pi/4) == -2+2i");
is-approx(4.unpolar($pi), -4, "4.unpolar(pi) == -4");
}
{
# Basic tests for unpolar()
my $s = 2 * sqrt(2);
is-approx(4.Rat.unpolar(0), 4, "4.Rat.unpolar(0) == 4");
is-approx(4.Rat.unpolar($pi/4), $s + ($s)i ,"4.Rat.unpolar(pi/4) == 2+2i");
is-approx(4.Rat.unpolar($pi/2), 4i, "4.Rat.unpolar(pi/2) == 4i");
is-approx(4.Rat.unpolar(3.Num*$pi/4), -$s + ($s)i,"4.Rat.unpolar(3*pi/4) == -2+2i");
is-approx(4.Rat.unpolar($pi), -4, "4.Rat.unpolar(pi) == -4");
}
{
# Basic tests for unpolar()
my $s = 2 * sqrt(2);
is-approx(4.Num.unpolar(0), 4, "4.Num.unpolar(0) == 4");
is-approx(4.Num.unpolar($pi/4), $s + ($s)i ,"4.Num.unpolar(pi/4) == 2+2i");
is-approx(4.Num.unpolar($pi/2), 4i, "4.Num.unpolar(pi/2) == 4i");
is-approx(4.Num.unpolar(3.Num*$pi/4), -$s + ($s)i,"4.Num.unpolar(3*pi/4) == -2+2i");
is-approx(4.Num.unpolar($pi), -4, "4.Num.unpolar(pi) == -4");
}
# vim: ft=perl6