astronomy_conversions.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: MIT
#pragma once
#ifndef _TIME_SHIELD_ASTRONOMY_CONVERSIONS_HPP_INCLUDED
#define _TIME_SHIELD_ASTRONOMY_CONVERSIONS_HPP_INCLUDED

 
#include "config.hpp"
#include "julian_conversions.hpp"
#include "MoonPhase.hpp"
 
#include <cmath>
 
namespace time_shield {

    struct MoonPhaseSineCosine {
        double phase_sin = 0.0; 
        double phase_cos = 0.0; 
        double phase_angle_rad = 0.0; 
        constexpr MoonPhaseSineCosine() = default;
        constexpr MoonPhaseSineCosine(double phase_sin_value, double phase_cos_value, double phase_angle_rad_value) noexcept
            : phase_sin(phase_sin_value),
              phase_cos(phase_cos_value),
              phase_angle_rad(phase_angle_rad_value) {}
    };

    inline double moon_phase_jd_approx(fts_t ts) noexcept {
        double temp = (static_cast<double>(fts_to_jd(ts)) - 2451550.1) / 29.530588853;
        temp = temp - std::floor(temp);
        if (temp < 0.0) temp += 1.0;
        return temp;
    }

    inline double moon_phase(fts_t ts) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.compute(static_cast<double>(ts)).phase;
    }

    inline MoonPhaseSineCosine moon_phase_sincos(fts_t ts) noexcept {
        static const astronomy::MoonPhase calculator{};
        const auto result = calculator.compute(static_cast<double>(ts));
        return MoonPhaseSineCosine{result.phase_sin, result.phase_cos, result.phase_angle_rad};
    }

    inline double moon_illumination(fts_t ts) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.compute(static_cast<double>(ts)).illumination;
    }

    inline double moon_age_days_jd_approx(fts_t ts) noexcept {
        return moon_phase_jd_approx(ts) * 29.530588853;
    }

    inline double moon_age_days(fts_t ts) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.compute(static_cast<double>(ts)).age_days;
    }

    inline astronomy::MoonQuarterInstants moon_quarters(fts_t ts) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.quarter_instants_unix(static_cast<double>(ts));
    }

    inline bool is_new_moon_window(fts_t ts, double window_seconds = astronomy::MoonPhase::kDefaultQuarterWindow_s) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.is_new_moon_window(static_cast<double>(ts), window_seconds);
    }

    inline bool is_full_moon_window(fts_t ts, double window_seconds = astronomy::MoonPhase::kDefaultQuarterWindow_s) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.is_full_moon_window(static_cast<double>(ts), window_seconds);
    }

    inline bool is_first_quarter_window(fts_t ts, double window_seconds = astronomy::MoonPhase::kDefaultQuarterWindow_s) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.is_first_quarter_window(static_cast<double>(ts), window_seconds);
    }

    inline bool is_last_quarter_window(fts_t ts, double window_seconds = astronomy::MoonPhase::kDefaultQuarterWindow_s) noexcept {
        static const astronomy::MoonPhase calculator{};
        return calculator.is_last_quarter_window(static_cast<double>(ts), window_seconds);
    }
 
} // namespace time_shield
 
#endif // _TIME_SHIELD_ASTRONOMY_CONVERSIONS_HPP_INCLUDED