equationoftime 4.0.0 → 4.1.1

Sign up to get free protection for your applications and to get access to all the features.
Files changed (265) hide show
  1. checksums.yaml +4 -4
  2. data/.gitignore +8 -3
  3. data/Gemfile.lock +48 -0
  4. data/README.md +22 -19
  5. data/Rakefile +14 -7
  6. data/equationoftime.gemspec +9 -6
  7. data/ext/ceot/ceot.c +52 -0
  8. data/ext/ceot/eot.c +79 -0
  9. data/ext/ceot/eot.h +8 -0
  10. data/ext/ceot/extconf.rb +5 -0
  11. data/lib/eot.rb +2 -1
  12. data/lib/eot/angles.rb +44 -43
  13. data/lib/eot/init.rb +13 -6
  14. data/lib/eot/version.rb +1 -1
  15. data/wiki.md +6 -7
  16. metadata +26 -255
  17. data/.settings/org.eclipse.ltk.core.refactoring.prefs +0 -2
  18. data/Makefile +0 -238
  19. data/analemma_data.yml +0 -2925
  20. data/examples/Equation_of_Time.jpg +0 -0
  21. data/examples/analemma_data_generator.rb +0 -53
  22. data/examples/celes_parts/.RUBYARCHDIR.time +0 -0
  23. data/examples/celes_parts/.cproject +0 -56
  24. data/examples/celes_parts/.project +0 -78
  25. data/examples/celes_parts/a2af.c +0 -120
  26. data/examples/celes_parts/a2tf.c +0 -116
  27. data/examples/celes_parts/af2a.c +0 -107
  28. data/examples/celes_parts/anp.c +0 -82
  29. data/examples/celes_parts/anpm.c +0 -82
  30. data/examples/celes_parts/bi00.c +0 -116
  31. data/examples/celes_parts/bp00.c +0 -173
  32. data/examples/celes_parts/bp06.c +0 -136
  33. data/examples/celes_parts/bpn2xy.c +0 -100
  34. data/examples/celes_parts/c2i00a.c +0 -139
  35. data/examples/celes_parts/c2i00b.c +0 -139
  36. data/examples/celes_parts/c2i06a.c +0 -136
  37. data/examples/celes_parts/c2ibpn.c +0 -142
  38. data/examples/celes_parts/c2ixy.c +0 -131
  39. data/examples/celes_parts/c2ixys.c +0 -123
  40. data/examples/celes_parts/c2s.c +0 -96
  41. data/examples/celes_parts/c2t00a.c +0 -154
  42. data/examples/celes_parts/c2t00b.c +0 -150
  43. data/examples/celes_parts/c2t06a.c +0 -152
  44. data/examples/celes_parts/c2tcio.c +0 -122
  45. data/examples/celes_parts/c2teqx.c +0 -122
  46. data/examples/celes_parts/c2tpe.c +0 -167
  47. data/examples/celes_parts/c2txy.c +0 -159
  48. data/examples/celes_parts/cal2jd.c +0 -139
  49. data/examples/celes_parts/celes_core.c +0 -2522
  50. data/examples/celes_parts/cp.c +0 -80
  51. data/examples/celes_parts/cpv.c +0 -82
  52. data/examples/celes_parts/cr.c +0 -83
  53. data/examples/celes_parts/d2dtf.c +0 -206
  54. data/examples/celes_parts/d2tf.c +0 -160
  55. data/examples/celes_parts/dat.c +0 -289
  56. data/examples/celes_parts/depend +0 -25
  57. data/examples/celes_parts/dtdb.c +0 -1213
  58. data/examples/celes_parts/dtf2d.c +0 -196
  59. data/examples/celes_parts/ee00.c +0 -128
  60. data/examples/celes_parts/ee00a.c +0 -135
  61. data/examples/celes_parts/ee00b.c +0 -141
  62. data/examples/celes_parts/ee06a.c +0 -122
  63. data/examples/celes_parts/eect00.c +0 -282
  64. data/examples/celes_parts/eform.c +0 -147
  65. data/examples/celes_parts/eo06a.c +0 -131
  66. data/examples/celes_parts/eors.c +0 -108
  67. data/examples/celes_parts/epb.c +0 -96
  68. data/examples/celes_parts/epb2jd.c +0 -91
  69. data/examples/celes_parts/epj.c +0 -93
  70. data/examples/celes_parts/epj2jd.c +0 -91
  71. data/examples/celes_parts/epv00.c +0 -2589
  72. data/examples/celes_parts/eqeq94.c +0 -131
  73. data/examples/celes_parts/era00.c +0 -136
  74. data/examples/celes_parts/extconf.rb +0 -4
  75. data/examples/celes_parts/fad03.c +0 -103
  76. data/examples/celes_parts/fae03.c +0 -102
  77. data/examples/celes_parts/faf03.c +0 -106
  78. data/examples/celes_parts/faju03.c +0 -102
  79. data/examples/celes_parts/fal03.c +0 -103
  80. data/examples/celes_parts/falp03.c +0 -103
  81. data/examples/celes_parts/fama03.c +0 -102
  82. data/examples/celes_parts/fame03.c +0 -102
  83. data/examples/celes_parts/fane03.c +0 -99
  84. data/examples/celes_parts/faom03.c +0 -104
  85. data/examples/celes_parts/fapa03.c +0 -103
  86. data/examples/celes_parts/fasa03.c +0 -102
  87. data/examples/celes_parts/faur03.c +0 -99
  88. data/examples/celes_parts/fave03.c +0 -102
  89. data/examples/celes_parts/fk52h.c +0 -143
  90. data/examples/celes_parts/fk5hip.c +0 -126
  91. data/examples/celes_parts/fk5hz.c +0 -160
  92. data/examples/celes_parts/fw2m.c +0 -134
  93. data/examples/celes_parts/fw2xy.c +0 -120
  94. data/examples/celes_parts/gc2gd.c +0 -134
  95. data/examples/celes_parts/gc2gde.c +0 -200
  96. data/examples/celes_parts/gd2gc.c +0 -134
  97. data/examples/celes_parts/gd2gce.c +0 -138
  98. data/examples/celes_parts/gmst00.c +0 -145
  99. data/examples/celes_parts/gmst06.c +0 -136
  100. data/examples/celes_parts/gmst82.c +0 -151
  101. data/examples/celes_parts/gst00a.c +0 -138
  102. data/examples/celes_parts/gst00b.c +0 -146
  103. data/examples/celes_parts/gst06.c +0 -140
  104. data/examples/celes_parts/gst06a.c +0 -131
  105. data/examples/celes_parts/gst94.c +0 -131
  106. data/examples/celes_parts/h2fk5.c +0 -148
  107. data/examples/celes_parts/hfk5z.c +0 -175
  108. data/examples/celes_parts/ir.c +0 -83
  109. data/examples/celes_parts/jd2cal.c +0 -155
  110. data/examples/celes_parts/jdcalf.c +0 -161
  111. data/examples/celes_parts/num00a.c +0 -121
  112. data/examples/celes_parts/num00b.c +0 -121
  113. data/examples/celes_parts/num06a.c +0 -125
  114. data/examples/celes_parts/numat.c +0 -109
  115. data/examples/celes_parts/nut00a.c +0 -2047
  116. data/examples/celes_parts/nut00b.c +0 -372
  117. data/examples/celes_parts/nut06a.c +0 -153
  118. data/examples/celes_parts/nut80.c +0 -325
  119. data/examples/celes_parts/nutm80.c +0 -117
  120. data/examples/celes_parts/obl06.c +0 -118
  121. data/examples/celes_parts/obl80.c +0 -118
  122. data/examples/celes_parts/p06e.c +0 -321
  123. data/examples/celes_parts/p2pv.c +0 -83
  124. data/examples/celes_parts/p2s.c +0 -91
  125. data/examples/celes_parts/pap.c +0 -139
  126. data/examples/celes_parts/pas.c +0 -96
  127. data/examples/celes_parts/pb06.c +0 -144
  128. data/examples/celes_parts/pdp.c +0 -84
  129. data/examples/celes_parts/pfw06.c +0 -165
  130. data/examples/celes_parts/plan94.c +0 -514
  131. data/examples/celes_parts/pm.c +0 -83
  132. data/examples/celes_parts/pmat00.c +0 -118
  133. data/examples/celes_parts/pmat06.c +0 -122
  134. data/examples/celes_parts/pmat76.c +0 -141
  135. data/examples/celes_parts/pmp.c +0 -85
  136. data/examples/celes_parts/pn.c +0 -109
  137. data/examples/celes_parts/pn00.c +0 -177
  138. data/examples/celes_parts/pn00a.c +0 -162
  139. data/examples/celes_parts/pn00b.c +0 -162
  140. data/examples/celes_parts/pn06.c +0 -187
  141. data/examples/celes_parts/pn06a.c +0 -152
  142. data/examples/celes_parts/pnm00a.c +0 -121
  143. data/examples/celes_parts/pnm00b.c +0 -121
  144. data/examples/celes_parts/pnm06a.c +0 -124
  145. data/examples/celes_parts/pnm80.c +0 -126
  146. data/examples/celes_parts/pom00.c +0 -115
  147. data/examples/celes_parts/ppp.c +0 -85
  148. data/examples/celes_parts/ppsp.c +0 -94
  149. data/examples/celes_parts/pr00.c +0 -142
  150. data/examples/celes_parts/prec76.c +0 -148
  151. data/examples/celes_parts/pv2p.c +0 -81
  152. data/examples/celes_parts/pv2s.c +0 -144
  153. data/examples/celes_parts/pvdpv.c +0 -102
  154. data/examples/celes_parts/pvm.c +0 -86
  155. data/examples/celes_parts/pvmpv.c +0 -87
  156. data/examples/celes_parts/pvppv.c +0 -87
  157. data/examples/celes_parts/pvstar.c +0 -207
  158. data/examples/celes_parts/pvu.c +0 -93
  159. data/examples/celes_parts/pvup.c +0 -88
  160. data/examples/celes_parts/pvxpv.c +0 -107
  161. data/examples/celes_parts/pxp.c +0 -94
  162. data/examples/celes_parts/rm2v.c +0 -111
  163. data/examples/celes_parts/rv2m.c +0 -118
  164. data/examples/celes_parts/rx.c +0 -110
  165. data/examples/celes_parts/rxp.c +0 -99
  166. data/examples/celes_parts/rxpv.c +0 -86
  167. data/examples/celes_parts/rxr.c +0 -99
  168. data/examples/celes_parts/ry.c +0 -110
  169. data/examples/celes_parts/rz.c +0 -110
  170. data/examples/celes_parts/s00.c +0 -371
  171. data/examples/celes_parts/s00a.c +0 -143
  172. data/examples/celes_parts/s00b.c +0 -143
  173. data/examples/celes_parts/s06.c +0 -368
  174. data/examples/celes_parts/s06a.c +0 -145
  175. data/examples/celes_parts/s2c.c +0 -85
  176. data/examples/celes_parts/s2p.c +0 -88
  177. data/examples/celes_parts/s2pv.c +0 -103
  178. data/examples/celes_parts/s2xpv.c +0 -87
  179. data/examples/celes_parts/sepp.c +0 -105
  180. data/examples/celes_parts/seps.c +0 -93
  181. data/examples/celes_parts/sofa.h +0 -379
  182. data/examples/celes_parts/sofam.h +0 -155
  183. data/examples/celes_parts/sp00.c +0 -118
  184. data/examples/celes_parts/starpm.c +0 -205
  185. data/examples/celes_parts/starpv.c +0 -264
  186. data/examples/celes_parts/sxp.c +0 -84
  187. data/examples/celes_parts/sxpv.c +0 -85
  188. data/examples/celes_parts/taitt.c +0 -110
  189. data/examples/celes_parts/taiut1.c +0 -112
  190. data/examples/celes_parts/taiutc.c +0 -182
  191. data/examples/celes_parts/tcbtdb.c +0 -132
  192. data/examples/celes_parts/tcgtt.c +0 -109
  193. data/examples/celes_parts/tdbtcb.c +0 -137
  194. data/examples/celes_parts/tdbtt.c +0 -122
  195. data/examples/celes_parts/test_celes.rb +0 -48
  196. data/examples/celes_parts/tf2a.c +0 -107
  197. data/examples/celes_parts/tf2d.c +0 -107
  198. data/examples/celes_parts/tr.c +0 -93
  199. data/examples/celes_parts/trxp.c +0 -93
  200. data/examples/celes_parts/trxpv.c +0 -93
  201. data/examples/celes_parts/tttai.c +0 -110
  202. data/examples/celes_parts/tttcg.c +0 -112
  203. data/examples/celes_parts/tttdb.c +0 -121
  204. data/examples/celes_parts/ttut1.c +0 -110
  205. data/examples/celes_parts/ut1tai.c +0 -111
  206. data/examples/celes_parts/ut1tt.c +0 -110
  207. data/examples/celes_parts/ut1utc.c +0 -193
  208. data/examples/celes_parts/utctai.c +0 -163
  209. data/examples/celes_parts/utcut1.c +0 -151
  210. data/examples/celes_parts/xy06.c +0 -2758
  211. data/examples/celes_parts/xys00a.c +0 -133
  212. data/examples/celes_parts/xys00b.c +0 -133
  213. data/examples/celes_parts/xys06a.c +0 -133
  214. data/examples/celes_parts/zp.c +0 -77
  215. data/examples/celes_parts/zpv.c +0 -79
  216. data/examples/celes_parts/zr.c +0 -83
  217. data/examples/check_date_type.rb +0 -55
  218. data/examples/compare_geoc_long_ra.rb +0 -39
  219. data/examples/data_table_for_astro_dog.rb +0 -40
  220. data/examples/earth_rotation.rb +0 -38
  221. data/examples/eot.c +0 -15
  222. data/examples/eot.h +0 -2
  223. data/examples/eot_methods_list.rb +0 -43
  224. data/examples/eot_plot.r +0 -57
  225. data/examples/eot_suntimes.rb +0 -144
  226. data/examples/equation_of_time.py +0 -186
  227. data/examples/extconf.rb +0 -2
  228. data/examples/figure_1.jpg +0 -0
  229. data/examples/file_converter.rb +0 -31
  230. data/examples/from_readme.rb +0 -10
  231. data/examples/geo_locator.rb +0 -12
  232. data/examples/getjd.rb +0 -45
  233. data/examples/gmst_gast_non_sofa.rb +0 -404
  234. data/examples/input_suntimes.rb +0 -21
  235. data/examples/julian_day_formula.rb +0 -29
  236. data/examples/julian_day_formula.txt +0 -12
  237. data/examples/my_time_conversion.rb +0 -21
  238. data/examples/nutation_series.txt +0 -678
  239. data/examples/nutation_series.yaml +0 -14239
  240. data/examples/nutation_table5_3a.txt +0 -682
  241. data/examples/nutation_table5_3a.yaml +0 -9532
  242. data/examples/ptime.rb +0 -162
  243. data/examples/rbeot.c +0 -24
  244. data/examples/suntimes.rb +0 -28
  245. data/examples/suntimes_test.rb +0 -47
  246. data/examples/test_eot.rb +0 -17
  247. data/examples/test_poly_eval.rb +0 -27
  248. data/examples/time_scales.rb +0 -26
  249. data/examples/times_year.rb +0 -50
  250. data/examples/usage_example.rb +0 -23
  251. data/examples/use_angles.rb +0 -219
  252. data/rise_set_data.yml +0 -1461
  253. data/tests/minitest/aliased_angles_spec.rb +0 -241
  254. data/tests/minitest/aliased_displays_spec.rb +0 -108
  255. data/tests/minitest/aliased_utilities_spec.rb +0 -38
  256. data/tests/minitest/angles_spec.rb +0 -266
  257. data/tests/minitest/constants_spec.rb +0 -21
  258. data/tests/minitest/delta_epsilon_spec.rb +0 -32
  259. data/tests/minitest/displays_spec.rb +0 -113
  260. data/tests/minitest/geo_spec.rb +0 -40
  261. data/tests/minitest/init_spec.rb +0 -46
  262. data/tests/minitest/nutation_spec.rb +0 -39
  263. data/tests/minitest/times_spec.rb +0 -139
  264. data/tests/minitest/utilities_spec.rb +0 -37
  265. data/tests/spec_config.rb +0 -3
