target-match 0.3.3

Given a telescope pointing and field of view, rank which catalogued sky objects fall on the frame. Matches by sky position; holds no catalogue data and does no I/O.
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

//! Property-based tests of the matcher invariants (SC-007, SC-008).
//!
//! The coordinate-math properties this suite once carried (separation
//! symmetry, sexagesimal and precession round-trips) moved to `skymath`,
//! which owns that geometry; only matcher behaviour is tested here.

use proptest::prelude::*;
use skymath::{Angle, Equatorial};
use target_match::{rank, Constraint, Matcher, SkyObject};

#[derive(Clone)]
struct Obj {
    id: usize,
    ra: f64,
    dec: f64,
}
impl SkyObject for Obj {
    fn position(&self) -> Equatorial {
        Equatorial::j2000(Angle::from_degrees(self.ra), Angle::from_degrees(self.dec)).unwrap()
    }
}

fn eq(ra: f64, dec: f64) -> Equatorial {
    Equatorial::j2000(Angle::from_degrees(ra), Angle::from_degrees(dec)).unwrap()
}

proptest! {
    /// SC-007 / SC-008: the prebuilt Matcher returns identical results to the stateless
    /// scan for the same objects, pointing, and constraint.
    #[test]
    fn matcher_equals_rank(
        seed in prop::collection::vec((0.0..360.0f64, -90.0..90.0f64), 1..40),
        pra in 0.0..360.0f64, pdec in -90.0..90.0f64, radius in 0.1..40.0f64,
    ) {
        let cat: Vec<Obj> = seed.iter().enumerate().map(|(i, &(ra, dec))| Obj { id: i, ra, dec }).collect();
        let pointing = eq(pra, pdec);

        for c in [
            Constraint::circular(Angle::from_degrees(radius)).all(),
            Constraint::circular(Angle::from_degrees(radius)).nearest_n(5),
            Constraint::circular(Angle::from_degrees(radius)).nearest_one(),
        ] {
            let via_rank: Vec<usize> = rank(pointing, &cat, c).iter().map(|m| m.object.id).collect();
            let matcher = Matcher::from_objects(cat.clone());
            let via_index: Vec<usize> = matcher.query(pointing, c).iter().map(|m| m.object.id).collect();
            prop_assert_eq!(via_rank, via_index);
        }
    }
}