equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/h2fk5.c

@@ -1,148 +0,0 @@
-#include "sofam.h"
-
-void iauH2fk5(double rh, double dh,
-              double drh, double ddh, double pxh, double rvh,
-              double *r5, double *d5,
-              double *dr5, double *dd5, double *px5, double *rv5)
-/*
-**  - - - - - - - - -
-**   i a u H 2 f k 5
-**  - - - - - - - - -
-**
-**  Transform Hipparcos star data into the FK5 (J2000.0) system.
-**
-**  Status:  support function.
-**
-**  Given (all Hipparcos, epoch J2000.0):
-**     rh      double    RA (radians)
-**     dh      double    Dec (radians)
-**     drh     double    proper motion in RA (dRA/dt, rad/Jyear)
-**     ddh     double    proper motion in Dec (dDec/dt, rad/Jyear)
-**     pxh     double    parallax (arcsec)
-**     rvh     double    radial velocity (km/s, positive = receding)
-**
-**  Returned (all FK5, equinox J2000.0, epoch J2000.0):
-**     r5      double    RA (radians)
-**     d5      double    Dec (radians)
-**     dr5     double    proper motion in RA (dRA/dt, rad/Jyear)
-**     dd5     double    proper motion in Dec (dDec/dt, rad/Jyear)
-**     px5     double    parallax (arcsec)
-**     rv5     double    radial velocity (km/s, positive = receding)
-**
-**  Notes:
-**
-**  1) This function transforms Hipparcos star positions and proper
-**     motions into FK5 J2000.0.
-**
-**  2) The proper motions in RA are dRA/dt rather than
-**     cos(Dec)*dRA/dt, and are per year rather than per century.
-**
-**  3) The FK5 to Hipparcos transformation is modeled as a pure
-**     rotation and spin;  zonal errors in the FK5 catalog are not
-**     taken into account.
-**
-**  4) See also iauFk52h, iauFk5hz, iauHfk5z.
-**
-**  Called:
-**     iauStarpv    star catalog data to space motion pv-vector
-**     iauFk5hip    FK5 to Hipparcos rotation and spin
-**     iauRv2m      r-vector to r-matrix
-**     iauRxp       product of r-matrix and p-vector
-**     iauTrxp      product of transpose of r-matrix and p-vector
-**     iauPxp       vector product of two p-vectors
-**     iauPmp       p-vector minus p-vector
-**     iauPvstar    space motion pv-vector to star catalog data
-**
-**  Reference:
-**
-**     F.Mignard & M.Froeschle, Astron. Astrophys. 354, 732-739 (2000).
-**
-**  This revision:  2009 December 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int i;
-   double pvh[2][3], r5h[3][3], s5h[3], sh[3], wxp[3], vv[3], pv5[2][3];
-
-
-/* Hipparcos barycentric position/velocity pv-vector (normalized). */
-   iauStarpv(rh, dh, drh, ddh, pxh, rvh, pvh);
-
-/* FK5 to Hipparcos orientation matrix and spin vector. */
-   iauFk5hip(r5h, s5h);
-
-/* Make spin units per day instead of per year. */
-   for ( i = 0; i < 3; s5h[i++] /= 365.25 );
-
-/* Orient the spin into the Hipparcos system. */
-   iauRxp(r5h, s5h, sh);
-
-/* De-orient the Hipparcos position into the FK5 system. */
-   iauTrxp(r5h, pvh[0], pv5[0]);
-
-/* Apply spin to the position giving an extra space motion component. */
-   iauPxp(pvh[0], sh, wxp);
-
-/* Subtract this component from the Hipparcos space motion. */
-   iauPmp(pvh[1], wxp, vv);
-
-/* De-orient the Hipparcos space motion into the FK5 system. */
-   iauTrxp(r5h, vv, pv5[1]);
-
-/* FK5 pv-vector to spherical. */
-   iauPvstar(pv5, r5, d5, dr5, dd5, px5, rv5);
-
-   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/hfk5z.c

