@openglobus/og 0.11.1 → 0.11.6

package/src/og/ajax.js

@@ -2,7 +2,7 @@
  * @module og/ajax
  */
 
-'use strict';
+"use strict";
 
 /**
  * Ajax parameters.
@@ -13,8 +13,7 @@ const ajax = {
      * Ajax ready state result.
      * @enum
      */
-    ReadyState:
-    {
+    ReadyState: {
         Uninitialized: 0,
         Loading: 1,
         Loaded: 2,
@@ -25,8 +24,7 @@ const ajax = {
      * Ajax status code.
      * @enum
      */
-    Status:
-    {
+    Status: {
         OK: 200,
         Created: 201,
         Accepted: 202,
@@ -88,24 +86,6 @@ const defaultParams = {
     responseType: "text"
 };
 
-function createXMLHttp() {
-    var xhr = null;
-    if (typeof XMLHttpRequest != "undefined") {
-        xhr = new XMLHttpRequest();
-        return xhr;
-    } else if (window.ActiveXObject) {
-        var ieXMLHttpVersions = ['MSXML2.XMLHttp.5.0', 'MSXML2.XMLHttp.4.0', 'MSXML2.XMLHttp.3.0', 'MSXML2.XMLHttp', 'Microsoft.XMLHttp'];
-        for (var i = 0; i < ieXMLHttpVersions.length; i++) {
-            try {
-                xhr = new ActiveXObject(ieXMLHttpVersions[i]);
-                return xhr;
-            } catch (e) {
-                throw new Error('og.ajax.createXMLHttp creation failed.');
-            }
-        }
-    }
-};
-
 /**
  * Send an ajax request.
  * @function
@@ -122,10 +102,10 @@ function createXMLHttp() {
  * @returns {ajax.Xhr} - Returns object that could be aborted.
  */
 ajax.request = function (url, params) {
-
     params = params || {};
 
-    var p = {}, i;
+    var p = {},
+        i;
 
     for (i in defaultParams) {
         p[i] = defaultParams[i];
@@ -137,18 +117,23 @@ ajax.request = function (url, params) {
 
     p.data = params.data;
 
-    var xhr = createXMLHttp();
+    var xhr = new XMLHttpRequest();
 
     var customXhr = new Xhr(xhr);
 
-    var body = null, d;
+    var body = null,
+        d;
 
     if (p.type === ajax.Method.Post) {
         if (p.data) {
             body = "";
             for (let key in p.data) {
                 d = p.data[key];
-                body += key + "=" + encodeURIComponent(d instanceof Object ? JSON.stringify(d) : d) + "&";
+                body +=
+                    key +
+                    "=" +
+                    encodeURIComponent(d instanceof Object ? JSON.stringify(d) : d) +
+                    "&";
             }
             body = body.slice(0, -1);
         }
@@ -158,7 +143,8 @@ ajax.request = function (url, params) {
         var tail = "?";
         for (let key in p.data) {
             d = p.data[key];
-            tail += key + "=" + encodeURIComponent(d instanceof Object ? JSON.stringify(d) : d) + "&";
+            tail +=
+                key + "=" + encodeURIComponent(d instanceof Object ? JSON.stringify(d) : d) + "&";
         }
         tail = tail.slice(0, -1);
         xhr.open(p.type, url + tail, p.async);
@@ -190,7 +176,8 @@ ajax.request = function (url, params) {
                  * @param {Object} Response data
                  * @param {Object} Status object
                  */
-                params.abort && params.abort.call(params.sender || customXhr, xhr.response, xhr.status);
+                params.abort &&
+                    params.abort.call(params.sender || customXhr, xhr.response, xhr.status);
             } else {
                 /**
                  * Error callback.
@@ -198,9 +185,10 @@ ajax.request = function (url, params) {
                  * @param {Object} Response data
                  * @param {Object} Status object
                  */
-                params.error && params.error.call(params.sender || customXhr, xhr.response, xhr.status);
+                params.error &&
+                    params.error.call(params.sender || customXhr, xhr.response, xhr.status);
             }
-            delete xhr['onreadystatechange'];
+            delete xhr["onreadystatechange"];
             xhr.onreadystatechange = null;
             xhr = null;
         } else {
@@ -213,4 +201,4 @@ ajax.request = function (url, params) {
     return customXhr;
 };
 
-export { ajax };
+export { ajax };

package/src/og/astro/astro.js

@@ -1,6 +1,6 @@
-'use strict';
+"use strict";
 
-import * as jd from './jd.js';
+import * as jd from "./jd.js";
 
 /**
  * Angle between J2000 mean equator and the ecliptic plane.
@@ -16,7 +16,7 @@ export const J2000_OBLIQUITY = 23.4392911;
  * @const
  * @type{Number}
  */
-export const AU_TO_METERS = 1.49597870e+11;
+export const AU_TO_METERS = 1.4959787e11;
 
 /**
  * Terestrial and atomic time difference.
@@ -45,7 +45,7 @@ export const SUN_GRAVITATIONAL_PARAMETER = 1.32712440018e20;
  * @returns {Number} - returns barycentric dynamical time.
  */
 export function TAItoTDB(tai) {
-    tai += og.astro.TDT_TAI * jd.ONE_BY_SECONDS_PER_DAY;
+    tai += TDT_TAI * jd.ONE_BY_SECONDS_PER_DAY;
     var g = 6.239996 + 0.0172019696544 * (tai - jd.J2000);
     return tai + 0.001658 * Math.sin(g + 1.671e-2 * Math.sin(g)) * jd.ONE_BY_SECONDS_PER_DAY;
-};
+}

package/src/og/astro/earth.js

@@ -2,18 +2,18 @@
  * @module og/astro/earth
  */
 
-'use strict';
+"use strict";
 
-import * as jd from './jd.js';
-import * as math from '../math.js';
-import * as astro from './astro.js';
-import { Quat } from '../math/Quat.js';
-import { Vec3 } from '../math/Vec3.js';
+import * as jd from "./jd.js";
+import * as math from "../math.js";
+import * as astro from "./astro.js";
+import { Quat } from "../math/Quat.js";
+import { Vec3 } from "../math/Vec3.js";
 
 /**
  * Returns Sun position in the geocentric coordinate system by the time.
  * @param {Number} jDate - Julian date time.
- * @returns {og.Vec3} - Sun geocentric coordinates.
+ * @returns {Vec3} - Sun geocentric coordinates.
  */
 export function getSunPosition(jDate) {
     // http://stjarnhimlen.se/comp/tutorial.html
@@ -21,7 +21,7 @@ export function getSunPosition(jDate) {
     // e  Eccentricity
     // T  Time at perihelion
 
-    // q  Perihelion distance  = a * (1 - e)    
+    // q  Perihelion distance  = a * (1 - e)
     // Q  Aphelion distance    = a * (1 + e)
 
     // i  Inclination, i.e. the "tilt" of the orbit relative to the
@@ -70,16 +70,17 @@ export function getSunPosition(jDate) {
 
     var d = jDate - jd.J2000;
 
-    var w = 282.9404 + 4.70935E-5 * d; // longitude of perihelion
+    var w = 282.9404 + 4.70935e-5 * d; // longitude of perihelion
     // var a = 1.000000; // mean distance, a.u.
-    var e = 0.016709 - 1.151E-9 * d; // eccentricity
-    var M = math.rev(356.0470 + 0.9856002585 * d); // mean anomaly
+    var e = 0.016709 - 1.151e-9 * d; // eccentricity
+    var M = math.rev(356.047 + 0.9856002585 * d); // mean anomaly
 
-    var oblecl = astro.J2000_OBLIQUITY - 3.563E-7 * d; // obliquity of the ecliptic
+    var oblecl = astro.J2000_OBLIQUITY - 3.563e-7 * d; // obliquity of the ecliptic
 
     // var L = math.rev(w + M); // Sun's mean longitude
 
-    var E = M + math.DEGREES * e * Math.sin(M * math.RADIANS) * (1 + e * Math.cos(M * math.RADIANS)); // eccentric anomaly
+    var E =
+        M + math.DEGREES * e * Math.sin(M * math.RADIANS) * (1 + e * Math.cos(M * math.RADIANS)); // eccentric anomaly
 
     // Sun rectangular coordiantes, where the X axis points towards the perihelion
     var x = Math.cos(E * math.RADIANS) - e;
@@ -99,12 +100,17 @@ export function getSunPosition(jDate) {
     var yequat = y * Math.cos(oblecl * math.RADIANS);
     var zequat = y * Math.sin(oblecl * math.RADIANS);
 
-    var theta = math.TWO_PI * (d * 24.0 / 23.9344694 - 259.853 / 360.0); // Siderial spin time
+    var theta = math.TWO_PI * ((d * 24.0) / 23.9344694 - 259.853 / 360.0); // Siderial spin time
 
-    return Quat.yRotation(-theta).mulVec3(new Vec3(-yequat * astro.AU_TO_METERS,
-        zequat * astro.AU_TO_METERS, -xequat * astro.AU_TO_METERS));
+    return Quat.yRotation(-theta).mulVec3(
+        new Vec3(
+            -yequat * astro.AU_TO_METERS,
+            zequat * astro.AU_TO_METERS,
+            -xequat * astro.AU_TO_METERS
+        )
+    );
 
     // Convert to RA and Decl
     // var RA = Math.atan2(yequat, xequat) * math.DEGREES;
     // var Decl = Math.atan2(zequat, Math.sqrt(xequat * xequat + yequat * yequat)) * math.DEGREES;
-};
+}

package/src/og/astro/jd.js

@@ -2,9 +2,9 @@
  * @module og/astro/jd
  */
 
-'use strict';
+"use strict";
 
-import { binarySearch } from '../utils/shared.js';
+import { binarySearch } from "../utils/shared.js";
 
 /**
  * Seconds in millisecond.
@@ -160,7 +160,7 @@ export const J2000 = 2451545.0;
  */
 export function T(jd) {
     return (jd - J2000) / DAYS_PER_JULIAN_CENTURY;
-};
+}
 
 /**
  * Gets the date's julian day.
@@ -172,9 +172,14 @@ export function T(jd) {
 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;
-};
+    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.
@@ -190,7 +195,8 @@ export function DateToUTC(date) {
     }
 
     var secondsOfDay =
-        date.getUTCSeconds() + hour * SECONDS_PER_HOUR +
+        date.getUTCSeconds() +
+        hour * SECONDS_PER_HOUR +
         date.getUTCMinutes() * SECONDS_PER_MINUTE +
         date.getUTCMilliseconds() * SECONDS_PER_MILLISECOND;
 
@@ -198,7 +204,7 @@ export function DateToUTC(date) {
         dayNumber--;
     }
 
-    var extraDays = secondsOfDay * ONE_BY_SECONDS_PER_DAY | 0;
+    var extraDays = (secondsOfDay * ONE_BY_SECONDS_PER_DAY) | 0;
     dayNumber += extraDays;
     secondsOfDay -= SECONDS_PER_DAY * extraDays;
 
@@ -208,7 +214,7 @@ export function DateToUTC(date) {
     }
 
     return dayNumber + secondsOfDay * ONE_BY_SECONDS_PER_DAY;
-};
+}
 
 /**
  * Converts javascript date to the atomic(TAI) julian date.
@@ -217,7 +223,7 @@ export function DateToUTC(date) {
  */
 export function DateToTAI(date) {
     return UTCtoTAI(DateToUTC(date));
-};
+}
 
 /**
  * Converts coordinated universal(UTC) julian date to atomic(TAI) julian date.
@@ -248,7 +254,7 @@ export function UTCtoTAI(jd) {
     }
 
     return jd + offset * ONE_BY_SECONDS_PER_DAY;
-};
+}
 
 /**
  * Converts atomic julian date(TAI) to the coordinated universal(UTC) julian date.
@@ -284,7 +290,7 @@ export function TAItoUTC(tai) {
     }
 
     return tai - leapSeconds[index - 1].leapSeconds * ONE_BY_SECONDS_PER_DAY;
-};
+}
 
 /**
  * Converts UTC julian date to the javascript date object.
@@ -300,7 +306,7 @@ export function UTCtoDate(utc) {
     }
 
     var L = (julianDayNumber + 68569) | 0;
-    var N = (4 * L / 146097) | 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;
@@ -310,12 +316,12 @@ export function UTCtoDate(utc) {
     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 hour = (secondsOfDay * ONE_BY_SECONDS_PER_HOUR) | 0;
     var remainingSeconds = secondsOfDay - hour * SECONDS_PER_HOUR;
-    var minute = remainingSeconds * ONE_BY_SECONDS_PER_MINUTE | 0;
+    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;
+    var millisecond = ((remainingSeconds - second) * MILLISECONDS_PER_SECOND) | 0;
 
     hour += 12;
     if (hour > 23) {
@@ -323,7 +329,7 @@ export function UTCtoDate(utc) {
     }
 
     return new Date(Date.UTC(year, month - 1, day, hour, minute, second, millisecond));
-};
+}
 
 /**
  * Converts TAI julian date to the javascript date object.
@@ -331,15 +337,14 @@ export function UTCtoDate(utc) {
  * @returns {Date} JavaScript Date object
  */
 export function TAItoDate(tai) {
-
     var utc = TAItoUTC(tai);
     if (!utc) {
         utc = TAItoUTC(addSeconds(tai, -1));
-        og.console.logWrn("TAItoDate:336 - can't conv utc.");
+        console.warn("TAItoDate:336 - can't conv utc.");
     }
 
     return UTCtoDate(utc);
-};
+}
 
 /**
  * Adds milliseconds to the julian date.
@@ -349,7 +354,7 @@ export function TAItoDate(tai) {
  */
 export function addMilliseconds(jd, milliseconds) {
     return jd + milliseconds * ONE_BY_MILLISECONDS_PER_DAY;
-};
+}
 
 /**
  * Adds seconds to the julian date.
@@ -359,7 +364,7 @@ export function addMilliseconds(jd, milliseconds) {
  */
 export function addSeconds(jd, seconds) {
     return jd + seconds * ONE_BY_SECONDS_PER_DAY;
-};
+}
 
 /**
  * Adds hours to the julian date.
@@ -369,7 +374,7 @@ export function addSeconds(jd, seconds) {
  */
 export function addHours(jd, hours) {
     return jd + hours * ONE_BY_HOURS_PER_DAY;
-};
+}
 
 /**
  * Adds minutes to the julian date.
@@ -379,7 +384,7 @@ export function addHours(jd, hours) {
  */
 export function addMinutes(jd, minutes) {
     return jd + minutes * MINUTES_PER_DAY;
-};
+}
 
 /**
  * Adds days to the julian date.
@@ -389,7 +394,7 @@ export function addMinutes(jd, minutes) {
  */
 export function addDays(jd, days) {
     return jd + days;
-};
+}
 
 /**
  * Gets milliseconds of a julian date.
@@ -399,8 +404,8 @@ export function addDays(jd, days) {
 export function getMilliseconds(jd) {
     var s = jd - (jd | 0);
     s *= SECONDS_PER_DAY;
-    return (s - (s | 0)) * MILLISECONDS_PER_SECOND | 0;
-};
+    return ((s - (s | 0)) * MILLISECONDS_PER_SECOND) | 0;
+}
 
 /**
  * Gets seconds of a julian date.
@@ -410,7 +415,7 @@ export function getMilliseconds(jd) {
 export function getSeconds(jd) {
     var s = jd - (jd | 0);
     return s * SECONDS_PER_DAY;
-};
+}
 
 /**
  * Gets hours of a julian date.
@@ -421,12 +426,12 @@ 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 hour = (secondsOfDay * ONE_BY_SECONDS_PER_HOUR) | 0;
     var remainingSeconds = secondsOfDay - hour * SECONDS_PER_HOUR;
-    var minute = remainingSeconds * ONE_BY_SECONDS_PER_MINUTE | 0;
+    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;
+    var millisecond = ((remainingSeconds - second) * MILLISECONDS_PER_SECOND) | 0;
 
     hour += 12 + minute / 60 + second / 3600 + millisecond / 1000;
     if (hour > 23) {
@@ -434,7 +439,7 @@ export function getHours(jd) {
     }
 
     return hour;
-};
+}
 
 /**
  * Gets minutes of a julian date.
@@ -443,8 +448,8 @@ export function getHours(jd) {
  */
 export function getMinutes(jd) {
     var s = jd - (jd | 0);
-    return s * MINUTES_PER_DAY | 0;
-};
+    return (s * MINUTES_PER_DAY) | 0;
+}
 
 /**
  * Gets days of a julian date.
@@ -453,7 +458,7 @@ export function getMinutes(jd) {
  */
 export function getDays(jd) {
     return jd | 0;
-};
+}
 
 /**
  * Returns days in seconds.
@@ -462,7 +467,7 @@ export function getDays(jd) {
  */
 export function secondsToDays(s) {
     return s * ONE_BY_SECONDS_PER_DAY;
-};
+}
 
 /**
  * Returns seconds in days.
@@ -471,14 +476,14 @@ export function secondsToDays(s) {
  */
 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
