@openglobus/og 0.17.2 → 0.18.1

package/src/og/astro/jd.js

@@ -1,514 +0,0 @@
-"use strict";
-
-import { binarySearch } from "../utils/shared.js";
-
-/**
- * Seconds in millisecond.
- * @const
- * @default
- */
-export const SECONDS_PER_MILLISECOND = 0.001;
-
-/**
- * Milliseconds in second.
- * @const
- * @default
- */
-export const MILLISECONDS_PER_SECOND = 1000.0;
-
-/**
- * Seconds in minute.
- * @const
- * @default
- */
-export const SECONDS_PER_MINUTE = 60.0;
-
-/**
- * One by seconds in minute.
- * @const
- * @default
- */
-export const ONE_BY_SECONDS_PER_MINUTE = 1.0 / SECONDS_PER_MINUTE;
-
-/**
- * Minutes in hour.
- * @const
- * @default
- */
-export const MINUTES_PER_HOUR = 60.0;
-
-/**
- * Hours in day.
- * @const
- * @default
- */
-export const HOURS_PER_DAY = 24.0;
-
-/**
- * One by hours in day.
- * @const
- * @default
- */
-export const ONE_BY_HOURS_PER_DAY = 1.0 / HOURS_PER_DAY;
-
-/**
- * Seconds in hour.
- * @const
- * @default
- */
-export const SECONDS_PER_HOUR = 3600.0;
-
-/**
- * One by seconds in hour.
- * @const
- * @default
- */
-export const ONE_BY_SECONDS_PER_HOUR = 1.0 / SECONDS_PER_HOUR;
-
-/**
- * Seconds in 12 hours.
- * @const
- * @default
- */
-export const SECONDS_PER_12_HOURS = 12.0 * SECONDS_PER_HOUR;
-
-/**
- * Minutes in day.
- * @const
- * @default
- */
-export const MINUTES_PER_DAY = 1440.0;
-
-/**
- * One by minutes in day.
- * @const
- * @default
- */
-export const ONE_BY_MINUTES_PER_DAY = 1.0 / MINUTES_PER_DAY;
-
-/**
- * Seconds in day.
- * @const
- * @default
- */
-export const SECONDS_PER_DAY = 86400.0;
-
-/**
- * Milliseconds in day.
- * @const
- * @default
- */
-export const MILLISECONDS_PER_DAY = 86400000.0;
-
-/**
- * One by milliseconds in day.
- * @const
- * @default
- */
-export const ONE_BY_MILLISECONDS_PER_DAY = 1.0 / MILLISECONDS_PER_DAY;
-
-/**
- * One by seconds in day.
- * @const
- * @default
- */
-export const ONE_BY_SECONDS_PER_DAY = 1.0 / SECONDS_PER_DAY;
-
-/**
- * Days in julian century.
- * @const
- * @default
- */
-export const DAYS_PER_JULIAN_CENTURY = 36525.0;
-
-/**
- * Days in julian year.
- * @const
- * @default
- */
-export const DAYS_PER_JULIAN_YEAR = 365.25;
-
-/**
- * Seconds in picosecond.
- * @const
- * @default
- */
-export const PICOSECOND = 0.000000001;
-
-/**
- * Modified julian date difference.
- * @const
- * @default
- */
-export const MODIFIED_JULIAN_DATE_DIFFERENCE = 2400000.5;
-
-/**
- * Julian date of 2000 year. Epoch.
- * @const
- * @default
- */
-export const J2000 = 2451545.0;
-
-/**
- * Returns julian days from Epoch.
- * @param {number} jd - Julian date.
- * @returns {number} Days from epoch
- */
-export function T(jd) {
-    return (jd - J2000) / DAYS_PER_JULIAN_CENTURY;
-}
-
-/**
- * Gets the date's julian day.
- * @param {number} year - Year.
- * @param {number} month - Month.
- * @param {number} day - Day.
- * @returns {Number} Day number
- */
-export function getDayNumber(year, month, day) {
-    var a = ((month - 14) / 12) | 0;
-    var b = year + 4800 + a;
-    return (
-        (((1461 * b) / 4) | 0) +
-        (((367 * (month - 2 - 12 * a)) / 12) | 0) -
-        (((3 * (((b + 100) / 100) | 0)) / 4) | 0) +
-        day -
-        32075
-    );
-}
-
-/**
- * Converts javascript date to the universal(UTC) julian date.
- * @param {Date} date - Date.
- * @returns {number} UTC julian date
- */
-export function DateToUTC(date) {
-    var dayNumber = getDayNumber(date.getUTCFullYear(), date.getUTCMonth() + 1, date.getUTCDate());
-
-    var hour = date.getUTCHours() - 12;
-    if (hour < 0) {
-        hour += 24;
-    }
-
-    var secondsOfDay =
-        date.getUTCSeconds() +
-        hour * SECONDS_PER_HOUR +
-        date.getUTCMinutes() * SECONDS_PER_MINUTE +
-        date.getUTCMilliseconds() * SECONDS_PER_MILLISECOND;
-
-    if (secondsOfDay >= SECONDS_PER_12_HOURS) {
-        dayNumber--;
-    }
-
-    var extraDays = (secondsOfDay * ONE_BY_SECONDS_PER_DAY) | 0;
-    dayNumber += extraDays;
-    secondsOfDay -= SECONDS_PER_DAY * extraDays;
-
-    if (secondsOfDay < 0) {
-        dayNumber--;
-        secondsOfDay += SECONDS_PER_DAY;
-    }
-
-    return dayNumber + secondsOfDay * ONE_BY_SECONDS_PER_DAY;
-}
-
-/**
- * Converts javascript date to the atomic(TAI) julian date.
- * @param {Date} date - Date.
- * @returns {number} TAI julian date
- */
-export function DateToTAI(date) {
-    return UTCtoTAI(DateToUTC(date));
-}
-
-/**
- * Converts coordinated universal(UTC) julian date to atomic(TAI) julian date.
- * @param {number} jd - UTC julian date.
- * @returns {number} TAI julian date
- */
-export function UTCtoTAI(jd) {
-    var leapSeconds = leapSecondsTable;
-
-    var index = binarySearch(leapSeconds, jd, function (a, b) {
-        return a - b.jd;
-    });
-
-    if (index < 0) {
-        index = ~index;
-    }
-
-    if (index >= leapSeconds.length) {
-        index = leapSeconds.length - 1;
-    }
-
-    var offset = leapSeconds[index].leapSeconds;
-
-    if (index !== 0) {
-        if ((leapSeconds[index].jd - jd) * SECONDS_PER_DAY > offset) {
-            offset = leapSeconds[index - 1].leapSeconds;
-        }
-    }
-
-    return jd + offset * ONE_BY_SECONDS_PER_DAY;
-}
-
-/**
- * Converts atomic julian date(TAI) to the coordinated universal(UTC) julian date.
- * @param {number} tai - TAI julian date.
- * @returns {number} UTC julian date
- */
-export function TAItoUTC(tai) {
-    var leapSeconds = leapSecondsTable;
-    var index = binarySearch(leapSeconds, tai, function (a, b) {
-        return a - b.jd;
-    });
-
-    if (index < 0) {
-        index = ~index;
-    }
-
-    if (index >= leapSeconds.length) {
-        return tai - leapSeconds[index - 1].leapSeconds * ONE_BY_SECONDS_PER_DAY;
-    }
-
-    if (index === 0) {
-        return tai - leapSeconds[0].leapSeconds * ONE_BY_SECONDS_PER_DAY;
-    }
-
-    var diff = (leapSeconds[index].jd - tai) * SECONDS_PER_DAY;
-
-    if (diff === 0) {
-        return tai - leapSeconds[index].leapSeconds * ONE_BY_SECONDS_PER_DAY;
-    }
-
-    if (diff <= 1.0) {
-        return null;
-    }
-
-    return tai - leapSeconds[index - 1].leapSeconds * ONE_BY_SECONDS_PER_DAY;
-}
-
-/**
- * Converts UTC julian date to the javascript date object.
- * @param {number} utc - UTC julian date.
- * @returns {Date} JavaScript Date object
- */
-export function UTCtoDate(utc) {
-    var julianDayNumber = utc | 0;
-    var secondsOfDay = (utc - julianDayNumber) * SECONDS_PER_DAY;
-
-    if (secondsOfDay >= SECONDS_PER_12_HOURS) {
-        julianDayNumber++;
-    }
-
-    var L = (julianDayNumber + 68569) | 0;
-    var N = ((4 * L) / 146097) | 0;
-    L = (L - (((146097 * N + 3) / 4) | 0)) | 0;
-    var I = ((4000 * (L + 1)) / 1461001) | 0;
-    L = (L - (((1461 * I) / 4) | 0) + 31) | 0;
-    var J = ((80 * L) / 2447) | 0;
-    var day = (L - (((2447 * J) / 80) | 0)) | 0;
-    L = (J / 11) | 0;
-    var month = (J + 2 - 12 * L) | 0;
-    var year = (100 * (N - 49) + I + L) | 0;
-
-    var hour = (secondsOfDay * ONE_BY_SECONDS_PER_HOUR) | 0;
-    var remainingSeconds = secondsOfDay - hour * SECONDS_PER_HOUR;
-    var minute = (remainingSeconds * ONE_BY_SECONDS_PER_MINUTE) | 0;
-    remainingSeconds = remainingSeconds - minute * SECONDS_PER_MINUTE;
-    var second = remainingSeconds | 0;
-    var millisecond = ((remainingSeconds - second) * MILLISECONDS_PER_SECOND) | 0;
-
-    hour += 12;
-    if (hour > 23) {
-        hour -= 24;
-    }
-
-    return new Date(Date.UTC(year, month - 1, day, hour, minute, second, millisecond));
-}
-
-/**
- * Converts TAI julian date to the javascript date object.
- * @param {number} tai - TAI julian date.
- * @returns {Date} JavaScript Date object
- */
-export function TAItoDate(tai) {
-    var utc = TAItoUTC(tai);
-    if (!utc) {
-        utc = TAItoUTC(addSeconds(tai, -1));
-        console.warn("TAItoDate:336 - can't conv utc.");
-    }
-
-    return UTCtoDate(utc);
-}
-
-/**
- * Adds milliseconds to the julian date.
- * @param {number} jd - Julian date.
- * @param {number} milliseconds - Milliseconds to add.
- * @returns {number} Julian date
- */
-export function addMilliseconds(jd, milliseconds) {
-    return jd + milliseconds * ONE_BY_MILLISECONDS_PER_DAY;
-}
-
-/**
- * Adds seconds to the julian date.
- * @param {number} jd - Julian date.
- * @param {number} seconds - Seconds to add.
- * @returns {number} Julian date
- */
-export function addSeconds(jd, seconds) {
-    return jd + seconds * ONE_BY_SECONDS_PER_DAY;
-}
-
-/**
- * Adds hours to the julian date.
- * @param {number} jd - Julian date.
- * @param {number} hours - Hours to add.
- * @returns {number} Julian date
- */
-export function addHours(jd, hours) {
-    return jd + hours * ONE_BY_HOURS_PER_DAY;
-}
-
-/**
- * Adds minutes to the julian date.
- * @param {number} jd - Julian date.
- * @param {number} minutes - Minutes to add.
- * @returns {number} Julian date
- */
-export function addMinutes(jd, minutes) {
-    return jd + minutes * MINUTES_PER_DAY;
-}
-
-/**
- * Adds days to the julian date.
- * @param {number} jd - Julian date.
- * @param {number} days - Days to add.
- * @returns {number} Julian date
- */
-export function addDays(jd, days) {
-    return jd + days;
-}
-
-/**
- * Gets milliseconds of a julian date.
- * @param {number} js - julian date.
- * @returns {number} Milliseconds
- */
-export function getMilliseconds(jd) {
-    var s = jd - (jd | 0);
-    s *= SECONDS_PER_DAY;
-    return ((s - (s | 0)) * MILLISECONDS_PER_SECOND) | 0;
-}
-
-/**
- * Gets seconds of a julian date.
- * @param {number} js - julian date.
- * @returns {number} Seconds
- */
-export function getSeconds(jd) {
-    var s = jd - (jd | 0);
-    return s * SECONDS_PER_DAY;
-}
-
-/**
- * Gets hours of a julian date.
- * @param {number} js - julian date.
- * @returns {number} Hours
- */
-export function getHours(jd) {
-    var julianDayNumber = jd | 0;
-    var secondsOfDay = (jd - julianDayNumber) * SECONDS_PER_DAY;
-
-    var hour = (secondsOfDay * ONE_BY_SECONDS_PER_HOUR) | 0;
-    var remainingSeconds = secondsOfDay - hour * SECONDS_PER_HOUR;
-    var minute = (remainingSeconds * ONE_BY_SECONDS_PER_MINUTE) | 0;
-    remainingSeconds = remainingSeconds - minute * SECONDS_PER_MINUTE;
-    var second = remainingSeconds | 0;
-    var millisecond = ((remainingSeconds - second) * MILLISECONDS_PER_SECOND) | 0;
-
-    hour += 12 + minute / 60 + second / 3600 + millisecond / 1000;
-    if (hour > 23) {
-        hour -= 24;
-    }
-
-    return hour;
-}
-
-/**
- * Gets minutes of a julian date.
- * @param {number} js - julian date.
- * @returns {number} Minutes
- */
-export function getMinutes(jd) {
-    var s = jd - (jd | 0);
-    return (s * MINUTES_PER_DAY) | 0;
-}
-
-/**
- * Gets days of a julian date.
- * @param {number} js - julian date.
- * @returns {number} Days
- */
-export function getDays(jd) {
-    return jd | 0;
-}
-
-/**
- * Returns days in seconds.
- * @param {number} s - Seconds.
- * @returns {number} Days
- */
-export function secondsToDays(s) {
-    return s * ONE_BY_SECONDS_PER_DAY;
-}
-
-/**
- * Returns seconds in days.
- * @param {number} d - Days.
- * @returns {number} Seconds
- */
-export function daysToSeconds(d) {
-    return d * SECONDS_PER_DAY;
-}
-
-function __ls(jd, leapSeconds) {
-    return {
-        jd: jd,
-        leapSeconds: leapSeconds
-    };
-}
-
-const leapSecondsTable = [
-    __ls(2441317.5, 10.0), // 1972-01-01T00:00:00.000Z
-    __ls(2441499.5, 11.0), // 1972-07-01T00:00:00.000Z
-    __ls(2441683.5, 12.0), // 1973-01-01T00:00:00.000Z
-    __ls(2442048.5, 13.0), // 1974-01-01T00:00:00.000Z
-    __ls(2442413.5, 14.0), // 1975-01-01T00:00:00.000Z
-    __ls(2442778.5, 15.0), // 1976-01-01T00:00:00.000Z
-    __ls(2443144.5, 16.0), // 1977-01-01T00:00:00.000Z
-    __ls(2443509.5, 17.0), // 1978-01-01T00:00:00.000Z
-    __ls(2443874.5, 18.0), // 1979-01-01T00:00:00.000Z
-    __ls(2444239.5, 19.0), // 1980-01-01T00:00:00.000Z
-    __ls(2444786.5, 20.0), // 1981-07-01T00:00:00.000Z
-    __ls(2445151.5, 21.0), // 1982-07-01T00:00:00.000Z
-    __ls(2445516.5, 22.0), // 1983-01-01T00:00:00.000Z
-    __ls(2446247.5, 23.0), // 1985-07-01T00:00:00.000Z
-    __ls(2447161.5, 24.0), // 1988-01-01T00:00:00.000Z
-    __ls(2447892.5, 25.0), // 1990-01-01T00:00:00.000Z
-    __ls(2448257.5, 26.0), // 1991-01-01T00:00:00.000Z
-    __ls(2448804.5, 27.0), // 1992-07-01T00:00:00.000Z
-    __ls(2449169.5, 28.0), // 1993-07-01T00:00:00.000Z
-    __ls(2449534.5, 29.0), // 1994-07-01T00:00:00.000Z
-    __ls(2450083.5, 30.0), // 1996-01-01T00:00:00.000Z
-    __ls(2450630.5, 31.0), // 1997-07-01T00:00:00.000Z
-    __ls(2451179.5, 32.0), // 1999-01-01T00:00:00.000Z
-    __ls(2453736.5, 33.0), // 2006-01-01T00:00:00.000Z
-    __ls(2454832.5, 34.0), // 2009-01-01T00:00:00.000Z
-    __ls(2456109.5, 35.0), // 2012-07-01T00:00:00.000Z
-    __ls(2457204.5, 36.0) // 2015-07-01T00:00:00.000Z
-];
-
-export const J2000TAI = UTCtoTAI(J2000);

