equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/fk5hz.c

@@ -1,160 +0,0 @@
-#include "sofam.h"
-
-void iauFk5hz(double r5, double d5, double date1, double date2,
-              double *rh, double *dh)
-/*
-**  - - - - - - - - -
-**   i a u F k 5 h z
-**  - - - - - - - - -
-**
-**  Transform an FK5 (J2000.0) star position into the system of the
-**  Hipparcos catalogue, assuming zero Hipparcos proper motion.
-**
-**  Status:  support function.
-**
-**  Given:
-**     r5           double   FK5 RA (radians), equinox J2000.0, at date
-**     d5           double   FK5 Dec (radians), equinox J2000.0, at date
-**     date1,date2  double   TDB date (Notes 1,2)
-**
-**  Returned:
-**     rh           double   Hipparcos RA (radians)
-**     dh           double   Hipparcos Dec (radians)
-**
-**  Notes:
-**
-**  1) This function converts a star position from the FK5 system to
-**     the Hipparcos system, in such a way that the Hipparcos proper
-**     motion is zero.  Because such a star has, in general, a non-zero
-**     proper motion in the FK5 system, the function requires the date
-**     at which the position in the FK5 system was determined.
-**
-**  2) The TT date date1+date2 is a Julian Date, apportioned in any
-**     convenient way between the two arguments.  For example,
-**     JD(TT)=2450123.7 could be expressed in any of these ways,
-**     among others:
-**
-**            date1          date2
-**
-**         2450123.7           0.0       (JD method)
-**         2451545.0       -1421.3       (J2000 method)
-**         2400000.5       50123.2       (MJD method)
-**         2450123.5           0.2       (date & time method)
-**
-**     The JD method is the most natural and convenient to use in
-**     cases where the loss of several decimal digits of resolution
-**     is acceptable.  The J2000 method is best matched to the way
-**     the argument is handled internally and will deliver the
-**     optimum resolution.  The MJD method and the date & time methods
-**     are both good compromises between resolution and convenience.
-**
-**  3) The FK5 to Hipparcos transformation is modeled as a pure
-**     rotation and spin;  zonal errors in the FK5 catalogue are not
-**     taken into account.
-**
-**  4) The position returned by this function is in the Hipparcos
-**     reference system but at date date1+date2.
-**
-**  5) See also iauFk52h, iauH2fk5, iauHfk5z.
-**
-**  Called:
-**     iauS2c       spherical coordinates to unit vector
-**     iauFk5hip    FK5 to Hipparcos rotation and spin
-**     iauSxp       multiply p-vector by scalar
-**     iauRv2m      r-vector to r-matrix
-**     iauTrxp      product of transpose of r-matrix and p-vector
-**     iauPxp       vector product of two p-vectors
-**     iauC2s       p-vector to spherical
-**     iauAnp       normalize angle into range 0 to 2pi
-**
-**  Reference:
-**
-**     F.Mignard & M.Froeschle, 2000, Astron.Astrophys. 354, 732-739.
-**
-**  This revision:  2009 December 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double t, p5e[3], r5h[3][3], s5h[3], vst[3], rst[3][3], p5[3],
-          ph[3], w;
-
-
-/* Interval from given date to fundamental epoch J2000.0 (JY). */
-   t = - ((date1 - DJ00) + date2) / DJY;
-
-/* FK5 barycentric position vector. */
-   iauS2c(r5, d5, p5e);
-
-/* FK5 to Hipparcos orientation matrix and spin vector. */
-   iauFk5hip(r5h, s5h);
-
-/* Accumulated Hipparcos wrt FK5 spin over that interval. */
-   iauSxp(t, s5h, vst);
-
-/* Express the accumulated spin as a rotation matrix. */
-   iauRv2m(vst, rst);
-
-/* Derotate the vector's FK5 axes back to date. */
-   iauTrxp(rst, p5e, p5);
-
-/* Rotate the vector into the Hipparcos system. */
-   iauRxp(r5h, p5, ph);
-
-/* Hipparcos vector to spherical. */
-   iauC2s(ph, &w, dh);
-   *rh = iauAnp(w);
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/fw2m.c

@@ -1,134 +0,0 @@
-#include "sofam.h"
-
-void iauFw2m(double gamb, double phib, double psi, double eps,
-             double r[3][3])
-/*
-**  - - - - - - - -
-**   i a u F w 2 m
-**  - - - - - - - -
-**
-**  Form rotation matrix given the Fukushima-Williams angles.
-**
-**  Status:  support function.
-**
-**  Given:
-**     gamb     double         F-W angle gamma_bar (radians)
-**     phib     double         F-W angle phi_bar (radians)
-**     psi      double         F-W angle psi (radians)
-**     eps      double         F-W angle epsilon (radians)
-**
-**  Returned:
-**     r        double[3][3]   rotation matrix
-**
-**  Notes:
-**
-**  1) Naming the following points:
-**
-**           e = J2000.0 ecliptic pole,
-**           p = GCRS pole,
-**           E = ecliptic pole of date,
-**     and   P = CIP,
-**
-**     the four Fukushima-Williams angles are as follows:
-**
-**        gamb = gamma = epE
-**        phib = phi = pE
-**        psi = psi = pEP
-**        eps = epsilon = EP
-**
-**  2) The matrix representing the combined effects of frame bias,
-**     precession and nutation is:
-**
-**        NxPxB = R_1(-eps).R_3(-psi).R_1(phib).R_3(gamb)
-**
-**  3) Three different matrices can be constructed, depending on the
-**     supplied angles:
-**
-**     o  To obtain the nutation x precession x frame bias matrix,
-**        generate the four precession angles, generate the nutation
-**        components and add them to the psi_bar and epsilon_A angles,
-**        and call the present function.
-**
-**     o  To obtain the precession x frame bias matrix, generate the
-**        four precession angles and call the present function.
-**
-**     o  To obtain the frame bias matrix, generate the four precession
-**        angles for date J2000.0 and call the present function.
-**
-**     The nutation-only and precession-only matrices can if necessary
-**     be obtained by combining these three appropriately.
-**
-**  Called:
-**     iauIr        initialize r-matrix to identity
-**     iauRz        rotate around Z-axis
-**     iauRx        rotate around X-axis
-**
-**  Reference:
-**
-**     Hilton, J. et al., 2006, Celest.Mech.Dyn.Astron. 94, 351
-**
-**  This revision:  2009 December 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-/* Construct the matrix. */
-   iauIr(r);
-   iauRz(gamb, r);
-   iauRx(phib, r);
-   iauRz(-psi, r);
-   iauRx(-eps, r);
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/fw2xy.c

