ksnumbers.h

Go to the documentation of this file.

/***************************************************************************
                          ksnumbers.h  -  description
                             -------------------
    begin                : Sun Jan 13 2002
    copyright            : (C) 2002-2005 by Jason Harris
    email                : kstars@30doradus.org
    copyright            : (C) 2004-2005 by Pablo de Vicente
    email                : p.devicente@wanadoo.es
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

#ifndef KSNUMBERS_H_
#define KSNUMBERS_H_

#define NUTTERMS 63

#include "dms.h"

class KSNumbers {
public:

    explicit KSNumbers( long double jd );

    ~KSNumbers();

    const dms* obliquity() const { return &Obliquity; }

    dms constAberr() const { return K; }

    dms sunMeanAnomaly() const { return M; }

    dms sunMeanLongitude() const { return L; }

    dms sunTrueAnomaly() const { return M0; }

    dms sunTrueLongitude() const { return L0; }

    dms earthPerihelionLongitude() const { return P; }

    double earthEccentricity() const { return e; }

    double dObliq() const { return deltaObliquity; }

    double dEcLong() const { return deltaEcLong; }

    double julianCenturies() const { return T; }

    long double julianDay() const { return days; }

    double julianMillenia() const { return jm; }

    double p1( int i1, int i2 ) const { return P1[i1][i2]; }

    double p2( int i1, int i2 ) const { return P2[i1][i2]; }

    double p1b( int i1, int i2 ) const { return P1B[i1][i2]; }

    double p2b( int i1, int i2 ) const { return P2B[i1][i2]; }

    void computeConstantValues();

    void updateValues( long double jd );

    inline long double getJD() const { return days; }

    double vEarth(int i) const {return vearth[i];}

private:
    dms Obliquity, K, L, L0, LM, M, M0, O, P, D, MM, F;
    dms XP, YP, ZP, XB, YB, ZB;
    double CX, SX, CY, SY, CZ, SZ;
    double CXB, SXB, CYB, SYB, CZB, SZB;
    double P1[3][3], P2[3][3], P1B[3][3], P2B[3][3];
    double deltaObliquity, deltaEcLong;
    double e, T, TB;
    long double days; // JD for which the last update was called
    double jm;
    static const int arguments[NUTTERMS][5];
    static const int amp[NUTTERMS][4];
    double vearth[3];
};

#endif