eyepiece 0.8.1

A crate to generate star fields as seen with different telescopes
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

use std::fmt::Display;

use crate::Observer;
use geotrans::{Conic, Segment, SegmentTrait, Transform, M1};

#[derive(Debug, Clone)]
/// Giant Magellan Telescope
///
/// <img src="https://raw.githubusercontent.com/rconan/eyepiece/main/examples/gmt/telescope_pupil.png" width="20%" alt="GMT pupil">
pub struct Gmt;
impl Gmt {
    /// Creates a Giant Magellan Telescope object
    pub fn new() -> Self {
        Self
    }
}
impl Display for Gmt {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        write!(
            f,
            "GMT: {}m diameter, {:.3}m^2 collection area",
            self.diameter(),
            self.area()
        )
    }
}
impl Observer for Gmt {
    fn diameter(&self) -> f64 {
        25.5
    }

    fn inside_pupil(&self, x: f64, y: f64) -> bool {
        let m1 = Conic::m1();
        let h = m1.height(x.hypot(y));
        let pt = [x, y, 3.9 + h];
        let r_outer = 8.365 * 0.5;
        let r_inner = 3.6 * 0.5;
        for i in 1..=7 {
            let xyz = pt.fro(Segment::<M1>::new(i)).unwrap();
            let r = xyz[0].hypot(xyz[1]);
            if i == 7 && r < r_inner {
                return false;
            }
            if r <= r_outer {
                return true;
            }
        }
        false
    }
}

use serde::ser::{Serialize, Serializer};
impl Serialize for Gmt {
    fn serialize<S>(&self, serializer: S) -> Result<S::Ok, S::Error>
    where
        S: Serializer,
    {
        "GMT".serialize(serializer)
    }
}