cherry-rs 0.2.0

Tools for designing sequential optical systems
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

use ndarray::{arr3, Array3};

use crate::{
    FieldSpec, GapSpec, ImagePlane, Pupil, PupilSampling, RealSpec, RefractiveIndexSpec,
    SequentialModel, SurfaceSpec, SurfaceType,
};

pub fn sequential_model() -> SequentialModel {
    let air = RefractiveIndexSpec {
        real: RealSpec::Constant(1.0),
        imag: None,
    };
    // Glass: N-BK7
    let nbk7 = RefractiveIndexSpec {
        real: RealSpec::Constant(1.515),
        imag: None,
    };

    let gap_0 = GapSpec {
        thickness: f64::INFINITY,
        refractive_index: air.clone(),
    };
    let gap_1 = GapSpec {
        thickness: 5.3,
        refractive_index: nbk7,
    };
    let gap_2 = GapSpec {
        thickness: 46.6,
        refractive_index: air,
    };
    let gaps = vec![gap_0, gap_1, gap_2];

    let surf_0 = SurfaceSpec::Object;
    let surf_1 = SurfaceSpec::Conic {
        semi_diameter: 12.5,
        radius_of_curvature: 25.8,
        conic_constant: 0.0,
        surf_type: SurfaceType::Refracting,
    };
    let surf_2 = SurfaceSpec::Conic {
        semi_diameter: 12.5,
        radius_of_curvature: f64::INFINITY,
        conic_constant: 0.0,
        surf_type: SurfaceType::Refracting,
    };
    let surf_3 = SurfaceSpec::Image;
    let surfaces = vec![surf_0, surf_1, surf_2, surf_3];

    let wavelengths: Vec<f64> = vec![0.567];

    SequentialModel::new(&gaps, &surfaces, &wavelengths).unwrap()
}

pub fn field_specs() -> Vec<FieldSpec> {
    vec![
        FieldSpec::Angle {
            angle: 0.0,
            pupil_sampling: PupilSampling::ChiefAndMarginalRays,
        },
        FieldSpec::Angle {
            angle: 5.0,
            pupil_sampling: PupilSampling::ChiefAndMarginalRays,
        },
    ]
}

// Paraxial View values
pub const APERTURE_STOP: usize = 1;

pub const BACK_FOCAL_DISTANCE: f64 = 46.5987;

pub const BACK_PRINCIPAL_PLANE: f64 = 1.8017;

pub const EFFECTIVE_FOCAL_LENGTH: f64 = 50.097;

pub const ENTRANCE_PUPIL: Pupil = Pupil {
    location: 0.0,
    semi_diameter: 12.5,
};

pub const EXIT_PUPIL: Pupil = Pupil {
    location: 1.8017,
    semi_diameter: 12.5,
};

pub const FRONT_FOCAL_DISTANCE: f64 = -EFFECTIVE_FOCAL_LENGTH;

pub const FRONT_PRINCIPAL_PLANE: f64 = 0.0;

// For a 5 degree field angle
pub const PARAXIAL_IMAGE_PLANE: ImagePlane = ImagePlane {
    location: 51.8987,
    semi_diameter: 4.3829,
};

// For a 5 degree field angle
// Paraxial angle = tan(field angle)
pub fn chief_ray() -> Array3<f64> {
    arr3(&[
        [[0.0], [0.087489]],
        [[0.0], [0.0577482]],
        [[0.306067], [0.087489]],
        [[4.382944], [0.087489]],
    ])
}

pub fn marginal_ray() -> Array3<f64> {
    arr3(&[
        [[12.5000], [0.0]],
        [[12.5000], [-0.1647]],
        [[11.6271], [-0.2495]],
        [[-0.0003], [-0.2495]],
    ])
}