equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/pvxpv.c

@@ -1,107 +0,0 @@
-#include "sofam.h"
-
-void iauPvxpv(double a[2][3], double b[2][3], double axb[2][3])
-/*
-**  - - - - - - - - -
-**   i a u P v x p v
-**  - - - - - - - - -
-**
-**  Outer (=vector=cross) product of two pv-vectors.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     a        double[2][3]      first pv-vector
-**     b        double[2][3]      second pv-vector
-**
-**  Returned:
-**     axb      double[2][3]      a x b
-**
-**  Notes:
-**
-**  1) If the position and velocity components of the two pv-vectors are
-**     ( ap, av ) and ( bp, bv ), the result, a x b, is the pair of
-**     vectors ( ap x bp, ap x bv + av x bp ).  The two vectors are the
-**     cross-product of the two p-vectors and its derivative.
-**
-**  2) It is permissible to re-use the same array for any of the
-**     arguments.
-**
-**  Called:
-**     iauCpv       copy pv-vector
-**     iauPxp       vector product of two p-vectors
-**     iauPpp       p-vector plus p-vector
-**
-**  This revision:  2008 November 18
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double wa[2][3], wb[2][3], axbd[3], adxb[3];
-
-
-/* Make copies of the inputs. */
-   iauCpv(a, wa);
-   iauCpv(b, wb);
-
-/* a x b = position part of result. */
-   iauPxp(wa[0], wb[0], axb[0]);
-
-/* a x bdot + adot x b = velocity part of result. */
-   iauPxp(wa[0], wb[1], axbd);
-   iauPxp(wa[1], wb[0], adxb);
-   iauPpp(axbd, adxb, axb[1]);
-
-   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/pxp.c

