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
/*
 * File: defined.rs
 * Project: src
 * Created Date: 05/12/2022
 * Author: Shun Suzuki
 * -----
 * Last Modified: 22/11/2023
 * Modified By: Shun Suzuki (suzuki@hapis.k.u-tokyo.ac.jp)
 * -----
 * Copyright (c) 2022-2023 Shun Suzuki. All rights reserved.
 *
 */

#[cfg(feature = "single_float")]
mod float_def {
    pub use f32 as float;
    pub use std::f32::consts::PI;
}
#[cfg(not(feature = "single_float"))]
mod float_def {
    pub use f64 as float;
    pub use std::f64::consts::PI;
}

pub use float_def::*;

#[cfg(feature = "use_meter")]
mod unit {
    pub const METER: super::float = 1.0;
}
#[cfg(not(feature = "use_meter"))]
mod unit {
    pub const METER: super::float = 1000.0;
}
pub use unit::*;

pub const MILLIMETER: float = METER / 1000.0;

pub const ULTRASOUND_FREQUENCY: float = 40e3;

pub type Complex = nalgebra::Complex<float>;

pub const ABSOLUTE_THRESHOLD_OF_HEARING: float = 20e-6;

pub const T4010A1_AMPLITUDE: float = 23.77004454874038 * 300.0 * MILLIMETER; // [Pa*mm] (121.5 dB SPL @ 30cm)

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn t4010a1_amplitude() {
        assert_approx_eq::assert_approx_eq!(
            121.5,
            20.0 * (T4010A1_AMPLITUDE * 1.0 / (300.0 * MILLIMETER) / ABSOLUTE_THRESHOLD_OF_HEARING)
                .log10()
        );
    }
}