time.h

Go to the documentation of this file.

#pragma once
 
#include "Math/constant.h"
 
#include <cmath>
#include <array>
#include <chrono>
#include <ratio>
 
namespace ph::math
{
 
template<typename TimeType, typename UnitType>
inline auto duration_to_HMS(const std::chrono::duration<TimeType, UnitType>& totalTime)
-> std::array<TimeType, 3>
{
    auto hours = std::chrono::duration_cast<std::chrono::hours>(totalTime);
    auto minutes = std::chrono::duration_cast<std::chrono::minutes>(totalTime - hours);
    auto seconds = std::chrono::duration_cast<std::chrono::seconds>(totalTime - hours - minutes);
 
    return {
        static_cast<TimeType>(hours.count()),
        static_cast<TimeType>(minutes.count()),
        static_cast<TimeType>(seconds.count())};
}
 
template<typename TimeType>
inline auto milliseconds_to_HMS(const TimeType& totalTime)
-> std::array<TimeType, 3>
{
    return duration_to_HMS(std::chrono::duration<TimeType, std::milli>(totalTime));
}
 
template<typename TimeType>
inline auto seconds_to_HMS(const TimeType& totalSeconds)
-> std::array<TimeType, 3>
{
    return duration_to_HMS(std::chrono::duration<TimeType>(totalSeconds));
}
 
template<typename T>
inline T mean_solar_time_to_solar_time_24H(
    const T standardTime24H,
    const T standardMeridianRadians,
    const T siteLongitudeRadians,
    const T julianDate)
{
    using namespace ph::math::constant;
 
    const T term1 = static_cast<T>(0.170) * std::sin(four_pi<T> * (julianDate - static_cast<T>(80)) / static_cast<T>(373));
    const T term2 = static_cast<T>(-0.129) * std::sin(two_pi<T> * (julianDate - static_cast<T>(8)) / static_cast<T>(355));
    const T term3 = static_cast<T>(12) * (standardMeridianRadians - siteLongitudeRadians) * rcp_pi<T>;
 
    return standardTime24H + term1 + term2 + term3;
}
 
}// end namespace ph::math