src/core/NEON/NEAsymm.inl

/*
 * Copyright (c) 2017-2020, 2023 Arm Limited.
 *
 * SPDX-License-Identifier: MIT
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to
 * deal in the Software without restriction, including without limitation the
 * rights to use, copy, modify, merge, publish, distribute, sublicense, and/or
 * sell copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in all
 * copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
 * SOFTWARE.
 */

#include "arm_compute/core/Rounding.h"

namespace arm_compute
{
template <RoundingPolicy round_policy>
inline qasymm8x16_t vmlaq_qasymm8(qasymm8x16_t vd, float32x4_t vs, float32x4_t vo)
{
    // Convert uint8 vectors to uint16 vectors
    const uint8x8_t vd_low        = vget_low_u8(vd);
    const uint8x8_t vd_high       = vget_high_u8(vd);
    uint16x8_t      vd_low_u16x8  = vmovl_u8(vd_low);
    uint16x8_t      vd_high_u16x8 = vmovl_u8(vd_high);
    // Convert uint16 vectors to uint32 vectors
    uint32x4_t A_u32x4 = vmovl_u16(vget_low_u16(vd_low_u16x8));
    uint32x4_t B_u32x4 = vmovl_u16(vget_high_u16(vd_low_u16x8));
    uint32x4_t C_u32x4 = vmovl_u16(vget_low_u16(vd_high_u16x8));
    uint32x4_t D_u32x4 = vmovl_u16(vget_high_u16(vd_high_u16x8));
    // Convert uint32 vectors to float32 vectors
    float32x4_t A_f32x4 = vcvtq_f32_u32(A_u32x4);
    float32x4_t B_f32x4 = vcvtq_f32_u32(B_u32x4);
    float32x4_t C_f32x4 = vcvtq_f32_u32(C_u32x4);
    float32x4_t D_f32x4 = vcvtq_f32_u32(D_u32x4);
    // vd = vd*vs + vo
    A_f32x4 = vmlaq_f32(vo, A_f32x4, vs);
    B_f32x4 = vmlaq_f32(vo, B_f32x4, vs);
    C_f32x4 = vmlaq_f32(vo, C_f32x4, vs);
    D_f32x4 = vmlaq_f32(vo, D_f32x4, vs);
    // Convert float32 vectors to uint32 vectors
#if __aarch64__
    if (round_policy == RoundingPolicy::TO_NEAREST_EVEN)
    {
        A_u32x4 = vcvtnq_u32_f32(A_f32x4);
        B_u32x4 = vcvtnq_u32_f32(B_f32x4);
        C_u32x4 = vcvtnq_u32_f32(C_f32x4);
        D_u32x4 = vcvtnq_u32_f32(D_f32x4);
    }
    else if (round_policy == RoundingPolicy::TO_NEAREST_UP)
    {
        A_u32x4 = vcvtaq_u32_f32(A_f32x4);
        B_u32x4 = vcvtaq_u32_f32(B_f32x4);
        C_u32x4 = vcvtaq_u32_f32(C_f32x4);
        D_u32x4 = vcvtaq_u32_f32(D_f32x4);
    }
    else
    {
        A_u32x4 = vcvtq_u32_f32(A_f32x4);
        B_u32x4 = vcvtq_u32_f32(B_f32x4);
        C_u32x4 = vcvtq_u32_f32(C_f32x4);
        D_u32x4 = vcvtq_u32_f32(D_f32x4);
    }
#else  // #if __aarch64__
    // rounding mode only supported in aarch64
    A_u32x4 = vcvtq_u32_f32(A_f32x4);
    B_u32x4 = vcvtq_u32_f32(B_f32x4);
    C_u32x4 = vcvtq_u32_f32(C_f32x4);
    D_u32x4 = vcvtq_u32_f32(D_f32x4);
#endif // #if __aarch64__
    // Convert uint32 vectors to uint16 vectors (with saturation)
    vd_low_u16x8  = vcombine_u16(vqmovn_u32(A_u32x4), vqmovn_u32(B_u32x4));
    vd_high_u16x8 = vcombine_u16(vqmovn_u32(C_u32x4), vqmovn_u32(D_u32x4));
    // convert uint16 vectors to uint8 vectors (with saturation)
    return vcombine_u8(vqmovn_u16(vd_low_u16x8), vqmovn_u16(vd_high_u16x8));
}

template <RoundingPolicy round_policy>
inline qasymm8x16_signed_t vmlaq_qasymm8_signed(qasymm8x16_signed_t vd, float32x4_t vs, float32x4_t vo)
{
    // Convert uint8 vectors to int16 vectors
    const int8x8_t vd_low        = vget_low_s8(vd);
    const int8x8_t vd_high       = vget_high_s8(vd);
    int16x8_t      vd_low_s16x8  = vmovl_s8(vd_low);
    int16x8_t      vd_high_s16x8 = vmovl_s8(vd_high);
    // Convert int16 vectors to int32 vectors
    int32x4_t A_s32x4 = vmovl_s16(vget_low_s16(vd_low_s16x8));
    int32x4_t B_s32x4 = vmovl_s16(vget_high_s16(vd_low_s16x8));
    int32x4_t C_s32x4 = vmovl_s16(vget_low_s16(vd_high_s16x8));
    int32x4_t D_s32x4 = vmovl_s16(vget_high_s16(vd_high_s16x8));
    // Convert int32 vectors to float32 vectors
    float32x4_t A_f32x4 = vcvtq_f32_s32(A_s32x4);
    float32x4_t B_f32x4 = vcvtq_f32_s32(B_s32x4);
    float32x4_t C_f32x4 = vcvtq_f32_s32(C_s32x4);
    float32x4_t D_f32x4 = vcvtq_f32_s32(D_s32x4);
    // vd = vd*vs + vo
    A_f32x4 = vmlaq_f32(vo, A_f32x4, vs);
    B_f32x4 = vmlaq_f32(vo, B_f32x4, vs);
    C_f32x4 = vmlaq_f32(vo, C_f32x4, vs);
    D_f32x4 = vmlaq_f32(vo, D_f32x4, vs);
#if __aarch64__
    if (round_policy == RoundingPolicy::TO_NEAREST_EVEN)
    {
        A_s32x4 = vcvtnq_s32_f32(A_f32x4);
        B_s32x4 = vcvtnq_s32_f32(B_f32x4);
        C_s32x4 = vcvtnq_s32_f32(C_f32x4);
        D_s32x4 = vcvtnq_s32_f32(D_f32x4);
    }
    else if (round_policy == RoundingPolicy::TO_NEAREST_UP)
    {
        A_s32x4 = vcvtaq_s32_f32(A_f32x4);
        B_s32x4 = vcvtaq_s32_f32(B_f32x4);
        C_s32x4 = vcvtaq_s32_f32(C_f32x4);
        D_s32x4 = vcvtaq_s32_f32(D_f32x4);
    }
    else
    {
        A_s32x4 = vcvtq_s32_f32(A_f32x4);
        B_s32x4 = vcvtq_s32_f32(B_f32x4);
        C_s32x4 = vcvtq_s32_f32(C_f32x4);
        D_s32x4 = vcvtq_s32_f32(D_f32x4);
    }
#else  // #if __aarch64__
    // rounding mode only supported in aarch64
    A_s32x4 = vcvtq_s32_f32(A_f32x4);
    B_s32x4 = vcvtq_s32_f32(B_f32x4);
    C_s32x4 = vcvtq_s32_f32(C_f32x4);
    D_s32x4 = vcvtq_s32_f32(D_f32x4);
#endif // #if __aarch64__

    // Convert int32 vectors to int16 vectors (with saturation)
    vd_low_s16x8  = vcombine_s16(vqmovn_s32(A_s32x4), vqmovn_s32(B_s32x4));
    vd_high_s16x8 = vcombine_s16(vqmovn_s32(C_s32x4), vqmovn_s32(D_s32x4));
    // convert int16 vectors to int8 vectors (with saturation)
    return vcombine_s8(vqmovn_s16(vd_low_s16x8), vqmovn_s16(vd_high_s16x8));
}
} // namespace arm_compute