magba 0.5.0

Magnetic computation library for Rust
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

/*
 * Magba is licensed under The 3-Clause BSD, see LICENSE.
 * Copyright 2025 Sira Pornsiriprasert <code@psira.me>
 */

use nalgebra::Vector3;

use crate::crate_utils::define_source;

define_source! {
    /// Cuboid magnet source.
    ///
    /// # References
    ///
    /// - Ortner, Michael, and Lucas Gabriel Coliado Bandeira. “Magpylib: A Free Python Package for Magnetic Field Computation.” SoftwareX 11 (January 1, 2020): 100466. <https://doi.org/10.1016/j.softx.2020.100466>.
    CuboidMagnet
    field_fn: cuboid_B
    args: {
        polarization:Vector3<T> = Vector3::z(),
        dimensions:Vector3<T> = Vector3::from_element(T::one());
            validate dimensions.iter().all(|&elem| elem >= T::zero());
            error "Dimensions must be non-negative."
    }
    arg_display: "pol={}, dim={}";
    arg_fmt: [format_vector3, format_vector3]
    docs: {
        new: {
            /// Construct a [CuboidMagnet].
            ///
            /// # Examples
            ///
            /// ```
            /// # use magba::magnets::CuboidMagnet;
            /// # use nalgebra::UnitQuaternion;
            /// let magnet = CuboidMagnet::new(
            ///     [0.0, 0.0, 0.0],              // position: Center of the cuboid (m)
            ///     UnitQuaternion::identity(),   // orientation as unit quaternion
            ///     [0.0, 0.0, 1.0],              // polarization (T)
            ///     [0.01, 0.01, 0.02],           // dimensions: Cuboid side lengths (m)
            /// );
            /// ```
        }
    }
}

#[cfg(all(test, feature = "std"))]
crate::testing_util::generate_tests! {
    CuboidMagnet
    filename: cuboid
    params: { polarization: vector![1.0, 2.0, 3.0], dimensions: vector![0.1, 0.2, 0.3]}
    rtols: {
        static: 2e-10,
        static_small: 2e-10,
        translate: 2e-10,
        rotate: 2e-10,
    }
    p95_rtols: {
        static: 2e-10,
        static_small: 2e-10,
        translate: 2e-10,
        rotate: 2e-10,
    }
    f32_rtols: {
        static: 1e-4,
        static_small: 2e-4,
        translate: 1e-4,
        rotate: 1e-4,
    }
    f32_p95_rtols: {
        static: 1e-4,
        static_small: 1e-4,
        translate: 1e-4,
        rotate: 1e-4,
    }
}

#[cfg(test)]
mod tests {
    use nalgebra::UnitQuaternion;

    use crate::magnets::CuboidMagnet;

    #[test]
    #[should_panic]
    fn test_init_validation() {
        CuboidMagnet::<f64>::new(
            [0.0; 3],
            UnitQuaternion::identity(),
            [0.0; 3],
            [0.0, -1.0, -2.0],
        );
    }

    #[test]
    #[should_panic]
    fn test_set_dimensions_validation() {
        let mut magnet = CuboidMagnet::default();
        magnet.set_dimensions([1.0, -1.0, 1.0]);
    }

    #[test]
    #[should_panic]
    fn test_with_dimensions_validation() {
        let _ = CuboidMagnet::default().with_dimensions([-1.0, 1.0, 1.0]);
    }
}