equationoftime 4.0.0 → 4.1.1

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
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,115 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPom00(double xp, double yp, double sp, double rpom[3][3])
4
- /*
5
- ** - - - - - - - - - -
6
- ** i a u P o m 0 0
7
- ** - - - - - - - - - -
8
- **
9
- ** Form the matrix of polar motion for a given date, IAU 2000.
10
- **
11
- ** Status: support function.
12
- **
13
- ** Given:
14
- ** xp,yp double coordinates of the pole (radians, Note 1)
15
- ** sp double the TIO locator s' (radians, Note 2)
16
- **
17
- ** Returned:
18
- ** rpom double[3][3] polar-motion matrix (Note 3)
19
- **
20
- ** Notes:
21
- **
22
- ** 1) The arguments xp and yp are the coordinates (in radians) of the
23
- ** Celestial Intermediate Pole with respect to the International
24
- ** Terrestrial Reference System (see IERS Conventions 2003),
25
- ** measured along the meridians to 0 and 90 deg west respectively.
26
- **
27
- ** 2) The argument sp is the TIO locator s', in radians, which
28
- ** positions the Terrestrial Intermediate Origin on the equator. It
29
- ** is obtained from polar motion observations by numerical
30
- ** integration, and so is in essence unpredictable. However, it is
31
- ** dominated by a secular drift of about 47 microarcseconds per
32
- ** century, and so can be taken into account by using s' = -47*t,
33
- ** where t is centuries since J2000.0. The function iauSp00
34
- ** implements this approximation.
35
- **
36
- ** 3) The matrix operates in the sense V(TRS) = rpom * V(CIP), meaning
37
- ** that it is the final rotation when computing the pointing
38
- ** direction to a celestial source.
39
- **
40
- ** Called:
41
- ** iauIr initialize r-matrix to identity
42
- ** iauRz rotate around Z-axis
43
- ** iauRy rotate around Y-axis
44
- ** iauRx rotate around X-axis
45
- **
46
- ** Reference:
47
- **
48
- ** McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
49
- ** IERS Technical Note No. 32, BKG (2004)
50
- **
51
- ** This revision: 2009 December 17
52
- **
53
- ** Original version 2012-03-01
54
- **
55
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
56
- */
57
- {
58
-
59
- /* Construct the matrix. */
60
- iauIr(rpom);
61
- iauRz(sp, rpom);
62
- iauRy(-xp, rpom);
63
- iauRx(-yp, rpom);
64
-
65
- return;
66
-
67
- /*----------------------------------------------------------------------
68
- **
69
- ** Celes is a wrapper of the SOFA Library for Ruby.
70
- **
71
- ** This file is redistributed and relicensed in accordance with
72
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
73
- **
74
- ** The original library is available from IAU Standards of
75
- ** Fundamental Astronomy (http://www.iausofa.org/).
76
- **
77
- **
78
- **
79
- **
80
- **
81
- ** Copyright (C) 2013, Naoki Arita
82
- ** All rights reserved.
83
- **
84
- ** Redistribution and use in source and binary forms, with or without
85
- ** modification, are permitted provided that the following conditions
86
- ** are met:
87
- **
88
- ** 1 Redistributions of source code must retain the above copyright
89
- ** notice, this list of conditions and the following disclaimer.
90
- **
91
- ** 2 Redistributions in binary form must reproduce the above copyright
92
- ** notice, this list of conditions and the following disclaimer in
93
- ** the documentation and/or other materials provided with the
94
- ** distribution.
95
- **
96
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
97
- ** the International Astronomical Union nor the names of its
98
- ** contributors may be used to endorse or promote products derived
99
- ** from this software without specific prior written permission.
100
- **
101
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
102
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
103
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
104
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
105
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
106
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
107
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
108
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
109
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
110
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
111
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
112
- ** POSSIBILITY OF SUCH DAMAGE.
113
- **
114
- **--------------------------------------------------------------------*/
115
- }
@@ -1,85 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPpp(double a[3], double b[3], double apb[3])
4
- /*
5
- ** - - - - - - -
6
- ** i a u P p p
7
- ** - - - - - - -
8
- **
9
- ** P-vector addition.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Given:
14
- ** a double[3] first p-vector
15
- ** b double[3] second p-vector
16
- **
17
- ** Returned:
18
- ** apb double[3] a + b
19
- **
20
- ** Note:
21
- ** It is permissible to re-use the same array for any of the
22
- ** arguments.
23
- **
24
- ** This revision: 2008 November 18
25
- **
26
- ** Original version 2012-03-01
27
- **
28
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
29
- */
30
- {
31
- apb[0] = a[0] + b[0];
32
- apb[1] = a[1] + b[1];
33
- apb[2] = a[2] + b[2];
34
-
35
- return;
36
-
37
- /*----------------------------------------------------------------------
38
- **
39
- ** Celes is a wrapper of the SOFA Library for Ruby.
40
- **
41
- ** This file is redistributed and relicensed in accordance with
42
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
43
- **
44
- ** The original library is available from IAU Standards of
45
- ** Fundamental Astronomy (http://www.iausofa.org/).
46
- **
47
- **
48
- **
49
- **
50
- **
51
- ** Copyright (C) 2013, Naoki Arita
52
- ** All rights reserved.
53
- **
54
- ** Redistribution and use in source and binary forms, with or without
55
- ** modification, are permitted provided that the following conditions
56
- ** are met:
57
- **
58
- ** 1 Redistributions of source code must retain the above copyright
59
- ** notice, this list of conditions and the following disclaimer.
60
- **
61
- ** 2 Redistributions in binary form must reproduce the above copyright
62
- ** notice, this list of conditions and the following disclaimer in
63
- ** the documentation and/or other materials provided with the
64
- ** distribution.
65
- **
66
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
67
- ** the International Astronomical Union nor the names of its
68
- ** contributors may be used to endorse or promote products derived
69
- ** from this software without specific prior written permission.
70
- **
71
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
72
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
73
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
74
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
75
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
76
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
77
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
78
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
79
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
80
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
81
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
82
- ** POSSIBILITY OF SUCH DAMAGE.
83
- **
84
- **--------------------------------------------------------------------*/
85
- }
@@ -1,94 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPpsp(double a[3], double s, double b[3], double apsb[3])
4
- /*
5
- ** - - - - - - - -
6
- ** i a u P p s p
7
- ** - - - - - - - -
8
- **
9
- ** P-vector plus scaled p-vector.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Given:
14
- ** a double[3] first p-vector
15
- ** s double scalar (multiplier for b)
16
- ** b double[3] second p-vector
17
- **
18
- ** Returned:
19
- ** apsb double[3] a + s*b
20
- **
21
- ** Note:
22
- ** It is permissible for any of a, b and apsb to be the same array.
23
- **
24
- ** Called:
25
- ** iauSxp multiply p-vector by scalar
26
- ** iauPpp p-vector plus p-vector
27
- **
28
- ** This revision: 2008 November 18
29
- **
30
- ** Original version 2012-03-01
31
- **
32
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
33
- */
34
- {
35
- double sb[3];
36
-
37
-
38
- /* s*b. */
39
- iauSxp(s, b, sb);
40
-
41
- /* a + s*b. */
42
- iauPpp(a, sb, apsb);
43
-
44
- return;
45
-
46
- /*----------------------------------------------------------------------
47
- **
48
- ** Celes is a wrapper of the SOFA Library for Ruby.
49
- **
50
- ** This file is redistributed and relicensed in accordance with
51
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
52
- **
53
- ** The original library is available from IAU Standards of
54
- ** Fundamental Astronomy (http://www.iausofa.org/).
55
- **
56
- **
57
- **
58
- **
59
- **
60
- ** Copyright (C) 2013, Naoki Arita
61
- ** All rights reserved.
62
- **
63
- ** Redistribution and use in source and binary forms, with or without
64
- ** modification, are permitted provided that the following conditions
65
- ** are met:
66
- **
67
- ** 1 Redistributions of source code must retain the above copyright
68
- ** notice, this list of conditions and the following disclaimer.
69
- **
70
- ** 2 Redistributions in binary form must reproduce the above copyright
71
- ** notice, this list of conditions and the following disclaimer in
72
- ** the documentation and/or other materials provided with the
73
- ** distribution.
74
- **
75
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
76
- ** the International Astronomical Union nor the names of its
77
- ** contributors may be used to endorse or promote products derived
78
- ** from this software without specific prior written permission.
79
- **
80
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
81
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
82
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
83
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
84
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
85
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
86
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
87
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
88
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
89
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
90
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
91
- ** POSSIBILITY OF SUCH DAMAGE.
92
- **
93
- **--------------------------------------------------------------------*/
94
- }
@@ -1,142 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPr00(double date1, double date2, double *dpsipr, double *depspr)
4
- /*
5
- ** - - - - - - - -
6
- ** i a u P r 0 0
7
- ** - - - - - - - -
8
- **
9
- ** Precession-rate part of the IAU 2000 precession-nutation models
10
- ** (part of MHB2000).
11
- **
12
- ** Status: canonical model.
13
- **
14
- ** Given:
15
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
16
- **
17
- ** Returned:
18
- ** dpsipr,depspr double precession corrections (Notes 2,3)
19
- **
20
- ** Notes:
21
- **
22
- ** 1) The TT date date1+date2 is a Julian Date, apportioned in any
23
- ** convenient way between the two arguments. For example,
24
- ** JD(TT)=2450123.7 could be expressed in any of these ways,
25
- ** among others:
26
- **
27
- ** date1 date2
28
- **
29
- ** 2450123.7 0.0 (JD method)
30
- ** 2451545.0 -1421.3 (J2000 method)
31
- ** 2400000.5 50123.2 (MJD method)
32
- ** 2450123.5 0.2 (date & time method)
33
- **
34
- ** The JD method is the most natural and convenient to use in
35
- ** cases where the loss of several decimal digits of resolution
36
- ** is acceptable. The J2000 method is best matched to the way
37
- ** the argument is handled internally and will deliver the
38
- ** optimum resolution. The MJD method and the date & time methods
39
- ** are both good compromises between resolution and convenience.
40
- **
41
- ** 2) The precession adjustments are expressed as "nutation
42
- ** components", corrections in longitude and obliquity with respect
43
- ** to the J2000.0 equinox and ecliptic.
44
- **
45
- ** 3) Although the precession adjustments are stated to be with respect
46
- ** to Lieske et al. (1977), the MHB2000 model does not specify which
47
- ** set of Euler angles are to be used and how the adjustments are to
48
- ** be applied. The most literal and straightforward procedure is to
49
- ** adopt the 4-rotation epsilon_0, psi_A, omega_A, xi_A option, and
50
- ** to add dpsipr to psi_A and depspr to both omega_A and eps_A.
51
- **
52
- ** 4) This is an implementation of one aspect of the IAU 2000A nutation
53
- ** model, formally adopted by the IAU General Assembly in 2000,
54
- ** namely MHB2000 (Mathews et al. 2002).
55
- **
56
- ** References:
57
- **
58
- ** Lieske, J.H., Lederle, T., Fricke, W. & Morando, B., "Expressions
59
- ** for the precession quantities based upon the IAU (1976) System of
60
- ** Astronomical Constants", Astron.Astrophys., 58, 1-16 (1977)
61
- **
62
- ** Mathews, P.M., Herring, T.A., Buffet, B.A., "Modeling of nutation
63
- ** and precession New nutation series for nonrigid Earth and
64
- ** insights into the Earth's interior", J.Geophys.Res., 107, B4,
65
- ** 2002. The MHB2000 code itself was obtained on 9th September 2002
66
- ** from ftp://maia.usno.navy.mil/conv2000/chapter5/IAU2000A.
67
- **
68
- ** Wallace, P.T., "Software for Implementing the IAU 2000
69
- ** Resolutions", in IERS Workshop 5.1 (2002).
70
- **
71
- ** This revision: 2009 December 17
72
- **
73
- ** Original version 2012-03-01
74
- **
75
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
76
- */
77
- {
78
- double t;
79
-
80
- /* Precession and obliquity corrections (radians per century) */
81
- static const double PRECOR = -0.29965 * DAS2R,
82
- OBLCOR = -0.02524 * DAS2R;
83
-
84
-
85
- /* Interval between fundamental epoch J2000.0 and given date (JC). */
86
- t = ((date1 - DJ00) + date2) / DJC;
87
-
88
- /* Precession rate contributions with respect to IAU 1976/80. */
89
- *dpsipr = PRECOR * t;
90
- *depspr = OBLCOR * t;
91
-
92
- return;
93
-
94
- /*----------------------------------------------------------------------
95
- **
96
- ** Celes is a wrapper of the SOFA Library for Ruby.
97
- **
98
- ** This file is redistributed and relicensed in accordance with
99
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
100
- **
101
- ** The original library is available from IAU Standards of
102
- ** Fundamental Astronomy (http://www.iausofa.org/).
103
- **
104
- **
105
- **
106
- **
107
- **
108
- ** Copyright (C) 2013, Naoki Arita
109
- ** All rights reserved.
110
- **
111
- ** Redistribution and use in source and binary forms, with or without
112
- ** modification, are permitted provided that the following conditions
113
- ** are met:
114
- **
115
- ** 1 Redistributions of source code must retain the above copyright
116
- ** notice, this list of conditions and the following disclaimer.
117
- **
118
- ** 2 Redistributions in binary form must reproduce the above copyright
119
- ** notice, this list of conditions and the following disclaimer in
120
- ** the documentation and/or other materials provided with the
121
- ** distribution.
122
- **
123
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
124
- ** the International Astronomical Union nor the names of its
125
- ** contributors may be used to endorse or promote products derived
126
- ** from this software without specific prior written permission.
127
- **
128
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
129
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
130
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
131
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
132
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
133
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
134
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
135
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
136
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
137
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
138
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
139
- ** POSSIBILITY OF SUCH DAMAGE.
140
- **
141
- **--------------------------------------------------------------------*/
142
- }
@@ -1,148 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPrec76(double ep01, double ep02, double ep11, double ep12,
4
- double *zeta, double *z, double *theta)
5
- /*
6
- ** - - - - - - - - - -
7
- ** i a u P r e c 7 6
8
- ** - - - - - - - - - -
9
- **
10
- ** IAU 1976 precession model.
11
- **
12
- ** This function forms the three Euler angles which implement general
13
- ** precession between two epochs, using the IAU 1976 model (as for
14
- ** the FK5 catalog).
15
- **
16
- ** Status: canonical model.
17
- **
18
- ** Given:
19
- ** ep01,ep02 double TDB starting epoch (Note 1)
20
- ** ep11,ep12 double TDB ending epoch (Note 1)
21
- **
22
- ** Returned:
23
- ** zeta double 1st rotation: radians cw around z
24
- ** z double 3rd rotation: radians cw around z
25
- ** theta double 2nd rotation: radians ccw around y
26
- **
27
- ** Notes:
28
- **
29
- ** 1) The epochs ep01+ep02 and ep11+ep12 are Julian Dates, apportioned
30
- ** in any convenient way between the arguments epn1 and epn2. For
31
- ** example, JD(TDB)=2450123.7 could be expressed in any of these
32
- ** ways, among others:
33
- **
34
- ** epn1 epn2
35
- **
36
- ** 2450123.7 0.0 (JD method)
37
- ** 2451545.0 -1421.3 (J2000 method)
38
- ** 2400000.5 50123.2 (MJD method)
39
- ** 2450123.5 0.2 (date & time method)
40
- **
41
- ** The JD method is the most natural and convenient to use in cases
42
- ** where the loss of several decimal digits of resolution is
43
- ** acceptable. The J2000 method is best matched to the way the
44
- ** argument is handled internally and will deliver the optimum
45
- ** optimum resolution. The MJD method and the date & time methods
46
- ** are both good compromises between resolution and convenience.
47
- ** The two epochs may be expressed using different methods, but at
48
- ** the risk of losing some resolution.
49
- **
50
- ** 2) The accumulated precession angles zeta, z, theta are expressed
51
- ** through canonical polynomials which are valid only for a limited
52
- ** time span. In addition, the IAU 1976 precession rate is known to
53
- ** be imperfect. The absolute accuracy of the present formulation
54
- ** is better than 0.1 arcsec from 1960AD to 2040AD, better than
55
- ** 1 arcsec from 1640AD to 2360AD, and remains below 3 arcsec for
56
- ** the whole of the period 500BC to 3000AD. The errors exceed
57
- ** 10 arcsec outside the range 1200BC to 3900AD, exceed 100 arcsec
58
- ** outside 4200BC to 5600AD and exceed 1000 arcsec outside 6800BC to
59
- ** 8200AD.
60
- **
61
- ** 3) The three angles are returned in the conventional order, which
62
- ** is not the same as the order of the corresponding Euler
63
- ** rotations. The precession matrix is
64
- ** R_3(-z) x R_2(+theta) x R_3(-zeta).
65
- **
66
- ** Reference:
67
- **
68
- ** Lieske, J.H., 1979, Astron.Astrophys. 73, 282, equations
69
- ** (6) & (7), p283.
70
- **
71
- ** This revision: 2009 December 17
72
- **
73
- ** Original version 2012-03-01
74
- **
75
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
76
- */
77
- {
78
- double t0, t, tas2r, w;
79
-
80
-
81
- /* Interval between fundamental epoch J2000.0 and start epoch (JC). */
82
- t0 = ((ep01 - DJ00) + ep02) / DJC;
83
-
84
- /* Interval over which precession required (JC). */
85
- t = ((ep11 - ep01) + (ep12 - ep02)) / DJC;
86
-
87
- /* Euler angles. */
88
- tas2r = t * DAS2R;
89
- w = 2306.2181 + (1.39656 - 0.000139 * t0) * t0;
90
-
91
- *zeta = (w + ((0.30188 - 0.000344 * t0) + 0.017998 * t) * t) * tas2r;
92
-
93
- *z = (w + ((1.09468 + 0.000066 * t0) + 0.018203 * t) * t) * tas2r;
94
-
95
- *theta = ((2004.3109 + (-0.85330 - 0.000217 * t0) * t0)
96
- + ((-0.42665 - 0.000217 * t0) - 0.041833 * t) * t) * tas2r;
97
-
98
- return;
99
-
100
- /*----------------------------------------------------------------------
101
- **
102
- ** Celes is a wrapper of the SOFA Library for Ruby.
103
- **
104
- ** This file is redistributed and relicensed in accordance with
105
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
106
- **
107
- ** The original library is available from IAU Standards of
108
- ** Fundamental Astronomy (http://www.iausofa.org/).
109
- **
110
- **
111
- **
112
- **
113
- **
114
- ** Copyright (C) 2013, Naoki Arita
115
- ** All rights reserved.
116
- **
117
- ** Redistribution and use in source and binary forms, with or without
118
- ** modification, are permitted provided that the following conditions
119
- ** are met:
120
- **
121
- ** 1 Redistributions of source code must retain the above copyright
122
- ** notice, this list of conditions and the following disclaimer.
123
- **
124
- ** 2 Redistributions in binary form must reproduce the above copyright
125
- ** notice, this list of conditions and the following disclaimer in
126
- ** the documentation and/or other materials provided with the
127
- ** distribution.
128
- **
129
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
130
- ** the International Astronomical Union nor the names of its
131
- ** contributors may be used to endorse or promote products derived
132
- ** from this software without specific prior written permission.
133
- **
134
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
135
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
136
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
137
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
138
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
139
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
140
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
141
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
142
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
143
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
144
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
145
- ** POSSIBILITY OF SUCH DAMAGE.
146
- **
147
- **--------------------------------------------------------------------*/
148
- }