equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/tdbtcb.c

@@ -1,137 +0,0 @@
-#include "sofam.h"
-
-int iauTdbtcb(double tdb1, double tdb2, double *tcb1, double *tcb2)
-/*
-**  - - - - - - - - - -
-**   i a u T d b t c b
-**  - - - - - - - - - -
-**
-**  Time scale transformation:  Barycentric Dynamical Time, TDB, to
-**  Barycentric Coordinate Time, TCB.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     tdb1,tdb2  double    TDB as a 2-part Julian Date
-**
-**  Returned:
-**     tcb1,tcb2  double    TCB as a 2-part Julian Date
-**
-**  Returned (function value):
-**                int       status:  0 = OK
-**
-**  Notes:
-**
-**  1) tdb1+tdb2 is Julian Date, apportioned in any convenient way
-**     between the two arguments, for example where tdb1 is the Julian
-**     Day Number and tdb2 is the fraction of a day.  The returned
-**     tcb1,tcb2 follow suit.
-**
-**  2) The 2006 IAU General Assembly introduced a conventional linear
-**     transformation between TDB and TCB.  This transformation
-**     compensates for the drift between TCB and terrestrial time TT,
-**     and keeps TDB approximately centered on TT.  Because the
-**     relationship between TT and TCB depends on the adopted solar
-**     system ephemeris, the degree of alignment between TDB and TT over
-**     long intervals will vary according to which ephemeris is used.
-**     Former definitions of TDB attempted to avoid this problem by
-**     stipulating that TDB and TT should differ only by periodic
-**     effects.  This is a good description of the nature of the
-**     relationship but eluded precise mathematical formulation.  The
-**     conventional linear relationship adopted in 2006 sidestepped
-**     these difficulties whilst delivering a TDB that in practice was
-**     consistent with values before that date.
-**
-**  3) TDB is essentially the same as Teph, the time argument for the
-**     JPL solar system ephemerides.
-**
-**  Reference:
-**
-**     IAU 2006 Resolution B3
-**
-**  This revision:  2011 May 14
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-
-/* 1977 Jan 1 00:00:32.184 TT, as two-part JD */
-   static const double t77td = DJM0 + DJM77;
-   static const double t77tf = TTMTAI/DAYSEC;
-
-/* TDB (days) at TAI 1977 Jan 1.0 */
-   static const double tdb0 = TDB0/DAYSEC;
-
-/* TDB to TCB rate */
-   static const double elbb = ELB/(1.0-ELB);
-
-   double d, f;
-
-
-/* Result, preserving date format but safeguarding precision. */
-   if ( tdb1 > tdb2 ) {
-      d = t77td - tdb1;
-      f  = tdb2 - tdb0;
-      *tcb1 = tdb1;
-      *tcb2 = f - ( d - ( f - t77tf ) ) * elbb;
-   } else {
-      d = t77td - tdb2;
-      f  = tdb1 - tdb0;
-      *tcb1 = f + ( d - ( f - t77tf ) ) * elbb;
-      *tcb2 = tdb2;
-   }
-
-/* Status (always OK). */
-   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.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/tdbtt.c

