eyepiece 0.1.0

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

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

/// 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 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
    }
}