@@ -1,120 +0,0 @@
-#include "sofam.h"
-
-void iauFw2xy(double gamb, double phib, double psi, double eps,
-              double *x, double *y)
-/*
-**  - - - - - - - - -
-**   i a u F w 2 x y
-**  - - - - - - - - -
-**
-**  CIP X,Y given Fukushima-Williams bias-precession-nutation angles.
-**
-**  Status:  support function.
-**
-**  Given:
-**     gamb     double    F-W angle gamma_bar (radians)
-**     phib     double    F-W angle phi_bar (radians)
-**     psi      double    F-W angle psi (radians)
-**     eps      double    F-W angle epsilon (radians)
-**
-**  Returned:
-**     x,y      double    CIP X,Y ("radians")
-**
-**  Notes:
-**
-**  1) Naming the following points:
-**
-**           e = J2000.0 ecliptic pole,
-**           p = GCRS pole
-**           E = ecliptic pole of date,
-**     and   P = CIP,
-**
-**     the four Fukushima-Williams angles are as follows:
-**
-**        gamb = gamma = epE
-**        phib = phi = pE
-**        psi = psi = pEP
-**        eps = epsilon = EP
-**
-**  2) The matrix representing the combined effects of frame bias,
-**     precession and nutation is:
-**
-**        NxPxB = R_1(-epsA).R_3(-psi).R_1(phib).R_3(gamb)
-**
-**     X,Y are elements (3,1) and (3,2) of the matrix.
-**
-**  Called:
-**     iauFw2m      F-W angles to r-matrix
-**     iauBpn2xy    extract CIP X,Y coordinates from NPB matrix
-**
-**  Reference:
-**
-**     Hilton, J. et al., 2006, Celest.Mech.Dyn.Astron. 94, 351
-**
-**  This revision:  2009 December 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double r[3][3];
-
-
-/* Form NxPxB matrix. */
-   iauFw2m(gamb, phib, psi, eps, r);
-
-/* Extract CIP X,Y. */
-   iauBpn2xy(r, x, y);
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/gc2gd.c

@@ -1,134 +0,0 @@
-#include "sofam.h"
-
-int iauGc2gd ( int n, double xyz[3],
-               double *elong, double *phi, double *height )
-/*
-**  - - - - - - - - -
-**   i a u G c 2 g d
-**  - - - - - - - - -
-**
-**  Transform geocentric coordinates to geodetic using the specified
-**  reference ellipsoid.
-**
-**  Status:  canonical transformation.
-**
-**  Given:
-**     n       int        ellipsoid identifier (Note 1)
-**     xyz     double[3]  geocentric vector (Note 2)
-**
-**  Returned:
-**     elong   double     longitude (radians, east +ve)
-**     phi     double     latitude (geodetic, radians, Note 3)
-**     height  double     height above ellipsoid (geodetic, Notes 2,3)
-**
-**  Returned (function value):
-**            int         status:  0 = OK
-**                                -1 = illegal identifier (Note 3)
-**                                -2 = internal error (Note 3)
-**
-**  Notes:
-**
-**  1) The identifier n is a number that specifies the choice of
-**     reference ellipsoid.  The following are supported:
-**
-**        n    ellipsoid
-**
-**        1     WGS84
-**        2     GRS80
-**        3     WGS72
-**
-**     The n value has no significance outside the SOFA software.  For
-**     convenience, symbols WGS84 etc. are defined in sofam.h.
-**
-**  2) The geocentric vector (xyz, given) and height (height, returned)
-**     are in meters.
-**
-**  3) An error status -1 means that the identifier n is illegal.  An
-**     error status -2 is theoretically impossible.  In all error cases,
-**     phi and height are both set to -1e9.
-**
-**  4) The inverse transformation is performed in the function iauGd2gc.
-**
-**  Called:
-**     iauEform     Earth reference ellipsoids
-**     iauGc2gde    geocentric to geodetic transformation, general
-**
-**  This revision:  2012 Febuary 23
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int j;
-   double a, f;
-
-
-/* Obtain reference ellipsoid parameters. */
-   j = iauEform ( n, &a, &f );
-
-/* If OK, transform x,y,z to longitude, geodetic latitude, height. */
-   if ( j == 0 ) {
-      j = iauGc2gde ( a, f, xyz, elong, phi, height );
-      if ( j < 0 ) j = -2;
-   }
-
-/* Deal with any errors. */
-   if ( j < 0 ) {
-      *phi = -1e9;
-      *height = -1e9;
-   }
-
-/* Return the status. */
-   return j;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-
-}