ootk 1.2.0 → 1.5.0

package/lib/utils/sun-math.d.ts

@@ -0,0 +1,230 @@
+/**
+ * @author Theodore Kruczek.
+ * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
+ * with satellites and other orbital objects.
+ *
+ * @file SunMath is a collection of functions for calculating the position of the sun
+ * and other stars. This was originally created by Vladimir Agafonkin. Robert Gester's
+ * update was referenced for documentation. There were a couple of bugs in both versions
+ * so there will be some differences if you are migrating from either to this library.
+ *
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2022 Theodore Kruczek
+ *
+ * @Copyright (c) 2011-2015, Vladimir Agafonkin
+ * SunCalc is a JavaScript library for calculating sun/moon position and light phases.
+ * https://github.com/mourner/suncalc
+ *
+ * Reworked and enhanced by Robert Gester
+ * @Copyright (c) 2022 Robert Gester
+ * https://github.com/hypnos3/suncalc3
+ *
+ * sun calculations are based on http://aa.quae.nl/en/reken/zonpositie.html formulas
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit 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 Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { AzEl, Meters, SunTime } from '../types/types';
+import { Degrees, Radians } from '../ootk';
+export declare class SunMath {
+    private static readonly J1970;
+    private static readonly J2000;
+    private static readonly J0;
+    private static readonly e;
+    private static times;
+    private static julianCycle;
+    /**
+     * calculates the obderver angle
+     * @param {number} alt  the observer altitude (in meters) relative to the horizon
+     * @returns {number} height for further calculations
+     */
+    private static observerAngle;
+    /**
+     * get solar mean anomaly
+     * @param {number} d - julian day
+     * @returns {number} solar mean anomaly
+     */
+    private static solarMeanAnomaly;
+    /**
+     * convert Julian calendar to date object
+     * @param {number} julian day number in Julian calendar to convert
+     * @return {number} result date as timestamp
+     */
+    static julian2date(julian: number): Date;
+    /**
+     * get number of days for a dateValue since 2000
+     *
+     * See: https://en.wikipedia.org/wiki/Epoch_(astronomy)
+     *
+     * @param {number} date date as timestamp to get days
+     * @return {number} count of days
+     */
+    static date2jSince2000(date: Date): number;
+    /**
+     * get right ascension
+     * @param {number} l - ecliptic longitude
+     * @param {number} b - ecliptic latitude
+     * @returns {number} right ascension
+     */
+    static rightAscension(l: number, b: number): Radians;
+    /**
+     * get declination
+     * @param {number} l - ecliptic longitude
+     * @param {number} b - ecliptic latitude
+     * @returns {number} declination
+     */
+    static declination(l: number, b: number): Radians;
+    /**
+     * get azimuth
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} azimuth in rad
+     */
+    static azimuth(H: number, phi: Radians, dec: Radians): Radians;
+    /**
+     * get elevation (sometimes called altitude)
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} elevation
+     */
+    static elevation(H: number, phi: Radians, dec: Radians): Radians;
+    /**
+     * side real time
+     * @param {number} d - julian day
+     * @param {Radians} lw - longitude of the observer
+     * @returns {number} sidereal time
+     */
+    static siderealTime(d: number, lw: Radians): number;
+    /**
+     * get astro refraction
+     * @param {Radians} h - elevation
+     * @returns {number} refraction
+     */
+    static astroRefraction(h: Degrees): Radians;
+    /**
+     * ecliptic longitude
+     * @param {number} M - solar mean anomaly
+     * @returns {number} ecliptic longitude
+     */
+    static eclipticLongitude(M: number): number;
+    /**
+     * ecliptic latitude
+     * @param {number} B - ?
+     * @returns {number} ecliptic latitude
+     */
+    static eclipticLatitude(B: number): number;
+    /**
+     * Julian cycle
+     * @param {number} d - number of days
+     * @param {number} lw - rad * -lng;
+     * @returns {number} julian cycle
+     */
+    static julianCyle(d: number, lw: number): number;
+    /**
+     * approx transit
+     * @param {number} Ht - hourAngle
+     * @param {number} lw - rad * -lng
+     * @param {number} n - Julian cycle
+     * @returns {number} approx transit
+     */
+    static approxTransit(Ht: number, lw: number, n: number): number;
+    /**
+     * solar transit in Julian
+     * @param {number} ds - approxTransit
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @returns {number} solar transit in Julian
+     */
+    static solarTransitJulian(ds: number, M: number, L: number): number;
+    /**
+     * hour angle
+     * @param {number} h - heigh at 0
+     * @param {number} phi -  rad * lat;
+     * @param {number} dec - declination
+     * @returns {number} hour angle
+     */
+    static hourAngle(h: number, phi: number, dec: number): number;
+    /**
+     * returns set time for the given sun altitude
+     * @param {number} h - heigh at 0
+     * @param {number} lw - rad * -lng
+     * @param {number} phi -  rad * lat;
+     * @param {number} dec - declination
+     * @param {number} n - Julian cycle
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @returns {number} set time
+     */
+    static getSetJulian(h: number, lw: number, phi: number, dec: number, n: number, M: number, L: number): number;
+    /**
+     * returns set time for the given sun altitude
+     * @param {Meters} alt - altitude at 0
+     * @param {Radians} lw - -lng
+     * @param {Radians} phi - lat;
+     * @param {Radians} dec - declination
+     * @param {number} n - Julian cycle
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @return {J2000Value} sunset time in days since 2000
+     */
+    private static getSetJ;
+    /**
+     * calculates sun times for a given date and latitude/longitude
+     * @param {Date|number} dateVal Date object for calculating sun-times
+     * @param {number} lat latitude for calculating sun-times
+     * @param {number} lon longitude for calculating sun-times
+     * @param {number} [alt=0]  the observer height (in meters) relative to the horizon
+     * @param {boolean} [isUtc=false] defines if the calculation should be in utc or local time (default is local)
+     * @return {SunTime} result object of sunTime
+     */
+    static getTimes(dateVal: Date, lat: Degrees, lon: Degrees, alt?: Meters, isUtc?: boolean): SunTime;
+    private static calculateJnoon;
+    /**
+     * calculates time for a given azimuth angle for a given date and latitude/longitude
+     * @param {number|Date} dateValue Date object or timestamp for calculating sun-time
+     * @param {Degrees} lat latitude for calculating sun-time
+     * @param {Degrees} lon longitude for calculating sun-time
+     * @param {Radians | Degrees} az azimuth for calculating sun-time
+     * @param {boolean} [isDegrees] true if the az angle is in degree and not in rad
+     * @return {Date} result time of sun-time
+     */
+    static getSunTimeByAzimuth(dateValue: any, lat: Degrees, lon: Degrees, az: Radians | Degrees, isDegrees?: boolean): Date;
+    /**
+     * get azimuth
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} azimuth
+     */
+    private static getAzimuth;
+    /**
+     * Convert Date to Solar Date/Time
+     *
+     * Based on https://www.pveducation.org/pvcdrom/properties-of-sunlight/solar-time
+     *
+     * See: https://github.com/mourner/suncalc/pull/156
+     *
+     * @param {date} date - Date to calculate for
+     * @param {number} utcOffset - UTC offset in hours
+     * @param {Degrees} lon - Longitude in degrees
+     * @returns {Date} - Date Object in Solar Time
+     */
+    static getSolarTime(date: Date, utcOffset: number, lon: number): Date;
+    static getStarAzEl(date: Date, lat: Degrees, lon: Degrees, ra: number, dec: number): AzEl;
+    static getSunAzEl(date: Date, lat: Degrees, lon: Degrees): AzEl;
+    static getSunRaDec(d: number): {
+        dec: number;
+        ra: number;
+    };
+}

