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,133 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauXys00a(double date1, double date2,
4
- double *x, double *y, double *s)
5
- /*
6
- ** - - - - - - - - - -
7
- ** i a u X y s 0 0 a
8
- ** - - - - - - - - - -
9
- **
10
- ** For a given TT date, compute the X,Y coordinates of the Celestial
11
- ** Intermediate Pole and the CIO locator s, using the IAU 2000A
12
- ** precession-nutation model.
13
- **
14
- ** Status: support function.
15
- **
16
- ** Given:
17
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
18
- **
19
- ** Returned:
20
- ** x,y double Celestial Intermediate Pole (Note 2)
21
- ** s double the CIO locator s (Note 2)
22
- **
23
- ** Notes:
24
- **
25
- ** 1) The TT date date1+date2 is a Julian Date, apportioned in any
26
- ** convenient way between the two arguments. For example,
27
- ** JD(TT)=2450123.7 could be expressed in any of these ways,
28
- ** among others:
29
- **
30
- ** date1 date2
31
- **
32
- ** 2450123.7 0.0 (JD method)
33
- ** 2451545.0 -1421.3 (J2000 method)
34
- ** 2400000.5 50123.2 (MJD method)
35
- ** 2450123.5 0.2 (date & time method)
36
- **
37
- ** The JD method is the most natural and convenient to use in
38
- ** cases where the loss of several decimal digits of resolution
39
- ** is acceptable. The J2000 method is best matched to the way
40
- ** the argument is handled internally and will deliver the
41
- ** optimum resolution. The MJD method and the date & time methods
42
- ** are both good compromises between resolution and convenience.
43
- **
44
- ** 2) The Celestial Intermediate Pole coordinates are the x,y
45
- ** components of the unit vector in the Geocentric Celestial
46
- ** Reference System.
47
- **
48
- ** 3) The CIO locator s (in radians) positions the Celestial
49
- ** Intermediate Origin on the equator of the CIP.
50
- **
51
- ** 4) A faster, but slightly less accurate result (about 1 mas for
52
- ** X,Y), can be obtained by using instead the iauXys00b function.
53
- **
54
- ** Called:
55
- ** iauPnm00a classical NPB matrix, IAU 2000A
56
- ** iauBpn2xy extract CIP X,Y coordinates from NPB matrix
57
- ** iauS00 the CIO locator s, given X,Y, IAU 2000A
58
- **
59
- ** Reference:
60
- **
61
- ** McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
62
- ** IERS Technical Note No. 32, BKG (2004)
63
- **
64
- ** This revision: 2008 May 12
65
- **
66
- ** Original version 2012-03-01
67
- **
68
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
69
- */
70
- {
71
- double rbpn[3][3];
72
-
73
-
74
- /* Form the bias-precession-nutation matrix, IAU 2000A. */
75
- iauPnm00a(date1, date2, rbpn);
76
-
77
- /* Extract X,Y. */
78
- iauBpn2xy(rbpn, x, y);
79
-
80
- /* Obtain s. */
81
- *s = iauS00(date1, date2, *x, *y);
82
-
83
- return;
84
-
85
- /*----------------------------------------------------------------------
86
- **
87
- ** Celes is a wrapper of the SOFA Library for Ruby.
88
- **
89
- ** This file is redistributed and relicensed in accordance with
90
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
91
- **
92
- ** The original library is available from IAU Standards of
93
- ** Fundamental Astronomy (http://www.iausofa.org/).
94
- **
95
- **
96
- **
97
- **
98
- **
99
- ** Copyright (C) 2013, Naoki Arita
100
- ** All rights reserved.
101
- **
102
- ** Redistribution and use in source and binary forms, with or without
103
- ** modification, are permitted provided that the following conditions
104
- ** are met:
105
- **
106
- ** 1 Redistributions of source code must retain the above copyright
107
- ** notice, this list of conditions and the following disclaimer.
108
- **
109
- ** 2 Redistributions in binary form must reproduce the above copyright
110
- ** notice, this list of conditions and the following disclaimer in
111
- ** the documentation and/or other materials provided with the
112
- ** distribution.
113
- **
114
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
115
- ** the International Astronomical Union nor the names of its
116
- ** contributors may be used to endorse or promote products derived
117
- ** from this software without specific prior written permission.
118
- **
119
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
120
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
121
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
122
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
123
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
124
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
125
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
126
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
127
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
128
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
129
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
130
- ** POSSIBILITY OF SUCH DAMAGE.
131
- **
132
- **--------------------------------------------------------------------*/
133
- }
@@ -1,133 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauXys00b(double date1, double date2,
4
- double *x, double *y, double *s)
5
- /*
6
- ** - - - - - - - - - -
7
- ** i a u X y s 0 0 b
8
- ** - - - - - - - - - -
9
- **
10
- ** For a given TT date, compute the X,Y coordinates of the Celestial
11
- ** Intermediate Pole and the CIO locator s, using the IAU 2000B
12
- ** precession-nutation model.
13
- **
14
- ** Status: support function.
15
- **
16
- ** Given:
17
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
18
- **
19
- ** Returned:
20
- ** x,y double Celestial Intermediate Pole (Note 2)
21
- ** s double the CIO locator s (Note 2)
22
- **
23
- ** Notes:
24
- **
25
- ** 1) The TT date date1+date2 is a Julian Date, apportioned in any
26
- ** convenient way between the two arguments. For example,
27
- ** JD(TT)=2450123.7 could be expressed in any of these ways,
28
- ** among others:
29
- **
30
- ** date1 date2
31
- **
32
- ** 2450123.7 0.0 (JD method)
33
- ** 2451545.0 -1421.3 (J2000 method)
34
- ** 2400000.5 50123.2 (MJD method)
35
- ** 2450123.5 0.2 (date & time method)
36
- **
37
- ** The JD method is the most natural and convenient to use in
38
- ** cases where the loss of several decimal digits of resolution
39
- ** is acceptable. The J2000 method is best matched to the way
40
- ** the argument is handled internally and will deliver the
41
- ** optimum resolution. The MJD method and the date & time methods
42
- ** are both good compromises between resolution and convenience.
43
- **
44
- ** 2) The Celestial Intermediate Pole coordinates are the x,y
45
- ** components of the unit vector in the Geocentric Celestial
46
- ** Reference System.
47
- **
48
- ** 3) The CIO locator s (in radians) positions the Celestial
49
- ** Intermediate Origin on the equator of the CIP.
50
- **
51
- ** 4) The present function is faster, but slightly less accurate (about
52
- ** 1 mas in X,Y), than the iauXys00a function.
53
- **
54
- ** Called:
55
- ** iauPnm00b classical NPB matrix, IAU 2000B
56
- ** iauBpn2xy extract CIP X,Y coordinates from NPB matrix
57
- ** iauS00 the CIO locator s, given X,Y, IAU 2000A
58
- **
59
- ** Reference:
60
- **
61
- ** McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
62
- ** IERS Technical Note No. 32, BKG (2004)
63
- **
64
- ** This revision: 2008 May 12
65
- **
66
- ** Original version 2012-03-01
67
- **
68
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
69
- */
70
- {
71
- double rbpn[3][3];
72
-
73
-
74
- /* Form the bias-precession-nutation matrix, IAU 2000A. */
75
- iauPnm00b(date1, date2, rbpn);
76
-
77
- /* Extract X,Y. */
78
- iauBpn2xy(rbpn, x, y);
79
-
80
- /* Obtain s. */
81
- *s = iauS00(date1, date2, *x, *y);
82
-
83
- return;
84
-
85
- /*----------------------------------------------------------------------
86
- **
87
- ** Celes is a wrapper of the SOFA Library for Ruby.
88
- **
89
- ** This file is redistributed and relicensed in accordance with
90
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
91
- **
92
- ** The original library is available from IAU Standards of
93
- ** Fundamental Astronomy (http://www.iausofa.org/).
94
- **
95
- **
96
- **
97
- **
98
- **
99
- ** Copyright (C) 2013, Naoki Arita
100
- ** All rights reserved.
101
- **
102
- ** Redistribution and use in source and binary forms, with or without
103
- ** modification, are permitted provided that the following conditions
104
- ** are met:
105
- **
106
- ** 1 Redistributions of source code must retain the above copyright
107
- ** notice, this list of conditions and the following disclaimer.
108
- **
109
- ** 2 Redistributions in binary form must reproduce the above copyright
110
- ** notice, this list of conditions and the following disclaimer in
111
- ** the documentation and/or other materials provided with the
112
- ** distribution.
113
- **
114
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
115
- ** the International Astronomical Union nor the names of its
116
- ** contributors may be used to endorse or promote products derived
117
- ** from this software without specific prior written permission.
118
- **
119
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
120
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
121
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
122
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
123
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
124
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
125
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
126
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
127
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
128
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
129
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
130
- ** POSSIBILITY OF SUCH DAMAGE.
131
- **
132
- **--------------------------------------------------------------------*/
133
- }
@@ -1,133 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauXys06a(double date1, double date2,
4
- double *x, double *y, double *s)
5
- /*
6
- ** - - - - - - - - - -
7
- ** i a u X y s 0 6 a
8
- ** - - - - - - - - - -
9
- **
10
- ** For a given TT date, compute the X,Y coordinates of the Celestial
11
- ** Intermediate Pole and the CIO locator s, using the IAU 2006
12
- ** precession and IAU 2000A nutation models.
13
- **
14
- ** Status: support function.
15
- **
16
- ** Given:
17
- ** date1,date2 double TT as a 2-part Julian Date (Note 1)
18
- **
19
- ** Returned:
20
- ** x,y double Celestial Intermediate Pole (Note 2)
21
- ** s double the CIO locator s (Note 2)
22
- **
23
- ** Notes:
24
- **
25
- ** 1) The TT date date1+date2 is a Julian Date, apportioned in any
26
- ** convenient way between the two arguments. For example,
27
- ** JD(TT)=2450123.7 could be expressed in any of these ways,
28
- ** among others:
29
- **
30
- ** date1 date2
31
- **
32
- ** 2450123.7 0.0 (JD method)
33
- ** 2451545.0 -1421.3 (J2000 method)
34
- ** 2400000.5 50123.2 (MJD method)
35
- ** 2450123.5 0.2 (date & time method)
36
- **
37
- ** The JD method is the most natural and convenient to use in
38
- ** cases where the loss of several decimal digits of resolution
39
- ** is acceptable. The J2000 method is best matched to the way
40
- ** the argument is handled internally and will deliver the
41
- ** optimum resolution. The MJD method and the date & time methods
42
- ** are both good compromises between resolution and convenience.
43
- **
44
- ** 2) The Celestial Intermediate Pole coordinates are the x,y components
45
- ** of the unit vector in the Geocentric Celestial Reference System.
46
- **
47
- ** 3) The CIO locator s (in radians) positions the Celestial
48
- ** Intermediate Origin on the equator of the CIP.
49
- **
50
- ** 4) Series-based solutions for generating X and Y are also available:
51
- ** see Capitaine & Wallace (2006) and iauXy06.
52
- **
53
- ** Called:
54
- ** iauPnm06a classical NPB matrix, IAU 2006/2000A
55
- ** iauBpn2xy extract CIP X,Y coordinates from NPB matrix
56
- ** iauS06 the CIO locator s, given X,Y, IAU 2006
57
- **
58
- ** References:
59
- **
60
- ** Capitaine, N. & Wallace, P.T., 2006, Astron.Astrophys. 450, 855
61
- **
62
- ** Wallace, P.T. & Capitaine, N., 2006, Astron.Astrophys. 459, 981
63
- **
64
- ** This revision: 2008 May 11
65
- **
66
- ** Original version 2012-03-01
67
- **
68
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
69
- */
70
- {
71
- double rbpn[3][3];
72
-
73
-
74
- /* Form the bias-precession-nutation matrix, IAU 2000A. */
75
- iauPnm06a(date1, date2, rbpn);
76
-
77
- /* Extract X,Y. */
78
- iauBpn2xy(rbpn, x, y);
79
-
80
- /* Obtain s. */
81
- *s = iauS06(date1, date2, *x, *y);
82
-
83
- return;
84
-
85
- /*----------------------------------------------------------------------
86
- **
87
- ** Celes is a wrapper of the SOFA Library for Ruby.
88
- **
89
- ** This file is redistributed and relicensed in accordance with
90
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
91
- **
92
- ** The original library is available from IAU Standards of
93
- ** Fundamental Astronomy (http://www.iausofa.org/).
94
- **
95
- **
96
- **
97
- **
98
- **
99
- ** Copyright (C) 2013, Naoki Arita
100
- ** All rights reserved.
101
- **
102
- ** Redistribution and use in source and binary forms, with or without
103
- ** modification, are permitted provided that the following conditions
104
- ** are met:
105
- **
106
- ** 1 Redistributions of source code must retain the above copyright
107
- ** notice, this list of conditions and the following disclaimer.
108
- **
109
- ** 2 Redistributions in binary form must reproduce the above copyright
110
- ** notice, this list of conditions and the following disclaimer in
111
- ** the documentation and/or other materials provided with the
112
- ** distribution.
113
- **
114
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
115
- ** the International Astronomical Union nor the names of its
116
- ** contributors may be used to endorse or promote products derived
117
- ** from this software without specific prior written permission.
118
- **
119
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
120
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
121
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
122
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
123
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
124
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
125
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
126
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
127
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
128
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
129
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
130
- ** POSSIBILITY OF SUCH DAMAGE.
131
- **
132
- **--------------------------------------------------------------------*/
133
- }
@@ -1,77 +0,0 @@
1
- #include "sofam.h"
2
-
3
- void iauZp(double p[3])
4
- /*
5
- ** - - - - - -
6
- ** i a u Z p
7
- ** - - - - - -
8
- **
9
- ** Zero a p-vector.
10
- **
11
- ** Status: vector/matrix support function.
12
- **
13
- ** Returned:
14
- ** p double[3] p-vector
15
- **
16
- ** This revision: 2008 May 11
17
- **
18
- ** Original version 2012-03-01
19
- **
20
- ** Copyright (C) 2013 Naoki Arita. See notes at end.
21
- */
22
- {
23
- p[0] = 0.0;
24
- p[1] = 0.0;
25
- p[2] = 0.0;
26
-
27
- return;
28
-
29
- /*----------------------------------------------------------------------
30
- **
31
- ** Celes is a wrapper of the SOFA Library for Ruby.
32
- **
33
- ** This file is redistributed and relicensed in accordance with
34
- ** the SOFA Software License (http://www.iausofa.org/tandc.html).
35
- **
36
- ** The original library is available from IAU Standards of
37
- ** Fundamental Astronomy (http://www.iausofa.org/).
38
- **
39
- **
40
- **
41
- **
42
- **
43
- ** Copyright (C) 2013, Naoki Arita
44
- ** All rights reserved.
45
- **
46
- ** Redistribution and use in source and binary forms, with or without
47
- ** modification, are permitted provided that the following conditions
48
- ** are met:
49
- **
50
- ** 1 Redistributions of source code must retain the above copyright
51
- ** notice, this list of conditions and the following disclaimer.
52
- **
53
- ** 2 Redistributions in binary form must reproduce the above copyright
54
- ** notice, this list of conditions and the following disclaimer in
55
- ** the documentation and/or other materials provided with the
56
- ** distribution.
57
- **
58
- ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
59
- ** the International Astronomical Union nor the names of its
60
- ** contributors may be used to endorse or promote products derived
61
- ** from this software without specific prior written permission.
62
- **
63
- ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
64
- ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
65
- ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
66
- ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
67
- ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
68
- ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
69
- ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
70
- ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
71
- ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
72
- ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
73
- ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
74
- ** POSSIBILITY OF SUCH DAMAGE.
75
- **
76
- **--------------------------------------------------------------------*/
77
- }