equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/p2s.c

@@ -1,91 +0,0 @@
-#include "sofam.h"
-
-void iauP2s(double p[3], double *theta, double *phi, double *r)
-/*
-**  - - - - - - -
-**   i a u P 2 s
-**  - - - - - - -
-**
-**  P-vector to spherical polar coordinates.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     p        double[3]    p-vector
-**
-**  Returned:
-**     theta    double       longitude angle (radians)
-**     phi      double       latitude angle (radians)
-**     r        double       radial distance
-**
-**  Notes:
-**
-**  1) If P is null, zero theta, phi and r are returned.
-**
-**  2) At either pole, zero theta is returned.
-**
-**  Called:
-**     iauC2s       p-vector to spherical
-**     iauPm        modulus of p-vector
-**
-**  This revision:  2008 May 22
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   iauC2s(p, theta, phi);
-   *r = iauPm(p);
-
-   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/pap.c

@@ -1,139 +0,0 @@
-#include "sofam.h"
-
-double iauPap(double a[3], double b[3])
-/*
-**  - - - - - - -
-**   i a u P a p
-**  - - - - - - -
-**
-**  Position-angle from two p-vectors.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     a      double[3]  direction of reference point
-**     b      double[3]  direction of point whose PA is required
-**
-**  Returned (function value):
-**            double     position angle of b with respect to a (radians)
-**
-**  Notes:
-**
-**  1) The result is the position angle, in radians, of direction b with
-**     respect to direction a.  It is in the range -pi to +pi.  The
-**     sense is such that if b is a small distance "north" of a the
-**     position angle is approximately zero, and if b is a small
-**     distance "east" of a the position angle is approximately +pi/2.
-**
-**  2) The vectors a and b need not be of unit length.
-**
-**  3) Zero is returned if the two directions are the same or if either
-**     vector is null.
-**
-**  4) If vector a is at a pole, the result is ill-defined.
-**
-**  Called:
-**     iauPn        decompose p-vector into modulus and direction
-**     iauPm        modulus of p-vector
-**     iauPxp       vector product of two p-vectors
-**     iauPmp       p-vector minus p-vector
-**     iauPdp       scalar product of two p-vectors
-**
-**  This revision:  2008 May 25
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double am, au[3], bm, st, ct, xa, ya, za, eta[3], xi[3], a2b[3], pa;
-
-
-/* Modulus and direction of the a vector. */
-   iauPn(a, &am, au);
-
-/* Modulus of the b vector. */
-   bm = iauPm(b);
-
-/* Deal with the case of a null vector. */
-   if ((am == 0.0) || (bm == 0.0)) {
-      st = 0.0;
-      ct = 1.0;
-   } else {
-
-   /* The "north" axis tangential from a (arbitrary length). */
-      xa = a[0];
-      ya = a[1];
-      za = a[2];
-      eta[0] = -xa * za;
-      eta[1] = -ya * za;
-      eta[2] =  xa*xa + ya*ya;
-
-   /* The "east" axis tangential from a (same length). */
-      iauPxp(eta, au, xi);
-
-   /* The vector from a to b. */
-      iauPmp(b, a, a2b);
-
-   /* Resolve into components along the north and east axes. */
-      st = iauPdp(a2b, xi);
-      ct = iauPdp(a2b, eta);
-
-   /* Deal with degenerate cases. */
-      if ((st == 0.0) && (ct == 0.0)) ct = 1.0;
-   }
-
-/* Position angle. */
-   pa = atan2(st, ct);
-
-   return pa;
-
-/*----------------------------------------------------------------------
-**
-**  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/pas.c

@@ -1,96 +0,0 @@
-#include "sofam.h"
-
-double iauPas(double al, double ap, double bl, double bp)
-/*
-**  - - - - - - -
-**   i a u P a s
-**  - - - - - - -
-**
-**  Position-angle from spherical coordinates.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     al     double     longitude of point A (e.g. RA) in radians
-**     ap     double     latitude of point A (e.g. Dec) in radians
-**     bl     double     longitude of point B
-**     bp     double     latitude of point B
-**
-**  Returned (function value):
-**            double     position angle of B with respect to A
-**
-**  Notes:
-**
-**  1) The result is the bearing (position angle), in radians, of point
-**     B with respect to point A.  It is in the range -pi to +pi.  The
-**     sense is such that if B is a small distance "east" of point A,
-**     the bearing is approximately +pi/2.
-**
-**  2) Zero is returned if the two points are coincident.
-**
-**  This revision:  2008 May 22
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double dl, x, y, pa;
-
-
-   dl = bl - al;
-   y = sin(dl) * cos(bp);
-   x = sin(bp) * cos(ap) - cos(bp) * sin(ap) * cos(dl);
-   pa = ((x != 0.0) || (y != 0.0)) ? atan2(y, x) : 0.0;
-
-   return pa;
-
-/*----------------------------------------------------------------------
-**
-**  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/pb06.c

@@ -1,144 +0,0 @@
-#include "sofam.h"
-
-void iauPb06(double date1, double date2,
-             double *bzeta, double *bz, double *btheta)
-/*
-**  - - - - - - - -
-**   i a u P b 0 6
-**  - - - - - - - -
-**
-**  This function forms three Euler angles which implement general
-**  precession from epoch J2000.0, using the IAU 2006 model.  Frame
-**  bias (the offset between ICRS and mean J2000.0) is included.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2  double   TT as a 2-part Julian Date (Note 1)
-**
-**  Returned:
-**     bzeta        double   1st rotation: radians cw around z
-**     bz           double   3rd rotation: radians cw around z
-**     btheta       double   2nd rotation: radians ccw around y
-**
-**  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 traditional accumulated precession angles zeta_A, z_A,
-**     theta_A cannot be obtained in the usual way, namely through
-**     polynomial expressions, because of the frame bias.  The latter
-**     means that two of the angles undergo rapid changes near this
-**     date.  They are instead the results of decomposing the
-**     precession-bias matrix obtained by using the Fukushima-Williams
-**     method, which does not suffer from the problem.  The
-**     decomposition returns values which can be used in the
-**     conventional formulation and which include frame bias.
-**
-**  3) The three angles are returned in the conventional order, which
-**     is not the same as the order of the corresponding Euler
-**     rotations.  The precession-bias matrix is
-**     R_3(-z) x R_2(+theta) x R_3(-zeta).
-**
-**  4) Should zeta_A, z_A, theta_A angles be required that do not
-**     contain frame bias, they are available by calling the SOFA
-**     function iauP06e.
-**
-**  Called:
-**     iauPmat06    PB matrix, IAU 2006
-**     iauRz        rotate around Z-axis
-**
-**  This revision:  2008 May 26
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double r[3][3], r31, r32;
-
-
-/* Precession matrix via Fukushima-Williams angles. */
-   iauPmat06(date1, date2, r);
-
-/* Solve for z. */
-   *bz = atan2(r[1][2], r[0][2]);
-
-/* Remove it from the matrix. */
-   iauRz(*bz, r);
-
-/* Solve for the remaining two angles. */
-   *bzeta = atan2 (r[1][0], r[1][1]);
-   r31 = r[2][0];
-   r32 = r[2][1];
-   *btheta = atan2(-dsign(sqrt(r31 * r31 + r32 * r32), r[0][2]),
-                   r[2][2]);
-
-   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.
-**
-**--------------------------------------------------------------------*/
-}