package/lib/utils/sun-math.js

@@ -0,0 +1,421 @@
+/**
+ * @author Theodore Kruczek.
+ * @description Orbital Object ToolKit (OOTK) is a collection of tools for working
+ * with satellites and other orbital objects.
+ *
+ * @file SunMath is a collection of functions for calculating the position of the sun
+ * and other stars. This was originally created by Vladimir Agafonkin. Robert Gester's
+ * update was referenced for documentation. There were a couple of bugs in both versions
+ * so there will be some differences if you are migrating from either to this library.
+ *
+ *
+ * @license AGPL-3.0-or-later
+ * @Copyright (c) 2020-2022 Theodore Kruczek
+ *
+ * @Copyright (c) 2011-2015, Vladimir Agafonkin
+ * SunCalc is a JavaScript library for calculating sun/moon position and light phases.
+ * https://github.com/mourner/suncalc
+ *
+ * Reworked and enhanced by Robert Gester
+ * @Copyright (c) 2022 Robert Gester
+ * https://github.com/hypnos3/suncalc3
+ *
+ * sun calculations are based on http://aa.quae.nl/en/reken/zonpositie.html formulas
+ *
+ * Orbital Object ToolKit is free software: you can redistribute it and/or modify it under the
+ * terms of the GNU Affero General Public License as published by the Free Software
+ * Foundation, either version 3 of the License, or (at your option) any later version.
+ *
+ * Orbital Object ToolKit 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 Affero General Public License for more details.
+ *
+ * You should have received a copy of the GNU Affero General Public License along with
+ * Orbital Object ToolKit. If not, see <http://www.gnu.org/licenses/>.
+ */
+import { DEG2RAD, MS_PER_DAY } from './constants';
+var TAU = Math.PI * 2; // https://tauday.com/tau-manifesto
+var SunMath = /** @class */ (function () {
+    function SunMath() {
+    }
+    // eslint-disable-next-line @typescript-eslint/no-unused-vars
+    SunMath.julianCycle = function (d, lw) {
+        var lonOffset = lw / TAU;
+        return Math.round(d - SunMath.J0 - lonOffset);
+    };
+    /**
+     * calculates the obderver angle
+     * @param {number} alt  the observer altitude (in meters) relative to the horizon
+     * @returns {number} height for further calculations
+     */
+    SunMath.observerAngle = function (alt) {
+        return (-2.076 * Math.sqrt(alt)) / 60;
+    };
+    /**
+     * get solar mean anomaly
+     * @param {number} d - julian day
+     * @returns {number} solar mean anomaly
+     */
+    SunMath.solarMeanAnomaly = function (d) {
+        return DEG2RAD * (357.5291 + 0.98560028 * d);
+    };
+    /**
+     * convert Julian calendar to date object
+     * @param {number} julian day number in Julian calendar to convert
+     * @return {number} result date as timestamp
+     */
+    SunMath.julian2date = function (julian) {
+        return new Date((julian - SunMath.J1970) * MS_PER_DAY);
+    };
+    /**
+     * get number of days for a dateValue since 2000
+     *
+     * See: https://en.wikipedia.org/wiki/Epoch_(astronomy)
+     *
+     * @param {number} date date as timestamp to get days
+     * @return {number} count of days
+     */
+    SunMath.date2jSince2000 = function (date) {
+        return date.getTime() / MS_PER_DAY + SunMath.J1970 - SunMath.J2000;
+    };
+    /**
+     * get right ascension
+     * @param {number} l - ecliptic longitude
+     * @param {number} b - ecliptic latitude
+     * @returns {number} right ascension
+     */
+    SunMath.rightAscension = function (l, b) {
+        return Math.atan2(Math.sin(l) * Math.cos(SunMath.e) - Math.tan(b) * Math.sin(SunMath.e), Math.cos(l));
+    };
+    /**
+     * get declination
+     * @param {number} l - ecliptic longitude
+     * @param {number} b - ecliptic latitude
+     * @returns {number} declination
+     */
+    SunMath.declination = function (l, b) {
+        return Math.asin(Math.sin(b) * Math.cos(SunMath.e) + Math.cos(b) * Math.sin(SunMath.e) * Math.sin(l));
+    };
+    /**
+     * get azimuth
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} azimuth in rad
+     */
+    SunMath.azimuth = function (H, phi, dec) {
+        return Math.PI + Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(dec) * Math.cos(phi));
+    };
+    /**
+     * get elevation (sometimes called altitude)
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} elevation
+     */
+    SunMath.elevation = function (H, phi, dec) {
+        return Math.asin(Math.sin(phi) * Math.sin(dec) + Math.cos(phi) * Math.cos(dec) * Math.cos(H));
+    };
+    /**
+     * side real time
+     * @param {number} d - julian day
+     * @param {Radians} lw - longitude of the observer
+     * @returns {number} sidereal time
+     */
+    SunMath.siderealTime = function (d, lw) {
+        return DEG2RAD * (280.16 + 360.9856235 * d) - lw;
+    };
+    /**
+     * get astro refraction
+     * @param {Radians} h - elevation
+     * @returns {number} refraction
+     */
+    SunMath.astroRefraction = function (h) {
+        if (h < 0) {
+            h = 0;
+        }
+        return 0.0002967 / Math.tan(h + 0.00312536 / (h + 0.08901179));
+    };
+    /**
+     * ecliptic longitude
+     * @param {number} M - solar mean anomaly
+     * @returns {number} ecliptic longitude
+     */
+    SunMath.eclipticLongitude = function (M) {
+        var C = DEG2RAD * (1.9148 * Math.sin(M) + 0.02 * Math.sin(2 * M) + 0.0003 * Math.sin(3 * M));
+        var P = DEG2RAD * 102.9372; // perihelion of Earth
+        return M + C + P + Math.PI; // Sun's mean longitude
+    };
+    /**
+     * ecliptic latitude
+     * @param {number} B - ?
+     * @returns {number} ecliptic latitude
+     */
+    SunMath.eclipticLatitude = function (B) {
+        var C = TAU / 360;
+        var L = B - 0.00569 - 0.00478 * Math.sin(C * B);
+        return TAU * (L + 0.0003 * Math.sin(C * 2 * L));
+    };
+    /**
+     * Julian cycle
+     * @param {number} d - number of days
+     * @param {number} lw - rad * -lng;
+     * @returns {number} julian cycle
+     */
+    SunMath.julianCyle = function (d, lw) {
+        return Math.round(d - SunMath.J0 - lw / ((2 * TAU) / 2));
+    };
+    /**
+     * approx transit
+     * @param {number} Ht - hourAngle
+     * @param {number} lw - rad * -lng
+     * @param {number} n - Julian cycle
+     * @returns {number} approx transit
+     */
+    SunMath.approxTransit = function (Ht, lw, n) {
+        return SunMath.J0 + (Ht + lw) / TAU + n;
+    };
+    /**
+     * solar transit in Julian
+     * @param {number} ds - approxTransit
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @returns {number} solar transit in Julian
+     */
+    SunMath.solarTransitJulian = function (ds, M, L) {
+        return SunMath.J2000 + ds + 0.0053 * Math.sin(M) - 0.0069 * Math.sin(2 * L);
+    };
+    /**
+     * hour angle
+     * @param {number} h - heigh at 0
+     * @param {number} phi -  rad * lat;
+     * @param {number} dec - declination
+     * @returns {number} hour angle
+     */
+    SunMath.hourAngle = function (h, phi, dec) {
+        return Math.acos((Math.sin(h) - Math.sin(phi) * Math.sin(dec)) / (Math.cos(phi) * Math.cos(dec)));
+    };
+    /**
+     * returns set time for the given sun altitude
+     * @param {number} h - heigh at 0
+     * @param {number} lw - rad * -lng
+     * @param {number} phi -  rad * lat;
+     * @param {number} dec - declination
+     * @param {number} n - Julian cycle
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @returns {number} set time
+     */
+    SunMath.getSetJulian = function (h, lw, phi, dec, n, M, L) {
+        var w = SunMath.hourAngle(h, phi, dec);
+        var a = SunMath.approxTransit(w, lw, n);
+        return SunMath.solarTransitJulian(a, M, L);
+    };
+    /**
+     * returns set time for the given sun altitude
+     * @param {Meters} alt - altitude at 0
+     * @param {Radians} lw - -lng
+     * @param {Radians} phi - lat;
+     * @param {Radians} dec - declination
+     * @param {number} n - Julian cycle
+     * @param {number} M - solar mean anomal
+     * @param {number} L - ecliptic longitude
+     * @return {J2000Value} sunset time in days since 2000
+     */
+    SunMath.getSetJ = function (alt, lw, phi, dec, n, M, L) {
+        var w = SunMath.hourAngle(alt, phi, dec);
+        var a = SunMath.approxTransit(w, lw, n);
+        return SunMath.solarTransitJulian(a, M, L);
+    };
+    /**
+     * calculates sun times for a given date and latitude/longitude
+     * @param {Date|number} dateVal Date object for calculating sun-times
+     * @param {number} lat latitude for calculating sun-times
+     * @param {number} lon longitude for calculating sun-times
+     * @param {number} [alt=0]  the observer height (in meters) relative to the horizon
+     * @param {boolean} [isUtc=false] defines if the calculation should be in utc or local time (default is local)
+     * @return {SunTime} result object of sunTime
+     */
+    SunMath.getTimes = function (dateVal, lat, lon, alt, isUtc) {
+        if (alt === void 0) { alt = 0; }
+        if (isUtc === void 0) { isUtc = false; }
+        if (isNaN(lat)) {
+            throw new Error('latitude missing');
+        }
+        if (isNaN(lon)) {
+            throw new Error('longitude missing');
+        }
+        var date = new Date(dateVal);
+        if (isUtc) {
+            date.setUTCHours(12, 0, 0, 0);
+        }
+        else {
+            date.setHours(12, 0, 0, 0);
+        }
+        var i = 0;
+        var len = 0;
+        var time = [];
+        var h0 = 0;
+        var Jset = 0;
+        var Jrise = 0;
+        var _a = SunMath.calculateJnoon(lon, lat, alt, date), Jnoon = _a.Jnoon, dh = _a.dh, lw = _a.lw, phi = _a.phi, dec = _a.dec, n = _a.n, M = _a.M, L = _a.L;
+        // Determine when the sun is at its highest and lowest (darkest) points.
+        var result = {
+            solarNoon: SunMath.julian2date(Jnoon),
+            nadir: SunMath.julian2date(Jnoon + 0.5), // https://github.com/mourner/suncalc/pull/125
+        };
+        // Add all other unique times using Jnoon as a reference
+        for (i = 0, len = SunMath.times.length; i < len; i += 1) {
+            time = SunMath.times[i];
+            h0 = (time[0] + dh) * DEG2RAD;
+            Jset = SunMath.getSetJ(h0, lw, phi, dec, n, M, L);
+            Jrise = Jnoon - (Jset - Jnoon);
+            result[time[1]] = SunMath.julian2date(Jrise);
+            result[time[2]] = SunMath.julian2date(Jset);
+        }
+        return result;
+    };
+    SunMath.calculateJnoon = function (lon, lat, alt, date) {
+        var lw = DEG2RAD * -lon;
+        var phi = DEG2RAD * lat;
+        var dh = SunMath.observerAngle(alt);
+        var d = SunMath.date2jSince2000(date);
+        var n = SunMath.julianCycle(d, lw);
+        var ds = SunMath.approxTransit(0, lw, n);
+        var M = SunMath.solarMeanAnomaly(ds);
+        var L = SunMath.eclipticLongitude(M);
+        var dec = SunMath.declination(L, 0);
+        var Jnoon = SunMath.solarTransitJulian(ds, M, L);
+        return { Jnoon: Jnoon, dh: dh, lw: lw, phi: phi, dec: dec, n: n, M: M, L: L };
+    };
+    /**
+     * calculates time for a given azimuth angle for a given date and latitude/longitude
+     * @param {number|Date} dateValue Date object or timestamp for calculating sun-time
+     * @param {Degrees} lat latitude for calculating sun-time
+     * @param {Degrees} lon longitude for calculating sun-time
+     * @param {Radians | Degrees} az azimuth for calculating sun-time
+     * @param {boolean} [isDegrees] true if the az angle is in degree and not in rad
+     * @return {Date} result time of sun-time
+     */
+    SunMath.getSunTimeByAzimuth = function (dateValue, lat, lon, az, isDegrees) {
+        if (isDegrees === void 0) { isDegrees = false; }
+        if (isNaN(az)) {
+            throw new Error('azimuth missing');
+        }
+        if (isNaN(lat)) {
+            throw new Error('latitude missing');
+        }
+        if (isNaN(lon)) {
+            throw new Error('longitude missing');
+        }
+        if (isDegrees) {
+            az *= DEG2RAD;
+        }
+        var date = new Date(dateValue);
+        var lw = DEG2RAD * -lon;
+        var phi = DEG2RAD * lat;
+        var dateVal = new Date(date.getFullYear(), date.getMonth(), date.getDate(), 0, 0, 0).getTime();
+        var addval = MS_PER_DAY; // / 2);
+        dateVal += addval;
+        while (addval > 200) {
+            // let nazi = this.getPosition(dateVal, lat, lng).azimuth;
+            var d = SunMath.date2jSince2000(new Date(dateVal));
+            var c = SunMath.getSunRaDec(d);
+            var H = SunMath.siderealTime(d, lw) - c.ra;
+            var newAz = SunMath.getAzimuth(H, phi, c.dec);
+            addval /= 2;
+            if (newAz < az) {
+                dateVal += addval;
+            }
+            else {
+                dateVal -= addval;
+            }
+        }
+        return new Date(Math.floor(dateVal));
+    };
+    /**
+     * get azimuth
+     * @param {number} H - siderealTime
+     * @param {Radians} phi - latitude
+     * @param {Radians} dec - The declination of the sun
+     * @returns {Radians} azimuth
+     */
+    SunMath.getAzimuth = function (H, phi, dec) {
+        return Math.atan2(Math.sin(H), Math.cos(H) * Math.sin(phi) - Math.tan(dec) * Math.cos(phi)) + Math.PI;
+    };
+    /**
+     * Convert Date to Solar Date/Time
+     *
+     * Based on https://www.pveducation.org/pvcdrom/properties-of-sunlight/solar-time
+     *
+     * See: https://github.com/mourner/suncalc/pull/156
+     *
+     * @param {date} date - Date to calculate for
+     * @param {number} utcOffset - UTC offset in hours
+     * @param {Degrees} lon - Longitude in degrees
+     * @returns {Date} - Date Object in Solar Time
+     */
+    SunMath.getSolarTime = function (date, utcOffset, lon) {
+        // calculate the day of year
+        var start = new Date();
+        start.setUTCFullYear(date.getUTCFullYear(), 0, 1);
+        start.setUTCHours(0, 0, 0, 0);
+        var diff = date.getTime() - start.getTime();
+        var dayOfYear = Math.floor(diff / MS_PER_DAY);
+        var b = (360 / 365) * (dayOfYear - 81) * DEG2RAD;
+        var equationOfTime = 9.87 * Math.sin(2 * b) - 7.53 * Math.cos(b) - 1.5 * Math.sin(b);
+        var localSolarTimeMeridian = 15 * utcOffset;
+        var timeCorrection = equationOfTime + 4 * (lon - localSolarTimeMeridian);
+        return new Date(date.getTime() + timeCorrection * 60 * 1000);
+    };
+    SunMath.getStarAzEl = function (date, lat, lon, ra, dec) {
+        var lw = -lon * DEG2RAD;
+        var phi = lat * DEG2RAD;
+        var d = SunMath.date2jSince2000(date);
+        var H = SunMath.siderealTime(d, lw) - (ra / 12) * Math.PI;
+        var h = SunMath.elevation(H, phi, (dec / 180) * Math.PI);
+        h += SunMath.astroRefraction(h); // elevation correction for refraction
+        return {
+            az: SunMath.azimuth(H, phi, (dec / 180) * Math.PI),
+            el: h,
+        };
+    };
+    SunMath.getSunAzEl = function (date, lat, lon) {
+        var lw = -lon * DEG2RAD;
+        var phi = lat * DEG2RAD;
+        var d = SunMath.date2jSince2000(date);
+        var c = SunMath.getSunRaDec(d);
+        var H = SunMath.siderealTime(d, lw) - c.ra;
+        return {
+            az: SunMath.azimuth(H, phi, c.dec),
+            el: SunMath.elevation(H, phi, c.dec),
+        };
+    };
+    SunMath.getSunRaDec = function (d) {
+        var M = SunMath.solarMeanAnomaly(d);
+        var L = SunMath.eclipticLongitude(M);
+        return {
+            dec: SunMath.declination(L, 0),
+            ra: SunMath.rightAscension(L, 0),
+        };
+    };
+    SunMath.J1970 = 2440587.5;
+    SunMath.J2000 = 2451545;
+    SunMath.J0 = 0.0009;
+    SunMath.e = DEG2RAD * 23.4397;
+    SunMath.times = [
+        [6, 'goldenHourDawnEnd', 'goldenHourDuskStart'],
+        [-0.3, 'sunriseEnd', 'sunsetStart'],
+        [-0.833, 'sunriseStart', 'sunsetEnd'],
+        [-1, 'goldenHourDawnStart', 'goldenHourDuskEnd'],
+        [-4, 'blueHourDawnEnd', 'blueHourDuskStart'],
+        [-6, 'civilDawn', 'civilDusk'],
+        [-8, 'blueHourDawnStart', 'blueHourDuskEnd'],
+        [-12, 'nauticalDawn', 'nauticalDusk'],
+        [-15, 'amateurDawn', 'amateurDusk'],
+        [-18, 'astronomicalDawn', 'astronomicalDusk'], // ASTRO_DAWN
+    ];
+    return SunMath;
+}());
+export { SunMath };
+//# sourceMappingURL=sun-math.js.map

