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
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
|
#!/usr/bin/perl -w
# -*- perl -*-
#
# $Id: base.t,v 0.1 2006/02/21 eserte Exp $
# Author: Michael R. Davis
#
=head1 NAME
base.t - Good examples concerning how to use this module
=cut
use strict;
use lib q{lib};
use lib q{../lib};
use constant NEAR_DEFAULT => 7;
sub near {
my $x=shift();
my $y=shift();
my $p=shift()||NEAR_DEFAULT;
if (($x-$y)/$y < 10**-$p) {
return 1;
} else {
return 0;
}
}
BEGIN {
if (!eval q{
use Test;
1;
}) {
print "1..0 # tests only works with installed Test module\n";
exit;
}
}
BEGIN { plan tests => 49 }
# just check that all modules can be compiled
ok(eval {require Geo::Ellipsoids; 1}, 1, $@);
my $o = Geo::Ellipsoids->new();
ok(ref $o, "Geo::Ellipsoids");
my $obj=Geo::Ellipsoids->new();
$obj->set({a=>1,b=>1}); # a 1 unit sphere
ok($obj->a, 1);
ok($obj->b, 1);
ok($obj->i, undef);
ok($obj->f, 0);
$obj->set({a=>1,b=>0.995});
ok($obj->a, 1);
ok($obj->b, 0.995);
ok($obj->i, 200);
ok($obj->f, 0.005);
$obj->set({a=>1,i=>200});
ok($obj->a, 1);
ok($obj->b, 0.995);
ok($obj->i, 200);
ok($obj->f, 0.005);
$obj->set({a=>1,f=>0.005});
ok($obj->a, 1);
ok($obj->b, 0.995);
ok($obj->i, 200);
ok($obj->f, 0.005);
$obj->set({a=>1});
ok($obj->a, 1);
ok($obj->b, 1);
ok($obj->i, undef());
ok($obj->f, 0);
ok($obj->e, 0);
ok near($obj->equatorial_circumference, 8 * atan2(1,1), 13);
ok near($obj->polar_circumference, 8 * atan2(1,1), 13);
$obj->set('WGS84');
ok($obj->a, 6378137);
ok near($obj->b, 6356752.31424518, 13);
ok($obj->i, 298.257223563);
ok near($obj->f, 0.00335281066474748, 12);
#For the Clarke 1866, e² is 0.006768658.
$obj->set('Clarke 1866');
ok near($obj->e2, 0.006768658, 7);
#For the GRS 80, e² is 0.0066943800.
$obj->set('GRS80');
ok near($obj->e2, 0.0066943800, 6);
$obj->set('WGS84');
ok($obj->shortname, "WGS84");
ok($obj->longname, "World Geodetic System of 1984");
ok near($obj->equatorial_circumference, 40075016.6855785, 13);
ok near($obj->polar_circumference, 39940652.7422451, 13);
#Just testing that n works. It is not my formula.
ok(near $obj->n(39.56789), 6386817.167912991, 13);
ok(near $obj->n_rad(0.690589958566939), 6386817.167912991, 13);
my $a=$obj->a;
my $b=$obj->b;
my $i=$obj->i;
my $f=$obj->f;
my $e=$obj->e;
my $e2=$obj->e2;
my $ep2=$obj->ep2;
#Run through a bunch of identities
ok(near $i, 1/$f, 13);
ok(near $i, $a/($a-$b), 13);
ok(near $f, ($a-$b)/$a, 13);
ok(near $b, $a*(1-$f), 13);
ok(near 1, (1-$e2)*(1+$ep2), 13);
ok(near $e2, $f*(2-$f), 13);
ok(near $e2, 2*$f-$f**2, 13);
ok(near $e2, 1-$b**2/$a**2, 13);
ok(near $ep2, ($e*$a/$b)**2, 13);
ok(near $ep2, $e2/(1-$e2), 13);
ok(near $ep2, $a**2/$b**2-1, 13);
ok(near $ep2, ($a**2-$b**2)/$b**2, 13);
|
