outfit 2.1.0

Orbit determination toolkit in Rust. Provides astrometric parsing, observer management, and initial orbit determination (Gauss method) with JPL ephemeris support.
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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96

use nom::{
    bytes::complete::{tag, take_until, take_while1},
    character::complete::{char, multispace0},
    combinator::{map, opt},
    number::complete::float,
    sequence::{preceded, separated_pair},
    IResult, Parser,
};

use crate::error_models::ParseResult;

fn is_word_char(c: char) -> bool {
    c.is_alphanumeric() || "*".contains(c)
}

fn parse_word(input: &str) -> IResult<&str, &str> {
    take_while1(is_word_char)(input)
}

fn parse_station(input: &str) -> IResult<&str, &str> {
    parse_word(input)
}

fn parse_catalog_codes(input: &str) -> IResult<&str, Vec<String>> {
    preceded(
        tag("c="),
        map(take_while1(is_word_char), |s: &str| {
            s.chars().map(|c| c.to_string()).collect()
        }),
    )
    .parse(input)
}

fn parse_rms_values(input: &str) -> IResult<&str, (f32, f32)> {
    preceded(
        take_until("@"),
        preceded(
            tag("@"),
            separated_pair(
                preceded(multispace0, float),
                preceded(multispace0, char(',')),
                preceded(multispace0, float),
            ),
        ),
    )
    .parse(input)
}

pub fn parse_vfcc17_line(input: &str) -> ParseResult {
    let (input, remain) = opt(take_until("!")).parse(input)?; // Ignore comments
    let input = remain.unwrap_or(input).trim();

    map(
        (
            parse_station,
            take_until("c="),
            parse_catalog_codes,
            parse_rms_values,
        ),
        |(station, _, catalogs, (rmsa, rmsd))| {
            catalogs
                .into_iter()
                .map(|cat| ((station.to_string(), cat), (rmsa, rmsd)))
                .collect()
        },
    )
    .parse(input)
}

#[cfg(test)]
mod test_vfcc17_parser {
    use super::*;

    #[test]
    fn test_vfcc17_parser() {
        let input = "ALL t=cBCVn c=*          p=    >            <            @  1.00,  1.00 ! Unknown catalog";
        let result = parse_vfcc17_line(input);
        assert!(result.is_ok());
        let (_, parsed) = result.unwrap();
        assert_eq!(parsed.len(), 1);
        assert_eq!(parsed[0].0 .0, "ALL");
        assert_eq!(parsed[0].0 .1, "*");
        assert_eq!(parsed[0].1 .0, 1.);
        assert_eq!(parsed[0].1 .1, 1.);

        let input = "568 t=cC    c=t          p=_   >            <            @  0.20,  0.20  ! Micheli updated ";
        let result = parse_vfcc17_line(input);
        assert!(result.is_ok());
        let (_, parsed) = result.unwrap();
        assert_eq!(parsed.len(), 1);
        assert_eq!(parsed[0].0 .0, "568");
        assert_eq!(parsed[0].0 .1, "t");
        assert_eq!(parsed[0].1 .0, 0.2);
        assert_eq!(parsed[0].1 .1, 0.2);
    }
}