package/src/og/astro/orbit.js

@@ -1,119 +0,0 @@
-"use strict";
-
-import * as math from "../math.js";
-import { Mat3 } from "../math/Mat3.js";
-
-export function getEccentricAnomaly(M, ecc) {
-    if (ecc == 0.0) {
-        // Circular orbit
-        return M;
-    } else if (ecc < 0.2) {
-        // Low eccentricity, so use the standard iteration technique
-        return math.solve_iteration_fixed(solveKeplerFunc1(ecc, M), M, 5);
-    } else if (ecc < 0.9) {
-        // Higher eccentricity elliptical orbit; use a more complex but
-        // much faster converging iteration.
-        return math.solve_iteration_fixed(solveKeplerFunc2(ecc, M), M, 6);
-    } else if (ecc < 1.0) {
-        // Extremely stable Laguerre-Conway method for solving Kepler's
-        // equation.  Only use this for high-eccentricity orbits, as it
-        // requires more calcuation.
-        let E = M + 0.85 * ecc * Math.sign(Math.sin(M));
-        return math.solve_iteration_fixed(solveKeplerLaguerreConway(ecc, M), E, 8);
-    } else if (ecc == 1.0) {
-        // TODO: Parabolic orbit
-        return M;
-    } else {
-        // Laguerre-Conway method for hyperbolic (ecc > 1) orbits.
-        let E = Math.log((2 * M) / ecc + 1.85);
-        return math.solve_iteration_fixed(solveKeplerLaguerreConwayHyp(ecc, M), E, 30);
-    }
-}
-
-// Standard iteration for solving Kepler's Equation
-function solveKeplerFunc1(ecc, M) {
-    return function (x) {
-        return M + ecc * Math.sin(x);
-    };
-}
-
-// Faster converging iteration for Kepler's Equation; more efficient
-// than above for orbits with eccentricities greater than 0.3.  This
-// is from Jean Meeus's _Astronomical Algorithms_ (2nd ed), p. 199
-function solveKeplerFunc2(ecc, M) {
-    return function (x) {
-        return x + (M + ecc * Math.sin(x) - x) / (1 - ecc * Math.cos(x));
-    };
-}
-
-function solveKeplerLaguerreConway(ecc, M) {
-    return function (x) {
-        var s = ecc * Math.sin(x);
-        var c = ecc * Math.cos(x);
-        var f = x - s - M;
-        var f1 = 1 - c;
-        var f2 = s;
-        x += (-5 * f) / (f1 + Math.sign(f1) * Math.sqrt(Math.abs(16 * f1 * f1 - 20 * f * f2)));
-        return x;
-    };
-}
-
-function solveKeplerLaguerreConwayHyp(ecc, M) {
-    return function (x) {
-        var s = ecc * Math.sinh(x);
-        var c = ecc * Math.cosh(x);
-        var f = s - x - M;
-        var f1 = c - 1;
-        var f2 = s;
-        x += (-5 * f) / (f1 + Math.sign(f1) * Math.sqrt(Math.abs(16 * f1 * f1 - 20 * f * f2)));
-        return x;
-    };
-}
-
-export function getEllipticalEccentricAnomaly(meanAnomaly, eccentricity) {
-    var tol = 0.00000001745;
-    var iterations = 20;
-    var e = meanAnomaly - 2.0 * Math.PI * ((meanAnomaly / (2.0 * Math.PI)) | 0);
-    var err = 1;
-    while (Math.abs(err) > tol && iterations > 0) {
-        err = e - eccentricity * Math.sin(e) - meanAnomaly;
-        var delta = err / (1 - eccentricity * Math.cos(e));
-        e -= delta;
-        iterations--;
-    }
-    return e;
-}
-
-export function getTrueAnomaly(eccentricAnomaly, eccentricity) {
-    var revs = Math.floor(eccentricAnomaly / math.TWO_PI);
-    eccentricAnomaly -= revs * math.TWO_PI;
-    var trueAnomaly = Math.atan2(
-        Math.sin(eccentricAnomaly) * Math.sqrt(1 - eccentricity * eccentricity),
-        Math.cos(eccentricAnomaly) - eccentricity
-    );
-    trueAnomaly = math.zeroTwoPI(trueAnomaly);
-    if (eccentricAnomaly < 0) {
-        trueAnomaly -= math.TWO_PI;
-    }
-    return trueAnomaly + revs * math.TWO_PI;
-}
-
-export function getPerifocalToCartesianMatrix(argumentOfPeriapsis, inclination, rightAscension) {
-    var res = new Mat3();
-    var cosap = Math.cos(argumentOfPeriapsis),
-        sinap = Math.sin(argumentOfPeriapsis),
-        cosi = Math.cos(inclination),
-        sini = Math.sin(inclination),
-        cosraan = Math.cos(rightAscension),
-        sinraan = Math.sin(rightAscension);
-    res._m[0] = cosraan * cosap - sinraan * sinap * cosi;
-    res._m[1] = sinraan * cosap + cosraan * sinap * cosi;
-    res._m[2] = sinap * sini;
-    res._m[3] = -cosraan * sinap - sinraan * cosap * cosi;
-    res._m[4] = -sinraan * sinap + cosraan * cosap * cosi;
-    res._m[5] = cosap * sini;
-    res._m[6] = sinraan * sini;
-    res._m[7] = -cosraan * sini;
-    res._m[8] = cosi;
-    return res;
-}

