dynarmic/src/frontend/A64/translate/impl/simd_three_same_extra.cpp

/* This file is part of the dynarmic project.
 * Copyright (c) 2018 MerryMage
 * This software may be used and distributed according to the terms of the GNU
 * General Public License version 2 or any later version.
 */

#include "frontend/A64/translate/impl/impl.h"

namespace Dynarmic::A64 {
namespace {

using ExtensionFunction = IR::U32 (IREmitter::*)(const IR::UAny&);

bool DotProduct(TranslatorVisitor& v, bool Q, Imm<2> size, Vec Vm, Vec Vn, Vec Vd,
                ExtensionFunction extension) {
    if (size != 0b10) {
        return v.ReservedValue();
    }

    const size_t esize = 8 << size.ZeroExtend();
    const size_t datasize = Q ? 128 : 64;
    const size_t elements = datasize / esize;

    const IR::U128 operand1 = v.V(datasize, Vn);
    const IR::U128 operand2 = v.V(datasize, Vm);
    IR::U128 result = v.V(datasize, Vd);

    for (size_t i = 0; i < elements; i++) {
        IR::U32 res_element = v.ir.Imm32(0);

        for (size_t j = 0; j < 4; j++) {
            const IR::U32 elem1 = (v.ir.*extension)(v.ir.VectorGetElement(8, operand1, 4 * i + j));
            const IR::U32 elem2 = (v.ir.*extension)(v.ir.VectorGetElement(8, operand2, 4 * i + j));

            res_element = v.ir.Add(res_element, v.ir.Mul(elem1, elem2));
        }

        res_element = v.ir.Add(v.ir.VectorGetElement(32, result, i), res_element);
        result = v.ir.VectorSetElement(32, result, i, res_element);
    }

    v.V(datasize, Vd, result);
    return true;
}

} // Anonymous namespace

bool TranslatorVisitor::SDOT_vec(bool Q, Imm<2> size, Vec Vm, Vec Vn, Vec Vd) {
    return DotProduct(*this, Q, size, Vm, Vn, Vd, &IREmitter::SignExtendToWord);
}

bool TranslatorVisitor::UDOT_vec(bool Q, Imm<2> size, Vec Vm, Vec Vn, Vec Vd) {
    return DotProduct(*this, Q, size, Vm, Vn, Vd, &IREmitter::ZeroExtendToWord);
}

bool TranslatorVisitor::FCMLA_vec(bool Q, Imm<2> size, Vec Vm, Imm<2> rot, Vec Vn, Vec Vd) {
    if (size == 0) {
        return ReservedValue();
    }

    if (!Q && size == 0b11) {
        return ReservedValue();
    }

    const size_t esize = 8U << size.ZeroExtend();

    // TODO: Currently we don't support half-precision floating point
    if (esize == 16) {
        return InterpretThisInstruction();
    }

    const size_t datasize = Q ? 128 : 64;
    const size_t num_elements = datasize / esize;
    const size_t num_iterations = num_elements / 2;

    const IR::U128 operand1 = V(datasize, Vn);
    const IR::U128 operand2 = V(datasize, Vm);
    const IR::U128 operand3 = V(datasize, Vd);
    IR::U128 result = ir.ZeroVector();

    IR::U32U64 element1;
    IR::U32U64 element2;
    IR::U32U64 element3;
    IR::U32U64 element4;
    for (size_t e = 0; e < num_iterations; ++e) {
        const size_t first = e * 2;
        const size_t second = first + 1;

        switch (rot.ZeroExtend()) {
        case 0b00: // 0 degrees
            element1 = ir.VectorGetElement(esize, operand2, first);
            element2 = ir.VectorGetElement(esize, operand1, first);
            element3 = ir.VectorGetElement(esize, operand2, second);
            element4 = ir.VectorGetElement(esize, operand1, first);
            break;
        case 0b01: // 90 degrees
            element1 = ir.FPNeg(ir.VectorGetElement(esize, operand2, second));
            element2 = ir.VectorGetElement(esize, operand1, second);
            element3 = ir.VectorGetElement(esize, operand2, first);
            element4 = ir.VectorGetElement(esize, operand1, second);
            break;
        case 0b10: // 180 degrees
            element1 = ir.FPNeg(ir.VectorGetElement(esize, operand2, first));
            element2 = ir.VectorGetElement(esize, operand1, first);
            element3 = ir.FPNeg(ir.VectorGetElement(esize, operand2, second));
            element4 = ir.VectorGetElement(esize, operand1, first);
            break;
        case 0b11: // 270 degrees
            element1 = ir.VectorGetElement(esize, operand2, second);
            element2 = ir.VectorGetElement(esize, operand1, second);
            element3 = ir.FPNeg(ir.VectorGetElement(esize, operand2, first));
            element4 = ir.VectorGetElement(esize, operand1, second);
            break;
        }

        const IR::U32U64 operand3_elem1 = ir.VectorGetElement(esize, operand3, first);
        const IR::U32U64 operand3_elem2 = ir.VectorGetElement(esize, operand3, second);

        result = ir.VectorSetElement(esize, result, first,
                                     ir.FPMulAdd(operand3_elem1, element2, element1, true));
        result = ir.VectorSetElement(esize, result, second,
                                     ir.FPMulAdd(operand3_elem2, element4, element3, true));
    }

    ir.SetQ(Vd, result);
    return true;
}

bool TranslatorVisitor::FCADD_vec(bool Q, Imm<2> size, Vec Vm, Imm<1> rot, Vec Vn, Vec Vd) {
    if (size == 0) {
        return ReservedValue();
    }

    if (!Q && size == 0b11) {
        return ReservedValue();
    }

    const size_t esize = 8U << size.ZeroExtend();

    // TODO: Currently we don't support half-precision floating point
    if (esize == 16) {
        return InterpretThisInstruction();
    }

    const size_t datasize = Q ? 128 : 64;
    const size_t num_elements = datasize / esize;
    const size_t num_iterations = num_elements / 2;

    const IR::U128 operand1 = V(datasize, Vn);
    const IR::U128 operand2 = V(datasize, Vm);
    IR::U128 result = ir.ZeroVector();

    IR::U32U64 element1;
    IR::U32U64 element3;
    for (size_t e = 0; e < num_iterations; ++e) {
        const size_t first = e * 2;
        const size_t second = first + 1;

        if (rot == 0) {
            element1 = ir.FPNeg(ir.VectorGetElement(esize, operand2, second));
            element3 = ir.VectorGetElement(esize, operand2, first);
        } else if (rot == 1) {
            element1 = ir.VectorGetElement(esize, operand2, second);
            element3 = ir.FPNeg(ir.VectorGetElement(esize, operand2, first));
        }

        const IR::U32U64 operand1_elem1 = ir.VectorGetElement(esize, operand1, first);
        const IR::U32U64 operand1_elem3 = ir.VectorGetElement(esize, operand1, second);

        result = ir.VectorSetElement(esize, result, first,
                                     ir.FPAdd(operand1_elem1, element1, true));
        result = ir.VectorSetElement(esize, result, second,
                                     ir.FPAdd(operand1_elem3, element3, true));
    }

    ir.SetQ(Vd, result);
    return true;
}

} // namespace Dynarmic::A64