@@ -507,7 +512,7 @@ const leapSecondsTable = [
     __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
+    __ls(2457204.5, 36.0) // 2015-07-01T00:00:00.000Z
 ];
 
 export const J2000TAI = UTCtoTAI(J2000);

package/src/og/astro/orbit.js

@@ -2,10 +2,10 @@
  * @module og/astro/orbit
  */
 
-'use strict';
+"use strict";
 
-import * as math from '../math.js';
-import { Mat3 } from '../math/Mat3.js';
+import * as math from "../math.js";
+import { Mat3 } from "../math/Mat3.js";
 
 export function getEccentricAnomaly(M, ecc) {
     if (ecc == 0.0) {
@@ -22,24 +22,24 @@ export function getEccentricAnomaly(M, ecc) {
         // 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 * sign(sin(M));
+        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 = log(2 * M / ecc + 1.85);
+        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
@@ -48,7 +48,7 @@ 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) {
@@ -57,10 +57,10 @@ function solveKeplerLaguerreConway(ecc, M) {
         var f = x - s - M;
         var f1 = 1 - c;
         var f2 = s;
-        x += -5 * f / (f1 + Math.sign(f1) * Math.sqrt(abs(16 * f1 * f1 - 20 * f * f2)));
+        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) {
@@ -69,15 +69,15 @@ function solveKeplerLaguerreConwayHyp(ecc, M) {
         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)));
+        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 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;
@@ -86,19 +86,21 @@ export function getEllipticalEccentricAnomaly(meanAnomaly, eccentricity) {
         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);
+    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();
@@ -118,4 +120,4 @@ export function getPerifocalToCartesianMatrix(argumentOfPeriapsis, inclination,
     res._m[7] = -cosraan * sini;
     res._m[8] = cosi;
     return res;
-};
+}

package/src/og/astro/rotation.js

@@ -2,20 +2,29 @@
  * @module og/astro/rotation
  */
 
-'use strict';
+"use strict";
 
-export function getRotationMatrix(rightAscension, declination) {
-    xAxis.x = Math.cos(rightAscension + og.math.PI_TWO);
-    xAxis.y = Math.sin(rightAscension + og.math.PI_TWO);
+import { Vec3 } from "../math/Vec3.js";
+import { Mat3 } from "../math/Mat3.js";
+import { PI_TWO } from "../math.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;
 
-    var cosDec = Math.cos(declination);
+    let cosDec = Math.cos(declination);
 
     zAxis.x = cosDec * Math.cos(rightAscension);
     zAxis.y = cosDec * Math.sin(rightAscension);
     zAxis.z = Math.sin(declination);
 
-    var yAxis = zAxis.cross(xAxis);
+    let yAxis = zAxis.cross(xAxis);
 
     res._m[0] = xAxis.x;
     res._m[1] = yAxis.x;
@@ -27,5 +36,5 @@ export function getRotationMatrix(rightAscension, declination) {
     res._m[7] = yAxis.z;
     res._m[8] = zAxis.z;
 
-    return result;
-};
+    return res;
+}