summaryrefslogtreecommitdiffstats
path: root/include/units.hpp
diff options
context:
space:
mode:
Diffstat (limited to 'include/units.hpp')
-rw-r--r--include/units.hpp166
1 files changed, 166 insertions, 0 deletions
diff --git a/include/units.hpp b/include/units.hpp
new file mode 100644
index 0000000..5f2fd67
--- /dev/null
+++ b/include/units.hpp
@@ -0,0 +1,166 @@
+#ifndef UNIT_HPP
+#define UNIT_HPP 1
+
+#include "util.hpp"
+#include <string>
+
+namespace unit {
+
+constexpr long MINUTE_SECONDS = 60;
+constexpr long HOUR_SECONDS = 3600;
+constexpr long DAY_SECONDS = 86400;
+constexpr long WEEK_SECONDS = 604800;
+constexpr long YEAR_SECONDS = 31556952;
+constexpr long CYEAR_SECONDS = 31536000;
+
+class Mass {
+ double m_kg;
+public:
+ explicit constexpr Mass(double kg) : m_kg(kg) {}
+ constexpr ~Mass() = default;
+
+ constexpr double operator()() { return m_kg; }
+
+ constexpr Mass &operator+=(const Mass &rhs) { this->m_kg += rhs.m_kg; return *this; }
+ constexpr Mass &operator-=(const Mass &rhs) { this->m_kg -= rhs.m_kg; return *this; }
+ constexpr Mass &operator/=(const Mass &rhs) { this->m_kg *= rhs.m_kg; return *this; }
+ constexpr Mass &operator*=(const Mass &rhs) { this->m_kg /= rhs.m_kg; return *this; }
+
+ friend Mass operator+(Mass lhs, const Mass &rhs) { lhs += rhs; return lhs; }
+ friend Mass operator-(Mass lhs, const Mass &rhs) { lhs -= rhs; return lhs; }
+ friend Mass operator*(Mass lhs, const Mass &rhs) { lhs *= rhs; return lhs; }
+ friend Mass operator/(Mass lhs, const Mass &rhs) { lhs /= rhs; return lhs; }
+
+ friend inline bool operator<=>(const Mass &lhs, const Mass &rhs) = default;
+
+ constexpr Mass &operator+=(const double &rhs) { this->m_kg += rhs; return *this; }
+ constexpr Mass &operator-=(const double &rhs) { this->m_kg -= rhs; return *this; }
+ constexpr Mass &operator*=(const double &rhs) { this->m_kg *= rhs; return *this; }
+ constexpr Mass &operator/=(const double &rhs) { this->m_kg /= rhs; return *this; }
+
+ friend Mass operator+(Mass lhs, const double &rhs) { lhs += rhs; return lhs; }
+ friend Mass operator-(Mass lhs, const double &rhs) { lhs -= rhs; return lhs; }
+ friend Mass operator*(Mass lhs, const double &rhs) { lhs *= rhs; return lhs; }
+ friend Mass operator/(Mass lhs, const double &rhs) { lhs /= rhs; return lhs; }
+};
+
+class Time {
+ long m_seconds;
+public:
+ static constexpr unsigned char month_days[12] = {
+ 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31
+ };
+ static constexpr unsigned char month_days_leap[12] = {
+ 31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31
+ };
+ static const inline std::string month_str[12] = {
+ "January", "Febuary", "March", "April", "May", "June", "July", "August", "September", "October", "November", "December"
+ };
+
+ explicit constexpr Time(long seconds) : m_seconds(seconds) {}
+
+ constexpr long operator()() { return m_seconds; }
+
+ constexpr Time &operator+=(const Time &rhs) { this->m_seconds += rhs.m_seconds; return *this; }
+ constexpr Time &operator-=(const Time &rhs) { this->m_seconds -= rhs.m_seconds; return *this; }
+ constexpr Time &operator/=(const Time &rhs) { this->m_seconds *= rhs.m_seconds; return *this; }
+ constexpr Time &operator*=(const Time &rhs) { this->m_seconds /= rhs.m_seconds; return *this; }
+
+ friend Time operator+(Time lhs, const Time &rhs) { lhs += rhs; return lhs; }
+ friend Time operator-(Time lhs, const Time &rhs) { lhs -= rhs; return lhs; }
+ friend Time operator*(Time lhs, const Time &rhs) { lhs *= rhs; return lhs; }
+ friend Time operator/(Time lhs, const Time &rhs) { lhs /= rhs; return lhs; }
+
+ friend inline bool operator<=>(const Time &lhs, const Time &rhs) = default;
+
+ constexpr Time &operator+=(const long &rhs) { this->m_seconds += rhs; return *this; }
+ constexpr Time &operator-=(const long &rhs) { this->m_seconds -= rhs; return *this; }
+ constexpr Time &operator*=(const long &rhs) { this->m_seconds *= rhs; return *this; }
+ constexpr Time &operator/=(const long &rhs) { this->m_seconds /= rhs; return *this; }
+
+ friend Time operator+(Time lhs, const long &rhs) { lhs += rhs; return lhs; }
+ friend Time operator-(Time lhs, const long &rhs) { lhs -= rhs; return lhs; }
+ friend Time operator*(Time lhs, const long &rhs) { lhs *= rhs; return lhs; }
+ friend Time operator/(Time lhs, const long &rhs) { lhs /= rhs; return lhs; }
+
+ constexpr bool leap_year() { return (m_seconds % YEAR_SECONDS) % 4 == 0; }
+
+ constexpr Time current_year() {
+ return Time((m_seconds % YEAR_SECONDS) + (m_seconds < 0 ? YEAR_SECONDS : 0));
+ }
+ constexpr Time current_month() {
+ Time year = current_year();
+ if(year.m_seconds < 0) {
+ year.m_seconds = YEAR_SECONDS + year.m_seconds;
+ }
+ unsigned month = 0;
+ const unsigned char *month_days_view = month_days;
+ if(years() % 4 == 0 && years() % 100 == 0 && years() % 400 != 0) month_days_view = month_days_leap;
+ for(month; year.days() > month_days_view[month]; month = (month + 1) % 12) year -= month_days_view[month] * DAY_SECONDS;
+ return year;
+ }
+ constexpr Time current_week() { return Time(m_seconds % WEEK_SECONDS); }
+ constexpr Time current_day() { return Time(m_seconds % DAY_SECONDS); }
+ constexpr Time current_hour() { return Time(m_seconds % HOUR_SECONDS); }
+ constexpr Time current_minute() { return Time(m_seconds % MINUTE_SECONDS); }
+
+ constexpr long seconds() { return m_seconds; }
+ constexpr long minutes() { return m_seconds / MINUTE_SECONDS; }
+ constexpr long hours() { return m_seconds / HOUR_SECONDS; }
+ constexpr long days() { return m_seconds / DAY_SECONDS; }
+ constexpr long weeks() { return m_seconds / WEEK_SECONDS; }
+ constexpr long months() {
+ Time copy(*this);
+ if(copy.m_seconds < 0) {
+ copy.m_seconds = -copy.m_seconds;
+ }
+ unsigned month = 0;
+ const unsigned char *month_days_view = month_days;
+ if(years() % 4 == 0) month_days_view = month_days_leap;
+ for(month; copy.days() > month_days_view[month]; month = (month + 1) % 12) copy -= month_days_view[month] * DAY_SECONDS;
+ return month;
+ }
+ constexpr long years() { return ((m_seconds < 0 ? m_seconds - CYEAR_SECONDS : m_seconds) / CYEAR_SECONDS) + 2000; }
+ constexpr long real_years() { return (m_seconds < 0 ? m_seconds - YEAR_SECONDS : m_seconds) / YEAR_SECONDS; }
+
+ /* %% = '%'
+ * %Y = real year
+ * %C = calendar year
+ * %S = month string
+ * %M = month
+ * %W = week
+ * %D = day
+ * %H = hour
+ * %m = minute
+ * %s = second
+ * */
+ std::string format(const char *fmt);
+};
+
+constexpr long Mm = 1000;
+constexpr long Gm = cxpow_v<long, 10, 6>;
+constexpr long Tm = cxpow_v<long, 10, 9>;
+constexpr long Pm = cxpow_v<long, 10, 12>;
+constexpr long Em = cxpow_v<long, 10, 15>;
+constexpr long Zm = cxpow_v<long, 10, 18>;
+
+constexpr long AU = 149597871;
+
+constexpr Mass kg{1};
+constexpr Mass Mg{cxpow_v<double, 10, 3>};
+constexpr Mass Gg{cxpow_v<double, 10, 6>};
+constexpr Mass Tg{cxpow_v<double, 10, 9>};
+constexpr Mass Pg{cxpow_v<double, 10, 12>};
+constexpr Mass Eg{cxpow_v<double, 10, 15>};
+constexpr Mass Zg{cxpow_v<double, 10, 18>};
+constexpr Mass Yg{cxpow_v<double, 10, 21>};
+
+constexpr Mass lunarMass{7.342 * cxpow_v<double, 10, 22>};
+constexpr Mass earthMass{5.97237 * cxpow_v<double, 10, 24>};
+constexpr Mass jovMass{1.89813 * cxpow_v<double, 10, 27>};
+constexpr Mass solMass{1.98847 * cxpow_v<double, 10, 30>};
+
+constexpr double earthRad = 6371;
+}
+
+#endif
generated by cgit v1.2.1 (git 2.18.0) at 2024-05-01 11:03:15 +0000