equationoftime 3.0.0 → 4.0.0

Sign up to get free protection for your applications and to get access to all the features.
Files changed (251) hide show
  1. checksums.yaml +4 -4
  2. data/.buildpath +5 -0
  3. data/.project +17 -0
  4. data/.rvmrc +1 -0
  5. data/.settings/org.eclipse.ltk.core.refactoring.prefs +2 -0
  6. data/Makefile +238 -0
  7. data/Rakefile +40 -3
  8. data/analemma_data.yml +2925 -0
  9. data/equationoftime.gemspec +6 -3
  10. data/examples/analemma_data_generator.rb +12 -12
  11. data/examples/celes_parts/.RUBYARCHDIR.time +0 -0
  12. data/examples/celes_parts/.cproject +56 -0
  13. data/examples/celes_parts/.project +78 -0
  14. data/examples/celes_parts/a2af.c +120 -0
  15. data/examples/celes_parts/a2tf.c +116 -0
  16. data/examples/celes_parts/af2a.c +107 -0
  17. data/examples/celes_parts/anp.c +82 -0
  18. data/examples/celes_parts/anpm.c +82 -0
  19. data/examples/celes_parts/bi00.c +116 -0
  20. data/examples/celes_parts/bp00.c +173 -0
  21. data/examples/celes_parts/bp06.c +136 -0
  22. data/examples/celes_parts/bpn2xy.c +100 -0
  23. data/examples/celes_parts/c2i00a.c +139 -0
  24. data/examples/celes_parts/c2i00b.c +139 -0
  25. data/examples/celes_parts/c2i06a.c +136 -0
  26. data/examples/celes_parts/c2ibpn.c +142 -0
  27. data/examples/celes_parts/c2ixy.c +131 -0
  28. data/examples/celes_parts/c2ixys.c +123 -0
  29. data/examples/celes_parts/c2s.c +96 -0
  30. data/examples/celes_parts/c2t00a.c +154 -0
  31. data/examples/celes_parts/c2t00b.c +150 -0
  32. data/examples/celes_parts/c2t06a.c +152 -0
  33. data/examples/celes_parts/c2tcio.c +122 -0
  34. data/examples/celes_parts/c2teqx.c +122 -0
  35. data/examples/celes_parts/c2tpe.c +167 -0
  36. data/examples/celes_parts/c2txy.c +159 -0
  37. data/examples/celes_parts/cal2jd.c +139 -0
  38. data/examples/celes_parts/celes_core.c +2522 -0
  39. data/examples/celes_parts/cp.c +80 -0
  40. data/examples/celes_parts/cpv.c +82 -0
  41. data/examples/celes_parts/cr.c +83 -0
  42. data/examples/celes_parts/d2dtf.c +206 -0
  43. data/examples/celes_parts/d2tf.c +160 -0
  44. data/examples/celes_parts/dat.c +289 -0
  45. data/examples/celes_parts/depend +25 -0
  46. data/examples/celes_parts/dtdb.c +1213 -0
  47. data/examples/celes_parts/dtf2d.c +196 -0
  48. data/examples/celes_parts/ee00.c +128 -0
  49. data/examples/celes_parts/ee00a.c +135 -0
  50. data/examples/celes_parts/ee00b.c +141 -0
  51. data/examples/celes_parts/ee06a.c +122 -0
  52. data/examples/celes_parts/eect00.c +282 -0
  53. data/examples/celes_parts/eform.c +147 -0
  54. data/examples/celes_parts/eo06a.c +131 -0
  55. data/examples/celes_parts/eors.c +108 -0
  56. data/examples/celes_parts/epb.c +96 -0
  57. data/examples/celes_parts/epb2jd.c +91 -0
  58. data/examples/celes_parts/epj.c +93 -0
  59. data/examples/celes_parts/epj2jd.c +91 -0
  60. data/examples/celes_parts/epv00.c +2589 -0
  61. data/examples/celes_parts/eqeq94.c +131 -0
  62. data/examples/celes_parts/era00.c +136 -0
  63. data/examples/celes_parts/extconf.rb +4 -0
  64. data/examples/celes_parts/fad03.c +103 -0
  65. data/examples/celes_parts/fae03.c +102 -0
  66. data/examples/celes_parts/faf03.c +106 -0
  67. data/examples/celes_parts/faju03.c +102 -0
  68. data/examples/celes_parts/fal03.c +103 -0
  69. data/examples/celes_parts/falp03.c +103 -0
  70. data/examples/celes_parts/fama03.c +102 -0
  71. data/examples/celes_parts/fame03.c +102 -0
  72. data/examples/celes_parts/fane03.c +99 -0
  73. data/examples/celes_parts/faom03.c +104 -0
  74. data/examples/celes_parts/fapa03.c +103 -0
  75. data/examples/celes_parts/fasa03.c +102 -0
  76. data/examples/celes_parts/faur03.c +99 -0
  77. data/examples/celes_parts/fave03.c +102 -0
  78. data/examples/celes_parts/fk52h.c +143 -0
  79. data/examples/celes_parts/fk5hip.c +126 -0
  80. data/examples/celes_parts/fk5hz.c +160 -0
  81. data/examples/celes_parts/fw2m.c +134 -0
  82. data/examples/celes_parts/fw2xy.c +120 -0
  83. data/examples/celes_parts/gc2gd.c +134 -0
  84. data/examples/celes_parts/gc2gde.c +200 -0
  85. data/examples/celes_parts/gd2gc.c +134 -0
  86. data/examples/celes_parts/gd2gce.c +138 -0
  87. data/examples/celes_parts/gmst00.c +145 -0
  88. data/examples/celes_parts/gmst06.c +136 -0
  89. data/examples/celes_parts/gmst82.c +151 -0
  90. data/examples/celes_parts/gst00a.c +138 -0
  91. data/examples/celes_parts/gst00b.c +146 -0
  92. data/examples/celes_parts/gst06.c +140 -0
  93. data/examples/celes_parts/gst06a.c +131 -0
  94. data/examples/celes_parts/gst94.c +131 -0
  95. data/examples/celes_parts/h2fk5.c +148 -0
  96. data/examples/celes_parts/hfk5z.c +175 -0
  97. data/examples/celes_parts/ir.c +83 -0
  98. data/examples/celes_parts/jd2cal.c +155 -0
  99. data/examples/celes_parts/jdcalf.c +161 -0
  100. data/examples/celes_parts/num00a.c +121 -0
  101. data/examples/celes_parts/num00b.c +121 -0
  102. data/examples/celes_parts/num06a.c +125 -0
  103. data/examples/celes_parts/numat.c +109 -0
  104. data/examples/celes_parts/nut00a.c +2047 -0
  105. data/examples/celes_parts/nut00b.c +372 -0
  106. data/examples/celes_parts/nut06a.c +153 -0
  107. data/examples/celes_parts/nut80.c +325 -0
  108. data/examples/celes_parts/nutm80.c +117 -0
  109. data/examples/celes_parts/obl06.c +118 -0
  110. data/examples/celes_parts/obl80.c +118 -0
  111. data/examples/celes_parts/p06e.c +321 -0
  112. data/examples/celes_parts/p2pv.c +83 -0
  113. data/examples/celes_parts/p2s.c +91 -0
  114. data/examples/celes_parts/pap.c +139 -0
  115. data/examples/celes_parts/pas.c +96 -0
  116. data/examples/celes_parts/pb06.c +144 -0
  117. data/examples/celes_parts/pdp.c +84 -0
  118. data/examples/celes_parts/pfw06.c +165 -0
  119. data/examples/celes_parts/plan94.c +514 -0
  120. data/examples/celes_parts/pm.c +83 -0
  121. data/examples/celes_parts/pmat00.c +118 -0
  122. data/examples/celes_parts/pmat06.c +122 -0
  123. data/examples/celes_parts/pmat76.c +141 -0
  124. data/examples/celes_parts/pmp.c +85 -0
  125. data/examples/celes_parts/pn.c +109 -0
  126. data/examples/celes_parts/pn00.c +177 -0
  127. data/examples/celes_parts/pn00a.c +162 -0
  128. data/examples/celes_parts/pn00b.c +162 -0
  129. data/examples/celes_parts/pn06.c +187 -0
  130. data/examples/celes_parts/pn06a.c +152 -0
  131. data/examples/celes_parts/pnm00a.c +121 -0
  132. data/examples/celes_parts/pnm00b.c +121 -0
  133. data/examples/celes_parts/pnm06a.c +124 -0
  134. data/examples/celes_parts/pnm80.c +126 -0
  135. data/examples/celes_parts/pom00.c +115 -0
  136. data/examples/celes_parts/ppp.c +85 -0
  137. data/examples/celes_parts/ppsp.c +94 -0
  138. data/examples/celes_parts/pr00.c +142 -0
  139. data/examples/celes_parts/prec76.c +148 -0
  140. data/examples/celes_parts/pv2p.c +81 -0
  141. data/examples/celes_parts/pv2s.c +144 -0
  142. data/examples/celes_parts/pvdpv.c +102 -0
  143. data/examples/celes_parts/pvm.c +86 -0
  144. data/examples/celes_parts/pvmpv.c +87 -0
  145. data/examples/celes_parts/pvppv.c +87 -0
  146. data/examples/celes_parts/pvstar.c +207 -0
  147. data/examples/celes_parts/pvu.c +93 -0
  148. data/examples/celes_parts/pvup.c +88 -0
  149. data/examples/celes_parts/pvxpv.c +107 -0
  150. data/examples/celes_parts/pxp.c +94 -0
  151. data/examples/celes_parts/rm2v.c +111 -0
  152. data/examples/celes_parts/rv2m.c +118 -0
  153. data/examples/celes_parts/rx.c +110 -0
  154. data/examples/celes_parts/rxp.c +99 -0
  155. data/examples/celes_parts/rxpv.c +86 -0
  156. data/examples/celes_parts/rxr.c +99 -0
  157. data/examples/celes_parts/ry.c +110 -0
  158. data/examples/celes_parts/rz.c +110 -0
  159. data/examples/celes_parts/s00.c +371 -0
  160. data/examples/celes_parts/s00a.c +143 -0
  161. data/examples/celes_parts/s00b.c +143 -0
  162. data/examples/celes_parts/s06.c +368 -0
  163. data/examples/celes_parts/s06a.c +145 -0
  164. data/examples/celes_parts/s2c.c +85 -0
  165. data/examples/celes_parts/s2p.c +88 -0
  166. data/examples/celes_parts/s2pv.c +103 -0
  167. data/examples/celes_parts/s2xpv.c +87 -0
  168. data/examples/celes_parts/sepp.c +105 -0
  169. data/examples/celes_parts/seps.c +93 -0
  170. data/examples/celes_parts/sofa.h +379 -0
  171. data/examples/celes_parts/sofam.h +155 -0
  172. data/examples/celes_parts/sp00.c +118 -0
  173. data/examples/celes_parts/starpm.c +205 -0
  174. data/examples/celes_parts/starpv.c +264 -0
  175. data/examples/celes_parts/sxp.c +84 -0
  176. data/examples/celes_parts/sxpv.c +85 -0
  177. data/examples/celes_parts/taitt.c +110 -0
  178. data/examples/celes_parts/taiut1.c +112 -0
  179. data/examples/celes_parts/taiutc.c +182 -0
  180. data/examples/celes_parts/tcbtdb.c +132 -0
  181. data/examples/celes_parts/tcgtt.c +109 -0
  182. data/examples/celes_parts/tdbtcb.c +137 -0
  183. data/examples/celes_parts/tdbtt.c +122 -0
  184. data/examples/celes_parts/test_celes.rb +48 -0
  185. data/examples/celes_parts/tf2a.c +107 -0
  186. data/examples/celes_parts/tf2d.c +107 -0
  187. data/examples/celes_parts/tr.c +93 -0
  188. data/examples/celes_parts/trxp.c +93 -0
  189. data/examples/celes_parts/trxpv.c +93 -0
  190. data/examples/celes_parts/tttai.c +110 -0
  191. data/examples/celes_parts/tttcg.c +112 -0
  192. data/examples/celes_parts/tttdb.c +121 -0
  193. data/examples/celes_parts/ttut1.c +110 -0
  194. data/examples/celes_parts/ut1tai.c +111 -0
  195. data/examples/celes_parts/ut1tt.c +110 -0
  196. data/examples/celes_parts/ut1utc.c +193 -0
  197. data/examples/celes_parts/utctai.c +163 -0
  198. data/examples/celes_parts/utcut1.c +151 -0
  199. data/examples/celes_parts/xy06.c +2758 -0
  200. data/examples/celes_parts/xys00a.c +133 -0
  201. data/examples/celes_parts/xys00b.c +133 -0
  202. data/examples/celes_parts/xys06a.c +133 -0
  203. data/examples/celes_parts/zp.c +77 -0
  204. data/examples/celes_parts/zpv.c +79 -0
  205. data/examples/celes_parts/zr.c +83 -0
  206. data/examples/check_date_type.rb +2 -4
  207. data/examples/compare_geoc_long_ra.rb +19 -11
  208. data/examples/data_table_for_astro_dog.rb +40 -0
  209. data/examples/earth_rotation.rb +28 -3
  210. data/examples/eot.c +15 -0
  211. data/examples/eot.h +2 -0
  212. data/examples/eot_methods_list.rb +30 -3
  213. data/examples/eot_suntimes.rb +99 -95
  214. data/examples/extconf.rb +2 -0
  215. data/examples/{read_nutation_data.rb → gmst_gast_non_sofa.rb} +103 -98
  216. data/examples/nutation_table5_3a.yaml +1658 -1658
  217. data/examples/rbeot.c +24 -0
  218. data/examples/test_eot.rb +17 -0
  219. data/examples/test_poly_eval.rb +6 -17
  220. data/examples/time_scales.rb +4 -7
  221. data/examples/times_year.rb +50 -0
  222. data/examples/usage_example.rb +15 -5
  223. data/examples/use_angles.rb +151 -87
  224. data/lib/eot/angles.rb +130 -177
  225. data/lib/eot/constants.rb +41 -44
  226. data/lib/eot/displays.rb +40 -37
  227. data/lib/eot/geo_lat_lng_smt.rb +7 -6
  228. data/lib/eot/init.rb +62 -39
  229. data/lib/eot/nutation.rb +53 -45
  230. data/lib/eot/times.rb +32 -41
  231. data/lib/eot/utilities.rb +11 -90
  232. data/lib/eot/version.rb +3 -2
  233. data/rise_set_data.yml +1461 -0
  234. data/run_tests_eclipse.rb +1 -0
  235. data/tests/minitest/aliased_angles_spec.rb +119 -165
  236. data/tests/minitest/aliased_displays_spec.rb +41 -39
  237. data/tests/minitest/aliased_utilities_spec.rb +11 -22
  238. data/tests/minitest/angles_spec.rb +134 -181
  239. data/tests/minitest/constants_spec.rb +4 -10
  240. data/tests/minitest/delta_epsilon_spec.rb +13 -16
  241. data/tests/minitest/displays_spec.rb +39 -37
  242. data/tests/minitest/geo_spec.rb +15 -11
  243. data/tests/minitest/init_spec.rb +30 -16
  244. data/tests/minitest/nutation_spec.rb +19 -13
  245. data/tests/minitest/times_spec.rb +77 -75
  246. data/tests/minitest/utilities_spec.rb +10 -94
  247. data/wiki.md +7 -9
  248. metadata +218 -26
  249. data/examples/data_table.rb +0 -26
  250. data/lib/eot/nutation_table5_3a.yaml +0 -9532
  251. data/tests/minitest/aliased_times_spec.rb +0 -36