@@ -1,80 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauCp(double p[3], double c[3])
4
- /*
5
- ** - - - - - -
6
- ** i a u C p
7
- ** - - - - - -
8
- **
9
- ** Copy a p-vector.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Given:
14
- ** p double[3] p-vector to be copied
15
- **
16
- ** Returned:
17
- ** c double[3] copy
18
- **
19
- ** This revision: 2008 May 11
20
- **
21
- ** Original version 2012-03-01
22
- **
23
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
24
- */
25
- {
26
- c[0] = p[0];
27
- c[1] = p[1];
28
- c[2] = p[2];
29
-
30
- return;
31
-
32
- /*----------------------------------------------------------------------
33
- **
34
- ** Celes is a wrapper of the SOFA Library for Ruby.
35
- **
36
- ** This file is redistributed and relicensed in accordance with
37
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
38
- **
39
- ** The original library is available from IAU Standards of
40
- ** Fundamental Astronomy (http://www.iausofa.org/).
41
- **
42
- **
43
- **
44
- **
45
- **
46
- ** Copyright (C) 2013, Naoki Arita
47
- ** All rights reserved.
48
- **
49
- ** Redistribution and use in source and binary forms, with or without
50
- ** modification, are permitted provided that the following conditions
51
- ** are met:
52
- **
53
- ** 1 Redistributions of source code must retain the above copyright
54
- ** notice, this list of conditions and the following disclaimer.
55
- **
56
- ** 2 Redistributions in binary form must reproduce the above copyright
57
- ** notice, this list of conditions and the following disclaimer in
58
- ** the documentation and/or other materials provided with the
59
- ** distribution.
60
- **
61
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
62
- ** the International Astronomical Union nor the names of its
63
- ** contributors may be used to endorse or promote products derived
64
- ** from this software without specific prior written permission.
65
- **
66
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
67
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
68
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
69
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
70
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
71
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
72
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
73
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
74
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
75
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
76
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
77
- ** POSSIBILITY OF SUCH DAMAGE.
78
- **
79
- **--------------------------------------------------------------------*/
80
- }
@@ -1,82 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauCpv(double pv[2][3], double c[2][3])
4
- /*
5
- ** - - - - - - -
6
- ** i a u C p v
7
- ** - - - - - - -
8
- **
9
- ** Copy a position/velocity vector.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Given:
14
- ** pv double[2][3] position/velocity vector to be copied
15
- **
16
- ** Returned:
17
- ** c double[2][3] copy
18
- **
19
- ** Called:
20
- ** iauCp copy p-vector
21
- **
22
- ** This revision: 2008 May 11
23
- **
24
- ** Original version 2012-03-01
25
- **
26
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
27
- */
28
- {
29
- iauCp(pv[0], c[0]);
30
- iauCp(pv[1], c[1]);
31
-
32
- return;
33
-
34
- /*----------------------------------------------------------------------
35
- **
36
- ** Celes is a wrapper of the SOFA Library for Ruby.
37
- **
38
- ** This file is redistributed and relicensed in accordance with
39
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
40
- **
41
- ** The original library is available from IAU Standards of
42
- ** Fundamental Astronomy (http://www.iausofa.org/).
43
- **
44
- **
45
- **
46
- **
47
- **
48
- ** Copyright (C) 2013, Naoki Arita
49
- ** All rights reserved.
50
- **
51
- ** Redistribution and use in source and binary forms, with or without
52
- ** modification, are permitted provided that the following conditions
53
- ** are met:
54
- **
55
- ** 1 Redistributions of source code must retain the above copyright
56
- ** notice, this list of conditions and the following disclaimer.
57
- **
58
- ** 2 Redistributions in binary form must reproduce the above copyright
59
- ** notice, this list of conditions and the following disclaimer in
60
- ** the documentation and/or other materials provided with the
61
- ** distribution.
62
- **
63
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
64
- ** the International Astronomical Union nor the names of its
65
- ** contributors may be used to endorse or promote products derived
66
- ** from this software without specific prior written permission.
67
- **
68
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
69
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
70
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
71
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
72
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
73
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
74
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
75
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
76
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
77
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
78
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
79
- ** POSSIBILITY OF SUCH DAMAGE.
80
- **
81
- **--------------------------------------------------------------------*/
82
- }
@@ -1,83 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauCr(double r[3][3], double c[3][3])
4
- /*
5
- ** - - - - - -
6
- ** i a u C r
7
- ** - - - - - -
8
- **
9
- ** Copy an r-matrix.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Given:
14
- ** r double[3][3] r-matrix to be copied
15
- **
16
- ** Returned:
17
- ** char[] double[3][3] copy
18
- **
19
- ** Called:
20
- ** iauCp copy p-vector
21
- **
22
- ** This revision: 2008 May 11
23
- **
24
- ** Original version 2012-03-01
25
- **
26
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
27
- */
28
- {
29
- iauCp(r[0], c[0]);
30
- iauCp(r[1], c[1]);
31
- iauCp(r[2], c[2]);
32
-
33
- return;
34
-
35
- /*----------------------------------------------------------------------
36
- **
37
- ** Celes is a wrapper of the SOFA Library for Ruby.
38
- **
39
- ** This file is redistributed and relicensed in accordance with
40
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
41
- **
42
- ** The original library is available from IAU Standards of
43
- ** Fundamental Astronomy (http://www.iausofa.org/).
44
- **
45
- **
46
- **
47
- **
48
- **
49
- ** Copyright (C) 2013, Naoki Arita
50
- ** All rights reserved.
51
- **
52
- ** Redistribution and use in source and binary forms, with or without
53
- ** modification, are permitted provided that the following conditions
54
- ** are met:
55
- **
56
- ** 1 Redistributions of source code must retain the above copyright
57
- ** notice, this list of conditions and the following disclaimer.
58
- **
59
- ** 2 Redistributions in binary form must reproduce the above copyright
60
- ** notice, this list of conditions and the following disclaimer in
61
- ** the documentation and/or other materials provided with the
62
- ** distribution.
63
- **
64
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
65
- ** the International Astronomical Union nor the names of its
66
- ** contributors may be used to endorse or promote products derived
67
- ** from this software without specific prior written permission.
68
- **
69
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
70
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
71
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
72
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
73
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
74
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
75
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
76
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
77
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
78
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
79
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
80
- ** POSSIBILITY OF SUCH DAMAGE.
81
- **
82
- **--------------------------------------------------------------------*/
83
- }
@@ -1,206 +0,0 @@
1
- #include "sofam.h"
2
- #include <string.h>
3
-
4
- int iauD2dtf(const char *scale, int ndp, double d1, double d2,
5
- int *iy, int *im, int *id, int ihmsf[4])
6
- /*
7
- ** - - - - - - - - -
8
- ** i a u D 2 d t f
9
- ** - - - - - - - - -
10
- **
11
- ** Format for output a 2-part Julian Date (or in the case of UTC a
12
- ** quasi-JD form that includes special provision for leap seconds).
13
- **
14
- ** Status: support function.
15
- **
16
- ** Given:
17
- ** scale char[] time scale ID (Note 1)
18
- ** ndp int resolution (Note 2)
19
- ** d1,d2 double time as a 2-part Julian Date (Notes 3,4)
20
- **
21
- ** Returned:
22
- ** iy,im,id int year, month, day in Gregorian calendar (Note 5)
23
- ** ihmsf int[4] hours, minutes, seconds, fraction (Note 1)
24
- **
25
- ** Returned (function value):
26
- ** int status: +1 = dubious year (Note 5)
27
- ** 0 = OK
28
- ** -1 = unacceptable date (Note 6)
29
- **
30
- ** Notes:
31
- **
32
- ** 1) scale identifies the time scale. Only the value "UTC" (in upper
33
- ** case) is significant, and enables handling of leap seconds (see
34
- ** Note 4).
35
- **
36
- ** 2) ndp is the number of decimal places in the seconds field, and can
37
- ** have negative as well as positive values, such as:
38
- **
39
- ** ndp resolution
40
- ** -4 1 00 00
41
- ** -3 0 10 00
42
- ** -2 0 01 00
43
- ** -1 0 00 10
44
- ** 0 0 00 01
45
- ** 1 0 00 00.1
46
- ** 2 0 00 00.01
47
- ** 3 0 00 00.001
48
- **
49
- ** The limits are platform dependent, but a safe range is -5 to +9.
50
- **
51
- ** 3) d1+d2 is Julian Date, apportioned in any convenient way between
52
- ** the two arguments, for example where d1 is the Julian Day Number
53
- ** and d2 is the fraction of a day. In the case of UTC, where the
54
- ** use of JD is problematical, special conventions apply: see the
55
- ** next note.
56
- **
57
- ** 4) JD cannot unambiguously represent UTC during a leap second unless
58
- ** special measures are taken. The SOFA internal convention is that
59
- ** the quasi-JD day represents UTC days whether the length is 86399,
60
- ** 86400 or 86401 SI seconds.
61
- **
62
- ** 5) The warning status "dubious year" flags UTCs that predate the
63
- ** introduction of the time scale and that are too far in the future
64
- ** to be trusted. See iauDat for further details.
65
- **
66
- ** 6) For calendar conventions and limitations, see iauCal2jd.
67
- **
68
- ** Called:
69
- ** iauJd2cal JD to Gregorian calendar
70
- ** iauD2tf decompose days to hms
71
- ** iauDat delta(AT) = TAI-UTC
72
- **
73
- ** This revision: 2012 February 12
74
- **
75
- ** Original version 2012-03-01
76
- **
77
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
78
- */
79
- {
80
- int leap;
81
- char s;
82
- int iy1, im1, id1, js, iy2, im2, id2, ihmsf1[4], i;
83
- double a1, b1, fd, dat1, w, dat2, ddt;
84
-
85
-
86
- /* The two-part JD. */
87
- a1 = d1;
88
- b1 = d2;
89
-
90
- /* Provisional calendar date. */
91
- js = iauJd2cal(a1, b1, &iy1, &im1, &id1, &fd);
92
- if ( js ) return js < 0 ? -1 : js;
93
-
94
- /* Is this a leap second day? */
95
- leap = 0;
96
- if ( ! strcmp(scale,"UTC") ) {
97
-
98
- /* TAI-UTC today. */
99
- js = iauDat(iy1, im1, id1, fd, &dat1);
100
- if ( js < 0 ) return -1;
101
-
102
- /* TAI-UTC tomorrow (at noon, to avoid rounding effects). */
103
- js = iauJd2cal(a1+1.5, b1-fd, &iy2, &im2, &id2, &w);
104
- js = iauDat(iy2, im2, id2, 0.0, &dat2);
105
- if ( js < 0 ) return -1;
106
-
107
- /* The change in TAI-UTC (seconds). */
108
- ddt = dat2 - dat1;
109
-
110
- /* If leap second day, scale the fraction of a day into SI. */
111
- leap = fabs(ddt) > 0.5;
112
- if (leap) fd += fd * ddt/DAYSEC;
113
- }
114
-
115
- /* Provisional time of day. */
116
- iauD2tf ( ndp, fd, &s, ihmsf1 );
117
-
118
- /* Is this a leap second day? */
119
- if ( ! leap ) {
120
-
121
- /* No. Has the time rounded up to 24h? */
122
- if ( ihmsf1[0] > 23 ) {
123
-
124
- /* Yes. We will need tomorrow's calendar date. */
125
- js = iauJd2cal(a1+1.5, b1-fd, &iy2, &im2, &id2, &w);
126
-
127
- /* Use 0h tomorrow. */
128
- iy1 = iy2;
129
- im1 = im2;
130
- id1 = id2;
131
- for ( i = 0; i < 4; i++ ) {
132
- ihmsf1[i] = 0;
133
- }
134
- }
135
- } else {
136
-
137
- /* This is a leap second day. Has the time reached or passed 24h? */
138
- if ( ihmsf1[0] > 23 ) {
139
-
140
- /* Yes. Use 23 59 60... */
141
- ihmsf1[0] = 23;
142
- ihmsf1[1] = 59;
143
- ihmsf1[2] = 60;
144
- }
145
- }
146
-
147
- /* Results. */
148
- *iy = iy1;
149
- *im = im1;
150
- *id = id1;
151
- for ( i = 0; i < 4; i++ ) {
152
- ihmsf[i] = ihmsf1[i];
153
- }
154
-
155
- /* Status. */
156
- return js < 0 ? -1 : js;
157
-
158
- /*----------------------------------------------------------------------
159
- **
160
- ** Celes is a wrapper of the SOFA Library for Ruby.
161
- **
162
- ** This file is redistributed and relicensed in accordance with
163
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
164
- **
165
- ** The original library is available from IAU Standards of
166
- ** Fundamental Astronomy (http://www.iausofa.org/).
167
- **
168
- **
169
- **
170
- **
171
- **
172
- ** Copyright (C) 2013, Naoki Arita
173
- ** All rights reserved.
174
- **
175
- ** Redistribution and use in source and binary forms, with or without
176
- ** modification, are permitted provided that the following conditions
177
- ** are met:
178
- **
179
- ** 1 Redistributions of source code must retain the above copyright
180
- ** notice, this list of conditions and the following disclaimer.
181
- **
182
- ** 2 Redistributions in binary form must reproduce the above copyright
183
- ** notice, this list of conditions and the following disclaimer in
184
- ** the documentation and/or other materials provided with the
185
- ** distribution.
186
- **
187
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
188
- ** the International Astronomical Union nor the names of its
189
- ** contributors may be used to endorse or promote products derived
190
- ** from this software without specific prior written permission.
191
- **
192
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
193
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
194
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
195
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
196
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
197
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
198
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
199
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
200
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
201
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
202
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
203
- ** POSSIBILITY OF SUCH DAMAGE.
204
- **
205
- **--------------------------------------------------------------------*/
206
- }