equationoftime 3.0.0 → 4.0.0

data/examples/celes_parts/cal2jd.c

@@ -0,0 +1,139 @@
+#include "sofam.h"
+
+int iauCal2jd(int iy, int im, int id, double *djm0, double *djm)
+/*
+**  - - - - - - - - - -
+**   i a u C a l 2 j d
+**  - - - - - - - - - -
+**
+**  Gregorian Calendar to Julian Date.
+**
+**  Status:  support function.
+**
+**  Given:
+**     iy,im,id  int     year, month, day in Gregorian calendar (Note 1)
+**
+**  Returned:
+**     djm0      double  MJD zero-point: always 2400000.5
+**     djm       double  Modified Julian Date for 0 hrs
+**
+**  Returned (function value):
+**               int     status:
+**                           0 = OK
+**                          -1 = bad year   (Note 3: JD not computed)
+**                          -2 = bad month  (JD not computed)
+**                          -3 = bad day    (JD computed)
+**
+**  Notes:
+**
+**  1) The algorithm used is valid from -4800 March 1, but this
+**     implementation rejects dates before -4799 January 1.
+**
+**  2) The Julian Date is returned in two pieces, in the usual SOFA
+**     manner, which is designed to preserve time resolution.  The
+**     Julian Date is available as a single number by adding djm0 and
+**     djm.
+**
+**  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:  2009 October 19
+**
+**  Original version 2012-03-01
+**
+**  Copyright (C) 2013 Naoki Arita.  See notes at end.
+*/
+{
+   int j, ly, my;
+   long iypmy;
+
+/* Earliest year allowed (4800BC) */
+   const int IYMIN = -4799;
+
+/* Month lengths in days */
+   static const int mtab[]
+                     = {31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
+
+
+/* Preset status. */
+   j = 0;
+
+/* Validate year and month. */
+   if (iy < IYMIN) return -1;
+   if (im < 1 || im > 12) return -2;
+
+/* If February in a leap year, 1, otherwise 0. */
+   ly = ((im == 2) && !(iy%4) && (iy%100 || !(iy%400)));
+
+/* Validate day, taking into account leap years. */
+   if ( (id < 1) || (id > (mtab[im-1] + ly))) j = -3;
+
+/* Return result. */
+   my = (im - 14) / 12;
+   iypmy = (long) (iy + my);
+   *djm0 = 2400000.5;
+   *djm = (double)((1461L * (iypmy + 4800L)) / 4L
+                 + (367L * (long) (im - 2 - 12 * my)) / 12L
+                 - (3L * ((iypmy + 4900L) / 100L)) / 4L
+                 + (long) id - 2432076L);
+
+/* Return status. */
+   return j;
+
+/*----------------------------------------------------------------------
+**
+**  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/celes_core.c

@@ -0,0 +1,2522 @@
+#include "sofa.h"
+#include "sofam.h"
+#include "ruby.h"
+
+#ifndef DBL2NUM
+# define DBL2NUM(dbl) rb_float_new(dbl)
+#endif
+
+VALUE mCeles;
+ID id_status;
+
+static VALUE
+celes_s_status(VALUE self){
+	return rb_ivar_get(self, id_status);
+}
+
+static void
+cvec_cp(VALUE vvector, double vector[3]){
+	int i;
+
+	Check_Type(vvector, T_ARRAY);
+	if(RARRAY_LEN(vvector) != 3)
+		rb_raise(rb_eArgError, "vector size is not 3");
+
+	for(i = 0; i < 3; i++)
+		vector[i] = NUM2DBL(RARRAY_PTR(vvector)[i]);
+}
+
+static void
+cmat_cp23(VALUE vmatrix, double matrix[2][3]){
+	int i, j;
+
+	Check_Type(vmatrix, T_ARRAY);
+	if(RARRAY_LEN(vmatrix) != 2)
+		rb_raise(rb_eArgError, "matrix size is not 2x3");
+
+	for(i = 0; i < 2; i++){
+		Check_Type(RARRAY_PTR(vmatrix)[i], T_ARRAY);
+		if(RARRAY_LEN(RARRAY_PTR(vmatrix)[i]) != 3)
+			rb_raise(rb_eArgError, "matrix size is not 2x3");
+		
+		for(j = 0; j < 3; j++)
+			matrix[i][j] = NUM2DBL(RARRAY_PTR(RARRAY_PTR(vmatrix)[i])[j]);
+	}
+}
+
+static void
+cmat_cp33(VALUE vmatrix, double matrix[3][3]){
+	int i, j;
+
+	Check_Type(vmatrix, T_ARRAY);
+	if(RARRAY_LEN(vmatrix) != 3)
+		rb_raise(rb_eArgError, "matrix size is not 3x3");
+
+	for(i = 0; i < 3; i++){
+		Check_Type(RARRAY_PTR(vmatrix)[i], T_ARRAY);
+		if(RARRAY_LEN(RARRAY_PTR(vmatrix)[i]) != 3)
+			rb_raise(rb_eArgError, "matrix size is not 3x3");
+		
+		for(j = 0; j < 3; j++)
+			matrix[i][j] = NUM2DBL(RARRAY_PTR(RARRAY_PTR(vmatrix)[i])[j]);
+	}
+}
+
+static VALUE
+vvec_cp(double vector[3], VALUE vvector){
+	int i;
+
+	Check_Type(vvector, T_ARRAY);
+	if(RARRAY_LEN(vvector) != 3)
+		rb_raise(rb_eArgError, "vector size is not 3");
+
+	for(i = 0; i < 3; i++)
+		RARRAY_PTR(vvector)[i] = DBL2NUM(vector[i]);
+
+	return vvector;
+}
+
+static VALUE
+vmat_cp23(double matrix[2][3], VALUE vmatrix){
+	int i, j;
+
+	Check_Type(vmatrix, T_ARRAY);
+	if(RARRAY_LEN(vmatrix) != 2)
+		rb_raise(rb_eArgError, "matrix size is not 2x3");
+
+	for(i = 0; i < 2; i++){
+		Check_Type(RARRAY_PTR(vmatrix)[i], T_ARRAY);
+		if(RARRAY_LEN(RARRAY_PTR(vmatrix)[i]) != 3)
+			rb_raise(rb_eArgError, "matrix size is not 2x3");
+		
+		for(j = 0; j < 3; j++)
+			RARRAY_PTR(RARRAY_PTR(vmatrix)[i])[j] = DBL2NUM(matrix[i][j]);
+	}
+
+	return vmatrix;
+}
+
+static VALUE
+vmat_cp33(double matrix[3][3], VALUE vmatrix){
+	int i, j;
+
+	Check_Type(vmatrix, T_ARRAY);
+	if(RARRAY_LEN(vmatrix) != 3)
+		rb_raise(rb_eArgError, "matrix size is not 3x3");
+
+	for(i = 0; i < 3; i++){
+		Check_Type(RARRAY_PTR(vmatrix)[i], T_ARRAY);
+		if(RARRAY_LEN(RARRAY_PTR(vmatrix)[i]) != 3)
+			rb_raise(rb_eArgError, "matrix size is not 3x3");
+		
+		for(j = 0; j < 3; j++)
+			RARRAY_PTR(RARRAY_PTR(vmatrix)[i])[j] = DBL2NUM(matrix[i][j]);
+	}
+
+	return vmatrix;
+}
+
+static VALUE
+vvec_new(double vector[3]){
+	return rb_ary_new3(3,
+			DBL2NUM(vector[0]),
+			DBL2NUM(vector[1]),
+			DBL2NUM(vector[2]));
+}
+
+static VALUE
+vmat_new23(double matrix[2][3]){
+	return rb_ary_new3(2,
+			rb_ary_new3(3,
+					DBL2NUM(matrix[0][0]),
+					DBL2NUM(matrix[0][1]),
+					DBL2NUM(matrix[0][2])),
+			rb_ary_new3(3,
+					DBL2NUM(matrix[1][0]),
+					DBL2NUM(matrix[1][1]),
+					DBL2NUM(matrix[1][2])));
+}
+
+static VALUE
+vmat_new33(double matrix[3][3]){
+	return rb_ary_new3(3,
+			rb_ary_new3(3,
+					DBL2NUM(matrix[0][0]),
+					DBL2NUM(matrix[0][1]),
+					DBL2NUM(matrix[0][2])),
+			rb_ary_new3(3,
+					DBL2NUM(matrix[1][0]),
+					DBL2NUM(matrix[1][1]),
+					DBL2NUM(matrix[1][2])),
+			rb_ary_new3(3,
+					DBL2NUM(matrix[2][0]),
+					DBL2NUM(matrix[2][1]),
+					DBL2NUM(matrix[2][2])));
+}
+
+
+
+static VALUE
+celes_s_zp_b(VALUE self, VALUE vp){
+	double p[3];
+
+	cvec_cp(vp, p);
+	iauZp(p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_cp(p, vp);
+}
+
+static VALUE
+celes_s_zp(VALUE self){
+	double p[3];
+
+	iauZp(p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(p);
+}
+
+static VALUE
+celes_s_zr_b(VALUE self, VALUE vr){
+	double r[3][3];
+
+	cmat_cp33(vr, r);
+	iauZr(r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp33(r, vr);
+}
+
+static VALUE
+celes_s_zr(VALUE self){
+	double r[3][3];
+
+	iauZr(r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_ir_b(VALUE self, VALUE vr){
+	double r[3][3];
+
+	cmat_cp33(vr, r);
+	iauIr(r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp33(r, vr);
+}
+
+static VALUE
+celes_s_ir(VALUE self){
+	double r[3][3];
+
+	iauIr(r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_cp(VALUE self, VALUE vp){
+	double p[3], c[3];
+
+	cvec_cp(vp, p);
+	iauCp(p, c);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(c);
+}
+
+static VALUE
+celes_s_cr(VALUE self, VALUE vr){
+	double r[3][3], c[3][3];
+
+	cmat_cp33(vr, r);
+	iauCr(r, c);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(c);
+}
+
+static VALUE
+celes_s_rx_b(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+	
+	cmat_cp33(vr, r);
+	iauRx(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp33(r, vr);
+}
+
+static VALUE
+celes_s_rx(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+	
+	cmat_cp33(vr, r);
+	iauRx(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_ry_b(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+
+	cmat_cp33(vr, r);
+	iauRy(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp33(r, vr);
+}
+
+static VALUE
+celes_s_ry(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+	
+	cmat_cp33(vr, r);
+	iauRy(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_rz_b(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+
+	cmat_cp33(vr, r);
+	iauRz(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp33(r, vr);
+}
+
+static VALUE
+celes_s_rz(VALUE self, VALUE vphi, VALUE vr){
+	double r[3][3];
+	
+	cmat_cp33(vr, r);
+	iauRz(NUM2DBL(vphi), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_s2c(VALUE self, VALUE vtheta, VALUE vphi){
+	double c[3];
+
+	iauS2c(NUM2DBL(vtheta), NUM2DBL(vphi), c);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(c);
+}
+
+static VALUE
+celes_s_c2s(VALUE self, VALUE vp){
+	double p[3], theta, phi;
+
+	cvec_cp(vp, p);
+	iauC2s(p, &theta, &phi);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(theta), DBL2NUM(phi));
+}
+
+static VALUE
+celes_s_s2p(VALUE self, VALUE vtheta, VALUE vphi, VALUE vr){
+	double p[3];
+
+	iauS2p(NUM2DBL(vtheta), NUM2DBL(vphi), NUM2DBL(vr), p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(p);
+}
+
+static VALUE
+celes_s_p2s(VALUE self, VALUE vp){
+	double p[3], theta, phi, r;
+
+	cvec_cp(vp, p);
+	iauP2s(p, &theta, &phi, &r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(theta), DBL2NUM(phi), DBL2NUM(r));
+}
+
+static VALUE
+celes_s_ppp(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3], apb[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	iauPpp(a, b, apb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(apb);
+}
+
+static VALUE
+celes_s_pmp(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3], amb[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	iauPmp(a, b, amb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(amb);
+}
+
+static VALUE
+celes_s_ppsp(VALUE self, VALUE va, VALUE vs, VALUE vb){
+	double a[3], b[3], apsb[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	iauPpsp(a, NUM2DBL(vs), b, apsb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(apsb);
+}
+
+static VALUE
+celes_s_pdp(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauPdp(a, b));
+}
+
+static VALUE
+celes_s_pxp(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3], axb[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	iauPxp(a, b, axb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(axb);
+}
+
+static VALUE
+celes_s_pm(VALUE self, VALUE vp){
+	double p[3];
+
+	cvec_cp(vp, p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauPm(p));
+}
+
+static VALUE
+celes_s_pn(VALUE self, VALUE vp){
+	double p[3], r, u[3];
+
+	cvec_cp(vp, p);
+	iauPn(p, &r, u);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(r), vvec_new(u));
+}
+
+static VALUE
+celes_s_sxp(VALUE self, VALUE vs, VALUE vp){
+	double p[3], sp[3];
+
+	cvec_cp(vp, p);
+	iauSxp(NUM2DBL(vs), p, sp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(sp);
+}
+
+
+static VALUE
+celes_s_rxr(VALUE self, VALUE va, VALUE vb){
+	double a[3][3], b[3][3], atb[3][3];
+
+	cmat_cp33(va, a);
+	cmat_cp33(vb, b);
+	iauRxr(a, b, atb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(atb);
+}
+
+static VALUE
+celes_s_tr(VALUE self, VALUE vr){
+	double r[3][3], rt[3][3];
+
+	cmat_cp33(vr, r);
+	iauTr(r, rt);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rt);
+}
+
+static VALUE
+celes_s_rxp(VALUE self, VALUE vr, VALUE vp){
+	double r[3][3], p[3], rp[3];
+
+	cmat_cp33(vr, r);
+	cvec_cp(vp, p);
+	iauRxp(r, p, rp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(rp);
+}
+
+static VALUE
+celes_s_trxp(VALUE self, VALUE vr, VALUE vp){
+	double r[3][3], p[3], trp[3];
+
+	cmat_cp33(vr, r);
+	cvec_cp(vp, p);
+	iauTrxp(r, p, trp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(trp);
+}
+
+static VALUE
+celes_s_sepp(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauSepp(a, b));
+}
+
+static VALUE
+celes_s_seps(VALUE self, VALUE val, VALUE vap, VALUE vbl, VALUE vbp){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauSeps(NUM2DBL(val), NUM2DBL(vap),
+			NUM2DBL(vbl), NUM2DBL(vbp)));
+}
+
+static VALUE
+celes_s_pap(VALUE self, VALUE va, VALUE vb){
+	double a[3], b[3];
+
+	cvec_cp(va, a);
+	cvec_cp(vb, b);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauPap(a, b));
+}
+
+static VALUE
+celes_s_pas(VALUE self, VALUE val, VALUE vap, VALUE vbl, VALUE vbp){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauPas(NUM2DBL(val), NUM2DBL(vap),
+			NUM2DBL(vbl), NUM2DBL(vbp)));
+}
+
+static VALUE
+celes_s_rv2m(VALUE self, VALUE vw){
+	double w[3], r[3][3];
+
+	cvec_cp(vw, w);
+	iauRv2m(w, r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_rm2v(VALUE self, VALUE vr){
+	double r[3][3], w[3];
+
+	cmat_cp33(vr, r);
+	iauRm2v(r, w);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(w);
+}
+
+static VALUE
+celes_s_zpv_b(VALUE self, VALUE vpv){
+	double pv[2][3];
+
+	cmat_cp23(vpv, pv);
+	iauZpv(pv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_cp23(pv, vpv);
+}
+
+static VALUE
+celes_s_zpv(VALUE self){
+	double pv[2][3];
+
+	iauZpv(pv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(pv);
+}
+
+static VALUE
+celes_s_cpv(VALUE self, VALUE vpv){
+	double pv[2][3], c[2][3];
+
+	cmat_cp23(vpv, pv);
+	iauCpv(pv, c);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(c);
+}
+
+static VALUE
+celes_s_p2pv(VALUE self, VALUE vp){
+	double p[3], pv[2][3];
+
+	cvec_cp(vp, p);
+	iauP2pv(p, pv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(pv);
+}
+
+static VALUE
+celes_s_pv2p(VALUE self, VALUE vpv){
+	double pv[2][3], p[3];
+
+	cmat_cp23(vpv, pv);
+	iauPv2p(pv, p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(p);
+}
+
+static VALUE
+celes_s_s2pv(VALUE self, VALUE vtheta, VALUE vphi, VALUE vr,
+		VALUE vtd, VALUE vpd, VALUE vrd){
+	double pv[2][3];
+
+	iauS2pv(NUM2DBL(vtheta), NUM2DBL(vphi), NUM2DBL(vr),
+			NUM2DBL(vtd), NUM2DBL(vpd), NUM2DBL(vrd), pv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(pv);
+}
+
+static VALUE
+celes_s_pv2s(VALUE self, VALUE vpv){
+	double pv[2][3], theta, phi, r, td, pd, rd;
+
+	cmat_cp23(vpv, pv);
+	iauPv2s(pv, &theta, &phi, &r, &td, &pd, &rd);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(6, DBL2NUM(theta), DBL2NUM(phi),
+			DBL2NUM(r), DBL2NUM(td), DBL2NUM(pd), DBL2NUM(rd));
+}
+
+static VALUE
+celes_s_pvppv(VALUE self, VALUE va, VALUE vb){
+	double a[2][3], b[2][3], apb[2][3];
+
+	cmat_cp23(va, a);
+	cmat_cp23(vb, b);
+	iauPvppv(a, b, apb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(apb);
+}
+
+static VALUE
+celes_s_pvmpv(VALUE self, VALUE va, VALUE vb){
+	double a[2][3], b[2][3], amb[2][3];
+
+	cmat_cp23(va, a);
+	cmat_cp23(vb, b);
+	iauPvmpv(a, b, amb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(amb);
+}
+
+static VALUE
+celes_s_pvdpv(VALUE self, VALUE va, VALUE vb){
+	double a[2][3], b[2][3], adb[2];
+
+	cmat_cp23(va, a);
+	cmat_cp23(vb, b);
+	iauPvdpv(a, b, adb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(adb[0]), DBL2NUM(adb[1]));
+}
+
+static VALUE
+celes_s_pvxpv(VALUE self, VALUE va, VALUE vb){
+	double a[2][3], b[2][3], axb[2][3];
+
+	cmat_cp23(va, a);
+	cmat_cp23(vb, b);
+	iauPvxpv(a, b, axb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(axb);
+}
+
+static VALUE
+celes_s_pvm(VALUE self, VALUE vpv){
+	double pv[2][3], r, s;
+
+	cmat_cp23(vpv, pv);
+	iauPvm(pv, &r, &s);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(r), DBL2NUM(s));
+}
+
+static VALUE
+celes_s_sxpv(VALUE self, VALUE vs, VALUE vpv){
+	double pv[2][3], spv[2][3];
+
+	cmat_cp23(vpv, pv);
+	iauSxpv(NUM2DBL(vs), pv, spv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(spv);
+}
+
+static VALUE
+celes_s_s2xpv(VALUE self, VALUE vs1, VALUE vs2, VALUE vpv){
+	double pv[2][3], spv[2][3];
+
+	cmat_cp23(vpv, pv);
+	iauS2xpv(NUM2DBL(vs1), NUM2DBL(vs2), pv, spv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(spv);
+}
+
+static VALUE
+celes_s_pvu(VALUE self, VALUE vdt, VALUE vpv){
+	double pv[2][3], upv[2][3];
+
+	cmat_cp23(vpv, pv);
+	iauPvu(NUM2DBL(vdt), pv, upv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(upv);
+}
+
+static VALUE
+celes_s_pvup(VALUE self, VALUE vdt, VALUE vpv){
+	double pv[2][3], p[3];
+
+	cmat_cp23(vpv, pv);
+	iauPvup(NUM2DBL(vdt), pv, p);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vvec_new(p);
+}
+
+static VALUE
+celes_s_rxpv(VALUE self, VALUE vr, VALUE vpv){
+	double r[3][3], pv[2][3], rpv[2][3];
+
+	cmat_cp33(vr, r);
+	cmat_cp23(vpv, pv);
+	iauRxpv(r, pv, rpv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(rpv);
+}
+
+static VALUE
+celes_s_trxpv(VALUE self, VALUE vr, VALUE vpv){
+	double r[3][3], pv[2][3], trpv[2][3];
+
+	cmat_cp33(vr, r);
+	cmat_cp23(vpv, pv);
+	iauTrxpv(r, pv, trpv);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new23(trpv);
+}
+
+static VALUE
+celes_s_anp(VALUE self, VALUE va){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauAnp(NUM2DBL(va)));
+}
+
+static VALUE
+celes_s_anpm(VALUE self, VALUE va){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauAnpm(NUM2DBL(va)));
+}
+
+static VALUE
+celes_s_a2tf(VALUE self, VALUE vndp, VALUE vangle){
+	char sign;
+	int ihmsf[4];
+
+	iauA2tf(NUM2INT(vndp), NUM2DBL(vangle), &sign, ihmsf);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, rb_str_new(&sign, 1),
+			rb_ary_new3(4, INT2NUM(ihmsf[0]), INT2NUM(ihmsf[1]),
+				INT2NUM(ihmsf[2]), INT2NUM(ihmsf[3])));
+}
+
+static VALUE
+celes_s_a2af(VALUE self, VALUE vndp, VALUE vangle){
+	char sign;
+	int idmsf[4];
+
+	iauA2af(NUM2INT(vndp), NUM2DBL(vangle), &sign, idmsf);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, rb_str_new(&sign, 1),
+			rb_ary_new3(4, INT2NUM(idmsf[0]), INT2NUM(idmsf[1]),
+				INT2NUM(idmsf[2]), INT2NUM(idmsf[3])));
+}
+
+static VALUE
+celes_s_af2a(VALUE self, VALUE vs, VALUE vdeg, VALUE vamin, VALUE vasec){
+	int ret;
+	double rad;
+
+	Check_Type(vs, T_STRING);
+
+	ret = iauAf2a(RSTRING_PTR(vs)[0], NUM2DBL(vdeg),
+			NUM2DBL(vamin), NUM2DBL(vasec), &rad);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	return DBL2NUM(rad);
+}
+
+static VALUE
+celes_s_d2tf(VALUE self, VALUE vndp, VALUE vdays){
+	char sign;
+	int ihmsf[4];
+
+	iauD2tf(NUM2INT(vndp), NUM2DBL(vdays), &sign, ihmsf);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, rb_str_new(&sign, 1),
+			rb_ary_new3(4, INT2NUM(ihmsf[0]), INT2NUM(ihmsf[1]),
+				INT2NUM(ihmsf[2]), INT2NUM(ihmsf[3])));
+}
+
+static VALUE
+celes_s_tf2a(VALUE self, VALUE vs, VALUE vhour, VALUE vmin, VALUE vsec){
+	int ret;
+	double rad;
+
+	Check_Type(vs, T_STRING);
+
+	ret = iauTf2a(RSTRING_PTR(vs)[0], NUM2DBL(vhour),
+			NUM2DBL(vmin), NUM2DBL(vsec), &rad);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	return DBL2NUM(rad);
+}
+
+static VALUE
+celes_s_tf2d(VALUE self, VALUE vs, VALUE vhour, VALUE vmin, VALUE vsec){
+	int ret;
+	double days;
+
+	Check_Type(vs, T_STRING);
+
+	ret = iauTf2d(RSTRING_PTR(vs)[0], NUM2DBL(vhour),
+			NUM2DBL(vmin), NUM2DBL(vsec), &days);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	return DBL2NUM(days);
+}
+
+static VALUE
+celes_s_cal2jd(VALUE self, VALUE vy, VALUE vm, VALUE vd){
+	int ret;
+	double djm0, djm;
+
+	ret = iauCal2jd(NUM2INT(vy), NUM2INT(vm), NUM2INT(vd), &djm0, &djm);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "bad year");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "bad month");
+
+	return rb_ary_new3(2, DBL2NUM(djm0), DBL2NUM(djm));
+}
+
+static VALUE
+celes_s_epb(VALUE self, VALUE vdj1, VALUE vdj2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEpb(NUM2DBL(vdj1), NUM2DBL(vdj2)));
+}
+
+static VALUE
+celes_s_epb2jd(VALUE self, VALUE vepb){
+	double djm0, djm;
+
+	iauEpb2jd(NUM2DBL(vepb), &djm0, &djm);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(djm0), DBL2NUM(djm));
+}
+
+static VALUE
+celes_s_epj(VALUE self, VALUE vdj1, VALUE vdj2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEpj(NUM2DBL(vdj1), NUM2DBL(vdj2)));
+}
+
+static VALUE
+celes_s_epj2jd(VALUE self, VALUE vepj){
+	double djm0, djm;
+
+	iauEpj2jd(NUM2DBL(vepj), &djm0, &djm);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(djm0), DBL2NUM(djm));
+}
+
+
+static VALUE
+celes_s_jd2cal(VALUE self, VALUE vdj1, VALUE vdj2){
+	int ret, iy, im, id;
+	double fd;
+
+	ret = iauJd2cal(NUM2DBL(vdj1), NUM2DBL(vdj2), &iy, &im, &id, &fd);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(4, INT2FIX(iy), INT2FIX(im),
+			INT2FIX(id), DBL2NUM(fd));
+}
+
+static VALUE
+celes_s_jdcalf(VALUE self, VALUE vndp, VALUE vdj1, VALUE vdj2){
+	int ret, iymdf[4];
+
+	ret = iauJdcalf(NUM2INT(vndp), NUM2DBL(vdj1), NUM2DBL(vdj2), iymdf);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "date out of range");
+
+	return rb_ary_new3(4, INT2FIX(iymdf[0]), INT2FIX(iymdf[1]),
+			INT2FIX(iymdf[2]), INT2NUM(iymdf[3]));
+}
+
+static VALUE
+celes_s_d2dtf(VALUE self, VALUE vscale, VALUE vndp, VALUE vd1, VALUE vd2){
+	int ret;
+	int iy, im, id, ihmsf[4];
+
+	if(NIL_P(vscale)) vscale = rb_str_new2("");
+	ret = iauD2dtf(RSTRING_PTR(vscale), NUM2INT(vndp),
+			NUM2DBL(vd1), NUM2DBL(vd2), &iy, &im, &id, ihmsf);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(4, INT2FIX(iy), INT2FIX(im), INT2FIX(id),
+			rb_ary_new3(4, INT2FIX(ihmsf[0]), INT2FIX(ihmsf[1]),
+			INT2FIX(ihmsf[2]), INT2NUM(ihmsf[3])));
+}
+
+static VALUE
+celes_s_dat(VALUE self, VALUE vy, VALUE vm, VALUE vd, VALUE vfd){
+	int ret;
+	double deltat;
+
+	ret = iauDat(NUM2INT(vy), NUM2INT(vm),
+			NUM2INT(vd), NUM2DBL(vfd), &deltat);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "bad year");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "bad month");
+	else if(ret == -3)
+		rb_raise(rb_eArgError, "bad day");
+	else if(ret == -4)
+		rb_raise(rb_eArgError, "bad fraction");
+
+	return DBL2NUM(deltat);
+}
+
+static VALUE
+celes_s_dtdb(VALUE self, VALUE vdate1, VALUE vdate2,
+		VALUE vut, VALUE velong, VALUE vu, VALUE vv){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauDtdb(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vut), NUM2DBL(velong), NUM2DBL(vu), NUM2DBL(vv)));
+}
+
+static VALUE
+celes_s_dtf2d(VALUE self, VALUE vscale, VALUE vy, VALUE vm, VALUE vd,
+		VALUE vhr, VALUE vmn, VALUE vsec){
+	int ret;
+	double d1, d2;
+
+	if(NIL_P(vscale)) vscale = rb_str_new2("");
+	ret = iauDtf2d(RSTRING_PTR(vscale), NUM2INT(vy), NUM2INT(vm),
+			NUM2INT(vd), NUM2INT(vhr), NUM2INT(vmn), NUM2DBL(vsec),
+			&d1, &d2);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "bad year");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "bad month");
+	else if(ret == -3)
+		rb_raise(rb_eArgError, "bad day");
+
+	return rb_ary_new3(2, DBL2NUM(d1), DBL2NUM(d2));
+}
+
+static VALUE
+celes_s_taitt(VALUE self, VALUE vtai1, VALUE vtai2){
+	double tt1, tt2;
+
+	iauTaitt(NUM2DBL(vtai1), NUM2DBL(vtai2), &tt1, &tt2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tt1), DBL2NUM(tt2));
+}
+
+static VALUE
+celes_s_taiut1(VALUE self, VALUE vtai1, VALUE vtai2, VALUE vdta){
+	double ut11, ut12;
+
+	iauTaiut1(NUM2DBL(vtai1), NUM2DBL(vtai2), NUM2DBL(vdta), &ut11, &ut12);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(ut11), DBL2NUM(ut12));
+}
+
+static VALUE
+celes_s_taiutc(VALUE self, VALUE vtai1, VALUE vtai2){
+	int ret;
+	double utc1, utc2;
+
+	ret = iauTaiutc(NUM2DBL(vtai1), NUM2DBL(vtai2), &utc1, &utc2);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(2, DBL2NUM(utc1), DBL2NUM(utc2));
+}
+
+static VALUE
+celes_s_tcbtdb(VALUE self, VALUE vtcb1, VALUE vtcb2){
+	double tdb1, tdb2;
+
+	iauTcbtdb(NUM2DBL(vtcb1), NUM2DBL(vtcb2), &tdb1, &tdb2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tdb1), DBL2NUM(tdb2));
+}
+
+static VALUE
+celes_s_tcgtt(VALUE self, VALUE vtcg1, VALUE vtcg2){
+	double tt1, tt2;
+
+	iauTcgtt(NUM2DBL(vtcg1), NUM2DBL(vtcg2), &tt1, &tt2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tt1), DBL2NUM(tt2));
+}
+
+static VALUE
+celes_s_tdbtcb(VALUE self, VALUE vtdb1, VALUE vtdb2){
+	double tcb1, tcb2;
+
+	iauTdbtcb(NUM2DBL(vtdb1), NUM2DBL(vtdb2), &tcb1, &tcb2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tcb1), DBL2NUM(tcb2));
+}
+
+static VALUE
+celes_s_tdbtt(VALUE self, VALUE vtdb1, VALUE vtdb2, VALUE vdtr){
+	double tt1, tt2;
+
+	iauTdbtt(NUM2DBL(vtdb1), NUM2DBL(vtdb2), NUM2DBL(vdtr), &tt1, &tt2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tt1), DBL2NUM(tt2));
+}
+
+static VALUE
+celes_s_tttai(VALUE self, VALUE vtt1, VALUE vtt2){
+	double tai1, tai2;
+
+	iauTttai(NUM2DBL(vtt1), NUM2DBL(vtt2), &tai1, &tai2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tai1), DBL2NUM(tai2));
+}
+
+static VALUE
+celes_s_tttcg(VALUE self, VALUE vtt1, VALUE vtt2){
+	double tcg1, tcg2;
+
+	iauTttcg(NUM2DBL(vtt1), NUM2DBL(vtt2), &tcg1, &tcg2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tcg1), DBL2NUM(tcg2));
+}
+
+static VALUE
+celes_s_tttdb(VALUE self, VALUE vtt1, VALUE vtt2, VALUE vdtr){
+	double tdb1, tdb2;
+
+	iauTttdb(NUM2DBL(vtt1), NUM2DBL(vtt2), NUM2DBL(vdtr), &tdb1, &tdb2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tdb1), DBL2NUM(tdb2));
+}
+
+static VALUE
+celes_s_ttut1(VALUE self, VALUE vtt1, VALUE vtt2, VALUE vdt){
+	double ut11, ut12;
+
+	iauTtut1(NUM2DBL(vtt1), NUM2DBL(vtt2), NUM2DBL(vdt), &ut11, &ut12);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(ut11), DBL2NUM(ut12));
+}
+
+static VALUE
+celes_s_ut1tai(VALUE self, VALUE vut11, VALUE vut12, VALUE vdta){
+	double tai1, tai2;
+
+	iauUt1tai(NUM2DBL(vut11), NUM2DBL(vut12), NUM2DBL(vdta), &tai1, &tai2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tai1), DBL2NUM(tai2));
+}
+
+static VALUE
+celes_s_ut1tt(VALUE self, VALUE vut11, VALUE vut12, VALUE vdt){
+	double tt1, tt2;
+
+	iauUt1tt(NUM2DBL(vut11), NUM2DBL(vut12), NUM2DBL(vdt), &tt1, &tt2);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(tt1), DBL2NUM(tt2));
+}
+
+static VALUE
+celes_s_ut1utc(VALUE self, VALUE vut11, VALUE vut12, VALUE vdut1){
+	int ret;
+	double utc1, utc2;
+
+	ret = iauUt1utc(NUM2DBL(vut11), NUM2DBL(vut12),
+			NUM2DBL(vdut1), &utc1, &utc2);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(2, DBL2NUM(utc1), DBL2NUM(utc2));
+}
+
+static VALUE
+celes_s_utctai(VALUE self, VALUE vutc1, VALUE vutc2){
+	int ret;
+	double tai1, tai2;
+
+	ret = iauUtctai(NUM2DBL(vutc1), NUM2DBL(vutc2), &tai1, &tai2);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(2, DBL2NUM(tai1), DBL2NUM(tai2));
+}
+
+static VALUE
+celes_s_utcut1(VALUE self, VALUE vutc1, VALUE vutc2, VALUE vdut1){
+	int ret;
+	double ut11, ut12;
+
+	ret = iauUtcut1(NUM2DBL(vutc1), NUM2DBL(vutc2),
+			NUM2DBL(vdut1), &ut11, &ut12);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "unacceptable date");
+
+	return rb_ary_new3(2, DBL2NUM(ut11), DBL2NUM(ut12));
+}
+
+static VALUE
+celes_s_ee00(VALUE self, VALUE vdate1, VALUE vdate2,
+		VALUE vepsa, VALUE vdpsi){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEe00(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vepsa), NUM2DBL(vdpsi)));
+}
+
+static VALUE
+celes_s_ee00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEe00a(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_ee00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEe00b(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_ee06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEe06a(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_eect00(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEect00(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_eqeq94(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEqeq94(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_era00(VALUE self, VALUE vdj1, VALUE vdj2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEra00(NUM2DBL(vdj1), NUM2DBL(vdj2)));
+}
+
+static VALUE
+celes_s_gmst00(VALUE self, VALUE vuta, VALUE vutb, VALUE vtta, VALUE vttb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGmst00(NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vtta), NUM2DBL(vttb)));
+}
+
+static VALUE
+celes_s_gmst06(VALUE self, VALUE vuta, VALUE vutb, VALUE vtta, VALUE vttb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGmst06(NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vtta), NUM2DBL(vttb)));
+}
+
+static VALUE
+celes_s_gmst82(VALUE self, VALUE vdj1, VALUE vdj2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGmst82(NUM2DBL(vdj1), NUM2DBL(vdj2)));
+}
+
+static VALUE
+celes_s_gst00a(VALUE self, VALUE vuta, VALUE vutb, VALUE vtta, VALUE vttb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGst00a(NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vtta), NUM2DBL(vttb)));
+}
+
+static VALUE
+celes_s_gst00b(VALUE self, VALUE vuta, VALUE vutb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGst00b(NUM2DBL(vuta), NUM2DBL(vutb)));
+}
+
+static VALUE
+celes_s_gst06(VALUE self, VALUE vuta, VALUE vutb,
+		VALUE vtta, VALUE vttb, VALUE vrnpb){
+	double rnpb[3][3];
+
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	cmat_cp33(vrnpb, rnpb);
+
+	return DBL2NUM(iauGst06(NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vtta), NUM2DBL(vttb), rnpb));
+}
+
+static VALUE
+celes_s_gst06a(VALUE self, VALUE vuta, VALUE vutb, VALUE vtta, VALUE vttb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGst06a(NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vtta), NUM2DBL(vttb)));
+}
+
+static VALUE
+celes_s_gst94(VALUE self, VALUE vuta, VALUE vutb){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauGst94(NUM2DBL(vuta), NUM2DBL(vutb)));
+}
+
+static VALUE
+celes_s_epv00(VALUE self, VALUE vdate1, VALUE vdate2){
+	int ret;
+	double pvh[2][3], pvb[2][3];
+
+	ret = iauEpv00(NUM2DBL(vdate1), NUM2DBL(vdate2), pvh, pvb);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	return rb_ary_new3(2, vmat_new23(pvh), vmat_new23(pvb));
+}
+
+static VALUE
+celes_s_plan94(VALUE self, VALUE vdate1, VALUE vdate2, VALUE vnp){
+	int ret;
+	double pv[2][3];
+
+	ret = iauPlan94(NUM2DBL(vdate1), NUM2DBL(vdate2), NUM2INT(vnp), pv);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal NP (outside 1-8)");
+
+	return vmat_new23(pv);
+}
+static VALUE
+celes_s_bi00(VALUE self){
+	double dpsibi, depsbi, dra;
+
+	iauBi00(&dpsibi, &depsbi, &dra);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(dpsibi), DBL2NUM(depsbi), DBL2NUM(dra));
+}
+
+static VALUE
+celes_s_bp00(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rb[3][3], rp[3][3], rbp[3][3];
+
+	iauBp00(NUM2DBL(vdate1), NUM2DBL(vdate2), rb, rp, rbp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, vmat_new33(rb), vmat_new33(rp), vmat_new33(rbp));
+}
+
+static VALUE
+celes_s_bp06(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rb[3][3], rp[3][3], rbp[3][3];
+
+	iauBp06(NUM2DBL(vdate1), NUM2DBL(vdate2), rb, rp, rbp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, vmat_new33(rb), vmat_new33(rp), vmat_new33(rbp));
+}
+
+static VALUE
+celes_s_bpn2xy(VALUE self, VALUE vrbpn){
+	double rbpn[3][3], x, y;
+
+	cmat_cp33(vrbpn, rbpn);
+	iauBpn2xy(rbpn, &x, &y);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(x), DBL2NUM(y));
+}
+
+static VALUE
+celes_s_c2i00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rc2i[3][3];
+
+	iauC2i00a(NUM2DBL(vdate1), NUM2DBL(vdate2), rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2i00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rc2i[3][3];
+
+	iauC2i00b(NUM2DBL(vdate1), NUM2DBL(vdate2), rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2i06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rc2i[3][3];
+
+	iauC2i06a(NUM2DBL(vdate1), NUM2DBL(vdate2), rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2ibpn(VALUE self, VALUE vdate1, VALUE vdate2, VALUE vrbpn){
+	double rbpn[3][3], rc2i[3][3];
+
+	cmat_cp33(vrbpn, rbpn);
+	iauC2ibpn(NUM2DBL(vdate1), NUM2DBL(vdate2), rbpn, rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2ixy(VALUE self, VALUE vdate1, VALUE vdate2, VALUE vx, VALUE vy){
+	double rc2i[3][3];
+
+	iauC2ixy(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vx), NUM2DBL(vy), rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2ixys(VALUE self, VALUE vx, VALUE vy, VALUE vs){
+	double rc2i[3][3];
+
+	iauC2ixys(NUM2DBL(vx), NUM2DBL(vy), NUM2DBL(vs), rc2i);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2i);
+}
+
+static VALUE
+celes_s_c2t00a(VALUE self, VALUE vtta, VALUE vttb,
+		VALUE vuta, VALUE vutb, VALUE vxp, VALUE vyp){
+	double rc2t[3][3];
+
+	iauC2t00a(NUM2DBL(vtta), NUM2DBL(vttb), NUM2DBL(vuta), NUM2DBL(vutb), 
+			NUM2DBL(vxp), NUM2DBL(vyp), rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2t00b(VALUE self, VALUE vtta, VALUE vttb,
+		VALUE vuta, VALUE vutb, VALUE vxp, VALUE vyp){
+	double rc2t[3][3];
+
+	iauC2t00b(NUM2DBL(vtta), NUM2DBL(vttb), NUM2DBL(vuta), NUM2DBL(vutb), 
+			NUM2DBL(vxp), NUM2DBL(vyp), rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2t06a(VALUE self, VALUE vtta, VALUE vttb,
+		VALUE vuta, VALUE vutb, VALUE vxp, VALUE vyp){
+	double rc2t[3][3];
+
+	iauC2t06a(NUM2DBL(vtta), NUM2DBL(vttb), NUM2DBL(vuta), NUM2DBL(vutb), 
+			NUM2DBL(vxp), NUM2DBL(vyp), rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2tcio(VALUE self, VALUE vrc2i, VALUE vera, VALUE vrpom){
+	double rc2i[3][3], rpom[3][3], rc2t[3][3];
+
+	cmat_cp33(vrc2i, rc2i);
+	cmat_cp33(vrpom, rpom);
+	iauC2tcio(rc2i, NUM2DBL(vera), rpom, rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2teqx(VALUE self, VALUE vrbpn, VALUE vgst, VALUE vrpom){
+	double rbpn[3][3], rpom[3][3], rc2t[3][3];
+
+	cmat_cp33(vrbpn, rbpn);
+	cmat_cp33(vrpom, rpom);
+	iauC2teqx(rbpn, NUM2DBL(vgst), rpom, rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2tpe(VALUE self, VALUE vtta, VALUE vttb, VALUE vuta, VALUE vutb,
+		VALUE vdpsi, VALUE vdeps, VALUE vxp, VALUE vyp){
+	double rc2t[3][3];
+
+	iauC2tpe(NUM2DBL(vtta), NUM2DBL(vttb), NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vdpsi), NUM2DBL(vdeps), NUM2DBL(vxp), NUM2DBL(vyp),
+			rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_c2txy(VALUE self, VALUE vtta, VALUE vttb, VALUE vuta, VALUE vutb,
+		VALUE vx, VALUE vy, VALUE vxp, VALUE vyp){
+	double rc2t[3][3];
+
+	iauC2txy(NUM2DBL(vtta), NUM2DBL(vttb), NUM2DBL(vuta), NUM2DBL(vutb),
+			NUM2DBL(vx), NUM2DBL(vy), NUM2DBL(vxp), NUM2DBL(vyp), rc2t);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rc2t);
+}
+
+static VALUE
+celes_s_eo06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEo06a(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_eors(VALUE self, VALUE vrnpb, VALUE vs){
+	double rnpb[3][3];
+
+	cmat_cp33(vrnpb, rnpb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauEors(rnpb, NUM2DBL(vs)));
+}
+
+static VALUE
+celes_s_fw2m(VALUE self, VALUE vgamb, VALUE vphib, VALUE vpsi, VALUE veps){
+	double r[3][3];
+
+	iauFw2m(NUM2DBL(vgamb), NUM2DBL(vphib),
+			NUM2DBL(vpsi), NUM2DBL(veps), r);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(r);
+}
+
+static VALUE
+celes_s_fw2xy(VALUE self, VALUE vgamb,
+		VALUE vphib, VALUE vpsi, VALUE veps){
+	double x, y;
+
+	iauFw2xy(NUM2DBL(vgamb), NUM2DBL(vphib),
+			NUM2DBL(vpsi), NUM2DBL(veps), &x, &y);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(x), DBL2NUM(y));
+}
+
+static VALUE
+celes_s_num00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatn[3][3];
+
+	iauNum00a(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatn);
+}
+
+static VALUE
+celes_s_num00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatn[3][3];
+
+	iauNum00b(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatn);
+}
+
+static VALUE
+celes_s_num06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatn[3][3];
+
+	iauNum06a(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatn);
+}
+
+static VALUE
+celes_s_numat(VALUE self, VALUE vepsa, VALUE vdpsi, VALUE vdeps){
+	double rmatn[3][3];
+
+	iauNumat(NUM2DBL(vepsa), NUM2DBL(vdpsi), NUM2DBL(vdeps), rmatn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatn);
+}
+
+static VALUE
+celes_s_nut00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps;
+
+	iauNut00a(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(dpsi), DBL2NUM(deps));
+}
+
+static VALUE
+celes_s_nut00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps;
+
+	iauNut00b(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(dpsi), DBL2NUM(deps));
+}
+
+static VALUE
+celes_s_nut06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps;
+
+	iauNut06a(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(dpsi), DBL2NUM(deps));
+}
+
+static VALUE
+celes_s_nut80(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps;
+
+	iauNut80(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(dpsi), DBL2NUM(deps));
+}
+
+static VALUE
+celes_s_nutm80(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatn[3][3];
+
+	iauNutm80(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatn);
+}
+
+static VALUE
+celes_s_obl06(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauObl06(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_obl80(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauObl80(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_pb06(VALUE self, VALUE vdate1, VALUE vdate2){
+	double bzeta, bz, btheta;
+
+	iauPb06(NUM2DBL(vdate1), NUM2DBL(vdate2), &bzeta, &bz, &btheta);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(bzeta), DBL2NUM(bz), DBL2NUM(btheta));
+}
+
+static VALUE
+celes_s_pfw06(VALUE self, VALUE vdate1, VALUE vdate2){
+	double gamb, phib, psib, epsa;
+
+	iauPfw06(NUM2DBL(vdate1), NUM2DBL(vdate2), &gamb, &phib, &psib, &epsa);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(4, DBL2NUM(gamb),
+			DBL2NUM(phib), DBL2NUM(psib), DBL2NUM(epsa));
+}
+
+static VALUE
+celes_s_pmat00(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rbp[3][3];
+
+	iauPmat00(NUM2DBL(vdate1), NUM2DBL(vdate2), rbp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rbp);
+}
+
+static VALUE
+celes_s_pmat06(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rbp[3][3];
+
+	iauPmat06(NUM2DBL(vdate1), NUM2DBL(vdate2), rbp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rbp);
+}
+
+static VALUE
+celes_s_pmat76(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatp[3][3];
+
+	iauPmat76(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatp);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatp);
+}
+
+static VALUE
+celes_s_pn00(VALUE self, VALUE vdate1, VALUE vdate2,
+		VALUE vdpsi, VALUE vdeps){
+	double epsa, rb[3][3], rp[3][3], rbp[3][3], rn[3][3], rbpn[3][3];
+
+	iauPn00(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vdpsi), NUM2DBL(vdeps), &epsa,
+			rb, rp, rbp, rn, rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(6, DBL2NUM(epsa), vmat_new33(rb),
+			vmat_new33(rp), vmat_new33(rbp), vmat_new33(rn), vmat_new33(rbpn));
+}
+
+static VALUE
+celes_s_pn00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps, epsa, rb[3][3],
+			rp[3][3], rbp[3][3], rn[3][3], rbpn[3][3];
+
+	iauPn00a(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps, &epsa,
+			rb, rp, rbp, rn, rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(8, DBL2NUM(dpsi), DBL2NUM(deps), DBL2NUM(epsa),
+			vmat_new33(rb), vmat_new33(rp), vmat_new33(rbp),
+			vmat_new33(rn), vmat_new33(rbpn));
+}
+
+static VALUE
+celes_s_pn00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps, epsa, rb[3][3],
+			rp[3][3], rbp[3][3], rn[3][3], rbpn[3][3];
+
+	iauPn00b(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps, &epsa,
+			rb, rp, rbp, rn, rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(8, DBL2NUM(dpsi), DBL2NUM(deps), DBL2NUM(epsa),
+			vmat_new33(rb), vmat_new33(rp), vmat_new33(rbp),
+			vmat_new33(rn), vmat_new33(rbpn));
+}
+
+static VALUE
+celes_s_pn06(VALUE self, VALUE vdate1, VALUE vdate2,
+		VALUE vdpsi, VALUE vdeps){
+	double epsa, rb[3][3], rp[3][3], rbp[3][3], rn[3][3], rbpn[3][3];
+
+	iauPn06(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vdpsi), NUM2DBL(vdeps), &epsa,
+			rb, rp, rbp, rn, rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(6, DBL2NUM(epsa), vmat_new33(rb),
+			vmat_new33(rp), vmat_new33(rbp), vmat_new33(rn), vmat_new33(rbpn));
+}
+
+static VALUE
+celes_s_pn06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsi, deps, epsa, rb[3][3],
+			rp[3][3], rbp[3][3], rn[3][3], rbpn[3][3];
+
+	iauPn06a(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsi, &deps, &epsa,
+			rb, rp, rbp, rn, rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(8, DBL2NUM(dpsi), DBL2NUM(deps), DBL2NUM(epsa),
+			vmat_new33(rb), vmat_new33(rp), vmat_new33(rbp),
+			vmat_new33(rn), vmat_new33(rbpn));
+}
+
+static VALUE
+celes_s_pnm00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rbpn[3][3];
+
+	iauPnm00a(NUM2DBL(vdate1), NUM2DBL(vdate2), rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rbpn);
+}
+
+static VALUE
+celes_s_pnm00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rbpn[3][3];
+
+	iauPnm00b(NUM2DBL(vdate1), NUM2DBL(vdate2), rbpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rbpn);
+}
+
+static VALUE
+celes_s_pnm06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rnpb[3][3];
+
+	iauPnm06a(NUM2DBL(vdate1), NUM2DBL(vdate2), rnpb);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rnpb);
+}
+
+static VALUE
+celes_s_pnm80(VALUE self, VALUE vdate1, VALUE vdate2){
+	double rmatpn[3][3];
+
+	iauPnm80(NUM2DBL(vdate1), NUM2DBL(vdate2), rmatpn);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rmatpn);
+}
+
+static VALUE
+celes_s_p06e(VALUE self, VALUE vdate1, VALUE vdate2){
+	double eps0, psia, oma, bpa, bqa, pia, bpia,
+			epsa, chia, za, zetaa, thetaa, pa, gam, phi, psi;
+	
+	iauP06e(NUM2DBL(vdate1), NUM2DBL(vdate2), &eps0, &psia,
+			&oma, &bpa, &bqa, &pia, &bpia, &epsa, &chia, &za,
+			&zetaa, &thetaa, &pa, &gam, &phi, &psi);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(16, DBL2NUM(eps0), DBL2NUM(psia),
+			DBL2NUM(oma), DBL2NUM(bpa), DBL2NUM(bqa), DBL2NUM(pia),
+			DBL2NUM(bpia), DBL2NUM(epsa), DBL2NUM(chia), DBL2NUM(za),
+			DBL2NUM(zetaa), DBL2NUM(thetaa), DBL2NUM(pa), DBL2NUM(gam),
+			DBL2NUM(phi), DBL2NUM(psi));
+}
+
+static VALUE
+celes_s_pom00(VALUE self, VALUE vxp, VALUE vyp, VALUE vsp){
+	double rpom[3][3];
+
+	iauPom00(NUM2DBL(vxp), NUM2DBL(vyp), NUM2DBL(vsp), rpom);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return vmat_new33(rpom);
+}
+
+static VALUE
+celes_s_pr00(VALUE self, VALUE vdate1, VALUE vdate2){
+	double dpsipr, depspr;
+
+	iauPr00(NUM2DBL(vdate1), NUM2DBL(vdate2), &dpsipr, &depspr);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(dpsipr), DBL2NUM(depspr));
+}
+
+static VALUE
+celes_s_prec76(VALUE self, VALUE vep01, VALUE vep02,
+		VALUE vep11, VALUE vep12){
+	double zeta, z, theta;
+
+	iauPrec76(NUM2DBL(vep01), NUM2DBL(vep02),
+			NUM2DBL(vep11), NUM2DBL(vep12), &zeta, &z, &theta);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(zeta), DBL2NUM(z), DBL2NUM(theta));
+}
+
+static VALUE
+celes_s_s00(VALUE self, VALUE vdate1, VALUE vdate2, VALUE vx, VALUE vy){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauS00(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vx), NUM2DBL(vy)));
+}
+
+static VALUE
+celes_s_s00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauS00a(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_s00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauS00b(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_s06(VALUE self, VALUE vdate1, VALUE vdate2, VALUE vx, VALUE vy){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauS06(NUM2DBL(vdate1), NUM2DBL(vdate2),
+			NUM2DBL(vx), NUM2DBL(vy)));
+}
+
+static VALUE
+celes_s_s06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauS06a(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_sp00(VALUE self, VALUE vdate1, VALUE vdate2){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauSp00(NUM2DBL(vdate1), NUM2DBL(vdate2)));
+}
+
+static VALUE
+celes_s_xy06(VALUE self, VALUE vdate1, VALUE vdate2){
+	double x, y;
+
+	iauXy06(NUM2DBL(vdate1), NUM2DBL(vdate2), &x, &y);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(x), DBL2NUM(y));
+}
+
+static VALUE
+celes_s_xys00a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double x, y, s;
+
+	iauXys00a(NUM2DBL(vdate1), NUM2DBL(vdate2), &x, &y, &s);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(x), DBL2NUM(y), DBL2NUM(s));
+}
+
+static VALUE
+celes_s_xys00b(VALUE self, VALUE vdate1, VALUE vdate2){
+	double x, y, s;
+
+	iauXys00b(NUM2DBL(vdate1), NUM2DBL(vdate2), &x, &y, &s);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(x), DBL2NUM(y), DBL2NUM(s));
+}
+
+static VALUE
+celes_s_xys06a(VALUE self, VALUE vdate1, VALUE vdate2){
+	double x, y, s;
+
+	iauXys06a(NUM2DBL(vdate1), NUM2DBL(vdate2), &x, &y, &s);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(3, DBL2NUM(x), DBL2NUM(y), DBL2NUM(s));
+}
+
+static VALUE
+celes_s_fad03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFad03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fae03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFae03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_faf03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFaf03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_faju03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFaju03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fal03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFal03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_falp03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFalp03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fama03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFama03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fame03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFame03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fane03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFane03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_faom03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFaom03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fapa03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFapa03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fasa03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFasa03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_faur03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFaur03(NUM2DBL(vt)));
+}
+
+static VALUE
+celes_s_fave03(VALUE self, VALUE vt){
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return DBL2NUM(iauFave03(NUM2DBL(vt)));
+}
+static VALUE
+celes_s_pvstar(VALUE self, VALUE vpv){
+	int ret;
+	double pv[2][3], ra, dec, pmr, pmd, px, rv;
+
+	cmat_cp23(vpv, pv);
+	ret = iauPvstar(pv, &ra, &dec, &pmr, &pmd, &px, &rv);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "superluminal speed");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "null position vector");
+	
+	return rb_ary_new3(6, DBL2NUM(ra), DBL2NUM(dec),
+			DBL2NUM(pmr), DBL2NUM(pmd), DBL2NUM(px), DBL2NUM(rv));
+}
+
+static VALUE
+celes_s_starpv(VALUE self, VALUE vra, VALUE vdec,
+		VALUE vpmr, VALUE vpmd, VALUE vpx, VALUE vrv){
+	int ret;
+	double pv[2][3];
+
+	ret = iauStarpv(NUM2DBL(vra), NUM2DBL(vdec), NUM2DBL(vpmr),
+			NUM2DBL(vpmd), NUM2DBL(vpx), NUM2DBL(vrv), pv);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	return vmat_new23(pv);
+}
+static VALUE
+celes_s_fk52h(VALUE self, VALUE vr5, VALUE vd5,
+		VALUE vdr5, VALUE vdd5, VALUE vpx5, VALUE vrv5){
+	double rh, dh, drh, ddh, pxh, rvh;
+
+	iauFk52h(NUM2DBL(vr5), NUM2DBL(vd5), NUM2DBL(vdr5),
+			NUM2DBL(vdd5), NUM2DBL(vpx5), NUM2DBL(vrv5),
+			&rh, &dh, &drh, &ddh, &pxh, &rvh);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(6, DBL2NUM(rh), DBL2NUM(dh),
+			DBL2NUM(drh), DBL2NUM(ddh), DBL2NUM(pxh), DBL2NUM(rvh));
+}
+
+static VALUE
+celes_s_fk5hip(VALUE self){
+	double r5h[3][3], s5h[3];
+
+	iauFk5hip(r5h, s5h);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, vmat_new33(r5h), vvec_new(s5h));
+}
+
+static VALUE
+celes_s_fk5hz(VALUE self, VALUE vr5, VALUE vd5,
+		VALUE vdate1, VALUE vdate2){
+	double rh, dh;
+
+	iauFk5hz(NUM2DBL(vr5), NUM2DBL(vd5),
+			NUM2DBL(vdate1), NUM2DBL(vdate2), &rh, &dh);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(2, DBL2NUM(rh), DBL2NUM(dh));
+}
+
+static VALUE
+celes_s_h2fk5(VALUE self, VALUE vrh, VALUE vdh,
+		VALUE vdrh, VALUE vddh, VALUE vpxh, VALUE vrvh){
+	double r5, d5, dr5, dd5, px5, rv5;
+
+	iauH2fk5(NUM2DBL(vrh), NUM2DBL(vdh), NUM2DBL(vdrh),
+			NUM2DBL(vddh), NUM2DBL(vpxh), NUM2DBL(vrvh),
+			&r5, &d5, &dr5, &dd5, &px5, &rv5);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(6, DBL2NUM(r5), DBL2NUM(d5),
+			DBL2NUM(dr5), DBL2NUM(dd5), DBL2NUM(px5), DBL2NUM(rv5));
+}
+
+static VALUE
+celes_s_hfk5z(VALUE self, VALUE vrh, VALUE vdh,
+		VALUE vdate1, VALUE vdate2){
+	double r5, d5, dr5, dd5;
+
+	iauHfk5z(NUM2DBL(vrh), NUM2DBL(vdh), NUM2DBL(vdate1), NUM2DBL(vdate2),
+			&r5, &d5, &dr5, &dd5);
+	rb_ivar_set(self, id_status, INT2FIX(0));
+
+	return rb_ary_new3(4, DBL2NUM(r5), DBL2NUM(d5),
+			DBL2NUM(dr5), DBL2NUM(dd5));
+}
+
+static VALUE
+celes_s_starpm(VALUE self, VALUE vra1, VALUE vdec1,
+		VALUE vpmr1, VALUE vpmd1, VALUE vpx1, VALUE vrv1,
+		VALUE vep1a, VALUE vep1b, VALUE vep2a, VALUE vep2b){
+	int ret;
+	double ra2, dec2, pmr2, pmd2, px2, rv2;
+
+	ret = iauStarpm(NUM2DBL(vra1), NUM2DBL(vdec1),
+			NUM2DBL(vpmr1), NUM2DBL(vpmd1), NUM2DBL(vpx1), NUM2DBL(vrv1),
+			NUM2DBL(vep1a), NUM2DBL(vep1b), NUM2DBL(vep2a), NUM2DBL(vep2b),
+			&ra2, &dec2, &pmr2, &pmd2, &px2, &rv2);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "system error");
+	
+	return rb_ary_new3(6,
+			DBL2NUM(ra2), DBL2NUM(dec2), DBL2NUM(pmr2),
+			DBL2NUM(pmd2), DBL2NUM(px2), DBL2NUM(rv2));
+}
+static VALUE
+celes_s_eform(VALUE self, VALUE vn){
+	int ret;
+	double a, f;
+
+	ret = iauEform(NUM2INT(vn), &a, &f);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal identifier");
+
+	return rb_ary_new3(2, DBL2NUM(a), DBL2NUM(f));
+}
+
+static VALUE
+celes_s_gc2gd(VALUE self, VALUE vn, VALUE vxyz){
+	int ret;
+	double xyz[3], elong, phi, height;
+
+	cvec_cp(vxyz, xyz);
+	ret = iauGc2gd(NUM2INT(vn), xyz, &elong, &phi, &height);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal identifier");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "internal error");
+
+	return rb_ary_new3(3, DBL2NUM(elong), DBL2NUM(phi), DBL2NUM(height));
+}
+
+static VALUE
+celes_s_gc2gde(VALUE self, VALUE va, VALUE vf, VALUE vxyz){
+	int ret;
+	double xyz[3], elong, phi, height;
+
+	cvec_cp(vxyz, xyz);
+	ret = iauGc2gde(NUM2DBL(va), NUM2DBL(vf), xyz, &elong, &phi, &height);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal f");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "illegal a");
+
+	return rb_ary_new3(3, DBL2NUM(elong), DBL2NUM(phi), DBL2NUM(height));
+}
+
+static VALUE
+celes_s_gd2gc(VALUE self, VALUE vn,
+		VALUE velong, VALUE vphi, VALUE vheight){
+	int ret;
+	double xyz[3];
+
+	ret = iauGd2gc(NUM2INT(vn), NUM2DBL(velong),
+			NUM2DBL(vphi), NUM2DBL(vheight), xyz);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal identifier");
+	else if(ret == -2)
+		rb_raise(rb_eArgError, "illegal case");
+
+	return vvec_new(xyz);
+}
+
+static VALUE
+celes_s_gd2gce(VALUE self, VALUE va, VALUE vf,
+		VALUE velong, VALUE vphi, VALUE vheight){
+	int ret;
+	double xyz[3];
+
+	ret = iauGd2gce(NUM2DBL(va), NUM2DBL(vf),
+			NUM2DBL(velong), NUM2DBL(vphi), NUM2DBL(vheight), xyz);
+	rb_ivar_set(self, id_status, INT2FIX(ret));
+
+	if(ret == -1)
+		rb_raise(rb_eArgError, "illegal case");
+
+	return vvec_new(xyz);
+}
+
+
+void
+Init_celes_core(void){
+	mCeles = rb_define_module("Celes");
+	id_status = rb_intern("@status");
+
+	rb_define_module_function(mCeles, "status", celes_s_status, 0);
+
+	rb_define_module_function(mCeles, "zp!",    celes_s_zp_b,   1);
+	rb_define_module_function(mCeles, "zp",     celes_s_zp,     0);
+	rb_define_module_function(mCeles, "zr!",    celes_s_zr_b,   1);
+	rb_define_module_function(mCeles, "zr",     celes_s_zr,     0);
+	rb_define_module_function(mCeles, "ir!",    celes_s_ir_b,   1);
+	rb_define_module_function(mCeles, "ir",     celes_s_ir,     0);
+	rb_define_module_function(mCeles, "cp",     celes_s_cp,     1);
+	rb_define_module_function(mCeles, "cr",     celes_s_cr,     1);
+	rb_define_module_function(mCeles, "rx!",    celes_s_rx_b,   2);
+	rb_define_module_function(mCeles, "rx",     celes_s_rx,     2);
+	rb_define_module_function(mCeles, "ry!",    celes_s_ry_b,   2);
+	rb_define_module_function(mCeles, "ry",     celes_s_ry,     2);
+	rb_define_module_function(mCeles, "rz!",    celes_s_rz_b,   2);
+	rb_define_module_function(mCeles, "rz",     celes_s_rz,     2);
+	rb_define_module_function(mCeles, "s2c",    celes_s_s2c,    2);
+	rb_define_module_function(mCeles, "c2s",    celes_s_c2s,    1);
+	rb_define_module_function(mCeles, "s2p",    celes_s_s2p,    3);
+	rb_define_module_function(mCeles, "p2s",    celes_s_p2s,    1);
+	rb_define_module_function(mCeles, "ppp",    celes_s_ppp,    2);
+	rb_define_module_function(mCeles, "pmp",    celes_s_pmp,    2);
+	rb_define_module_function(mCeles, "ppsp",   celes_s_ppsp,   3);
+	rb_define_module_function(mCeles, "pdp",    celes_s_pdp,    2);
+	rb_define_module_function(mCeles, "pxp",    celes_s_pxp,    2);
+	rb_define_module_function(mCeles, "pm",     celes_s_pm,     1);
+	rb_define_module_function(mCeles, "pn",     celes_s_pn,     1);
+	rb_define_module_function(mCeles, "sxp",    celes_s_sxp,    2);
+	rb_define_module_function(mCeles, "rxr",    celes_s_rxr,    2);
+	rb_define_module_function(mCeles, "tr",     celes_s_tr,     1);
+	rb_define_module_function(mCeles, "rxp",    celes_s_rxp,    2);
+	rb_define_module_function(mCeles, "trxp",   celes_s_trxp,   2);
+	rb_define_module_function(mCeles, "sepp",   celes_s_sepp,   2);
+	rb_define_module_function(mCeles, "seps",   celes_s_seps,   4);
+	rb_define_module_function(mCeles, "pap",    celes_s_pap,    2);
+	rb_define_module_function(mCeles, "pas",    celes_s_pas,    4);
+	rb_define_module_function(mCeles, "rv2m",   celes_s_rv2m,   1);
+	rb_define_module_function(mCeles, "rm2v",   celes_s_rm2v,   1);
+	rb_define_module_function(mCeles, "zpv!",   celes_s_zpv_b,  1);
+	rb_define_module_function(mCeles, "zpv",    celes_s_zpv,    0);
+	rb_define_module_function(mCeles, "cpv",    celes_s_cpv,    1);
+	rb_define_module_function(mCeles, "p2pv",   celes_s_p2pv,   1);
+	rb_define_module_function(mCeles, "pv2p",   celes_s_pv2p,   1);
+	rb_define_module_function(mCeles, "s2pv",   celes_s_s2pv,   6);
+	rb_define_module_function(mCeles, "pv2s",   celes_s_pv2s,   1);
+	rb_define_module_function(mCeles, "pvppv",  celes_s_pvppv,  2);
+	rb_define_module_function(mCeles, "pvmpv",  celes_s_pvmpv,  2);
+	rb_define_module_function(mCeles, "pvdpv",  celes_s_pvdpv,  2);
+	rb_define_module_function(mCeles, "pvxpv",  celes_s_pvxpv,  2);
+	rb_define_module_function(mCeles, "pvm",    celes_s_pvm,    1);
+	rb_define_module_function(mCeles, "sxpv",   celes_s_sxpv,   2);
+	rb_define_module_function(mCeles, "s2xpv",  celes_s_s2xpv,  3);
+	rb_define_module_function(mCeles, "pvu",    celes_s_pvu,    2);
+	rb_define_module_function(mCeles, "pvup",   celes_s_pvup,   2);
+	rb_define_module_function(mCeles, "rxpv",   celes_s_rxpv,   2);
+	rb_define_module_function(mCeles, "trxpv",  celes_s_trxpv,  2);
+
+	rb_define_module_function(mCeles, "anp",    celes_s_anp,    1);
+	rb_define_module_function(mCeles, "anpm",   celes_s_anpm,   1);
+	rb_define_module_function(mCeles, "a2tf",   celes_s_a2tf,   2);
+	rb_define_module_function(mCeles, "a2af",   celes_s_a2af,   2);
+	rb_define_module_function(mCeles, "af2a",   celes_s_af2a,   4);
+	rb_define_module_function(mCeles, "d2tf",   celes_s_d2tf,   2);
+	rb_define_module_function(mCeles, "tf2a",   celes_s_tf2a,   4);
+	rb_define_module_function(mCeles, "tf2d",   celes_s_tf2d,   4);
+
+	rb_define_module_function(mCeles, "cal2jd", celes_s_cal2jd, 3);
+	rb_define_module_function(mCeles, "epb",    celes_s_epb,    2);
+	rb_define_module_function(mCeles, "epb2jd", celes_s_epb2jd, 1);
+	rb_define_module_function(mCeles, "epj",    celes_s_epj,    2);
+	rb_define_module_function(mCeles, "epj2jd", celes_s_epj2jd, 1);
+	rb_define_module_function(mCeles, "jd2cal", celes_s_jd2cal, 2);
+	rb_define_module_function(mCeles, "jdcalf", celes_s_jdcalf, 3);
+
+	rb_define_module_function(mCeles, "d2dtf",  celes_s_d2dtf,  4);
+	rb_define_module_function(mCeles, "dat",    celes_s_dat,    4);
+	rb_define_module_function(mCeles, "dtdb",   celes_s_dtdb,   6);
+	rb_define_module_function(mCeles, "dtf2d",  celes_s_dtf2d,  7);
+	rb_define_module_function(mCeles, "taitt",  celes_s_taitt,  2);
+	rb_define_module_function(mCeles, "taiut1", celes_s_taiut1, 3);
+	rb_define_module_function(mCeles, "taiutc", celes_s_taiutc, 2);
+	rb_define_module_function(mCeles, "tcbtdb", celes_s_tcbtdb, 2);
+	rb_define_module_function(mCeles, "tcgtt" , celes_s_tcgtt,  2);
+	rb_define_module_function(mCeles, "tdbtcb", celes_s_tdbtcb, 2);
+	rb_define_module_function(mCeles, "tdbtt" , celes_s_tdbtt,  3);
+	rb_define_module_function(mCeles, "tttai" , celes_s_tttai,  2);
+	rb_define_module_function(mCeles, "tttcg" , celes_s_tttcg,  2);
+	rb_define_module_function(mCeles, "tttdb" , celes_s_tttdb,  3);
+	rb_define_module_function(mCeles, "ttut1" , celes_s_ttut1,  3);
+	rb_define_module_function(mCeles, "ut1tai", celes_s_ut1tai, 3);
+	rb_define_module_function(mCeles, "ut1tt",  celes_s_ut1tt,  3);
+	rb_define_module_function(mCeles, "ut1utc", celes_s_ut1utc, 3);
+	rb_define_module_function(mCeles, "utctai", celes_s_utctai, 2);
+	rb_define_module_function(mCeles, "utcut1", celes_s_utcut1, 3);
+
+	rb_define_module_function(mCeles, "ee00",   celes_s_ee00,   4);
+	rb_define_module_function(mCeles, "ee00a",  celes_s_ee00a,  2);
+	rb_define_module_function(mCeles, "ee00b",  celes_s_ee00b,  2);
+	rb_define_module_function(mCeles, "ee06a",  celes_s_ee06a,  2);
+	rb_define_module_function(mCeles, "eect00", celes_s_eect00, 2);
+	rb_define_module_function(mCeles, "eqeq94", celes_s_eqeq94, 2);
+	rb_define_module_function(mCeles, "era00",  celes_s_era00,  2);
+	rb_define_module_function(mCeles, "gmst00", celes_s_gmst00, 4);
+	rb_define_module_function(mCeles, "gmst06", celes_s_gmst06, 4);
+	rb_define_module_function(mCeles, "gmst82", celes_s_gmst82, 2);
+	rb_define_module_function(mCeles, "gst00a", celes_s_gst00a, 4);
+	rb_define_module_function(mCeles, "gst00b", celes_s_gst00b, 2);
+	rb_define_module_function(mCeles, "gst06",  celes_s_gst06,  5);
+	rb_define_module_function(mCeles, "gst06a", celes_s_gst06a, 4);
+	rb_define_module_function(mCeles, "gst94",  celes_s_gst94,  2);
+
+	rb_define_module_function(mCeles, "epv00",  celes_s_epv00,  2);
+	rb_define_module_function(mCeles, "plan94", celes_s_plan94, 3);
+
+	rb_define_module_function(mCeles, "bi00",   celes_s_bi00,   0);
+	rb_define_module_function(mCeles, "bp00",   celes_s_bp00,   2);
+	rb_define_module_function(mCeles, "bp06",   celes_s_bp06,   2);
+	rb_define_module_function(mCeles, "bpn2xy", celes_s_bpn2xy, 1);
+	rb_define_module_function(mCeles, "c2i00a", celes_s_c2i00a, 2);
+	rb_define_module_function(mCeles, "c2i00b", celes_s_c2i00b, 2);
+	rb_define_module_function(mCeles, "c2i06a", celes_s_c2i06a, 2);
+	rb_define_module_function(mCeles, "c2ibpn", celes_s_c2ibpn, 3);
+	rb_define_module_function(mCeles, "c2ixy",  celes_s_c2ixy,  4);
+	rb_define_module_function(mCeles, "c2ixys", celes_s_c2ixys, 3);
+	rb_define_module_function(mCeles, "c2t00a", celes_s_c2t00a, 6);
+	rb_define_module_function(mCeles, "c2t00b", celes_s_c2t00b, 6);
+	rb_define_module_function(mCeles, "c2t06a", celes_s_c2t06a, 6);
+	rb_define_module_function(mCeles, "c2tcio", celes_s_c2tcio, 3);
+	rb_define_module_function(mCeles, "c2teqx", celes_s_c2teqx, 3);
+	rb_define_module_function(mCeles, "c2tpe",  celes_s_c2tpe,  8);
+	rb_define_module_function(mCeles, "c2txy",  celes_s_c2txy,  8);
+	rb_define_module_function(mCeles, "eo06a",  celes_s_eo06a,  2);
+	rb_define_module_function(mCeles, "eors",   celes_s_eors,   2);
+	rb_define_module_function(mCeles, "fw2m",   celes_s_fw2m,   4);
+	rb_define_module_function(mCeles, "fw2xy",  celes_s_fw2xy,  4);
+	rb_define_module_function(mCeles, "num00a", celes_s_num00a, 2);
+	rb_define_module_function(mCeles, "num00b", celes_s_num00b, 2);
+	rb_define_module_function(mCeles, "num06a", celes_s_num06a, 2);
+	rb_define_module_function(mCeles, "numat",  celes_s_numat,  3);
+	rb_define_module_function(mCeles, "nut00a", celes_s_nut00a, 2);
+	rb_define_module_function(mCeles, "nut00b", celes_s_nut00b, 2);
+	rb_define_module_function(mCeles, "nut06a", celes_s_nut06a, 2);
+	rb_define_module_function(mCeles, "nut80",  celes_s_nut80,  2);
+	rb_define_module_function(mCeles, "nutm80", celes_s_nutm80, 2);
+	rb_define_module_function(mCeles, "obl06",  celes_s_obl06,  2);
+	rb_define_module_function(mCeles, "obl80",  celes_s_obl80,  2);
+	rb_define_module_function(mCeles, "pb06",   celes_s_pb06,   2);
+	rb_define_module_function(mCeles, "pfw06",  celes_s_pfw06,  2);
+	rb_define_module_function(mCeles, "pmat00", celes_s_pmat00, 2);
+	rb_define_module_function(mCeles, "pmat06", celes_s_pmat06, 2);
+	rb_define_module_function(mCeles, "pmat76", celes_s_pmat76, 2);
+	rb_define_module_function(mCeles, "pn00",   celes_s_pn00,   4);
+	rb_define_module_function(mCeles, "pn00a",  celes_s_pn00a,  2);
+	rb_define_module_function(mCeles, "pn00b",  celes_s_pn00b,  2);
+	rb_define_module_function(mCeles, "pn06",   celes_s_pn06,   4);
+	rb_define_module_function(mCeles, "pn06a",  celes_s_pn06a,  2);
+	rb_define_module_function(mCeles, "pnm00a", celes_s_pnm00a, 2);
+	rb_define_module_function(mCeles, "pnm00b", celes_s_pnm00b, 2);
+	rb_define_module_function(mCeles, "pnm06a", celes_s_pnm06a, 2);
+	rb_define_module_function(mCeles, "pnm80",  celes_s_pnm80,  2);
+	rb_define_module_function(mCeles, "p06e",   celes_s_p06e,   2);
+	rb_define_module_function(mCeles, "pom00",  celes_s_pom00,  3);
+	rb_define_module_function(mCeles, "pr00",   celes_s_pr00,   2);
+	rb_define_module_function(mCeles, "prec76", celes_s_prec76, 4);
+	rb_define_module_function(mCeles, "s00",    celes_s_s00,    4);
+	rb_define_module_function(mCeles, "s00a",   celes_s_s00a,   2);
+	rb_define_module_function(mCeles, "s00b",   celes_s_s00b,   2);
+	rb_define_module_function(mCeles, "s06",    celes_s_s06,    4);
+	rb_define_module_function(mCeles, "s06a",   celes_s_s06a,   2);
+	rb_define_module_function(mCeles, "sp00",   celes_s_sp00,   2);
+	rb_define_module_function(mCeles, "xy06",   celes_s_xy06,   2);
+	rb_define_module_function(mCeles, "xys00a", celes_s_xys00a, 2);
+	rb_define_module_function(mCeles, "xys00b", celes_s_xys00b, 2);
+	rb_define_module_function(mCeles, "xys06a", celes_s_xys06a, 2);
+
+	rb_define_module_function(mCeles, "fad03",  celes_s_fad03,  1);
+	rb_define_module_function(mCeles, "fae03",  celes_s_fae03,  1);
+	rb_define_module_function(mCeles, "faf03",  celes_s_faf03,  1);
+	rb_define_module_function(mCeles, "faju03", celes_s_faju03, 1);
+	rb_define_module_function(mCeles, "fal03",  celes_s_fal03,  1);
+	rb_define_module_function(mCeles, "falp03", celes_s_falp03, 1);
+	rb_define_module_function(mCeles, "fama03", celes_s_fama03, 1);
+	rb_define_module_function(mCeles, "fame03", celes_s_fame03, 1);
+	rb_define_module_function(mCeles, "fane03", celes_s_fane03, 1);
+	rb_define_module_function(mCeles, "faom03", celes_s_faom03, 1);
+	rb_define_module_function(mCeles, "fapa03", celes_s_fapa03, 1);
+	rb_define_module_function(mCeles, "fasa03", celes_s_fasa03, 1);
+	rb_define_module_function(mCeles, "faur03", celes_s_faur03, 1);
+	rb_define_module_function(mCeles, "fave03", celes_s_fave03, 1);
+
+	rb_define_module_function(mCeles, "pvstar", celes_s_pvstar, 1);
+	rb_define_module_function(mCeles, "starpv", celes_s_starpv, 6);
+
+	rb_define_module_function(mCeles, "fk52h",  celes_s_fk52h,  6);
+	rb_define_module_function(mCeles, "fk5hip", celes_s_fk5hip, 0);
+	rb_define_module_function(mCeles, "fk5hz",  celes_s_fk5hz,  4);
+	rb_define_module_function(mCeles, "h2fk5",  celes_s_h2fk5,  6);
+	rb_define_module_function(mCeles, "hfk5z",  celes_s_hfk5z,  4);
+	rb_define_module_function(mCeles, "starpm", celes_s_starpm, 10);
+
+	rb_define_module_function(mCeles, "eform",  celes_s_eform,  1);
+	rb_define_module_function(mCeles, "gc2gd",  celes_s_gc2gd,  2);
+	rb_define_module_function(mCeles, "gc2gde", celes_s_gc2gde, 3);
+	rb_define_module_function(mCeles, "gd2gc",  celes_s_gd2gc,  4);
+	rb_define_module_function(mCeles, "gd2gce", celes_s_gd2gce, 5);
+
+	rb_define_const(mCeles, "DPI",    DBL2NUM(DPI));
+	rb_define_const(mCeles, "D2PI",   DBL2NUM(D2PI));
+	rb_define_const(mCeles, "DD2R",   DBL2NUM(DD2R));
+	rb_define_const(mCeles, "DR2AS",  DBL2NUM(DR2AS));
+	rb_define_const(mCeles, "DAS2R",  DBL2NUM(DAS2R));
+	rb_define_const(mCeles, "DS2R",   DBL2NUM(DS2R));
+	rb_define_const(mCeles, "TURNAS", DBL2NUM(TURNAS));
+	rb_define_const(mCeles, "DMAS2R", DBL2NUM(DMAS2R));
+	rb_define_const(mCeles, "DTY",    DBL2NUM(DTY));
+	rb_define_const(mCeles, "DAYSEC", DBL2NUM(DAYSEC));
+	rb_define_const(mCeles, "DJY",    DBL2NUM(DJY));
+	rb_define_const(mCeles, "DJC",    DBL2NUM(DJC));
+	rb_define_const(mCeles, "DJM",    DBL2NUM(DJM));
+	rb_define_const(mCeles, "DJ00",   DBL2NUM(DJ00));
+	rb_define_const(mCeles, "DJM0",   DBL2NUM(DJM0));
+	rb_define_const(mCeles, "DJM00",  DBL2NUM(DJM00));
+	rb_define_const(mCeles, "DJM77",  DBL2NUM(DJM77));
+	rb_define_const(mCeles, "TTMTAI", DBL2NUM(TTMTAI));
+	rb_define_const(mCeles, "DAU",    DBL2NUM(DAU));
+	rb_define_const(mCeles, "DC",     DBL2NUM(DC));
+	rb_define_const(mCeles, "ELG",    DBL2NUM(ELG));
+	rb_define_const(mCeles, "ELB",    DBL2NUM(ELB));
+	rb_define_const(mCeles, "TDB0",   DBL2NUM(TDB0));
+	rb_define_const(mCeles, "WGS84",  INT2NUM(WGS84));
+	rb_define_const(mCeles, "GRS80",  INT2NUM(GRS80));
+	rb_define_const(mCeles, "WGS72",  INT2NUM(WGS72));
+}
+
+
+
+
+/*----------------------------------------------------------------------
+**
+**  Celes is a wrapper of the SOFA Library for Ruby.
+**
+**
+**
+**
+**
+**  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.
+**
+**--------------------------------------------------------------------*/