@@ -0,0 +1,107 @@
1
+ #include "sofam.h"
2
+
3
+ void iauPvxpv(double a[2][3], double b[2][3], double axb[2][3])
4
+ /*
5
+ ** - - - - - - - - -
6
+ ** i a u P v x p v
7
+ ** - - - - - - - - -
8
+ **
9
+ ** Outer (=vector=cross) product of two pv-vectors.
10
+ **
11
+ ** Status: vector/matrix support function.
12
+ **
13
+ ** Given:
14
+ ** a double[2][3] first pv-vector
15
+ ** b double[2][3] second pv-vector
16
+ **
17
+ ** Returned:
18
+ ** axb double[2][3] a x b
19
+ **
20
+ ** Notes:
21
+ **
22
+ ** 1) If the position and velocity components of the two pv-vectors are
23
+ ** ( ap, av ) and ( bp, bv ), the result, a x b, is the pair of
24
+ ** vectors ( ap x bp, ap x bv + av x bp ). The two vectors are the
25
+ ** cross-product of the two p-vectors and its derivative.
26
+ **
27
+ ** 2) It is permissible to re-use the same array for any of the
28
+ ** arguments.
29
+ **
30
+ ** Called:
31
+ ** iauCpv copy pv-vector
32
+ ** iauPxp vector product of two p-vectors
33
+ ** iauPpp p-vector plus p-vector
34
+ **
35
+ ** This revision: 2008 November 18
36
+ **
37
+ ** Original version 2012-03-01
38
+ **
39
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
40
+ */
41
+ {
42
+ double wa[2][3], wb[2][3], axbd[3], adxb[3];
43
+
44
+
45
+ /* Make copies of the inputs. */
46
+ iauCpv(a, wa);
47
+ iauCpv(b, wb);
48
+
49
+ /* a x b = position part of result. */
50
+ iauPxp(wa[0], wb[0], axb[0]);
51
+
52
+ /* a x bdot + adot x b = velocity part of result. */
53
+ iauPxp(wa[0], wb[1], axbd);
54
+ iauPxp(wa[1], wb[0], adxb);
55
+ iauPpp(axbd, adxb, axb[1]);
56
+
57
+ return;
58
+
59
+ /*----------------------------------------------------------------------
60
+ **
61
+ ** Celes is a wrapper of the SOFA Library for Ruby.
62
+ **
63
+ ** This file is redistributed and relicensed in accordance with
64
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
65
+ **
66
+ ** The original library is available from IAU Standards of
67
+ ** Fundamental Astronomy (http://www.iausofa.org/).
68
+ **
69
+ **
70
+ **
71
+ **
72
+ **
73
+ ** Copyright (C) 2013, Naoki Arita
74
+ ** All rights reserved.
75
+ **
76
+ ** Redistribution and use in source and binary forms, with or without
77
+ ** modification, are permitted provided that the following conditions
78
+ ** are met:
79
+ **
80
+ ** 1 Redistributions of source code must retain the above copyright
81
+ ** notice, this list of conditions and the following disclaimer.
82
+ **
83
+ ** 2 Redistributions in binary form must reproduce the above copyright
84
+ ** notice, this list of conditions and the following disclaimer in
85
+ ** the documentation and/or other materials provided with the
86
+ ** distribution.
87
+ **
88
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
89
+ ** the International Astronomical Union nor the names of its
90
+ ** contributors may be used to endorse or promote products derived
91
+ ** from this software without specific prior written permission.
92
+ **
93
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
94
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
95
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
96
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
97
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
98
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
99
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
100
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
101
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
102
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
103
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
104
+ ** POSSIBILITY OF SUCH DAMAGE.
105
+ **
106
+ **--------------------------------------------------------------------*/
107
+ }
@@ -0,0 +1,94 @@
1
+ #include "sofam.h"
2
+
3
+ void iauPxp(double a[3], double b[3], double axb[3])
4
+ /*
5
+ ** - - - - - - -
6
+ ** i a u P x p
7
+ ** - - - - - - -
8
+ **
9
+ ** p-vector outer (=vector=cross) product.
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
+ ** axb double[3] a x 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
+ double xa, ya, za, xb, yb, zb;
32
+
33
+
34
+ xa = a[0];
35
+ ya = a[1];
36
+ za = a[2];
37
+ xb = b[0];
38
+ yb = b[1];
39
+ zb = b[2];
40
+ axb[0] = ya*zb - za*yb;
41
+ axb[1] = za*xb - xa*zb;
42
+ axb[2] = xa*yb - ya*xb;
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
+ }
@@ -0,0 +1,111 @@
1
+ #include "sofam.h"
2
+
3
+ void iauRm2v(double r[3][3], double w[3])
4
+ /*
5
+ ** - - - - - - - -
6
+ ** i a u R m 2 v
7
+ ** - - - - - - - -
8
+ **
9
+ ** Express an r-matrix as an r-vector.
10
+ **
11
+ ** Status: vector/matrix support function.
12
+ **
13
+ ** Given:
14
+ ** r double[3][3] rotation matrix
15
+ **
16
+ ** Returned:
17
+ ** w double[3] rotation vector (Note 1)
18
+ **
19
+ ** Notes:
20
+ **
21
+ ** 1) A rotation matrix describes a rotation through some angle about
22
+ ** some arbitrary axis called the Euler axis. The "rotation vector"
23
+ ** returned by this function has the same direction as the Euler axis,
24
+ ** and its magnitude is the angle in radians. (The magnitude and
25
+ ** direction can be separated by means of the function iauPn.)
26
+ **
27
+ ** 2) If r is null, so is the result. If r is not a rotation matrix
28
+ ** the result is undefined; r must be proper (i.e. have a positive
29
+ ** determinant) and real orthogonal (inverse = transpose).
30
+ **
31
+ ** 3) The reference frame rotates clockwise as seen looking along
32
+ ** the rotation vector from the origin.
33
+ **
34
+ ** This revision: 2010 August 27
35
+ **
36
+ ** Original version 2012-03-01
37
+ **
38
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
39
+ */
40
+ {
41
+ double x, y, z, s2, c2, phi, f;
42
+
43
+
44
+ x = r[1][2] - r[2][1];
45
+ y = r[2][0] - r[0][2];
46
+ z = r[0][1] - r[1][0];
47
+ s2 = sqrt(x*x + y*y + z*z);
48
+ if (s2 != 0) {
49
+ c2 = r[0][0] + r[1][1] + r[2][2] - 1.0;
50
+ phi = atan2(s2, c2);
51
+ f = phi / s2;
52
+ w[0] = x * f;
53
+ w[1] = y * f;
54
+ w[2] = z * f;
55
+ } else {
56
+ w[0] = 0.0;
57
+ w[1] = 0.0;
58
+ w[2] = 0.0;
59
+ }
60
+
61
+ return;
62
+
63
+ /*----------------------------------------------------------------------
64
+ **
65
+ ** Celes is a wrapper of the SOFA Library for Ruby.
66
+ **
67
+ ** This file is redistributed and relicensed in accordance with
68
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
69
+ **
70
+ ** The original library is available from IAU Standards of
71
+ ** Fundamental Astronomy (http://www.iausofa.org/).
72
+ **
73
+ **
74
+ **
75
+ **
76
+ **
77
+ ** Copyright (C) 2013, Naoki Arita
78
+ ** All rights reserved.
79
+ **
80
+ ** Redistribution and use in source and binary forms, with or without
81
+ ** modification, are permitted provided that the following conditions
82
+ ** are met:
83
+ **
84
+ ** 1 Redistributions of source code must retain the above copyright
85
+ ** notice, this list of conditions and the following disclaimer.
86
+ **
87
+ ** 2 Redistributions in binary form must reproduce the above copyright
88
+ ** notice, this list of conditions and the following disclaimer in
89
+ ** the documentation and/or other materials provided with the
90
+ ** distribution.
91
+ **
92
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
93
+ ** the International Astronomical Union nor the names of its
94
+ ** contributors may be used to endorse or promote products derived
95
+ ** from this software without specific prior written permission.
96
+ **
97
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
98
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
99
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
100
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
101
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
102
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
103
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
104
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
105
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
106
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
107
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
108
+ ** POSSIBILITY OF SUCH DAMAGE.
109
+ **
110
+ **--------------------------------------------------------------------*/
111
+ }
@@ -0,0 +1,118 @@
1
+ #include "sofam.h"
2
+
3
+ void iauRv2m(double w[3], double r[3][3])
4
+ /*
5
+ ** - - - - - - - -
6
+ ** i a u R v 2 m
7
+ ** - - - - - - - -
8
+ **
9
+ ** Form the r-matrix corresponding to a given r-vector.
10
+ **
11
+ ** Status: vector/matrix support function.
12
+ **
13
+ ** Given:
14
+ ** w double[3] rotation vector (Note 1)
15
+ **
16
+ ** Returned:
17
+ ** r double[3][3] rotation matrix
18
+ **
19
+ ** Notes:
20
+ **
21
+ ** 1) A rotation matrix describes a rotation through some angle about
22
+ ** some arbitrary axis called the Euler axis. The "rotation vector"
23
+ ** supplied to This function has the same direction as the Euler
24
+ ** axis, and its magnitude is the angle in radians.
25
+ **
26
+ ** 2) If w is null, the unit matrix is returned.
27
+ **
28
+ ** 3) The reference frame rotates clockwise as seen looking along the
29
+ ** rotation vector from the origin.
30
+ **
31
+ ** This revision: 2008 May 11
32
+ **
33
+ ** Original version 2012-03-01
34
+ **
35
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
36
+ */
37
+ {
38
+ double x, y, z, phi, s, c, f;
39
+
40
+
41
+ /* Euler angle (magnitude of rotation vector) and functions. */
42
+ x = w[0];
43
+ y = w[1];
44
+ z = w[2];
45
+ phi = sqrt(x*x + y*y + z*z);
46
+ s = sin(phi);
47
+ c = cos(phi);
48
+ f = 1.0 - c;
49
+
50
+ /* Euler axis (direction of rotation vector), perhaps null. */
51
+ if (phi != 0.0) {
52
+ x /= phi;
53
+ y /= phi;
54
+ z /= phi;
55
+ }
56
+
57
+ /* Form the rotation matrix. */
58
+ r[0][0] = x*x*f + c;
59
+ r[0][1] = x*y*f + z*s;
60
+ r[0][2] = x*z*f - y*s;
61
+ r[1][0] = y*x*f - z*s;
62
+ r[1][1] = y*y*f + c;
63
+ r[1][2] = y*z*f + x*s;
64
+ r[2][0] = z*x*f + y*s;
65
+ r[2][1] = z*y*f - x*s;
66
+ r[2][2] = z*z*f + c;
67
+
68
+ return;
69
+
70
+ /*----------------------------------------------------------------------
71
+ **
72
+ ** Celes is a wrapper of the SOFA Library for Ruby.
73
+ **
74
+ ** This file is redistributed and relicensed in accordance with
75
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
76
+ **
77
+ ** The original library is available from IAU Standards of
78
+ ** Fundamental Astronomy (http://www.iausofa.org/).
79
+ **
80
+ **
81
+ **
82
+ **
83
+ **
84
+ ** Copyright (C) 2013, Naoki Arita
85
+ ** All rights reserved.
86
+ **
87
+ ** Redistribution and use in source and binary forms, with or without
88
+ ** modification, are permitted provided that the following conditions
89
+ ** are met:
90
+ **
91
+ ** 1 Redistributions of source code must retain the above copyright
92
+ ** notice, this list of conditions and the following disclaimer.
93
+ **
94
+ ** 2 Redistributions in binary form must reproduce the above copyright
95
+ ** notice, this list of conditions and the following disclaimer in
96
+ ** the documentation and/or other materials provided with the
97
+ ** distribution.
98
+ **
99
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
100
+ ** the International Astronomical Union nor the names of its
101
+ ** contributors may be used to endorse or promote products derived
102
+ ** from this software without specific prior written permission.
103
+ **
104
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
105
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
106
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
107
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
108
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
109
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
110
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
111
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
112
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
113
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
114
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
115
+ ** POSSIBILITY OF SUCH DAMAGE.
116
+ **
117
+ **--------------------------------------------------------------------*/
118
+ }
@@ -0,0 +1,110 @@
1
+ #include "sofam.h"
2
+
3
+ void iauRx(double phi, double r[3][3])
4
+ /*
5
+ ** - - - - - -
6
+ ** i a u R x
7
+ ** - - - - - -
8
+ **
9
+ ** Rotate an r-matrix about the x-axis.
10
+ **
11
+ ** Status: vector/matrix support function.
12
+ **
13
+ ** Given:
14
+ ** phi double angle (radians)
15
+ **
16
+ ** Given and returned:
17
+ ** r double[3][3] r-matrix, rotated
18
+ **
19
+ ** Notes:
20
+ **
21
+ ** 1) Calling this function with positive phi incorporates in the
22
+ ** supplied r-matrix r an additional rotation, about the x-axis,
23
+ ** anticlockwise as seen looking towards the origin from positive x.
24
+ **
25
+ ** 2) The additional rotation can be represented by this matrix:
26
+ **
27
+ ** ( 1 0 0 )
28
+ ** ( )
29
+ ** ( 0 + cos(phi) + sin(phi) )
30
+ ** ( )
31
+ ** ( 0 - sin(phi) + cos(phi) )
32
+ **
33
+ ** This revision: 2012 April 3
34
+ **
35
+ ** Original version 2012-03-01
36
+ **
37
+ ** Copyright (C) 2013 Naoki Arita. See notes at end.
38
+ */
39
+ {
40
+ double s, c, a10, a11, a12, a20, a21, a22;
41
+
42
+
43
+ s = sin(phi);
44
+ c = cos(phi);
45
+
46
+ a10 = c*r[1][0] + s*r[2][0];
47
+ a11 = c*r[1][1] + s*r[2][1];
48
+ a12 = c*r[1][2] + s*r[2][2];
49
+ a20 = - s*r[1][0] + c*r[2][0];
50
+ a21 = - s*r[1][1] + c*r[2][1];
51
+ a22 = - s*r[1][2] + c*r[2][2];
52
+
53
+ r[1][0] = a10;
54
+ r[1][1] = a11;
55
+ r[1][2] = a12;
56
+ r[2][0] = a20;
57
+ r[2][1] = a21;
58
+ r[2][2] = a22;
59
+
60
+ return;
61
+
62
+ /*----------------------------------------------------------------------
63
+ **
64
+ ** Celes is a wrapper of the SOFA Library for Ruby.
65
+ **
66
+ ** This file is redistributed and relicensed in accordance with
67
+ ** the SOFA Software License (http://www.iausofa.org/tandc.html).
68
+ **
69
+ ** The original library is available from IAU Standards of
70
+ ** Fundamental Astronomy (http://www.iausofa.org/).
71
+ **
72
+ **
73
+ **
74
+ **
75
+ **
76
+ ** Copyright (C) 2013, Naoki Arita
77
+ ** All rights reserved.
78
+ **
79
+ ** Redistribution and use in source and binary forms, with or without
80
+ ** modification, are permitted provided that the following conditions
81
+ ** are met:
82
+ **
83
+ ** 1 Redistributions of source code must retain the above copyright
84
+ ** notice, this list of conditions and the following disclaimer.
85
+ **
86
+ ** 2 Redistributions in binary form must reproduce the above copyright
87
+ ** notice, this list of conditions and the following disclaimer in
88
+ ** the documentation and/or other materials provided with the
89
+ ** distribution.
90
+ **
91
+ ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
92
+ ** the International Astronomical Union nor the names of its
93
+ ** contributors may be used to endorse or promote products derived
94
+ ** from this software without specific prior written permission.
95
+ **
96
+ ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
97
+ ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
98
+ ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
99
+ ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
100
+ ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
101
+ ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
102
+ ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
103
+ ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
104
+ ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
105
+ ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
106
+ ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
107
+ ** POSSIBILITY OF SUCH DAMAGE.
108
+ **
109
+ **--------------------------------------------------------------------*/
110
+ }