panchang-ts 0.4.0 → 0.4.1

package/dist/index.cjs

@@ -76,6 +76,13 @@ var JUPITER_GM = 2825345909524226e-22;
 var SATURN_GM = 8459715185680659e-23;
 var URANUS_GM = 1292024916781969e-23;
 var NEPTUNE_GM = 1524358900784276e-23;
+function VerifyBoolean(b) {
+  if (b !== true && b !== false) {
+    console.trace();
+    throw `Value is not boolean: ${b}`;
+  }
+  return b;
+}
 function VerifyNumber(x) {
   if (!Number.isFinite(x)) {
     console.trace();
@@ -1598,6 +1605,8 @@ function SunPosition(date) {
 }
 function Equator(body, date, observer, ofdate, aberration) {
   VerifyObserver(observer);
+  VerifyBoolean(ofdate);
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   const gc_observer = geo_pos(time, observer);
   const gc = GeoVector(body, time, aberration);
@@ -2103,24 +2112,30 @@ var BodyPosition = class {
   }
 };
 function BackdatePosition(date, observerBody, targetBody, aberration) {
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   if (UserDefinedStar(targetBody)) {
     const tvec = HelioVector(targetBody, time);
-    {
+    if (aberration) {
       const ostate = HelioState(observerBody, time);
       const rvec = new Vector(tvec.x - ostate.x, tvec.y - ostate.y, tvec.z - ostate.z, time);
       const s = C_AUDAY / rvec.Length();
       return new Vector(rvec.x + ostate.vx / s, rvec.y + ostate.vy / s, rvec.z + ostate.vz / s, time);
     }
+    const ovec = HelioVector(observerBody, time);
+    return new Vector(tvec.x - ovec.x, tvec.y - ovec.y, tvec.z - ovec.z, time);
   }
   let observerPos;
-  {
+  if (aberration) {
     observerPos = new Vector(0, 0, 0, time);
+  } else {
+    observerPos = HelioVector(observerBody, time);
   }
   const bpos = new BodyPosition(observerBody, targetBody, aberration, observerPos);
   return CorrectLightTravel((t) => bpos.Position(t), time);
 }
 function GeoVector(body, date, aberration) {
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   switch (body) {
     case Body.Earth:
@@ -2266,13 +2281,6 @@ function Search(f, t1, t2, options) {
     return null;
   }
 }
-function EclipticLongitude(body, date) {
-  if (body === Body.Sun)
-    throw "Cannot calculate heliocentric longitude of the Sun.";
-  const hv = HelioVector(body, date);
-  const eclip = Ecliptic(hv);
-  return eclip.elon;
-}
 var AtmosphereInfo = class {
   constructor(pressure, temperature, density) {
     this.pressure = pressure;
@@ -4028,12 +4036,13 @@ function meanObliquity(T) {
   return 23.439291111 - 0.013004167 * T - 164e-9 * T * T + 504e-9 * T * T * T;
 }
 function getTropicalPlanetLongitude(body, date) {
-  return EclipticLongitude(body, MakeTime(date));
+  const vec = GeoVector(body, MakeTime(date), true);
+  return Ecliptic(vec).elon;
 }
 function isRetrograde(body, date) {
   const dt = 36e5;
-  const lon0 = EclipticLongitude(body, MakeTime(new Date(date.getTime() - dt)));
-  const lon1 = EclipticLongitude(body, MakeTime(new Date(date.getTime() + dt)));
+  const lon0 = Ecliptic(GeoVector(body, MakeTime(new Date(date.getTime() - dt)), false)).elon;
+  const lon1 = Ecliptic(GeoVector(body, MakeTime(new Date(date.getTime() + dt)), false)).elon;
   let delta = lon1 - lon0;
   if (delta > 180) delta -= 360;
   if (delta < -180) delta += 360;

package/dist/index.js

@@ -74,6 +74,13 @@ var JUPITER_GM = 2825345909524226e-22;
 var SATURN_GM = 8459715185680659e-23;
 var URANUS_GM = 1292024916781969e-23;
 var NEPTUNE_GM = 1524358900784276e-23;
+function VerifyBoolean(b) {
+  if (b !== true && b !== false) {
+    console.trace();
+    throw `Value is not boolean: ${b}`;
+  }
+  return b;
+}
 function VerifyNumber(x) {
   if (!Number.isFinite(x)) {
     console.trace();
@@ -1596,6 +1603,8 @@ function SunPosition(date) {
 }
 function Equator(body, date, observer, ofdate, aberration) {
   VerifyObserver(observer);
+  VerifyBoolean(ofdate);
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   const gc_observer = geo_pos(time, observer);
   const gc = GeoVector(body, time, aberration);
@@ -2101,24 +2110,30 @@ var BodyPosition = class {
   }
 };
 function BackdatePosition(date, observerBody, targetBody, aberration) {
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   if (UserDefinedStar(targetBody)) {
     const tvec = HelioVector(targetBody, time);
-    {
+    if (aberration) {
       const ostate = HelioState(observerBody, time);
       const rvec = new Vector(tvec.x - ostate.x, tvec.y - ostate.y, tvec.z - ostate.z, time);
       const s = C_AUDAY / rvec.Length();
       return new Vector(rvec.x + ostate.vx / s, rvec.y + ostate.vy / s, rvec.z + ostate.vz / s, time);
     }
+    const ovec = HelioVector(observerBody, time);
+    return new Vector(tvec.x - ovec.x, tvec.y - ovec.y, tvec.z - ovec.z, time);
   }
   let observerPos;
-  {
+  if (aberration) {
     observerPos = new Vector(0, 0, 0, time);
+  } else {
+    observerPos = HelioVector(observerBody, time);
   }
   const bpos = new BodyPosition(observerBody, targetBody, aberration, observerPos);
   return CorrectLightTravel((t) => bpos.Position(t), time);
 }
 function GeoVector(body, date, aberration) {
+  VerifyBoolean(aberration);
   const time = MakeTime(date);
   switch (body) {
     case Body.Earth:
@@ -2264,13 +2279,6 @@ function Search(f, t1, t2, options) {
     return null;
   }
 }
-function EclipticLongitude(body, date) {
-  if (body === Body.Sun)
-    throw "Cannot calculate heliocentric longitude of the Sun.";
-  const hv = HelioVector(body, date);
-  const eclip = Ecliptic(hv);
-  return eclip.elon;
-}
 var AtmosphereInfo = class {
   constructor(pressure, temperature, density) {
     this.pressure = pressure;
@@ -4026,12 +4034,13 @@ function meanObliquity(T) {
   return 23.439291111 - 0.013004167 * T - 164e-9 * T * T + 504e-9 * T * T * T;
 }
 function getTropicalPlanetLongitude(body, date) {
-  return EclipticLongitude(body, MakeTime(date));
+  const vec = GeoVector(body, MakeTime(date), true);
+  return Ecliptic(vec).elon;
 }
 function isRetrograde(body, date) {
   const dt = 36e5;
-  const lon0 = EclipticLongitude(body, MakeTime(new Date(date.getTime() - dt)));
-  const lon1 = EclipticLongitude(body, MakeTime(new Date(date.getTime() + dt)));
+  const lon0 = Ecliptic(GeoVector(body, MakeTime(new Date(date.getTime() - dt)), false)).elon;
+  const lon1 = Ecliptic(GeoVector(body, MakeTime(new Date(date.getTime() + dt)), false)).elon;
   let delta = lon1 - lon0;
   if (delta > 180) delta -= 360;
   if (delta < -180) delta += 360;

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "panchang-ts",
-  "version": "0.4.0",
+  "version": "0.4.1",
   "description": "Pure TypeScript Hindu Panchang calculations. Tithi, Nakshatra, Yoga, Karana, Vara, and more. Offline-first, React Native compatible.",
   "author": "Ishank",
   "license": "MIT",