Method _decode_siho

src/Module/Decoder/RA/Decoder_RA.cpp:81–148

template<typename B, typename R>
int
Decoder_RA<B, R>::_decode_siho(const R* Y_N, B* V_K, const size_t frame_id)
{
    //	auto t_load = std::chrono::steady_clock::now(); // -----------------------------------------------------------
    // LOAD set F, B and Td at 0
    for (auto i = 0; i < this->N; i++)
    {
        Fw[i] = 0;
        Bw[i] = 0;
        Td[i] = 0;
    }
    //	auto d_load = std::chrono::steady_clock::now() - t_load;

    //	auto t_decod = std::chrono::steady_clock::now(); // --------------------------------------------------------
    // DECODE
    for (auto iter = 0; iter < max_iter; iter++)
    {
        ///////////////////
        // SISO decoding //
        ///////////////////

        // Forward
        Fw[0] = Td[0];
        for (auto i = 1; i < this->N; i++)
            Fw[i] = check_node(Fw[i - 1] + Y_N[i - 1], Td[i]);

        // Backward
        Bw[this->N - 2] = check_node(Y_N[this->N - 1], Td[this->N - 1]);
        for (auto i = this->N - 3; i >= 0; i--)
            Bw[i] = check_node(Bw[i + 1] + Y_N[i + 1], Td[i + 1]);

        // Extrinsic
        Tu[0] = Bw[0] + Y_N[0];
        Tu[this->N - 1] = check_node(Fw[this->N - 2] + Y_N[this->N - 2], Y_N[this->N - 1]);
        for (auto i = 1; i < this->N - 1; i++)
            Tu[i] = check_node(Fw[i - 1] + Y_N[i - 1], Bw[i] + Y_N[i]);

        // Deinterleave
        interleaver->deinterleave(Tu.data(), Wu.data(), frame_id, false);

        // U computation
        R tmp;
        for (auto i = 0; i < this->K; i++)
        {
            tmp = 0;
            for (auto j = 0; j < rep_count; j++)
                tmp += Wu[i * rep_count + j];
            for (auto j = 0; j < rep_count; j++)
                Wd[i * rep_count + j] = tmp - Wu[i * rep_count + j];
            U[i] = tmp;
        }

        // Interleaving
        interleaver->interleave(Wd.data(), Td.data(), frame_id, false);
    }
    //	auto d_decod = std::chrono::steady_clock::now() - t_decod;

    //	auto t_store = std::chrono::steady_clock::now(); // ---------------------------------------------------------