autd3-link-visualizer 27.0.0

Link for visualizing the output of the AUTD3
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
115
116
117
118
119
120
121

use std::path::Path;

use anyhow::Result;

use autd3::prelude::*;
use autd3_link_visualizer::*;

#[cfg(not(feature = "python"))]
use PlotConfig as Config;
#[cfg(feature = "python")]
use PyPlotConfig as Config;

#[tokio::main]
async fn main() -> Result<()> {
    #[cfg(feature = "python")]
    let link = Visualizer::python();
    #[cfg(not(feature = "python"))]
    let link = Visualizer::plotters();

    #[cfg(feature = "gpu")]
    let link = link.with_gpu(-1);

    let mut autd = Controller::builder([AUTD3::new(Vector3::zeros())])
        .open(link)
        .await?;

    let center = autd.geometry().center() + Vector3::new(0., 0., 150.0 * mm);

    let g = Focus::new(center);
    let m = Square::new(150. * Hz);

    autd.send((m, g)).await?;

    autd.link().plot_phase(
        Config {
            fname: Path::new("phase.png").into(),
            ..Config::default()
        },
        Segment::S0,
        0,
    )?;
    autd.link().plot_field(
        Config {
            fname: Path::new("x.png").into(),
            ..Config::default()
        },
        PlotRange {
            x_range: center.x - 50.0..center.x + 50.0,
            y_range: center.y..center.y,
            z_range: center.z..center.z,
            resolution: 1.,
        },
        Segment::S0,
        0,
    )?;
    autd.link().plot_field(
        Config {
            fname: Path::new("xy.png").into(),
            ..Config::default()
        },
        PlotRange {
            x_range: center.x - 20.0..center.x + 20.0,
            y_range: center.y - 30.0..center.y + 30.0,
            z_range: center.z..center.z,
            resolution: 1.,
        },
        Segment::S0,
        0,
    )?;
    autd.link().plot_field(
        Config {
            fname: Path::new("yz.png").into(),
            ..Config::default()
        },
        PlotRange {
            x_range: center.x..center.x,
            y_range: center.y - 30.0..center.y + 30.0,
            z_range: 0.0..center.z + 50.0,
            resolution: 2.,
        },
        Segment::S0,
        0,
    )?;
    autd.link().plot_field(
        Config {
            fname: Path::new("zx.png").into(),
            ticks_step: 20.,
            ..Config::default()
        },
        PlotRange {
            x_range: center.x - 30.0..center.x + 30.0,
            y_range: center.y..center.y,
            z_range: 0.0..center.z + 50.0,
            resolution: 2.,
        },
        Segment::S0,
        0,
    )?;

    autd.link().plot_modulation(
        Config {
            fname: Path::new("mod.png").into(),
            ..Config::default()
        },
        Segment::S0,
    )?;

    // Calculate acoustic pressure without plotting
    let p = autd.link().calc_field(&[center], Segment::S0, 0)?;
    println!(
        "Acoustic pressure at ({}, {}, {}) = {} [Pa]",
        center.x,
        center.y,
        center.z,
        p[0].norm()
    );

    autd.close().await?;

    Ok(())
}