equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/pm.c

@@ -1,83 +0,0 @@
-#include "sofam.h"
-
-double iauPm(double p[3])
-/*
-**  - - - - - -
-**   i a u P m
-**  - - - - - -
-**
-**  Modulus of p-vector.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     p      double[3]     p-vector
-**
-**  Returned (function value):
-**            double        modulus
-**
-**  This revision:  2008 May 22
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double w;
-
-
-   w  = sqrt( p[0] * p[0]
-            + p[1] * p[1]
-            + p[2] * p[2] );
-
-   return w;
-
-/*----------------------------------------------------------------------
-**
-**  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/pmat00.c

@@ -1,118 +0,0 @@
-#include "sofam.h"
-
-void iauPmat00(double date1, double date2, double rbp[3][3])
-/*
-**  - - - - - - - - - -
-**   i a u P m a t 0 0
-**  - - - - - - - - - -
-**
-**  Precession matrix (including frame bias) from GCRS to a specified
-**  date, IAU 2000 model.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2  double          TT as a 2-part Julian Date (Note 1)
-**
-**  Returned:
-**     rbp          double[3][3]    bias-precession matrix (Note 2)
-**
-**  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 matrix operates in the sense V(date) = rbp * V(GCRS), where
-**     the p-vector V(GCRS) is with respect to the Geocentric Celestial
-**     Reference System (IAU, 2000) and the p-vector V(date) is with
-**     respect to the mean equatorial triad of the given date.
-**
-**  Called:
-**     iauBp00      frame bias and precession matrices, IAU 2000
-**
-**  Reference:
-**
-**     IAU: Trans. International Astronomical Union, Vol. XXIVB;  Proc.
-**     24th General Assembly, Manchester, UK.  Resolutions B1.3, B1.6.
-**     (2000)
-**
-**  This revision:  2009 December 21
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double rb[3][3], rp[3][3];
-
-
-/* Obtain the required matrix (discarding others). */
-   iauBp00(date1, date2, rb, rp, rbp);
-
-   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/pmat06.c

@@ -1,122 +0,0 @@
-#include "sofam.h"
-
-void iauPmat06(double date1, double date2, double rbp[3][3])
-/*
-**  - - - - - - - - - -
-**   i a u P m a t 0 6
-**  - - - - - - - - - -
-**
-**  Precession matrix (including frame bias) from GCRS to a specified
-**  date, IAU 2006 model.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2  double          TT as a 2-part Julian Date (Note 1)
-**
-**  Returned:
-**     rbp          double[3][3]    bias-precession matrix (Note 2)
-**
-**  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 matrix operates in the sense V(date) = rbp * V(GCRS), where
-**     the p-vector V(GCRS) is with respect to the Geocentric Celestial
-**     Reference System (IAU, 2000) and the p-vector V(date) is with
-**     respect to the mean equatorial triad of the given date.
-**
-**  Called:
-**     iauPfw06     bias-precession F-W angles, IAU 2006
-**     iauFw2m      F-W angles to r-matrix
-**
-**  References:
-**
-**     Capitaine, N. & Wallace, P.T., 2006, Astron.Astrophys. 450, 855
-**
-**     Wallace, P.T. & Capitaine, N., 2006, Astron.Astrophys. 459, 981
-**
-**  This revision:  2009 December 21
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double gamb, phib, psib, epsa;
-
-
-/* Bias-precession Fukushima-Williams angles. */
-   iauPfw06(date1, date2, &gamb, &phib, &psib, &epsa);
-
-/* Form the matrix. */
-   iauFw2m(gamb, phib, psib, epsa, rbp);
-
-   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/pmat76.c

@@ -1,141 +0,0 @@
-#include "sofam.h"
-
-void iauPmat76(double date1, double date2, double rmatp[3][3])
-/*
-**  - - - - - - - - - -
-**   i a u P m a t 7 6
-**  - - - - - - - - - -
-**
-**  Precession matrix from J2000.0 to a specified date, IAU 1976 model.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2 double       ending date, TT (Note 1)
-**
-**  Returned:
-**     rmatp       double[3][3] precession matrix, J2000.0 -> date1+date2
-**
-**  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 matrix operates in the sense V(date) = RMATP * V(J2000),
-**     where the p-vector V(J2000) is with respect to the mean
-**     equatorial triad of epoch J2000.0 and the p-vector V(date)
-**     is with respect to the mean equatorial triad of the given
-**     date.
-**
-**  3) Though the matrix method itself is rigorous, the precession
-**     angles are expressed through canonical polynomials which are
-**     valid only for a limited time span.  In addition, the IAU 1976
-**     precession rate is known to be imperfect.  The absolute accuracy
-**     of the present formulation is better than 0.1 arcsec from
-**     1960AD to 2040AD, better than 1 arcsec from 1640AD to 2360AD,
-**     and remains below 3 arcsec for the whole of the period
-**     500BC to 3000AD.  The errors exceed 10 arcsec outside the
-**     range 1200BC to 3900AD, exceed 100 arcsec outside 4200BC to
-**     5600AD and exceed 1000 arcsec outside 6800BC to 8200AD.
-**
-**  Called:
-**     iauPrec76    accumulated precession angles, IAU 1976
-**     iauIr        initialize r-matrix to identity
-**     iauRz        rotate around Z-axis
-**     iauRy        rotate around Y-axis
-**     iauCr        copy r-matrix
-**
-**  References:
-**
-**     Lieske, J.H., 1979, Astron.Astrophys. 73, 282.
-**      equations (6) & (7), p283.
-**
-**     Kaplan,G.H., 1981. USNO circular no. 163, pA2.
-**
-**  This revision:  2009 December 18
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double zeta, z, theta, wmat[3][3];
-
-
-/* Precession Euler angles, J2000.0 to specified date. */
-   iauPrec76(DJ00, 0.0, date1, date2, &zeta, &z, &theta);
-
-/* Form the rotation matrix. */
-   iauIr(  wmat);
-   iauRz( -zeta, wmat);
-   iauRy(  theta, wmat);
-   iauRz( -z, wmat);
-   iauCr( wmat, rmatp);
-
-   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.
-**
-**--------------------------------------------------------------------*/
-}