1 #include "erfa.h"
2
eraS2pv(double theta,double phi,double r,double td,double pd,double rd,double pv[2][3])3 void eraS2pv(double theta, double phi, double r,
4 double td, double pd, double rd,
5 double pv[2][3])
6 /*
7 ** - - - - - - - -
8 ** e r a S 2 p v
9 ** - - - - - - - -
10 **
11 ** Convert position/velocity from spherical to Cartesian coordinates.
12 **
13 ** Given:
14 ** theta double longitude angle (radians)
15 ** phi double latitude angle (radians)
16 ** r double radial distance
17 ** td double rate of change of theta
18 ** pd double rate of change of phi
19 ** rd double rate of change of r
20 **
21 ** Returned:
22 ** pv double[2][3] pv-vector
23 **
24 ** Copyright (C) 2013-2014, NumFOCUS Foundation.
25 ** Derived, with permission, from the SOFA library. See notes at end of file.
26 */
27 {
28 double st, ct, sp, cp, rcp, x, y, rpd, w;
29
30
31 st = sin(theta);
32 ct = cos(theta);
33 sp = sin(phi);
34 cp = cos(phi);
35 rcp = r * cp;
36 x = rcp * ct;
37 y = rcp * st;
38 rpd = r * pd;
39 w = rpd*sp - cp*rd;
40
41 pv[0][0] = x;
42 pv[0][1] = y;
43 pv[0][2] = r * sp;
44 pv[1][0] = -y*td - w*ct;
45 pv[1][1] = x*td - w*st;
46 pv[1][2] = rpd*cp + sp*rd;
47
48 return;
49
50 }
51 /*----------------------------------------------------------------------
52 **
53 **
54 ** Copyright (C) 2013-2014, NumFOCUS Foundation.
55 ** All rights reserved.
56 **
57 ** This library is derived, with permission, from the International
58 ** Astronomical Union's "Standards of Fundamental Astronomy" library,
59 ** available from http://www.iausofa.org.
60 **
61 ** The ERFA version is intended to retain identical functionality to
62 ** the SOFA library, but made distinct through different function and
63 ** file names, as set out in the SOFA license conditions. The SOFA
64 ** original has a role as a reference standard for the IAU and IERS,
65 ** and consequently redistribution is permitted only in its unaltered
66 ** state. The ERFA version is not subject to this restriction and
67 ** therefore can be included in distributions which do not support the
68 ** concept of "read only" software.
69 **
70 ** Although the intent is to replicate the SOFA API (other than
71 ** replacement of prefix names) and results (with the exception of
72 ** bugs; any that are discovered will be fixed), SOFA is not
73 ** responsible for any errors found in this version of the library.
74 **
75 ** If you wish to acknowledge the SOFA heritage, please acknowledge
76 ** that you are using a library derived from SOFA, rather than SOFA
77 ** itself.
78 **
79 **
80 ** TERMS AND CONDITIONS
81 **
82 ** Redistribution and use in source and binary forms, with or without
83 ** modification, are permitted provided that the following conditions
84 ** are met:
85 **
86 ** 1 Redistributions of source code must retain the above copyright
87 ** notice, this list of conditions and the following disclaimer.
88 **
89 ** 2 Redistributions in binary form must reproduce the above copyright
90 ** notice, this list of conditions and the following disclaimer in
91 ** the documentation and/or other materials provided with the
92 ** distribution.
93 **
94 ** 3 Neither the name of the Standards Of Fundamental Astronomy Board,
95 ** the International Astronomical Union nor the names of its
96 ** contributors may be used to endorse or promote products derived
97 ** from this software without specific prior written permission.
98 **
99 ** THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
100 ** "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
101 ** LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
102 ** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
103 ** COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
104 ** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
105 ** BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
106 ** LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
107 ** CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
108 ** LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
109 ** ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
110 ** POSSIBILITY OF SUCH DAMAGE.
111 **
112 */
113