1 /*
2 * - - - - - - - - - - - - - -
3 * g a l _ f w 2 m _ t e s t
4 * - - - - - - - - - - - - - -
5 *
6 * This routine is part of the General Astrodynamics Library
7 *
8 * Description:
9 *
10 * Test gal_fw2m routine.
11 *
12 * Status:
13 *
14 * Internal test routine
15 *
16 * Called:
17 *
18 * gal_fw2m
19 * gal_vdv
20 *
21 * This Revision:
22 *
23 * 2008 January 11 ( c version 2008 June 14 )
24 *
25 * Copyright (C) 2008 Paul C. L. Willmott. See notes at end.
26 */
27
28 #include "gal_fw2m.h"
29 #include "gal_fw2m_test.h"
30 #include "gal_test.h"
31
32 void
gal_fw2m_test()33 gal_fw2m_test
34 (
35 )
36 {
37
38 double gamb, phib, psi, eps, r[3][3] ;
39
40 extern int gal_tfunc ;
41
42 /*
43 * - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
44 */
45
46 gamb = -0.2243387670997992368e-05 ;
47 phib = 0.4091014602391312982 ;
48 psi = -0.9501954178013015092e-03 ;
49 eps = 0.4091014316587367472 ;
50
51 gal_fw2m ( gamb, phib, psi, eps, r ) ;
52 gal_vdv ( r[0][0], 0.9999995505176006994, 1e-12, "gal_fw2m", "r[0][0]" ) ;
53 gal_vdv ( r[0][1], 0.8695404617348192425e-03, 1e-12, "gal_fw2m", "r[0][1]" ) ;
54 gal_vdv ( r[0][2], 0.3779735201865582514e-03, 1e-12, "gal_fw2m", "r[0][2]" ) ;
55 gal_vdv ( r[1][0], -0.8695404723772016172e-03, 1e-12, "gal_fw2m", "r[1][0]" ) ;
56 gal_vdv ( r[1][1], 0.9999996219496026795, 1e-12, "gal_fw2m", "r[1][1]" ) ;
57 gal_vdv ( r[1][2], -0.1361752496022920582e-06, 1e-12, "gal_fw2m", "r[1][2]" ) ;
58 gal_vdv ( r[2][0], -0.3779734957034082605e-03, 1e-12, "gal_fw2m", "r[2][0]" ) ;
59 gal_vdv ( r[2][1], -0.1924880848614416492e-06, 1e-12, "gal_fw2m", "r[2][1]" ) ;
60 gal_vdv ( r[2][2], 0.9999999285679972427, 1e-12, "gal_fw2m", "r[2][2]" ) ;
61
62 gal_tfunc++ ;
63
64 /*
65 * Finished.
66 */
67
68 }
69
70 /*
71 * gal - General Astrodynamics Library
72 * Copyright (C) 2008 Paul C. L. Willmott
73 *
74 * This program is free software; you can redistribute it and/or modify
75 * it under the terms of the GNU General Public License as published by
76 * the Free Software Foundation; either version 2 of the License, or
77 * (at your option) any later version.
78 *
79 * This program is distributed in the hope that it will be useful,
80 * but WITHOUT ANY WARRANTY; without even the implied warranty of
81 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
82 * GNU General Public License for more details.
83 *
84 * You should have received a copy of the GNU General Public License along
85 * with this program; if not, write to the Free Software Foundation, Inc.,
86 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
87 *
88 * Contact:
89 *
90 * Paul Willmott
91 * vp9mu@amsat.org
92 */
93