Improve performance (#230)

* Heavily reduce PPP processing time (by about x50)
* Improve -g -o Observation plotting performances

* reduce time to render this plot (slowest to render)
 by iterating the whole set only once, at the expense of more
 complex inner data types

* Improve doc
* Remove secondary cfg print: print only once
* Position: orbit interpolation improvements
* Code based PPP 

---------

Signed-off-by: Guillaume W. Bres <[email protected]>
Co-authored-by: Laurențiu Nicola <[email protected]>

Cargo.toml

@@ -8,7 +8,6 @@ members = [
     "rinex-qc", 
     "rinex-cli", 
     "rnx2crx", 
-    "sinex",
     "sp3",
     "ublox-rnx", 
 ]

rinex-cli/Cargo.toml

@@ -42,8 +42,9 @@ plotly =  "0.8.4"
 # plotly = { git = "https://github.com/gwbres/plotly", branch = "density-mapbox" }
 
 # solver
-gnss-rtk = { version = "0.4.3", features = ["serde"] }
+# gnss-rtk = { version = "0.4.3", features = ["serde"] }
 # gnss-rtk = { path = "../../rtk-rs/gnss-rtk", features = ["serde"] }
+gnss-rtk = { git = "https://github.com/rtk-rs/gnss-rtk", branch = "main", features = ["serde"] }
 
 # cggtts
 cggtts = { version = "4.1.4", features = ["serde", "scheduler"] }

rinex-cli/config/rtk/gpst_ppp_basic.json

@@ -0,0 +1,18 @@
+{
+    "method": "CodePPP",
+    "timescale": "GPST",
+    "interp_order": 13,
+    "min_sv_elevation": 1.0,
+    "solver": {
+        "filter": "LSQ",
+        "gdop_threshold": 10.0
+    },
+    "modeling": {
+        "iono_delay": true,
+        "tropo_delay": true,
+        "earth_rotation": true,
+        "sv_clock_bias": true,
+        "sv_total_group_delay": true,
+        "relativistic_clock_bias": true
+    } 
+}

rinex-cli/config/rtk/gpst_10sv_basic.json → rinex-cli/config/rtk/gpst_spp_basic.json

@@ -1,18 +1,18 @@
 {
+    "method": "SPP",
     "timescale": "GPST",
     "interp_order": 11,
-    "max_sv": 10,
-    "min_sv_elev": 20.0,
-    "min_snr": 20.0,
+    "min_sv_elevation": 10.0,
     "solver": {
-        "gdop_threshold": 3.0
+        "filter": "LSQ",
+        "gdop_threshold": 10.0
     },
     "modeling": {
+        "iono_delay": true,
+        "tropo_delay": true,
+        "earth_rotation": true,
         "sv_clock_bias": true,
         "sv_total_group_delay": true,
-        "relativistic_clock_bias": true,
-        "tropo_delay": true,
-        "iono_delay": true,
-        "earth_rotation": true
+        "relativistic_clock_bias": true
     }
 }

rinex-cli/src/cli/positioning.rs

@@ -17,13 +17,8 @@ You should provide Observations from a unique receiver.")
             .action(ArgAction::Append)
             .help("Pass a Position Solver configuration file (JSON).
 [https://docs.rs/gnss-rtk/latest/gnss_rtk/prelude/struct.Config.html] is the structure to represent in JSON.
+Refer to [https://docs.rs/gnss-rtk/latest/gnss_rtk/prelude/enum.Method.html] for solving strategies.
 See [] for meaningful examples."))
-        .arg(Arg::new("spp")
-            .long("spp")
-            .action(ArgAction::SetTrue)
-            .help("Force resolution method to Single Point Positioning (SPP).
-Otherwise, the Default method is used.
-Refer to [https://docs.rs/gnss-rtk/latest/gnss_rtk/prelude/enum.Method.html]."))
         .arg(Arg::new("gpx")
             .long("gpx")
             .action(ArgAction::SetTrue)