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
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163

mod cache;
mod radiation_pressure;
mod transform;

use std::collections::HashMap;
use std::time::Duration;

pub use cache::Cache as ModulationCache;
pub use cache::IntoCache as IntoModulationCache;
pub use radiation_pressure::IntoRadiationPressure;
pub use radiation_pressure::RadiationPressure;
pub use transform::IntoTransform as IntoModulationTransform;
pub use transform::Transform as ModulationTransform;

use crate::defined::DEFAULT_TIMEOUT;
use crate::derive::Device;
use crate::derive::Geometry;
use crate::{
    error::AUTDInternalError,
    firmware::{
        fpga::{LoopBehavior, SamplingConfig, Segment, TransitionMode},
        operation::{ModulationOp, NullOp},
    },
};

use super::DatagramST;

use rayon::prelude::*;

const PARALLEL_THRESHOLD: usize = 4;

pub trait ModulationProperty {
    fn sampling_config(&self) -> SamplingConfig;
    fn loop_behavior(&self) -> LoopBehavior;
}

/// Modulation controls the amplitude modulation data.
///
/// Modulation has following restrictions:
/// * The buffer size is up to 65536.
/// * The sampling rate is [crate::firmware::fpga::fpga_clk_freq()]/N, where N is a 32-bit unsigned integer and must be at least [crate::fpga::SAMPLING_FREQ_DIV_MIN].
#[allow(clippy::len_without_is_empty)]
pub trait Modulation: ModulationProperty {
    fn calc(&self, geometry: &Geometry) -> Result<HashMap<usize, Vec<u8>>, AUTDInternalError>;
    fn transform<F: Fn(&Device) -> Result<Vec<u8>, AUTDInternalError> + Sync>(
        geometry: &Geometry,
        f: F,
    ) -> Result<HashMap<usize, Vec<u8>>, AUTDInternalError>
    where
        Self: Sized,
    {
        #[cfg(all(feature = "force_parallel", feature = "force_serial"))]
        compile_error!("Cannot specify both force_parallel and force_serial");
        #[cfg(all(feature = "force_parallel", not(feature = "force_serial")))]
        let n = usize::MAX;
        #[cfg(all(not(feature = "force_parallel"), feature = "force_serial"))]
        let n = 0;
        #[cfg(all(not(feature = "force_parallel"), not(feature = "force_serial")))]
        let n = geometry.devices().count();

        if n > PARALLEL_THRESHOLD {
            geometry
                .devices()
                .par_bridge()
                .map(|dev| Ok((dev.idx(), f(dev)?)))
                .collect::<Result<HashMap<usize, Vec<u8>>, AUTDInternalError>>()
        } else {
            geometry
                .devices()
                .map(|dev| Ok((dev.idx(), f(dev)?)))
                .collect::<Result<HashMap<usize, Vec<u8>>, AUTDInternalError>>()
        }
    }
}

// GRCOV_EXCL_START
impl ModulationProperty for Box<dyn Modulation> {
    fn sampling_config(&self) -> SamplingConfig {
        self.as_ref().sampling_config()
    }

    fn loop_behavior(&self) -> LoopBehavior {
        self.as_ref().loop_behavior()
    }
}

impl Modulation for Box<dyn Modulation> {
    fn calc(&self, geometry: &Geometry) -> Result<HashMap<usize, Vec<u8>>, AUTDInternalError> {
        self.as_ref().calc(geometry)
    }
}

impl DatagramST for Box<dyn Modulation> {
    type O1 = ModulationOp<Self>;
    type O2 = NullOp;

    fn operation_with_segment(
        self,
        segment: Segment,
        transition_mode: Option<TransitionMode>,
    ) -> (Self::O1, Self::O2) {
        (
            Self::O1::new(self, segment, transition_mode),
            Self::O2::default(),
        )
    }

    fn timeout(&self) -> Option<Duration> {
        Some(DEFAULT_TIMEOUT)
    }
}
// GRCOV_EXCL_STOP

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

    use crate::{defined::kHz, derive::*};

    #[derive(Modulation, Clone, PartialEq, Debug)]
    pub struct TestModulation {
        pub buf: Vec<u8>,
        pub config: SamplingConfig,
        pub loop_behavior: LoopBehavior,
    }

    impl Modulation for TestModulation {
        fn calc(&self, geometry: &Geometry) -> Result<HashMap<usize, Vec<u8>>, AUTDInternalError> {
            Self::transform(geometry, |_| Ok(self.buf.clone()))
        }
    }

    #[rstest::rstest]
    #[test]
    #[case(SamplingConfig::Freq(4 * kHz))]
    fn test_sampling_config(#[case] config: SamplingConfig) {
        assert_eq!(
            config,
            TestModulation {
                config,
                buf: vec![],
                loop_behavior: LoopBehavior::infinite(),
            }
            .sampling_config()
        );
    }

    #[rstest::rstest]
    #[test]
    #[case::infinite(LoopBehavior::infinite())]
    #[case::once(LoopBehavior::once())]
    fn test_loop_behavior(#[case] loop_behavior: LoopBehavior) {
        assert_eq!(
            loop_behavior,
            TestModulation {
                config: SamplingConfig::Freq(4 * kHz),
                buf: vec![],
                loop_behavior,
            }
            .loop_behavior()
        );
    }
}