spdcalc 2.0.1

SPDCalc, the Spontaneous Parametric Downconversion Calculator
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

//! Sellmeier standard form
use super::*;
use dim::{
  f64prefixes::MICRO,
  ucum::{self, M},
};

/// Standard form:
/// > n² = a + b1 * λ² / (λ² - c1) + b2 * λ² / (λ² - c2) + b3 * λ² / (λ² - c3)
#[derive(Debug, Clone, Copy, Serialize, Deserialize)]
pub struct SellmeierStandard {
  /// a
  pub a: [f64; 3],
  /// b1
  pub b1: [f64; 3],
  /// b2
  pub b2: [f64; 3],
  /// b3
  pub b3: [f64; 3],

  /// c1
  pub c1: [f64; 3],
  /// c2
  pub c2: [f64; 3],
  /// c3
  pub c3: [f64; 3],
}

impl SellmeierEquation for SellmeierStandard {
  fn get_indices(&self, wavelength: Wavelength) -> Indices {
    let a = na::Vector3::from_column_slice(&self.a);

    let b1 = na::Vector3::from_column_slice(&self.b1);
    let b2 = na::Vector3::from_column_slice(&self.b2);
    let b3 = na::Vector3::from_column_slice(&self.b3);

    let c1 = na::Vector3::from_column_slice(&self.c1);
    let c2 = na::Vector3::from_column_slice(&self.c2);
    let c3 = na::Vector3::from_column_slice(&self.c3);

    // convert to micrometers
    let l = *(wavelength / (MICRO * M));
    let l_sq = l * l;
    let one_by_l_sq = Vector3::repeat(l_sq);

    let n = a
      + (b1.component_div(&(one_by_l_sq - c1))
        + b2.component_div(&(one_by_l_sq - c2))
        + b3.component_div(&(one_by_l_sq - c3)))
        * l_sq;

    ucum::Unitless::new(n.map(|i| i.sqrt()))
  }
}