package/lib/utils/sun-math.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"sun-math.js","sourceRoot":"","sources":["../../src/utils/sun-math.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAkCG;AAKH,OAAO,EAAE,OAAO,EAAE,UAAU,EAAE,MAAM,aAAa,CAAC;AAIlD,IAAM,GAAG,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,mCAAmC;AAE5D;IAAA;IAubA,CAAC;IApaC,6DAA6D;IAC9C,mBAAW,GAA1B,UAA2B,CAAS,EAAE,EAAU;QAC9C,IAAM,SAAS,GAAG,EAAE,GAAG,GAAG,CAAC;QAE3B,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,SAAS,CAAC,CAAC;IAChD,CAAC;IAED;;;;OAIG;IACY,qBAAa,GAA5B,UAA6B,GAAW;QACtC,OAAO,CAAC,CAAC,KAAK,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IACxC,CAAC;IAED;;;;OAIG;IACY,wBAAgB,GAA/B,UAAgC,CAAS;QACvC,OAAO,OAAO,GAAG,CAAC,QAAQ,GAAG,UAAU,GAAG,CAAC,CAAC,CAAC;IAC/C,CAAC;IAED;;;;OAIG;IACI,mBAAW,GAAlB,UAAmB,MAAc;QAC/B,OAAO,IAAI,IAAI,CAAC,CAAC,MAAM,GAAG,OAAO,CAAC,KAAK,CAAC,GAAG,UAAU,CAAC,CAAC;IACzD,CAAC;IAED;;;;;;;OAOG;IACI,uBAAe,GAAtB,UAAuB,IAAU;QAC/B,OAAO,IAAI,CAAC,OAAO,EAAE,GAAG,UAAU,GAAG,OAAO,CAAC,KAAK,GAAG,OAAO,CAAC,KAAK,CAAC;IACrE,CAAC;IAED;;;;;OAKG;IACI,sBAAc,GAArB,UAAsB,CAAS,EAAE,CAAS;QACxC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACxG,CAAC;IAED;;;;;OAKG;IACI,mBAAW,GAAlB,UAAmB,CAAS,EAAE,CAAS;QACrC,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACxG,CAAC;IAED;;;;;;OAMG;IACI,eAAO,GAAd,UAAe,CAAS,EAAE,GAAY,EAAE,GAAY;QAClD,OAAO,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;IACxG,CAAC;IAED;;;;;;OAMG;IACI,iBAAS,GAAhB,UAAiB,CAAS,EAAE,GAAY,EAAE,GAAY;QACpD,OAAO,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAChG,CAAC;IAED;;;;;OAKG;IACI,oBAAY,GAAnB,UAAoB,CAAS,EAAE,EAAW;QACxC,OAAO,OAAO,GAAG,CAAC,MAAM,GAAG,WAAW,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IACnD,CAAC;IAED;;;;OAIG;IACI,uBAAe,GAAtB,UAAuB,CAAU;QAC/B,IAAI,CAAC,GAAG,CAAC,EAAE;YACT,CAAC,GAAG,CAAC,CAAC;SACP;QAED,OAAO,SAAS,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,UAAU,GAAG,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC;IACjE,CAAC;IAED;;;;OAIG;IACI,yBAAiB,GAAxB,UAAyB,CAAS;QAChC,IAAM,CAAC,GAAG,OAAO,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QAC/F,IAAM,CAAC,GAAG,OAAO,GAAG,QAAQ,CAAC,CAAC,sBAAsB;QAEpD,OAAO,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC,uBAAuB;IACrD,CAAC;IAED;;;;OAIG;IACI,wBAAgB,GAAvB,UAAwB,CAAS;QAC/B,IAAM,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC;QACpB,IAAM,CAAC,GAAG,CAAC,GAAG,OAAO,GAAG,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAElD,OAAO,GAAG,GAAG,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;IAED;;;;;OAKG;IACI,kBAAU,GAAjB,UAAkB,CAAS,EAAE,EAAU;QACrC,OAAO,IAAI,CAAC,KAAK,CAAC,CAAC,GAAG,OAAO,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC3D,CAAC;IAED;;;;;;OAMG;IACI,qBAAa,GAApB,UAAqB,EAAU,EAAE,EAAU,EAAE,CAAS;QACpD,OAAO,OAAO,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,CAAC;IAC1C,CAAC;IAED;;;;;;OAMG;IACI,0BAAkB,GAAzB,UAA0B,EAAU,EAAE,CAAS,EAAE,CAAS;QACxD,OAAO,OAAO,CAAC,KAAK,GAAG,EAAE,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;IAC9E,CAAC;IAED;;;;;;OAMG;IACI,iBAAS,GAAhB,UAAiB,CAAS,EAAE,GAAW,EAAE,GAAW;QAClD,OAAO,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IACpG,CAAC;IAED;;;;;;;;;;OAUG;IACI,oBAAY,GAAnB,UAAoB,CAAS,EAAE,EAAU,EAAE,GAAW,EAAE,GAAW,EAAE,CAAS,EAAE,CAAS,EAAE,CAAS;QAClG,IAAM,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QACzC,IAAM,CAAC,GAAG,OAAO,CAAC,aAAa,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAE1C,OAAO,OAAO,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,CAAC;IAED;;;;;;;;;;OAUG;IACY,eAAO,GAAtB,UACE,GAAW,EACX,EAAW,EACX,GAAY,EACZ,GAAY,EACZ,CAAS,EACT,CAAS,EACT,CAAU;QAEV,IAAM,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;QAC3C,IAAM,CAAC,GAAG,OAAO,CAAC,aAAa,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAE1C,OAAO,OAAO,CAAC,kBAAkB,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;IAC7C,CAAC;IAED;;;;;;;;OAQG;IACI,gBAAQ,GAAf,UAAgB,OAAa,EAAE,GAAY,EAAE,GAAY,EAAE,GAAe,EAAE,KAAa;QAA9B,oBAAA,EAAA,OAAe;QAAE,sBAAA,EAAA,aAAa;QACvF,IAAI,KAAK,CAAC,GAAG,CAAC,EAAE;YACd,MAAM,IAAI,KAAK,CAAC,kBAAkB,CAAC,CAAC;SACrC;QACD,IAAI,KAAK,CAAC,GAAG,CAAC,EAAE;YACd,MAAM,IAAI,KAAK,CAAC,mBAAmB,CAAC,CAAC;SACtC;QAED,IAAM,IAAI,GAAG,IAAI,IAAI,CAAC,OAAO,CAAC,CAAC;QAE/B,IAAI,KAAK,EAAE;YACT,IAAI,CAAC,WAAW,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;SAC/B;aAAM;YACL,IAAI,CAAC,QAAQ,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;SAC5B;QAED,IAAI,CAAC,GAAG,CAAC,CAAC;QACV,IAAI,GAAG,GAAG,CAAC,CAAC;QACZ,IAAI,IAAI,GAAG,EAAE,CAAC;QACd,IAAI,EAAE,GAAG,CAAC,CAAC;QACX,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,IAAI,KAAK,GAAG,CAAC,CAAC;QAER,IAAA,KAAuC,OAAO,CAAC,cAAc,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,IAAI,CAAC,EAAhF,KAAK,WAAA,EAAE,EAAE,QAAA,EAAE,EAAE,QAAA,EAAE,GAAG,SAAA,EAAE,GAAG,SAAA,EAAE,CAAC,OAAA,EAAE,CAAC,OAAA,EAAE,CAAC,OAAgD,CAAC;QAEzF,wEAAwE;QACxE,IAAM,MAAM,GAAG;YACb,SAAS,EAAE,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC;YACrC,KAAK,EAAE,OAAO,CAAC,WAAW,CAAC,KAAK,GAAG,GAAG,CAAC,EAAE,8CAA8C;SACxF,CAAC;QAEF,wDAAwD;QACxD,KAAK,CAAC,GAAG,CAAC,EAAE,GAAG,GAAG,OAAO,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE;YACvD,IAAI,GAAG,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;YACxB,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC;YAE9B,IAAI,GAAG,OAAO,CAAC,OAAO,CAAC,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;YAClD,KAAK,GAAG,KAAK,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC,CAAC;YAE/B,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,KAAK,CAAC,CAAC;YAC7C,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,IAAI,CAAC,CAAC;SAC7C;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAEc,sBAAc,GAA7B,UAA8B,GAAW,EAAE,GAAW,EAAE,GAAW,EAAE,IAAU;QAC7E,IAAM,EAAE,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;QAC1B,IAAM,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;QAC1B,IAAM,EAAE,GAAG,OAAO,CAAC,aAAa,CAAC,GAAG,CAAC,CAAC;QACtC,IAAM,CAAC,GAAG,OAAO,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC;QAExC,IAAM,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QACrC,IAAM,EAAE,GAAG,OAAO,CAAC,aAAa,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;QAC3C,IAAM,CAAC,GAAG,OAAO,CAAC,gBAAgB,CAAC,EAAE,CAAC,CAAC;QACvC,IAAM,CAAC,GAAG,OAAO,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;QACvC,IAAM,GAAG,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACtC,IAAM,KAAK,GAAG,OAAO,CAAC,kBAAkB,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAEnD,OAAO,EAAE,KAAK,OAAA,EAAE,EAAE,IAAA,EAAE,EAAE,IAAA,EAAE,GAAG,KAAA,EAAE,GAAG,KAAA,EAAE,CAAC,GAAA,EAAE,CAAC,GAAA,EAAE,CAAC,GAAA,EAAE,CAAC;IAC9C,CAAC;IAED;;;;;;;;OAQG;IACI,2BAAmB,GAA1B,UAA2B,SAAS,EAAE,GAAY,EAAE,GAAY,EAAE,EAAqB,EAAE,SAAiB;QAAjB,0BAAA,EAAA,iBAAiB;QACxG,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE;YACb,MAAM,IAAI,KAAK,CAAC,iBAAiB,CAAC,CAAC;SACpC;QACD,IAAI,KAAK,CAAC,GAAG,CAAC,EAAE;YACd,MAAM,IAAI,KAAK,CAAC,kBAAkB,CAAC,CAAC;SACrC;QACD,IAAI,KAAK,CAAC,GAAG,CAAC,EAAE;YACd,MAAM,IAAI,KAAK,CAAC,mBAAmB,CAAC,CAAC;SACtC;QAED,IAAI,SAAS,EAAE;YACb,EAAE,IAAI,OAAO,CAAC;SACf;QACD,IAAM,IAAI,GAAG,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC;QACjC,IAAM,EAAE,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC;QAC1B,IAAM,GAAG,GAAG,OAAO,GAAG,GAAG,CAAC;QAE1B,IAAI,OAAO,GAAG,IAAI,IAAI,CAAC,IAAI,CAAC,WAAW,EAAE,EAAE,IAAI,CAAC,QAAQ,EAAE,EAAE,IAAI,CAAC,OAAO,EAAE,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,OAAO,EAAE,CAAC;QAC/F,IAAI,MAAM,GAAG,UAAU,CAAC,CAAC,QAAQ;QAEjC,OAAO,IAAI,MAAM,CAAC;QAElB,OAAO,MAAM,GAAG,GAAG,EAAE;YACnB,0DAA0D;YAC1D,IAAM,CAAC,GAAG,OAAO,CAAC,eAAe,CAAC,IAAI,IAAI,CAAC,OAAO,CAAC,CAAC,CAAC;YACrD,IAAM,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;YACjC,IAAM,CAAC,GAAG,OAAO,CAAC,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YAC7C,IAAM,KAAK,GAAG,OAAO,CAAC,UAAU,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YAEhD,MAAM,IAAI,CAAC,CAAC;YACZ,IAAI,KAAK,GAAG,EAAE,EAAE;gBACd,OAAO,IAAI,MAAM,CAAC;aACnB;iBAAM;gBACL,OAAO,IAAI,MAAM,CAAC;aACnB;SACF;QAED,OAAO,IAAI,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,OAAO,CAAC,CAAC,CAAC;IACvC,CAAC;IAED;;;;;;OAMG;IACY,kBAAU,GAAzB,UAA0B,CAAC,EAAE,GAAG,EAAE,GAAG;QACnC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;IACxG,CAAC;IAED;;;;;;;;;;;OAWG;IACI,oBAAY,GAAnB,UAAoB,IAAU,EAAE,SAAiB,EAAE,GAAW;QAC5D,4BAA4B;QAC5B,IAAM,KAAK,GAAG,IAAI,IAAI,EAAE,CAAC;QAEzB,KAAK,CAAC,cAAc,CAAC,IAAI,CAAC,cAAc,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAClD,KAAK,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QAC9B,IAAM,IAAI,GAAG,IAAI,CAAC,OAAO,EAAE,GAAG,KAAK,CAAC,OAAO,EAAE,CAAC;QAC9C,IAAM,SAAS,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,GAAG,UAAU,CAAC,CAAC;QAEhD,IAAM,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,SAAS,GAAG,EAAE,CAAC,GAAG,OAAO,CAAC;QACnD,IAAM,cAAc,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACvF,IAAM,sBAAsB,GAAG,EAAE,GAAG,SAAS,CAAC;QAC9C,IAAM,cAAc,GAAG,cAAc,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,sBAAsB,CAAC,CAAC;QAE3E,OAAO,IAAI,IAAI,CAAC,IAAI,CAAC,OAAO,EAAE,GAAG,cAAc,GAAG,EAAE,GAAG,IAAI,CAAC,CAAC;IAC/D,CAAC;IAEM,mBAAW,GAAlB,UAAmB,IAAU,EAAE,GAAY,EAAE,GAAY,EAAE,EAAU,EAAE,GAAW;QAChF,IAAM,EAAE,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC;QAC1B,IAAM,GAAG,GAAG,GAAG,GAAG,OAAO,CAAC;QAC1B,IAAM,CAAC,GAAG,OAAO,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC;QACxC,IAAM,CAAC,GAAG,OAAO,CAAC,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QAC5D,IAAI,CAAC,GAAG,OAAO,CAAC,SAAS,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QAEzD,CAAC,IAAI,OAAO,CAAC,eAAe,CAAC,CAAC,CAAC,CAAC,CAAC,sCAAsC;QAEvE,OAAO;YACL,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;YAClD,EAAE,EAAE,CAAC;SACN,CAAC;IACJ,CAAC;IAEM,kBAAU,GAAjB,UAAkB,IAAU,EAAE,GAAY,EAAE,GAAY;QACtD,IAAM,EAAE,GAAG,CAAC,GAAG,GAAG,OAAO,CAAC;QAC1B,IAAM,GAAG,GAAG,GAAG,GAAG,OAAO,CAAC;QAC1B,IAAM,CAAC,GAAG,OAAO,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC;QACxC,IAAM,CAAC,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAC;QACjC,IAAM,CAAC,GAAG,OAAO,CAAC,YAAY,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QAE7C,OAAO;YACL,EAAE,EAAE,OAAO,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC;YAClC,EAAE,EAAE,OAAO,CAAC,SAAS,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC;SACrC,CAAC;IACJ,CAAC;IAEM,mBAAW,GAAlB,UAAmB,CAAS;QAC1B,IAAM,CAAC,GAAG,OAAO,CAAC,gBAAgB,CAAC,CAAC,CAAC,CAAC;QACtC,IAAM,CAAC,GAAG,OAAO,CAAC,iBAAiB,CAAC,CAAC,CAAC,CAAC;QAEvC,OAAO;YACL,GAAG,EAAE,OAAO,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9B,EAAE,EAAE,OAAO,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,CAAC;SACjC,CAAC;IACJ,CAAC;IArbuB,aAAK,GAAG,SAAS,CAAC;IAClB,aAAK,GAAG,OAAO,CAAC;IAChB,UAAE,GAAG,MAAM,CAAC;IACZ,SAAC,GAAG,OAAO,GAAG,OAAO,CAAC;IAE/B,aAAK,GAAG;QACrB,CAAC,CAAC,EAAE,mBAAmB,EAAE,qBAAqB,CAAC;QAC/C,CAAC,CAAC,GAAG,EAAE,YAAY,EAAE,aAAa,CAAC;QACnC,CAAC,CAAC,KAAK,EAAE,cAAc,EAAE,WAAW,CAAC;QACrC,CAAC,CAAC,CAAC,EAAE,qBAAqB,EAAE,mBAAmB,CAAC;QAChD,CAAC,CAAC,CAAC,EAAE,iBAAiB,EAAE,mBAAmB,CAAC;QAC5C,CAAC,CAAC,CAAC,EAAE,WAAW,EAAE,WAAW,CAAC;QAC9B,CAAC,CAAC,CAAC,EAAE,mBAAmB,EAAE,iBAAiB,CAAC;QAC5C,CAAC,CAAC,EAAE,EAAE,cAAc,EAAE,cAAc,CAAC;QACrC,CAAC,CAAC,EAAE,EAAE,aAAa,EAAE,aAAa,CAAC;QACnC,CAAC,CAAC,EAAE,EAAE,kBAAkB,EAAE,kBAAkB,CAAC,EAAE,aAAa;KAC7D,CAAC;IAsaJ,cAAC;CAAA,AAvbD,IAubC;SAvbY,OAAO"}