vec3 0.2.1

use core::ops::{Div, Mul};
use num_traits::real::Real;

use super::{len, sdiv};

/// # Example
/// ```
/// let mut v = vec3::new_one::<f32>();
/// assert_eq!(vec3::norm(&mut v, &vec3::new_one()), 3_f32.sqrt());
/// assert_eq!(v, [1_f32 / 3_f32.sqrt(); 3]);
/// ```
#[inline]
pub fn norm<'out, T>(out: &'out mut [T; 3], b: &[T; 3]) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T> + Div<&'b T, Output = T>,
{
    let s = len(b);
    sdiv(out, b, &s);
    s
}

/// # Example
/// ```
/// let mut v = vec3::new_one::<f32>();
/// assert_eq!(vec3::norm_mut(&mut v), 3_f32.sqrt());
/// assert_eq!(v, [1_f32 / 3_f32.sqrt(); 3]);
/// ```
#[inline]
pub fn norm_mut<'out, T>(out: &'out mut [T; 3]) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T> + Div<&'b T, Output = T>,
{
    let tmp = out.clone();
    norm(out, &tmp)
}