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,121 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPnm00a(double date1, double date2, double rbpn[3][3])
4
- /*
5
- ** - - - - - - - - - -
6
- ** i a u P n m 0 0 a
7
- ** - - - - - - - - - -
8
- **
9
- ** Form the matrix of precession-nutation for a given date (including
10
- ** frame bias), equinox-based, IAU 2000A model.
11
- **
12
- ** Status: support function.
13
- **
14
- ** Given:
15
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
16
- **
17
- ** Returned:
18
- ** rbpn double[3][3] classical NPB matrix (Note 2)
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 matrix operates in the sense V(date) = rbpn * V(GCRS), where
42
- ** the p-vector V(date) is with respect to the true equatorial triad
43
- ** of date date1+date2 and the p-vector V(GCRS) is with respect to
44
- ** the Geocentric Celestial Reference System (IAU, 2000).
45
- **
46
- ** 3) A faster, but slightly less accurate result (about 1 mas), can be
47
- ** obtained by using instead the iauPnm00b function.
48
- **
49
- ** Called:
50
- ** iauPn00a bias/precession/nutation, IAU 2000A
51
- **
52
- ** Reference:
53
- **
54
- ** IAU: Trans. International Astronomical Union, Vol. XXIVB; Proc.
55
- ** 24th General Assembly, Manchester, UK. Resolutions B1.3, B1.6.
56
- ** (2000)
57
- **
58
- ** This revision: 2009 December 21
59
- **
60
- ** Original version 2012-03-01
61
- **
62
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
63
- */
64
- {
65
- double dpsi, deps, epsa, rb[3][3], rp[3][3], rbp[3][3], rn[3][3];
66
-
67
-
68
- /* Obtain the required matrix (discarding other results). */
69
- iauPn00a(date1, date2, &dpsi, &deps, &epsa, rb, rp, rbp, rn, rbpn);
70
-
71
- return;
72
-
73
- /*----------------------------------------------------------------------
74
- **
75
- ** Celes is a wrapper of the SOFA Library for Ruby.
76
- **
77
- ** This file is redistributed and relicensed in accordance with
78
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
79
- **
80
- ** The original library is available from IAU Standards of
81
- ** Fundamental Astronomy (http://www.iausofa.org/).
82
- **
83
- **
84
- **
85
- **
86
- **
87
- ** Copyright (C) 2013, Naoki Arita
88
- ** All rights reserved.
89
- **
90
- ** Redistribution and use in source and binary forms, with or without
91
- ** modification, are permitted provided that the following conditions
92
- ** are met:
93
- **
94
- ** 1 Redistributions of source code must retain the above copyright
95
- ** notice, this list of conditions and the following disclaimer.
96
- **
97
- ** 2 Redistributions in binary form must reproduce the above copyright
98
- ** notice, this list of conditions and the following disclaimer in
99
- ** the documentation and/or other materials provided with the
100
- ** distribution.
101
- **
102
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
103
- ** the International Astronomical Union nor the names of its
104
- ** contributors may be used to endorse or promote products derived
105
- ** from this software without specific prior written permission.
106
- **
107
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
108
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
109
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
110
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
111
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
112
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
113
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
114
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
115
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
116
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
117
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
118
- ** POSSIBILITY OF SUCH DAMAGE.
119
- **
120
- **--------------------------------------------------------------------*/
121
- }
@@ -1,121 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPnm00b(double date1, double date2, double rbpn[3][3])
4
- /*
5
- ** - - - - - - - - - -
6
- ** i a u P n m 0 0 b
7
- ** - - - - - - - - - -
8
- **
9
- ** Form the matrix of precession-nutation for a given date (including
10
- ** frame bias), equinox-based, IAU 2000B model.
11
- **
12
- ** Status: support function.
13
- **
14
- ** Given:
15
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
16
- **
17
- ** Returned:
18
- ** rbpn double[3][3] bias-precession-nutation matrix (Note 2)
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 matrix operates in the sense V(date) = rbpn * V(GCRS), where
42
- ** the p-vector V(date) is with respect to the true equatorial triad
43
- ** of date date1+date2 and the p-vector V(GCRS) is with respect to
44
- ** the Geocentric Celestial Reference System (IAU, 2000).
45
- **
46
- ** 3) The present function is faster, but slightly less accurate (about
47
- ** 1 mas), than the iauPnm00a function.
48
- **
49
- ** Called:
50
- ** iauPn00b bias/precession/nutation, IAU 2000B
51
- **
52
- ** Reference:
53
- **
54
- ** IAU: Trans. International Astronomical Union, Vol. XXIVB; Proc.
55
- ** 24th General Assembly, Manchester, UK. Resolutions B1.3, B1.6.
56
- ** (2000)
57
- **
58
- ** This revision: 2009 December 21
59
- **
60
- ** Original version 2012-03-01
61
- **
62
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
63
- */
64
- {
65
- double dpsi, deps, epsa, rb[3][3], rp[3][3], rbp[3][3], rn[3][3];
66
-
67
-
68
- /* Obtain the required matrix (discarding other results). */
69
- iauPn00b(date1, date2, &dpsi, &deps, &epsa, rb, rp, rbp, rn, rbpn);
70
-
71
- return;
72
-
73
- /*----------------------------------------------------------------------
74
- **
75
- ** Celes is a wrapper of the SOFA Library for Ruby.
76
- **
77
- ** This file is redistributed and relicensed in accordance with
78
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
79
- **
80
- ** The original library is available from IAU Standards of
81
- ** Fundamental Astronomy (http://www.iausofa.org/).
82
- **
83
- **
84
- **
85
- **
86
- **
87
- ** Copyright (C) 2013, Naoki Arita
88
- ** All rights reserved.
89
- **
90
- ** Redistribution and use in source and binary forms, with or without
91
- ** modification, are permitted provided that the following conditions
92
- ** are met:
93
- **
94
- ** 1 Redistributions of source code must retain the above copyright
95
- ** notice, this list of conditions and the following disclaimer.
96
- **
97
- ** 2 Redistributions in binary form must reproduce the above copyright
98
- ** notice, this list of conditions and the following disclaimer in
99
- ** the documentation and/or other materials provided with the
100
- ** distribution.
101
- **
102
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
103
- ** the International Astronomical Union nor the names of its
104
- ** contributors may be used to endorse or promote products derived
105
- ** from this software without specific prior written permission.
106
- **
107
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
108
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
109
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
110
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
111
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
112
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
113
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
114
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
115
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
116
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
117
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
118
- ** POSSIBILITY OF SUCH DAMAGE.
119
- **
120
- **--------------------------------------------------------------------*/
121
- }
@@ -1,124 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPnm06a(double date1, double date2, double rnpb[3][3])
4
- /*
5
- ** - - - - - - - - - -
6
- ** i a u P n m 0 6 a
7
- ** - - - - - - - - - -
8
- **
9
- ** Form the matrix of precession-nutation for a given date (including
10
- ** frame bias), IAU 2006 precession and IAU 2000A nutation models.
11
- **
12
- ** Status: support function.
13
- **
14
- ** Given:
15
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
16
- **
17
- ** Returned:
18
- ** rnpb double[3][3] bias-precession-nutation matrix (Note 2)
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 matrix operates in the sense V(date) = rnpb * V(GCRS), where
42
- ** the p-vector V(date) is with respect to the true equatorial triad
43
- ** of date date1+date2 and the p-vector V(GCRS) is with respect to
44
- ** the Geocentric Celestial Reference System (IAU, 2000).
45
- **
46
- ** Called:
47
- ** iauPfw06 bias-precession F-W angles, IAU 2006
48
- ** iauNut06a nutation, IAU 2006/2000A
49
- ** iauFw2m F-W angles to r-matrix
50
- **
51
- ** Reference:
52
- **
53
- ** Capitaine, N. & Wallace, P.T., 2006, Astron.Astrophys. 450, 855.
54
- **
55
- ** This revision: 2009 December 21
56
- **
57
- ** Original version 2012-03-01
58
- **
59
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
60
- */
61
- {
62
- double gamb, phib, psib, epsa, dp, de;
63
-
64
-
65
- /* Fukushima-Williams angles for frame bias and precession. */
66
- iauPfw06(date1, date2, &gamb, &phib, &psib, &epsa);
67
-
68
- /* Nutation components. */
69
- iauNut06a(date1, date2, &dp, &de);
70
-
71
- /* Equinox based nutation x precession x bias matrix. */
72
- iauFw2m(gamb, phib, psib + dp, epsa + de, rnpb);
73
-
74
- return;
75
-
76
- /*----------------------------------------------------------------------
77
- **
78
- ** Celes is a wrapper of the SOFA Library for Ruby.
79
- **
80
- ** This file is redistributed and relicensed in accordance with
81
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
82
- **
83
- ** The original library is available from IAU Standards of
84
- ** Fundamental Astronomy (http://www.iausofa.org/).
85
- **
86
- **
87
- **
88
- **
89
- **
90
- ** Copyright (C) 2013, Naoki Arita
91
- ** All rights reserved.
92
- **
93
- ** Redistribution and use in source and binary forms, with or without
94
- ** modification, are permitted provided that the following conditions
95
- ** are met:
96
- **
97
- ** 1 Redistributions of source code must retain the above copyright
98
- ** notice, this list of conditions and the following disclaimer.
99
- **
100
- ** 2 Redistributions in binary form must reproduce the above copyright
101
- ** notice, this list of conditions and the following disclaimer in
102
- ** the documentation and/or other materials provided with the
103
- ** distribution.
104
- **
105
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
106
- ** the International Astronomical Union nor the names of its
107
- ** contributors may be used to endorse or promote products derived
108
- ** from this software without specific prior written permission.
109
- **
110
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
111
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
112
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
113
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
114
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
115
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
116
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
117
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
118
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
119
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
120
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
121
- ** POSSIBILITY OF SUCH DAMAGE.
122
- **
123
- **--------------------------------------------------------------------*/
124
- }
@@ -1,126 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauPnm80(double date1, double date2, double rmatpn[3][3])
4
- /*
5
- ** - - - - - - - - -
6
- ** i a u P n m 8 0
7
- ** - - - - - - - - -
8
- **
9
- ** Form the matrix of precession/nutation for a given date, IAU 1976
10
- ** precession model, IAU 1980 nutation model.
11
- **
12
- ** Status: support function.
13
- **
14
- ** Given:
15
- ** date1,date2 double TDB date (Note 1)
16
- **
17
- ** Returned:
18
- ** rmatpn double[3][3] combined precession/nutation matrix
19
- **
20
- ** Notes:
21
- **
22
- ** 1) The TDB date date1+date2 is a Julian Date, apportioned in any
23
- ** convenient way between the two arguments. For example,
24
- ** JD(TDB)=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 matrix operates in the sense V(date) = rmatpn * V(J2000),
42
- ** where the p-vector V(date) is with respect to the true equatorial
43
- ** triad of date date1+date2 and the p-vector V(J2000) is with
44
- ** respect to the mean equatorial triad of epoch J2000.0.
45
- **
46
- ** Called:
47
- ** iauPmat76 precession matrix, IAU 1976
48
- ** iauNutm80 nutation matrix, IAU 1980
49
- ** iauRxr product of two r-matrices
50
- **
51
- ** Reference:
52
- **
53
- ** Explanatory Supplement to the Astronomical Almanac,
54
- ** P. Kenneth Seidelmann (ed), University Science Books (1992),
55
- ** Section 3.3 (p145).
56
- **
57
- ** This revision: 2010 January 23
58
- **
59
- ** Original version 2012-03-01
60
- **
61
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
62
- */
63
- {
64
- double rmatp[3][3], rmatn[3][3];
65
-
66
-
67
- /* Precession matrix, J2000.0 to date. */
68
- iauPmat76(date1, date2, rmatp);
69
-
70
- /* Nutation matrix. */
71
- iauNutm80(date1, date2, rmatn);
72
-
73
- /* Combine the matrices: PN = N x P. */
74
- iauRxr(rmatn, rmatp, rmatpn);
75
-
76
- return;
77
-
78
- /*----------------------------------------------------------------------
79
- **
80
- ** Celes is a wrapper of the SOFA Library for Ruby.
81
- **
82
- ** This file is redistributed and relicensed in accordance with
83
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
84
- **
85
- ** The original library is available from IAU Standards of
86
- ** Fundamental Astronomy (http://www.iausofa.org/).
87
- **
88
- **
89
- **
90
- **
91
- **
92
- ** Copyright (C) 2013, Naoki Arita
93
- ** All rights reserved.
94
- **
95
- ** Redistribution and use in source and binary forms, with or without
96
- ** modification, are permitted provided that the following conditions
97
- ** are met:
98
- **
99
- ** 1 Redistributions of source code must retain the above copyright
100
- ** notice, this list of conditions and the following disclaimer.
101
- **
102
- ** 2 Redistributions in binary form must reproduce the above copyright
103
- ** notice, this list of conditions and the following disclaimer in
104
- ** the documentation and/or other materials provided with the
105
- ** distribution.
106
- **
107
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
108
- ** the International Astronomical Union nor the names of its
109
- ** contributors may be used to endorse or promote products derived
110
- ** from this software without specific prior written permission.
111
- **
112
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
113
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
114
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
115
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
116
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
117
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
118
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
119
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
120
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
121
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
122
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
123
- ** POSSIBILITY OF SUCH DAMAGE.
124
- **
125
- **--------------------------------------------------------------------*/
126
- }