File Coverage

palsrc/palPrebn.c
Criterion Covered Total %
statement 11 11 100.0
branch n/a
condition n/a
subroutine n/a
pod n/a
total 11 11 100.0


line stmt bran cond sub pod time code
1             /*
2             *+
3             * Name:
4             * palPrebn
5              
6             * Purpose:
7             * Generate the matrix of precession between two objects (old)
8              
9             * Language:
10             * Starlink ANSI C
11              
12             * Type of Module:
13             * Library routine
14              
15             * Invocation:
16             * void palPrebn ( double bep0, double bep1, double rmatp[3][3] );
17              
18             * Arguments:
19             * bep0 = double (Given)
20             * Beginning Besselian epoch.
21             * bep1 = double (Given)
22             * Ending Besselian epoch
23             * rmatp = double[3][3] (Returned)
24             * precession matrix in the sense V(BEP1) = RMATP * V(BEP0)
25              
26             * Description:
27             * Generate the matrix of precession between two epochs,
28             * using the old, pre-IAU1976, Bessel-Newcomb model, using
29             * Kinoshita's formulation
30              
31             * Authors:
32             * PTW: Pat Wallace (STFC)
33             * TIMJ: Tim Jenness (JAC, Hawaii)
34             * {enter_new_authors_here}
35              
36             * See Also:
37             * Kinoshita, H. (1975) 'Formulas for precession', SAO Special
38             * Report No. 364, Smithsonian Institution Astrophysical
39             * Observatory, Cambridge, Massachusetts.
40              
41             * History:
42             * 2012-02-12(TIMJ):
43             * Initial version with documentation taken from Fortran SLA
44             * Adapted with permission from the Fortran SLALIB library.
45             * {enter_further_changes_here}
46              
47             * Copyright:
48             * Copyright (C) 1996 Rutherford Appleton Laboratory
49             * Copyright (C) 2012 Science and Technology Facilities Council.
50             * All Rights Reserved.
51              
52             * Licence:
53             * This program is free software: you can redistribute it and/or
54             * modify it under the terms of the GNU Lesser General Public
55             * License as published by the Free Software Foundation, either
56             * version 3 of the License, or (at your option) any later
57             * version.
58             *
59             * This program is distributed in the hope that it will be useful,
60             * but WITHOUT ANY WARRANTY; without even the implied warranty of
61             * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
62             * GNU Lesser General Public License for more details.
63             *
64             * You should have received a copy of the GNU Lesser General
65             * License along with this program. If not, see
66             * .
67              
68             * Bugs:
69             * {note_any_bugs_here}
70             *-
71             */
72              
73             #include "pal.h"
74             #include "palmac.h"
75              
76 3           void palPrebn ( double bep0, double bep1, double rmatp[3][3] ) {
77              
78             double t,bigt, zeta, theta, z, tas2r, w;
79              
80             /* Interval between basic epoch B1850.0 and beginning epoch in TC */
81 3           bigt = (bep0-1850)/100.;
82              
83             /* Interval over which precession required, in tropical centuries */
84 3           t = (bep1-bep0)/100.;
85              
86             /* Euler angles */
87 3           tas2r = t * PAL__DAS2R;
88 3           w = 2303.5548 + ( 1.39720 + 0.000059 * bigt) * bigt;
89              
90 3           zeta = ( w + ( 0.30242 - 0.000269 * bigt + 0.017996 * t ) * t ) * tas2r;
91 3           z = ( w + ( 1.09478 + 0.000387 * bigt + 0.018324 * t ) * t ) * tas2r;
92 6           theta = ( 2005.1125 + ( -0.85294 - 0.000365 * bigt ) * bigt +
93 3           (-0.42647 - 0.000365 * bigt - 0.041802 * t ) * t ) * tas2r;
94              
95             /* Rotation matrix */
96 3           palDeuler("ZYZ", -zeta, theta, -z, rmatp);
97              
98 3           }