lodestone_core 0.2.3

User friendly magnetic field calculations in 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

use crate::magnets::magnet3d::bulirsch::cel;
use crate::points::Point2;
use crate::{MagnetError, PI};
// use crate::magnets::

/// Calculate the magnetic field due to a solenoid at any point
/// returns Bz,Br as a Point2 structure
pub fn solenoid_field(a: f64, b: f64, rho: f64, z: f64, b_0: f64) -> Result<Point2, MagnetError> {
    let zp = z + b;
    let zn = z - b;

    let rho_a = rho + a;
    let rho_na = rho - a;

    let alpha_p = a / (zp.powi(2) + (rho_a).powi(2)).sqrt();
    let alpha_n = a / (zn.powi(2) + (rho_a).powi(2)).sqrt();

    let beta_p = zp / (zp.powi(2) + (rho_a).powi(2)).sqrt();
    let beta_n = zn / (zn.powi(2) + (rho_a).powi(2)).sqrt();

    let gamma = (rho_na) / (rho_a);

    let kp = ((zp.powi(2) + (rho_na).powi(2)) / (zp.powi(2) + (rho_a).powi(2))).sqrt();

    let kn = ((zn.powi(2) + (rho_na).powi(2)) / (zn.powi(2) + (rho_a).powi(2))).sqrt();

    let b_r = (b_0 / PI) * (alpha_p * cel(kp, 1.0, 1.0, -1.0) - alpha_n * cel(kn, 1.0, 1.0, -1.0));

    let b_z = ((b_0 / PI) * a / (rho_a))
        * (beta_p * cel(kp, gamma.powi(2), 1.0, gamma)
            - beta_n * cel(kn, gamma.powi(2), 1.0, gamma));

    Ok(Point2::new(b_z, b_r))
}