surge-math 0.2.11-alpha.0

crate::ix!();

///# Safety 
///
///input must point to a contiguous range of (NQUADS * 4) values
pub unsafe fn clear_block<const NQUADS: usize>(input: *mut f32) 
{
    let zero: __m128 = _mm_set1_ps(0.0);

    for i in (0..(NQUADS << 2)).step_by(4)
    {
        _mm_store_ps(input.add(i), zero);
    }
}

///# Safety 
///
///input must point to a contiguous range of (NQUADS * 4) values
pub unsafe fn clear_block_antidenormalnoise<const NQUADS: usize>( input: *mut f32) 
{
    let smallvalue: __m128 = 
        _mm_set_ps(
            0.000000000000001, 
            0.000000000000001, 
            -0.000000000000001, 
            -0.000000000000001);

    for i in (0..(NQUADS << 2)).step_by(8)
    {
        _mm_store_ps(input.add(i),     smallvalue);
        _mm_store_ps(input.add(i + 4), smallvalue);
    }
}