rinex 0.22.0

RINEX file parsing, analysis and production
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

use crate::meteo::Record as MetRecord;
use crate::observation::Record as ObsRecord;
use crate::Rinex;

/*
 * Test: panic if given RINEX content is not equal to given constant
 */
pub fn is_constant_rinex(rnx: &Rinex, constant: f64, tolerance: f64) {
    if let Some(record) = rnx.record.as_obs() {
        is_constant_obs_record(record, constant, tolerance)
    } else if let Some(record) = rnx.record.as_meteo() {
        is_constant_meteo_record(record, constant, tolerance)
    } else {
        unimplemented!("is_constant_rinex({})", rnx.header.rinex_type);
    }
}

pub fn is_null_rinex(rnx: &Rinex, tolerance: f64) {
    is_constant_rinex(rnx, 0.0_f64, tolerance)
}

fn is_constant_obs_record(record: &ObsRecord, constant: f64, tolerance: f64) {
    for (_, (clk, svnn)) in record {
        if let Some(clk) = clk {
            let err = (clk - constant).abs();
            if err > tolerance {
                panic!("rcvr clock {} != {}", clk, constant);
            }
        }
        for (_, observables) in svnn {
            for (observable, observation) in observables {
                let err = (observation.obs - constant).abs();
                if err > tolerance {
                    panic!(
                        "{} observation {} != {}",
                        observable, observation.obs, constant
                    );
                }
            }
        }
    }
}

fn is_constant_meteo_record(record: &MetRecord, constant: f64, tolerance: f64) {
    for (_, observables) in record {
        for (observable, observation) in observables {
            let err = (observation - constant).abs();
            if err > tolerance {
                panic!("{} observation {} != {}", observable, observation, constant);
            }
        }
    }
}