/*
   Copyright (C) 2005, 2006, 2007, 2008, 2009, 2010, 2011 Her Majesty
   the Queen in Right of Canada (Communications Research Center Canada)

   Copyright (C) 2023
   Matthias P. Braendli, matthias.braendli@mpb.li

    http://opendigitalradio.org

    This flowgraph block converts complexf to signed integer.
 */
/*
   This file is part of ODR-DabMod.

   ODR-DabMod is free software: you can redistribute it and/or modify
   it under the terms of the GNU General Public License as
   published by the Free Software Foundation, either version 3 of the
   License, or (at your option) any later version.

   ODR-DabMod is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU General Public License for more details.

   You should have received a copy of the GNU General Public License
   along with ODR-DabMod.  If not, see <http://www.gnu.org/licenses/>.
 */

#include "FormatConverter.h"
#include "PcDebug.h"

#include <sys/types.h>
#include <string.h>
#include <stdexcept>
#include <assert.h>

FormatConverter::FormatConverter(const std::string& format) :
    ModCodec(),
    m_format(format)
{ }

FormatConverter::~FormatConverter()
{
    etiLog.level(debug) << "FormatConverter: "
        << m_num_clipped_samples.load() <<
        " clipped samples";
}


/* Expect the input samples to be in the correct range for the required format */
int FormatConverter::process(Buffer* const dataIn, Buffer* dataOut)
{
    PDEBUG("FormatConverter::process(dataIn: %p, dataOut: %p)\n",
            dataIn, dataOut);

    size_t num_clipped_samples = 0;

    bool source_is_complexf = m_format != "fixedpoint";

    if (source_is_complexf) {
        size_t sizeIn = dataIn->getLength() / sizeof(float);
        float* in = reinterpret_cast<float*>(dataIn->getData());

        if (m_format == "s16") {
            dataOut->setLength(sizeIn * sizeof(int16_t));
            int16_t* out = reinterpret_cast<int16_t*>(dataOut->getData());

            for (size_t i = 0; i < sizeIn; i++) {
                if (in[i] < INT16_MIN) {
                    out[i] = INT16_MIN;
                    num_clipped_samples++;
                }
                else if (in[i] > INT16_MAX) {
                    out[i] = INT16_MAX;
                    num_clipped_samples++;
                }
                else {
                    out[i] = in[i];
                }
            }
        }
        else if (m_format == "u8") {
            dataOut->setLength(sizeIn * sizeof(int8_t));
            uint8_t* out = reinterpret_cast<uint8_t*>(dataOut->getData());

            for (size_t i = 0; i < sizeIn; i++) {
                const auto samp = in[i] + 128.0f;
                if (samp < 0) {
                    out[i] = 0;
                    num_clipped_samples++;
                }
                else if (samp > UINT8_MAX) {
                    out[i] = UINT8_MAX;
                    num_clipped_samples++;
                }
                else {
                    out[i] = samp;
                }

            }
        }
        else if (m_format == "s8") {
            dataOut->setLength(sizeIn * sizeof(int8_t));
            int8_t* out = reinterpret_cast<int8_t*>(dataOut->getData());

            for (size_t i = 0; i < sizeIn; i++) {
                if (in[i] < INT8_MIN) {
                    out[i] = INT8_MIN;
                    num_clipped_samples++;
                }
                else if (in[i] > INT8_MAX) {
                    out[i] = INT8_MAX;
                    num_clipped_samples++;
                }
                else {
                    out[i] = in[i];
                }
            }
        }
        else {
            throw std::runtime_error("FormatConverter: Invalid format " + m_format);
        }

    }
    else {
        // Output is always sc16, because that's what UHD accepts

        using fixed_t = complexfix::value_type;
        size_t sizeIn = dataIn->getLength() / sizeof(fixed_t);
        fixed_t* in = reinterpret_cast<fixed_t*>(dataIn->getData());

        dataOut->setLength(sizeIn * sizeof(int16_t));
        int16_t* out = reinterpret_cast<int16_t*>(dataOut->getData());

        for (size_t i = 0; i < sizeIn; i++) {
            const auto v = (in[i] * 2).raw_value();

            if (v < INT16_MIN) {
                out[i] = INT16_MIN;
                num_clipped_samples++;
            }
            else if (v > INT16_MAX) {
                out[i] = INT16_MAX;
                num_clipped_samples++;
            }
            else {
                out[i] = (int16_t)v;
            }
        }
    }

    m_num_clipped_samples.store(num_clipped_samples);
    return dataOut->getLength();
}

const char* FormatConverter::name()
{
    return "FormatConverter";
}

size_t FormatConverter::get_num_clipped_samples() const
{
    return m_num_clipped_samples.load();
}


size_t FormatConverter::get_format_size(const std::string& format)
{
    // Returns 2*sizeof(SAMPLE_TYPE) because we have I + Q
    if (format == "fixedpoint") {
        return 4;
    }
    else if (format == "s16") {
        return 4;
    }
    else if (format == "u8") {
        return 2;
    }
    else if (format == "s8") {
        return 2;
    }
    else {
        throw std::runtime_error("FormatConverter: Invalid format " + format);
    }
}