Documentation
use core::ops::{Add, Mul};
use num_traits::real::Real;
/// # Example
/// ```
/// assert_eq!(vec4::dot_values(&1_f32, &1_f32, &1_f32, &1_f32, &1_f32, &1_f32, &1_f32, &1_f32), 4_f32);
/// ```
#[inline(always)]
pub fn dot_values<T>(ax: &T, ay: &T, az: &T, aw: &T, bx: &T, by: &T, bz: &T, bw: &T) -> T
where
    T: Add<T, Output = T>,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    ax * bx + ay * by + az * bz + aw * bw
}
/// # Example
/// ```
/// assert_eq!(vec4::len_values_sq(&1_f32, &1_f32, &1_f32, &1_f32), 4_f32);
/// ```
#[inline]
pub fn len_values_sq<T>(x: &T, y: &T, z: &T, w: &T) -> T
where
    T: Add<T, Output = T>,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    dot_values(x, y, z, w, x, y, z, w)
}
/// # Example
/// ```
/// assert_eq!(vec4::inv_len_values(&1_f32, &1_f32, &1_f32, &1_f32), 1_f32 / 4_f32.sqrt());
/// ```
#[inline]
pub fn inv_len_values<T>(x: &T, y: &T, z: &T, w: &T) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    let len_sq = len_values_sq(x, y, z, w);

    if len_sq.is_zero() {
        T::zero()
    } else {
        T::one() / len_sq.sqrt()
    }
}
/// # Example
/// ```
/// assert_eq!(vec4::len_values(&1_f32, &1_f32, &1_f32, &1_f32), 4_f32.sqrt());
/// ```
#[inline]
pub fn len_values<T>(x: &T, y: &T, z: &T, w: &T) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    let len_sq = len_values_sq(x, y, z, w);

    if len_sq.is_zero() {
        T::zero()
    } else {
        len_sq.sqrt()
    }
}
/// # Example
/// ```
/// assert_eq!(vec4::len_sq::<f32>(&vec4::new_one()), 4_f32);
/// ```
#[inline(always)]
pub fn len_sq<T>(v: &[T; 4]) -> T
where
    T: Add<T, Output = T>,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    len_values_sq(&v[0], &v[1], &v[2], &v[3])
}
/// # Example
/// ```
/// assert_eq!(vec4::dot::<f32>(&vec4::new_one(), &vec4::new_one()), 4_f32);
/// ```
#[inline(always)]
pub fn dot<T>(a: &[T; 4], b: &[T; 4]) -> T
where
    T: Add<T, Output = T>,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    dot_values(&a[0], &a[1], &a[2], &a[3], &b[0], &b[1], &b[2], &b[3])
}
/// # Example
/// ```
/// assert_eq!(vec4::len::<f32>(&vec4::new_one()), 4_f32.sqrt());
/// ```
#[inline(always)]
pub fn len<T>(v: &[T; 4]) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    len_values(&v[0], &v[1], &v[2], &v[3])
}
/// # Example
/// ```
/// assert_eq!(vec4::inv_len::<f32>(&vec4::new_one()), 1_f32 / 4_f32.sqrt());
/// ```
#[inline(always)]
pub fn inv_len<T>(v: &[T; 4]) -> T
where
    T: Real,
    for<'a, 'b> &'a T: Mul<&'b T, Output = T>,
{
    inv_len_values(&v[0], &v[1], &v[2], &v[3])
}