@@ -1,175 +0,0 @@
-#include "sofam.h"
-
-void iauHfk5z(double rh, double dh, double date1, double date2,
-              double *r5, double *d5, double *dr5, double *dd5)
-/*
-**  - - - - - - - - -
-**   i a u H f k 5 z
-**  - - - - - - - - -
-**
-**  Transform a Hipparcos star position into FK5 J2000.0, assuming
-**  zero Hipparcos proper motion.
-**
-**  Status:  support function.
-**
-**  Given:
-**     rh            double    Hipparcos RA (radians)
-**     dh            double    Hipparcos Dec (radians)
-**     date1,date2   double    TDB date (Note 1)
-**
-**  Returned (all FK5, equinox J2000.0, date date1+date2):
-**     r5            double    RA (radians)
-**     d5            double    Dec (radians)
-**     dr5           double    FK5 RA proper motion (rad/year, Note 4)
-**     dd5           double    Dec proper motion (rad/year, Note 4)
-**
-**  Notes:
-**
-**  1) 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.
-**
-**  2) The proper motion in RA is dRA/dt rather than cos(Dec)*dRA/dt.
-**
-**  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) It was the intention that Hipparcos should be a close
-**     approximation to an inertial frame, so that distant objects have
-**     zero proper motion;  such objects have (in general) non-zero
-**     proper motion in FK5, and this function returns those fictitious
-**     proper motions.
-**
-**  5) The position returned by this function is in the FK5 J2000.0
-**     reference system but at date date1+date2.
-**
-**  6) See also iauFk52h, iauH2fk5, iauFk5zhz.
-**
-**  Called:
-**     iauS2c       spherical coordinates to unit vector
-**     iauFk5hip    FK5 to Hipparcos rotation and spin
-**     iauRxp       product of r-matrix and p-vector
-**     iauSxp       multiply p-vector by scalar
-**     iauRxr       product of two r-matrices
-**     iauTrxp      product of transpose of r-matrix and p-vector
-**     iauPxp       vector product of two p-vectors
-**     iauPv2s      pv-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, ph[3], r5h[3][3], s5h[3], sh[3], vst[3],
-   rst[3][3], r5ht[3][3], pv5e[2][3], vv[3],
-   w, r, v;
-
-
-/* Time interval from fundamental epoch J2000.0 to given date (JY). */
-   t = ((date1 - DJ00) + date2) / DJY;
-
-/* Hipparcos barycentric position vector (normalized). */
-   iauS2c(rh, dh, ph);
-
-/* FK5 to Hipparcos orientation matrix and spin vector. */
-   iauFk5hip(r5h, s5h);
-
-/* Rotate the spin into the Hipparcos system. */
-   iauRxp(r5h, s5h, sh);
-
-/* Accumulated Hipparcos wrt FK5 spin over that interval. */
-   iauSxp(t, s5h, vst);
-
-/* Express the accumulated spin as a rotation matrix. */
-   iauRv2m(vst, rst);
-
-/* Rotation matrix:  accumulated spin, then FK5 to Hipparcos. */
-   iauRxr(r5h, rst, r5ht);
-
-/* De-orient & de-spin the Hipparcos position into FK5 J2000.0. */
-   iauTrxp(r5ht, ph, pv5e[0]);
-
-/* Apply spin to the position giving a space motion. */
-   iauPxp(sh, ph, vv);
-
-/* De-orient & de-spin the Hipparcos space motion into FK5 J2000.0. */
-   iauTrxp(r5ht, vv, pv5e[1]);
-
-/* FK5 position/velocity pv-vector to spherical. */
-   iauPv2s(pv5e, &w, d5, &r, dr5, dd5, &v);
-   *r5 = 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/ir.c

