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
|
use strict;
use warnings;
use Test::More;
use DateTime;
# Borrowed from Matt Sergeant's Time::Piece
# A table of MJD and components
my @mjd = (
'51603.524' => {
year => 2000,
month => 2,
day => 29,
hour => 12,
minute => 34,
second => 56,
},
'40598.574' => {
year => 1970,
month => 1,
day => 12,
hour => 13,
minute => 46,
second => 51,
},
'52411.140' => {
year => 2002,
month => 5,
day => 17,
hour => 3,
minute => 21,
second => 43,
},
'53568.547' => {
year => 2005,
month => 7,
day => 17,
hour => 13,
minute => 8,
second => 23,
},
'52295.218' => {
year => 2002,
month => 1,
day => 21,
hour => 5,
minute => 13,
second => 20,
},
'52295.399' => {
year => 2002,
month => 1,
day => 21,
hour => 9,
minute => 35,
second => 3,
},
# beginning of MJD
'0.000' => {
year => 1858,
month => 11,
day => 17,
hour => 0,
minute => 0,
second => 0,
},
# beginning of JD
'-2400000.500' => {
year => -4713,
month => 11,
day => 24,
hour => 12,
minute => 0,
second => 0,
},
);
while ( my ( $mjd, $comps ) = splice @mjd, 0, 2 ) {
my $dt = DateTime->new(
%$comps,
time_zone => 'UTC',
);
is( sprintf( '%.3f', $dt->mjd ), $mjd, "MJD should be $mjd" );
my $jd = sprintf( '%.3f', $mjd + 2_400_000.5 );
is( sprintf( '%.3f', $dt->jd ), $jd, "JD should be $jd" );
}
done_testing();
|
