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,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
- }