@@ -1,122 +0,0 @@
-#include "sofam.h"
-
-int iauTdbtt(double tdb1, double tdb2, double dtr,
-             double *tt1, double *tt2 )
-/*
-**  - - - - - - - - -
-**   i a u T d b t t
-**  - - - - - - - - -
-**
-**  Time scale transformation:  Barycentric Dynamical Time, TDB, to
-**  Terrestrial Time, TT.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     tdb1,tdb2  double    TDB as a 2-part Julian Date
-**     dtr        double    TDB-TT in seconds
-**
-**  Returned:
-**     tt1,tt2    double    TT as a 2-part Julian Date
-**
-**  Returned (function value):
-**                int       status:  0 = OK
-**
-**  Notes:
-**
-**  1) tdb1+tdb2 is Julian Date, apportioned in any convenient way
-**     between the two arguments, for example where tdb1 is the Julian
-**     Day Number and tdb2 is the fraction of a day.  The returned
-**     tt1,tt2 follow suit.
-**
-**  2) The argument dtr represents the quasi-periodic component of the
-**     GR transformation between TT and TCB.  It is dependent upon the
-**     adopted solar-system ephemeris, and can be obtained by numerical
-**     integration, by interrogating a precomputed time ephemeris or by
-**     evaluating a model such as that implemented in the SOFA function
-**     iauDtdb.   The quantity is dominated by an annual term of 1.7 ms
-**     amplitude.
-**
-**  3) TDB is essentially the same as Teph, the time argument for the
-**     JPL solar system ephemerides.
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     IAU 2006 Resolution 3
-**
-**  This revision:  2011 May 14
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-**
-*/
-{
-   double dtrd;
-
-
-/* Result, safeguarding precision. */
-   dtrd = dtr / DAYSEC;
-   if ( tdb1 > tdb2 ) {
-      *tt1 = tdb1;
-      *tt2 = tdb2 - dtrd;
-   } else {
-      *tt1 = tdb1 - dtrd;
-      *tt2 = tdb2;
-   }
-
-/* Status (always OK). */
-   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.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/test_celes.rb

@@ -1,48 +0,0 @@
-#require_relative 'celes_core'
-require 'celes'
-
-lib = File.expand_path('../../../lib', __FILE__)
-$LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
-
-require 'eot'
-
-eot = Eot.new
-
-# Celes.methods.sort.each {|m| p m }
-
-dt = DateTime.now.to_time.utc.to_datetime
-jd = Celes.cal2jd(dt.year, dt.month, dt.day)
-
-ajd = jd[0] + jd[1] + dt.day_fraction
-
-p eot.ajd = ajd
-
-p eot.delta_equinox()[0]
-p Celes.nut06a(ajd,0)[0]
-
-p eot.delta_equinox()[1]
-p Celes.nut06a(ajd,0)[1]
-
-p eot.delta_equinox()[2]
-p Celes.falp03(eot.ta)
-
-p eot.delta_equinox()[3]
-p Celes.faom03(eot.ta)
-
-require 'benchmark'
-
-n = 5_000
-Benchmark.bm do |x|
-  #~ x.report("celes")  { n.times { Celes.cal2jd(dt.year, dt.month, dt.day) } }
-  #~ x.report("ruby")  { n.times { DateTime.new(dt.year, dt.month, dt.day).jd } }
-  x.report("celes")  { n.times { Celes.nut06a(ajd,0)[0] }}
-  x.report("celes")  { n.times { Celes.nut06a(ajd,0)[1] }}
-  x.report("celes")  { n.times { Celes.falp03(eot.ta) }}
-  x.report("celes")  { n.times { Celes.faom03(eot.ta) }}
-
-  x.report("ruby")  { n.times { eot.delta_equinox()[0] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[1] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[2] }}
-  x.report("ruby")  { n.times { eot.delta_equinox()[3] }}
-end
-

data/examples/celes_parts/tf2a.c

@@ -1,107 +0,0 @@
-#include "sofam.h"
-#include <stdlib.h>
-
-int iauTf2a(char s, int ihour, int imin, double sec, double *rad)
-/*
-**  - - - - - - - -
-**   i a u T f 2 a
-**  - - - - - - - -
-**
-**  Convert hours, minutes, seconds to radians.
-**
-**  Status:  support function.
-**
-**  Given:
-**     s         char    sign:  '-' = negative, otherwise positive
-**     ihour     int     hours
-**     imin      int     minutes
-**     sec       double  seconds
-**
-**  Returned:
-**     rad       double  angle in radians
-**
-**  Returned (function value):
-**               int     status:  0 = OK
-**                                1 = ihour outside range 0-23
-**                                2 = imin outside range 0-59
-**                                3 = sec outside range 0-59.999...
-**
-**  Notes:
-**
-**  1)  The result is computed even if any of the range checks fail.
-**
-**  2)  Negative ihour, imin and/or sec produce a warning status, but
-**      the absolute value is used in the conversion.
-**
-**  3)  If there are multiple errors, the status value reflects only the
-**      first, the smallest taking precedence.
-**
-**  This revision:  2012 February 13
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-
-/* Compute the interval. */
-   *rad  = ( s == '-' ? -1.0 : 1.0 ) *
-           ( 60.0 * ( 60.0 * ( (double) abs(ihour) ) +
-                             ( (double) abs(imin) ) ) +
-                                        fabs(sec) ) * DS2R;
-
-/* Validate arguments and return status. */
-   if ( ihour < 0 || ihour > 23 ) return 1;
-   if ( imin < 0 || imin > 59 ) return 2;
-   if ( sec < 0.0 || sec >= 60.0 ) return 3;
-   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.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/tf2d.c

@@ -1,107 +0,0 @@
-#include "sofam.h"
-#include <stdlib.h>
-
-int iauTf2d(char s, int ihour, int imin, double sec, double *days)
-/*
-**  - - - - - - - -
-**   i a u T f 2 d
-**  - - - - - - - -
-**
-**  Convert hours, minutes, seconds to days.
-**
-**  Status:  support function.
-**
-**  Given:
-**     s         char    sign:  '-' = negative, otherwise positive
-**     ihour     int     hours
-**     imin      int     minutes
-**     sec       double  seconds
-**
-**  Returned:
-**     days      double  interval in days
-**
-**  Returned (function value):
-**               int     status:  0 = OK
-**                                1 = ihour outside range 0-23
-**                                2 = imin outside range 0-59
-**                                3 = sec outside range 0-59.999...
-**
-**  Notes:
-**
-**  1)  The result is computed even if any of the range checks fail.
-**
-**  2)  Negative ihour, imin and/or sec produce a warning status, but
-**      the absolute value is used in the conversion.
-**
-**  3)  If there are multiple errors, the status value reflects only the
-**      first, the smallest taking precedence.
-**
-**  This revision:  2012 February 13
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-
-/* Compute the interval. */
-   *days  = ( s == '-' ? -1.0 : 1.0 ) *
-            ( 60.0 * ( 60.0 * ( (double) abs(ihour) ) +
-                              ( (double) abs(imin) ) ) +
-                                         fabs(sec) ) / DAYSEC;
-
-/* Validate arguments and return status. */
-   if ( ihour < 0 || ihour > 23 ) return 1;
-   if ( imin < 0 || imin > 59 ) return 2;
-   if ( sec < 0.0 || sec >= 60.0 ) return 3;
-   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.
-**
-**--------------------------------------------------------------------*/
-}