gnss-js 1.32.0 → 1.33.0

package/dist/{chunk-HOH27N42.js → chunk-PVOF4IGX.js}

@@ -1454,7 +1454,7 @@ function applyCorrections(corrections, geom, ctx) {
   return { rangeM: 0, phaseWindupM: 0 };
 }
 function createPppCorrections(cfg) {
-  const useSat = cfg.satPco ?? false;
+  const useSat = cfg.satPco ?? true;
   const useRcv = cfg.rcvPco ?? true;
   const useTides = cfg.tides ?? true;
   const useWindup = cfg.windup ?? true;
@@ -1788,8 +1788,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -1840,15 +1838,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);

package/dist/index.cjs

@@ -12850,7 +12850,7 @@ function applyCorrections(corrections, geom, ctx) {
   return { rangeM: 0, phaseWindupM: 0 };
 }
 function createPppCorrections(cfg) {
-  const useSat = cfg.satPco ?? false;
+  const useSat = cfg.satPco ?? true;
   const useRcv = cfg.rcvPco ?? true;
   const useTides = cfg.tides ?? true;
   const useWindup = cfg.windup ?? true;
@@ -13184,8 +13184,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -13236,15 +13234,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);

package/dist/index.js

@@ -167,7 +167,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-HOH27N42.js";
+} from "./chunk-PVOF4IGX.js";
 import {
   almanacEpochMs,
   almanacSatPosition,

package/dist/positioning.cjs

@@ -1978,7 +1978,7 @@ function applyCorrections(corrections, geom, ctx) {
   return { rangeM: 0, phaseWindupM: 0 };
 }
 function createPppCorrections(cfg) {
-  const useSat = cfg.satPco ?? false;
+  const useSat = cfg.satPco ?? true;
   const useRcv = cfg.rcvPco ?? true;
   const useTides = cfg.tides ?? true;
   const useWindup = cfg.windup ?? true;
@@ -2312,8 +2312,6 @@ function solvePpp(epochs, sp3, opts) {
       P[ci][ci] = CLK_VAR;
     }
     let nSats = 0;
-    let sumPhaseResSq = 0;
-    let nPhaseRes = 0;
     for (const v of vis) {
       const ci = clockOf(v.prn[0]);
       const common = v.range + x[ci] - v.satClkM + v.tropoFixed + zwd0 * v.mw;
@@ -2364,15 +2362,37 @@ function solvePpp(epochs, sp3, opts) {
           h[ci] = 1;
           h[IDX_ZWD] = v.mw;
           h[ai] = 1;
-          if (ekfUpdateScalar(x, P, h, innov, rPhase, 4)) {
-            sumPhaseResSq += innov * innov;
-            nPhaseRes++;
-          }
+          ekfUpdateScalar(x, P, h, innov, rPhase, 4);
         }
       }
       lastSeen.set(v.prn, ei);
       nSats++;
     }
+    let sumPhaseResSq = 0;
+    let nPhaseRes = 0;
+    {
+      const bySysRes = /* @__PURE__ */ new Map();
+      for (const v of vis) {
+        const ai = ambIdx.get(v.prn);
+        if (ai === void 0) continue;
+        const sys = v.prn[0];
+        const ci = clockOf(sys);
+        const commonPh = v.range + x[ci] - v.satClkM + v.tropoFixed + x[IDX_ZWD] * v.mw;
+        const innov = v.lIf - (commonPh + v.windupM + x[ai]);
+        if (Math.abs(innov) >= 0.15) continue;
+        const arr = bySysRes.get(sys);
+        if (arr) arr.push(innov);
+        else bySysRes.set(sys, [innov]);
+      }
+      for (const res of bySysRes.values()) {
+        const mean = res.reduce((a, b) => a + b, 0) / res.length;
+        for (const r of res) {
+          const d = r - mean;
+          sumPhaseResSq += d * d;
+          nPhaseRes++;
+        }
+      }
+    }
     for (const [prn, seen] of [...lastSeen]) {
       if (ei - seen > 20) {
         const idx = ambIdx.get(prn);

package/dist/positioning.js

@@ -24,7 +24,7 @@ import {
   sunEcef,
   toRtkEpoch,
   wlWavelength
-} from "./chunk-HOH27N42.js";
+} from "./chunk-PVOF4IGX.js";
 import "./chunk-MHGNKSSH.js";
 import "./chunk-37QNKGTC.js";
 import "./chunk-6FAL6P4G.js";

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "gnss-js",
-  "version": "1.32.0",
+  "version": "1.33.0",
   "description": "Comprehensive GNSS library for JavaScript — time scales, coordinates, RINEX parsing, RTCM3 decoding, orbit computation, and signal analysis",
   "type": "module",
   "sideEffects": false,