Function checkQComplement

test/QAlg/ArithmeticUnit_V2.test.cpp:168–232

// use DDT test Complement

void checkQComplement(int bit_len, int value)
{
    auto qvm = initQuantumMachine(QMachineType::CPU);

    QVec a = qvm->qAllocMany(bit_len);
    int ancil_capacity = int(2 * bit_len + 1 - std::floor(std::log2(bit_len)));
    QVec ancil = qvm->qAllocMany(ancil_capacity);

    // qcla complement
    QProg qcla_qprog = createEmptyQProg();
    if (0 != value)
    {
        qcla_qprog << bind_data(value, a);
    }

    qcla_qprog << QComplement_V2(a, ancil);

    auto result_qcla = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qcla_qprog, a);
    auto result_qcla_anc = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qcla_qprog, ancil);

    // qft complement
    QProg qft_qprog = createEmptyQProg();
    if (0 != value)
    {
        qft_qprog << bind_data(value, a);
    }

    qft_qprog << QComplement_V2(a, ancil[0]);

    auto result_qft = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qft_qprog, a);
    auto result_qft_anc = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qft_qprog, ancil);

    // original complement
    QProg qprog = createEmptyQProg();
    if (0 != value)
    {
        qprog << bind_data(value, a);
    }
    qprog << QComplement(a, ancil);

    auto result = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qprog, a);
    auto result_anc = dynamic_cast<CPUQVM *>(qvm)->probRunDict(qprog, ancil);

    // ASSERT_EQ(result.size(), result_qcla.size());

    for (auto it : result)
    {
        double total = it.second + result_qcla.at(it.first) + result_qft.at(it.first);
        if (total <= 0.01)
            continue;
        ASSERT_TRUE(euqal_in_tolerance(it.second, result_qcla.at(it.first)));
        ASSERT_TRUE(euqal_in_tolerance(it.second, result_qft.at(it.first)));
    }

    for (auto it : result_anc)
    {
        double total = it.second + result_qcla_anc.at(it.first) + result_qft_anc.at(it.first);
        if (total <= 0.01)