@@ -1,94 +0,0 @@
-#include "sofam.h"
-
-void iauPxp(double a[3], double b[3], double axb[3])
-/*
-**  - - - - - - -
-**   i a u P x p
-**  - - - - - - -
-**
-**  p-vector outer (=vector=cross) product.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     a        double[3]      first p-vector
-**     b        double[3]      second p-vector
-**
-**  Returned:
-**     axb      double[3]      a x b
-**
-**  Note:
-**     It is permissible to re-use the same array for any of the
-**     arguments.
-**
-**  This revision:  2008 November 18
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double xa, ya, za, xb, yb, zb;
-
-
-   xa = a[0];
-   ya = a[1];
-   za = a[2];
-   xb = b[0];
-   yb = b[1];
-   zb = b[2];
-   axb[0] = ya*zb - za*yb;
-   axb[1] = za*xb - xa*zb;
-   axb[2] = xa*yb - ya*xb;
-
-   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/rm2v.c

@@ -1,111 +0,0 @@
-#include "sofam.h"
-
-void iauRm2v(double r[3][3], double w[3])
-/*
-**  - - - - - - - -
-**   i a u R m 2 v
-**  - - - - - - - -
-**
-**  Express an r-matrix as an r-vector.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     r        double[3][3]    rotation matrix
-**
-**  Returned:
-**     w        double[3]       rotation vector (Note 1)
-**
-**  Notes:
-**
-**  1) A rotation matrix describes a rotation through some angle about
-**     some arbitrary axis called the Euler axis.  The "rotation vector"
-**     returned by this function has the same direction as the Euler axis,
-**     and its magnitude is the angle in radians.  (The magnitude and
-**     direction can be separated by means of the function iauPn.)
-**
-**  2) If r is null, so is the result.  If r is not a rotation matrix
-**     the result is undefined;  r must be proper (i.e. have a positive
-**     determinant) and real orthogonal (inverse = transpose).
-**
-**  3) The reference frame rotates clockwise as seen looking along
-**     the rotation vector from the origin.
-**
-**  This revision:  2010 August 27
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double x, y, z, s2, c2, phi, f;
-
-
-   x = r[1][2] - r[2][1];
-   y = r[2][0] - r[0][2];
-   z = r[0][1] - r[1][0];
-   s2 = sqrt(x*x + y*y + z*z);
-   if (s2 != 0) {
-      c2 = r[0][0] + r[1][1] + r[2][2] - 1.0;
-      phi = atan2(s2, c2);
-      f =  phi / s2;
-      w[0] = x * f;
-      w[1] = y * f;
-      w[2] = z * f;
-   } else {
-      w[0] = 0.0;
-      w[1] = 0.0;
-      w[2] = 0.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/rv2m.c

@@ -1,118 +0,0 @@
-#include "sofam.h"
-
-void iauRv2m(double w[3], double r[3][3])
-/*
-**  - - - - - - - -
-**   i a u R v 2 m
-**  - - - - - - - -
-**
-**  Form the r-matrix corresponding to a given r-vector.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     w        double[3]      rotation vector (Note 1)
-**
-**  Returned:
-**     r        double[3][3]    rotation matrix
-**
-**  Notes:
-**
-**  1) A rotation matrix describes a rotation through some angle about
-**     some arbitrary axis called the Euler axis.  The "rotation vector"
-**     supplied to This function has the same direction as the Euler
-**     axis, and its magnitude is the angle in radians.
-**
-**  2) If w is null, the unit matrix is returned.
-**
-**  3) The reference frame rotates clockwise as seen looking along the
-**     rotation vector from the origin.
-**
-**  This revision:  2008 May 11
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double x, y, z, phi, s, c, f;
-
-
-/* Euler angle (magnitude of rotation vector) and functions. */
-   x = w[0];
-   y = w[1];
-   z = w[2];
-   phi = sqrt(x*x + y*y + z*z);
-   s = sin(phi);
-   c = cos(phi);
-   f = 1.0 - c;
-
-/* Euler axis (direction of rotation vector), perhaps null. */
-   if (phi != 0.0) {
-       x /= phi;
-       y /= phi;
-       z /= phi;
-   }
-
-/* Form the rotation matrix. */
-   r[0][0] = x*x*f + c;
-   r[0][1] = x*y*f + z*s;
-   r[0][2] = x*z*f - y*s;
-   r[1][0] = y*x*f - z*s;
-   r[1][1] = y*y*f + c;
-   r[1][2] = y*z*f + x*s;
-   r[2][0] = z*x*f + y*s;
-   r[2][1] = z*y*f - x*s;
-   r[2][2] = z*z*f + c;
-
-   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/rx.c

@@ -1,110 +0,0 @@
-#include "sofam.h"
-
-void iauRx(double phi, double r[3][3])
-/*
-**  - - - - - -
-**   i a u R x
-**  - - - - - -
-**
-**  Rotate an r-matrix about the x-axis.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     phi    double          angle (radians)
-**
-**  Given and returned:
-**     r      double[3][3]    r-matrix, rotated
-**
-**  Notes:
-**
-**  1) Calling this function with positive phi incorporates in the
-**     supplied r-matrix r an additional rotation, about the x-axis,
-**     anticlockwise as seen looking towards the origin from positive x.
-**
-**  2) The additional rotation can be represented by this matrix:
-**
-**         (  1        0            0      )
-**         (                               )
-**         (  0   + cos(phi)   + sin(phi)  )
-**         (                               )
-**         (  0   - sin(phi)   + cos(phi)  )
-**
-**  This revision:  2012 April 3
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double s, c, a10, a11, a12, a20, a21, a22;
-
-
-   s = sin(phi);
-   c = cos(phi);
-
-   a10 =   c*r[1][0] + s*r[2][0];
-   a11 =   c*r[1][1] + s*r[2][1];
-   a12 =   c*r[1][2] + s*r[2][2];
-   a20 = - s*r[1][0] + c*r[2][0];
-   a21 = - s*r[1][1] + c*r[2][1];
-   a22 = - s*r[1][2] + c*r[2][2];
-
-   r[1][0] = a10;
-   r[1][1] = a11;
-   r[1][2] = a12;
-   r[2][0] = a20;
-   r[2][1] = a21;
-   r[2][2] = a22;
-
-   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.
-**
-**--------------------------------------------------------------------*/
-}