core/RefractionExtinction.hpp

/*
 * Stellarium
 * Copyright (C) 2010 Fabien Chereau
 *
 * 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.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
 *
 * Refraction and extinction computations.
 * Implementation: 2010-03-23 GZ=Georg Zotti, Georg.Zotti@univie.ac.at
 * 2010-12 FC split into 2 classes, implemented Refraction */

#ifndef _REFRACTIONEXTINCTION_HPP_
#define _REFRACTIONEXTINCTION_HPP_
// USABILITY: added 3 more flags/switches in GUI:
// Temperature [C] [influences refraction]
// Pressure [mbar]  [influences refraction]
// extinction Coeff. k=0...(0.01)...1, [if k=0, no ext. effect]
// SUGGESTION: Allow Temperature/Pressure/ex.Coeff./LightPollution set in the landscape files

#include "VecMath.hpp"
#include "StelProjector.hpp"

class Extinction
{
public:
    Extinction();
    void forward(const Vec3d *altAzPos, float *mag, const int num) const;
    void forward(const Vec3f *altAzPos, float *mag, const int num) const;
    void forward(const double *sinAlt,  float *mag, const int num) const;
    void forward(const float  *sinAlt,  float *mag, const int num) const;

    void backward(const Vec3d *altAzPos, float *mag, const int num) const;
    void backward(const Vec3f *altAzPos, float *mag, const int num) const;
    void backward(const double *sinAlt,  float *mag, const int num) const;
    void backward(const float  *sinAlt,  float *mag, const int num) const;

    void setExtinctionCoefficient(float k) { ext_coeff=k; }
    float getExtinctionCoefficient() const {return ext_coeff;}

private:
    float airmass(const float cosZ, const bool apparent_z=true) const;

    float ext_coeff;
    static const float SUBHORIZONTAL_AIRMASS;
};


class Refraction: public StelProjector::ModelViewTranform
{
public:
    Refraction();

    void forward(Vec3d& altAzPos) const;

    void backward(Vec3d& altAzPos) const;

    void forward(Vec3f& altAzPos) const;

    void backward(Vec3f& altAzPos) const;

    void combine(const Mat4d& m)
    {
        setPreTransfoMat(preTransfoMat*m);
    }

    Mat4d getApproximateLinearTransfo() const {return postTransfoMat*preTransfoMat;}

    StelProjector::ModelViewTranformP clone() const {Refraction* refr = new Refraction(); *refr=*this; return StelProjector::ModelViewTranformP(refr);}

    void setPressure(float p_mbar);
    float getPressure() const {return pressure;}

    void setTemperature(float t_C);
    float getTemperature() const {return temperature;}

    void setPreTransfoMat(const Mat4d& m);
    void setPostTransfoMat(const Mat4d& m);

private:
    void updatePrecomputed();

    float pressure;
    float temperature;
    float press_temp_corr_Saemundson;
    float press_temp_corr_Bennett;

    static const double MIN_GEO_ALTITUDE_DEG;
    static const double MIN_GEO_ALTITUDE_RAD;
    static const double MIN_GEO_ALTITUDE_SIN;
    static const double MIN_APP_ALTITUDE_DEG;
    static const double MIN_APP_ALTITUDE_RAD;
    static const double MIN_APP_ALTITUDE_SIN;
    static const float MIN_GEO_ALTITUDE_DEG_F;
    static const float MIN_GEO_ALTITUDE_RAD_F;
    static const float MIN_GEO_ALTITUDE_SIN_F;
    static const float MIN_APP_ALTITUDE_DEG_F;
    static const float MIN_APP_ALTITUDE_RAD_F;
    static const float MIN_APP_ALTITUDE_SIN_F;
    static const double TRANSITION_WIDTH_GEO_DEG;
    static const double TRANSITION_WIDTH_GEO_DEG_F;
    static const double TRANSITION_WIDTH_APP_DEG;
    static const double TRANSITION_WIDTH_APP_DEG_F;

    Mat4d preTransfoMat;
    Mat4d invertPreTransfoMat;
    Mat4f preTransfoMatf;
    Mat4f invertPreTransfoMatf;

    Mat4d postTransfoMat;
    Mat4d invertPostTransfoMat;
    Mat4f postTransfoMatf;
    Mat4f invertPostTransfoMatf;
};

#endif  // _REFRACTIONEXTINCTION_HPP_