package/src/og/astro/rotation.js

@@ -1,36 +0,0 @@
-"use strict";
-
-import { PI_TWO } from "../math.js";
-import { Mat3 } from "../math/Mat3.js";
-import { Vec3 } from "../math/Vec3.js";
-
-export function getRotationMatrix(rightAscension, declination, res) {
-    let xAxis = new Vec3(),
-        zAxis = new Vec3();
-
-    res = res || new Mat3();
-
-    xAxis.x = Math.cos(rightAscension + PI_TWO);
-    xAxis.y = Math.sin(rightAscension + PI_TWO);
-    xAxis.z = 0.0;
-
-    let cosDec = Math.cos(declination);
-
-    zAxis.x = cosDec * Math.cos(rightAscension);
-    zAxis.y = cosDec * Math.sin(rightAscension);
-    zAxis.z = Math.sin(declination);
-
-    let yAxis = zAxis.cross(xAxis);
-
-    res._m[0] = xAxis.x;
-    res._m[1] = yAxis.x;
-    res._m[2] = zAxis.x;
-    res._m[3] = xAxis.y;
-    res._m[4] = yAxis.y;
-    res._m[5] = zAxis.y;
-    res._m[6] = xAxis.z;
-    res._m[7] = yAxis.z;
-    res._m[8] = zAxis.z;
-
-    return res;
-}