autd3-rs 0.4.0

Core async client library for the AUTD3 phased-array kit.
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

use core::f32::consts::PI;

use crate::Length;
use crate::geometry::{Point3, UnitVector3, Vector3};
use crate::value::Intensity;
use crate::value::{ControlPoint, ControlPoints};

pub fn circle(
    center: Point3<f32>,
    radius: Length,
    num_points: usize,
    normal: UnitVector3<f32>,
    intensity: Intensity,
    dst: &mut Vec<ControlPoints<1>>,
) {
    let z = Vector3::z();
    let v0 = if normal.dot(&z).abs() < 0.9 {
        z
    } else {
        Vector3::y()
    };
    let u = normal.cross(&v0).normalize();
    let v = normal.cross(&u).normalize();
    let radius = radius.mm();
    dst.clear();
    dst.extend((0..num_points).map(|i| {
        let theta = 2.0 * PI * i as f32 / num_points as f32;
        let point = center + (u * theta.cos() + v * theta.sin()) * radius;
        ControlPoints::new([ControlPoint::from(point)], intensity)
    }));
}

pub fn line(
    start: Point3<f32>,
    end: Point3<f32>,
    num_points: usize,
    intensity: Intensity,
    dst: &mut Vec<ControlPoints<1>>,
) {
    let dir = end - start;
    let denom = (num_points.max(2) - 1) as f32;
    dst.clear();
    dst.extend((0..num_points).map(|i| {
        let point = start + dir * (i as f32 / denom);
        ControlPoints::new([ControlPoint::from(point)], intensity)
    }));
}

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

    fn approx(a: Point3<f32>, b: Point3<f32>) {
        assert!((a - b).norm() < 1e-3, "{a:?} != {b:?}");
    }

    #[test]
    fn circle_in_xy_plane_walks_around_center() {
        let mut pts = Vec::new();
        circle(
            Point3::origin(),
            30.0 * mm,
            4,
            Vector3::z_axis(),
            Intensity::MAX,
            &mut pts,
        );
        assert_eq!(pts.len(), 4);
        approx(pts[0].points[0].point, Point3::new(-30.0, 0.0, 0.0));
        approx(pts[1].points[0].point, Point3::new(0.0, -30.0, 0.0));
        approx(pts[2].points[0].point, Point3::new(30.0, 0.0, 0.0));
        approx(pts[3].points[0].point, Point3::new(0.0, 30.0, 0.0));
        assert!(pts.iter().all(|p| p.intensity == Intensity::MAX));
    }

    #[test]
    fn line_interpolates_endpoints_inclusive() {
        let mut pts = Vec::new();
        line(
            Point3::new(0.0, -15.0, 0.0),
            Point3::new(0.0, 15.0, 0.0),
            3,
            Intensity(0x40),
            &mut pts,
        );
        assert_eq!(pts.len(), 3);
        approx(pts[0].points[0].point, Point3::new(0.0, -15.0, 0.0));
        approx(pts[1].points[0].point, Point3::new(0.0, 0.0, 0.0));
        approx(pts[2].points[0].point, Point3::new(0.0, 15.0, 0.0));
        assert!(pts.iter().all(|p| p.intensity == Intensity(0x40)));
    }
}