@@ -1,83 +0,0 @@
-#include "sofam.h"
-
-void iauIr(double r[3][3])
-/*
-**  - - - - - -
-**   i a u I r
-**  - - - - - -
-**
-**  Initialize an r-matrix to the identity matrix.
-**
-**  Status:  vector/matrix support function.
-**
-**  Returned:
-**     r       double[3][3]    r-matrix
-**
-**  This revision:  2012 April 3
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   r[0][0] = 1.0;
-   r[0][1] = 0.0;
-   r[0][2] = 0.0;
-   r[1][0] = 0.0;
-   r[1][1] = 1.0;
-   r[1][2] = 0.0;
-   r[2][0] = 0.0;
-   r[2][1] = 0.0;
-   r[2][2] = 1.0;
-
-   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/jd2cal.c

@@ -1,155 +0,0 @@
-#include "sofam.h"
-
-int iauJd2cal(double dj1, double dj2,
-              int *iy, int *im, int *id, double *fd)
-/*
-**  - - - - - - - - - -
-**   i a u J d 2 c a l
-**  - - - - - - - - - -
-**
-**  Julian Date to Gregorian year, month, day, and fraction of a day.
-**
-**  Status:  support function.
-**
-**  Given:
-**     dj1,dj2   double   Julian Date (Notes 1, 2)
-**
-**  Returned (arguments):
-**     iy        int      year
-**     im        int      month
-**     id        int      day
-**     fd        double   fraction of day
-**
-**  Returned (function value):
-**               int      status:
-**                           0 = OK
-**                          -1 = unacceptable date (Note 3)
-**
-**  Notes:
-**
-**  1) The earliest valid date is -68569.5 (-4900 March 1).  The
-**     largest value accepted is 10^9.
-**
-**  2) The Julian Date is apportioned in any convenient way between
-**     the arguments dj1 and dj2.  For example, JD=2450123.7 could
-**     be expressed in any of these ways, among others:
-**
-**            dj1             dj2
-**
-**         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)
-**
-**  3) In early eras the conversion is from the "proleptic Gregorian
-**     calendar";  no account is taken of the date(s) of adoption of
-**     the Gregorian calendar, nor is the AD/BC numbering convention
-**     observed.
-**
-**  Reference:
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992),
-**     Section 12.92 (p604).
-**
-**  This revision:  2008 May 26
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-/* Minimum and maximum allowed JD */
-   static const double djmin = -68569.5;
-   static const double djmax = 1e9;
-
-   long jd, l, n, i, k;
-   double dj, d1, d2, f1, f2, f, d;
-
-
-/* Verify date is acceptable. */
-   dj = dj1 + dj2;
-   if (dj < djmin || dj > djmax) return -1;
-
-/* Copy the date, big then small, and re-align to midnight. */
-   if (dj1 >= dj2) {
-      d1 = dj1;
-      d2 = dj2;
-   } else {
-      d1 = dj2;
-      d2 = dj1;
-   }
-   d2 -= 0.5;
-
-/* Separate day and fraction. */
-   f1 = fmod(d1, 1.0);
-   f2 = fmod(d2, 1.0);
-   f = fmod(f1 + f2, 1.0);
-   if (f < 0.0) f += 1.0;
-   d = floor(d1 - f1) + floor(d2 - f2) + floor(f1 + f2 - f);
-   jd = (long) floor(d) + 1L;
-
-/* Express day in Gregorian calendar. */
-   l = jd + 68569L;
-   n = (4L * l) / 146097L;
-   l -= (146097L * n + 3L) / 4L;
-   i = (4000L * (l + 1L)) / 1461001L;
-   l -= (1461L * i) / 4L - 31L;
-   k = (80L * l) / 2447L;
-   *id = (int) (l - (2447L * k) / 80L);
-   l = k / 11L;
-   *im = (int) (k + 2L - 12L * l);
-   *iy = (int) (100L * (n - 49L) + i + l);
-   *fd = f;
-
-   return 0;
-
-/*----------------------------------------------------------------------
-**
-**  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.
-**
-**--------------------------------------------------------------------*/
-}