softfloat 1.0.0

use crate::abs_diff;
use crate::soft_f32::F32;

type F = F32;

pub(crate) const fn pow(a: F, b: i32) -> F {
    let mut a = a;
    let recip = b < 0;
    let mut pow = abs_diff(b, 0);
    let mut mul = F::ONE;
    loop {
        if (pow & 1) != 0 {
            mul = mul.mul(a);
        }
        pow >>= 1;
        if pow == 0 {
            break;
        }
        a = a.mul(a);
    }

    if recip {
        F::ONE.div(mul)
    } else {
        mul
    }
}

#[cfg(test)]
mod test {
    #[test]
    fn sanity_check() {
        assert_eq!(f32!(2.0).powi(2), f32!(4.0))
    }
}