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
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
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119

use std::{fmt::Display, ops::Deref};

use num_complex::Complex;

use crate::Observer;

/// Hexagonal pupil
#[derive(Debug, Clone, serde::Serialize)]
pub struct Hexagon {
    origin: (f64, f64),
    flat_to_flat: f64,
}
impl Hexagon {
    /// Creates a new hexagon
    ///
    /// The hexagon is set by its origin and flat-to-flat size
    pub fn new(origin: (f64, f64), flat_to_flat: f64) -> Self {
        Self {
            origin,
            flat_to_flat,
        }
    }
}

impl Observer for Hexagon {
    fn diameter(&self) -> f64 {
        let (cx, cy) = self.origin;
        2. * (cx.hypot(cy) + 0.5 * self.flat_to_flat / 30f64.to_radians().cos())
    }

    fn inside_pupil(&self, x: f64, y: f64) -> bool {
        let (cx, cy) = self.origin;
        let d = self.flat_to_flat * 0.5;
        for o in [-30f64.to_radians(), 30f64.to_radians()] {
            let (so, co) = o.sin_cos();
            let xp = (x - cx) * co + (y - cy) * so;
            if xp.abs() > d {
                return false;
            }
        }
        if (y - cy).abs() > d {
            return false;
        }
        true
    }
}

/// James Webb Telescope
///
/// <img src="https://raw.githubusercontent.com/rconan/eyepiece/main/examples/jwst/telescope_pupil.png" width="20%" alt="JWST pupil">
#[derive(Debug, Clone)]
pub struct Jwst(Vec<Hexagon>);
impl Deref for Jwst {
    type Target = Vec<Hexagon>;

    fn deref(&self) -> &Self::Target {
        &self.0
    }
}
impl Display for Jwst {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        write!(
            f,
            "JWST: {}m diameter, {:.3}m^2 collection area",
            self.diameter(),
            self.area()
        )
    }
}
impl Jwst {
    /// Creates a James Webb Telescope object
    pub fn new() -> Self {
        let f2f = 1.32;
        Self(
            (0..6)
                .map(|i| {
                    let o = (30. + i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(f2f, o);
                    Hexagon::new((z.re, z.im), f2f)
                })
                .chain((0..6).map(|i| {
                    let o = (i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(3. * f2f / 3f64.sqrt(), o);
                    Hexagon::new((z.re, z.im), f2f)
                }))
                .chain((0..6).map(|i| {
                    let o = (30. + i as f64 * 60.).to_radians();
                    let z = Complex::from_polar(2. * f2f, o);
                    Hexagon::new((z.re, z.im), f2f)
                }))
                .collect(),
        )
    }
}

impl Observer for Jwst {
    fn diameter(&self) -> f64 {
        6.6
    }

    fn inside_pupil(&self, x: f64, y: f64) -> bool {
        for hex in self.iter() {
            if hex.inside_pupil(x, y) {
                return true;
            }
        }
        false
    }
}

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