/*
 * Copyright (C) 2025 Linux Studio Plugins Project <https://lsp-plug.in/>
 *           (C) 2025 Vladimir Sadovnikov <sadko4u@gmail.com>
 *
 * This file is part of lsp-plugins
 * Created on: 03 авг. 2016 г.
 *
 * lsp-plugins is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * any later version.
 *
 * lsp-plugins 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 Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public License
 * along with lsp-plugins. If not, see <https://www.gnu.org/licenses/>.
 */

#include <lsp-plug.in/common/alloc.h>
#include <lsp-plug.in/common/debug.h>
#include <lsp-plug.in/dsp-units/util/Crossover.h>
#include <lsp-plug.in/dsp-units/units.h>

namespace lsp
{
    namespace dspu
    {
        Crossover::Crossover()
        {
            construct();
        }

        Crossover::~Crossover()
        {
            destroy();
        }

        void Crossover::construct()
        {
            nReconfigure    = R_ALL;
            nSplits         = 0;
            nBufSize        = 0;
            nSampleRate     = LSP_DSP_UNITS_DEFAULT_SAMPLE_RATE;
            nPlanSize       = 0;

            vBands          = NULL;
            vSplit          = NULL;
            vPlan           = NULL;

            vLpfBuf         = NULL;
            vHpfBuf         = NULL;
            pData           = NULL;
        }

        void Crossover::destroy()
        {
            if (vSplit != NULL)
            {
                for (size_t i=0; i<nSplits; ++i)
                {
                    split_t *sp     = &vSplit[i];
                    sp->sLPF.destroy();
                    sp->sHPF.destroy();
                }
            }

            free_aligned(pData);
            construct();
        }

        bool Crossover::init(size_t bands, size_t buf_size)
        {
            if (bands < 1)
                return false;

            size_t splits       = bands - 1;
            size_t xbuf_size    = align_size(buf_size * sizeof(float), DEFAULT_ALIGN);
            size_t band_size    = align_size(bands * sizeof(band_t), DEFAULT_ALIGN);
            size_t split_size   = align_size(splits * sizeof(split_t), DEFAULT_ALIGN);
            size_t plan_size    = align_size(splits * sizeof(split_t *), DEFAULT_ALIGN);
            size_t to_alloc     = band_size +
                                  split_size +
                                  plan_size +
                                  xbuf_size * 2;

            // Allocate buffers
            uint8_t *data       = NULL;
            uint8_t *ptr        = alloc_aligned<uint8_t>(data, to_alloc);
            if (ptr == NULL)
                return false;

            // Distribute the allocated space
            vBands              = advance_ptr_bytes<band_t>(ptr, band_size);
            vSplit              = advance_ptr_bytes<split_t>(ptr, split_size);
            vPlan               = advance_ptr_bytes<split_t *>(ptr, plan_size);
            vLpfBuf             = advance_ptr_bytes<float>(ptr, xbuf_size);
            vHpfBuf             = advance_ptr_bytes<float>(ptr, xbuf_size);

            // Initialize fields, keep sample_rate unchanged
            nReconfigure        = R_ALL;
            nSplits             = uint32_t(splits);
            nBufSize            = uint32_t(buf_size);
            nPlanSize           = 0;

            // Store allocated data pointer
            pData               = data;

            // Construct all splits
            float step          = logf(LSP_DSP_UNITS_SPEC_FREQ_MAX / LSP_DSP_UNITS_SPEC_FREQ_MIN) / bands;

            for (size_t i=0; i<nSplits; ++i)
            {
                split_t *sp         = &vSplit[i];

                // Initialize filters
                sp->sLPF.construct();
                sp->sHPF.construct();

                if (!sp->sLPF.init(splits, 0))
                {
                    destroy();
                    return false;
                }
                sp->sLPF.set_sample_rate(nSampleRate);

                if (!sp->sHPF.init(NULL))
                {
                    destroy();
                    return false;
                }
                sp->sHPF.set_sample_rate(nSampleRate);

                // Set IIR mode for each filter
                sp->sLPF.set_mode(EQM_IIR);

                // Initialize split point parameters
                sp->nBandId         = i + 1; // Band N+1 is attached to split point N
                sp->nSlope          = 0;
                sp->fFreq           = LSP_DSP_UNITS_SPEC_FREQ_MIN * expf((i+1) * step);
                sp->nMode           = CROSS_MODE_BT;
            }

            // Construct all bands
            for (size_t i=0; i<=nSplits; ++i)
            {
                band_t *sb          = &vBands[i];

                sb->fGain           = GAIN_AMP_0_DB;
                sb->fStart          = (i == 0) ? LSP_DSP_UNITS_SPEC_FREQ_MIN : vSplit[i-1].fFreq;
                sb->fEnd            = (i < nSplits) ? vSplit[i].fFreq : nSampleRate >> 1;
                sb->bEnabled        = false;
                sb->pStart          = NULL;
                sb->pEnd            = NULL;

                sb->pFunc           = NULL;
                sb->pObject         = NULL;
                sb->pSubject        = NULL;
                sb->nId             = i;
            }

            return true;
        }

        filter_type_t Crossover::select_filter(xover_type_t type, crossover_mode_t mode, size_t slope)
        {
            if (slope == CROSS_SLOPE_LR2)
            {
                switch (type)
                {
                    case FILTER_LPF: return (mode == CROSS_MODE_BT) ? FLT_BT_RLC_LOPASS  : FLT_MT_RLC_LOPASS;
                    case FILTER_HPF: return (mode == CROSS_MODE_BT) ? FLT_BT_RLC_HIPASS  : FLT_MT_RLC_HIPASS;
                    case FILTER_APF: return (mode == CROSS_MODE_BT) ? FLT_BT_RLC_ALLPASS : FLT_MT_RLC_ALLPASS;
                    default:
                        return FLT_NONE;
                }
            }

            switch (type)
            {
                case FILTER_LPF: return (mode == CROSS_MODE_BT) ? FLT_BT_LRX_LOPASS  : FLT_MT_LRX_LOPASS;
                case FILTER_HPF: return (mode == CROSS_MODE_BT) ? FLT_BT_LRX_HIPASS  : FLT_MT_LRX_HIPASS;
                case FILTER_APF: return (mode == CROSS_MODE_BT) ? FLT_BT_LRX_ALLPASS : FLT_MT_LRX_ALLPASS;
                default:
                    return FLT_NONE;
            }
        }

        uint32_t Crossover::select_slope(xover_type_t type, size_t slope)
        {
            if (slope == CROSS_SLOPE_LR2)
                return (type == FILTER_APF) ? 1 : 2;
            return uint32_t(slope - 1);
        }

        void Crossover::set_slope(size_t sp, size_t slope)
        {
            if (sp >= nSplits)
                return;
            if (slope == vSplit[sp].nSlope)
                return;

            vSplit[sp].nSlope   = slope;
            nReconfigure       |= R_SPLIT;
        }

        ssize_t Crossover::get_slope(size_t sp) const
        {
            return (sp < nSplits) ? vSplit[sp].nSlope : -1;
        }

        void Crossover::set_frequency(size_t sp, float freq)
        {
            if (sp >= nSplits)
                return;
            if (freq == vSplit[sp].fFreq)
                return;

            vSplit[sp].fFreq    = freq;
            nReconfigure       |= R_SPLIT;
        }

        float Crossover::get_frequency(size_t sp) const
        {
            return (sp < nSplits) ? vSplit[sp].fFreq : -1.0f;
        }

        void Crossover::set_mode(size_t sp, crossover_mode_t mode)
        {
            if (sp >= nSplits)
                return;
            if (mode == vSplit[sp].nMode)
                return;

            vSplit[sp].nMode    = mode;
            nReconfigure       |= R_SPLIT;
        }

        ssize_t Crossover::get_mode(size_t sp) const
        {
            return (sp < nSplits) ? vSplit[sp].nMode : -1.0f;
        }

        void Crossover::set_gain(size_t band, float gain)
        {
            if (band > nSplits)
                return;
            if (gain == vBands[band].fGain)
                return;

            vBands[band].fGain = gain;
            nReconfigure       |= R_GAIN;
        }

        float Crossover::get_gain(size_t band) const
        {
            return (band <= nSplits) ? vBands[band].fGain: -1.0f;
        }

        float Crossover::get_band_start(size_t band)
        {
            reconfigure();
            return (band <= nSplits) ? vBands[band].fStart : -1.0f;
        }

        float Crossover::get_band_end(size_t band)
        {
            reconfigure();
            return (band <= nSplits) ? vBands[band].fEnd : -1.0f;
        }

        bool Crossover::set_handler(size_t band, crossover_func_t func, void *object, void *subject)
        {
            if (band > nSplits)
                return false;

            band_t *b       = &vBands[band];
            b->pFunc        = func;
            b->pObject      = object;
            b->pSubject     = subject;

            return true;
        }

        bool Crossover::unset_handler(size_t band)
        {
            if (band > nSplits)
                return false;

            band_t *b       = &vBands[band];
            b->pFunc        = NULL;
            b->pObject      = NULL;
            b->pSubject     = NULL;
            return true;
        }

        bool Crossover::band_active(size_t band)
        {
            if (band > nSplits)
                return false;
            else if (band == 0)
                return true;

            reconfigure();
            return vBands[band].bEnabled;
        }

        void Crossover::set_sample_rate(size_t sr)
        {
            if (nSampleRate == sr)
                return;

            nSampleRate     = uint32_t(sr);
            for (size_t i=0; i<nSplits; ++i)
            {
                vSplit[i].sLPF.set_sample_rate(sr);
                vSplit[i].sHPF.set_sample_rate(sr);
            }
            vBands[nSplits].fEnd = sr >> 1;

            nReconfigure   |= R_ALL;
        }

        void Crossover::reconfigure()
        {
            if (!nReconfigure)
                return;

            // Form the plan and reset band state
            nPlanSize       = 0;
            for (size_t i=0; i<nSplits; ++i)
            {
                if (vSplit[i].nSlope != CROSS_SLOPE_OFF)
                    vPlan[nPlanSize++]  = &vSplit[i];
            }
            for (size_t i=0; i<=nSplits; ++i)
                vBands[i].bEnabled  = false;

            // Sort split bands in ascending order
            for (ssize_t si=0, n=nPlanSize; si < n-1; ++si)
                for (ssize_t sj=si+1; sj < n; ++sj)
                    if (vPlan[sj]->fFreq < vPlan[si]->fFreq)
                        swap(vPlan[si], vPlan[sj]);

            band_t *left        = &vBands[0];
            left->fStart        = LSP_DSP_UNITS_SPEC_FREQ_MIN;
            left->bEnabled      = true;
            left->pStart        = NULL;

            // Configure LPF and HPF bands
            for (size_t i=0; i<nPlanSize; ++i)
            {
                split_t *sp         = vPlan[i];
                band_t *right       = &vBands[sp->nBandId];

                left->fEnd          = sp->fFreq;
                left->pEnd          = sp;
                right->fStart       = sp->fFreq;
                right->pStart       = sp;
                right->bEnabled     = true;

                // Set LPF parameters
                size_t filter_id    = 0;
                filter_params_t fp;

                fp.nType            = select_filter(FILTER_LPF, sp->nMode, sp->nSlope);
                fp.fFreq            = sp->fFreq;
                fp.fFreq2           = sp->fFreq;
                fp.fGain            = left->fGain;
                fp.nSlope           = select_slope(FILTER_LPF, sp->nSlope);
                fp.fQuality         = 0.0f;

                sp->sLPF.set_params(filter_id++, &fp);

                // Append APF filters
                for (size_t j=i+1; j<nPlanSize; ++j)
                {
                    split_t *xsp        = vPlan[j];

                    fp.nType            = select_filter(FILTER_APF, xsp->nMode, xsp->nSlope);
                    fp.fFreq            = xsp->fFreq;
                    fp.fFreq2           = xsp->fFreq;
                    fp.fGain            = GAIN_AMP_0_DB;
                    fp.nSlope           = select_slope(FILTER_APF, xsp->nSlope);
                    fp.fQuality         = 0.0f;

                    sp->sLPF.set_params(filter_id++, &fp);
                }

                // Disable all other filters in the chain
                while (filter_id < nSplits)
                {
                    fp.nType            = FLT_NONE;
                    fp.fFreq            = 0;
                    fp.fFreq2           = 0;
                    fp.fGain            = GAIN_AMP_0_DB;
                    fp.nSlope           = 0;
                    fp.fQuality         = 0.0f;

                    sp->sLPF.set_params(filter_id++, &fp);
                }

                // Set HPF parameters
                fp.nType            = select_filter(FILTER_HPF, sp->nMode, sp->nSlope);
                fp.fFreq            = sp->fFreq;
                fp.fFreq2           = sp->fFreq;
                fp.fGain            = (i < (nPlanSize-1)) ? GAIN_AMP_0_DB : right->fGain;
                if (sp->nSlope == CROSS_SLOPE_LR2)
                    fp.fGain            = -fp.fGain;
                fp.nSlope           = select_slope(FILTER_HPF, sp->nSlope);
                fp.fQuality         = 0.0f;

                sp->sHPF.update(nSampleRate, &fp);
                sp->sHPF.rebuild();

                // Move to next band
                left                = right;
            }

            // Update frequency of the last band
            left->fEnd          = nSampleRate * 0.5f;
            left->pEnd          = NULL;

            // DEBUG BEGIN
        #ifdef LSP_TRACE
            lsp_trace("Execution plan:");
            for (size_t i=0; i<nPlanSize; ++i)
            {
                split_t *sp         = vPlan[i];
                lsp_trace("  split point #%d: this=%p, band=%d, freq=%.2f, slope=%d",
                    int(i), sp, int(sp->nBandId), sp->fFreq, int(sp->nSlope)
                );
            }
            lsp_trace("Bands:");
            for (size_t i=0; i<=nSplits; ++i)
            {
                band_t *b           = &vBands[i];
                lsp_trace("  band #%d: this=%p, enabled=%s, gain=%f, start=%.2f, end=%.2f, start=%p, end=%p",
                    int(i), b, (b->bEnabled) ? "true " : "false",
                    b->fGain, b->fStart, b->fEnd, b->pStart, b->pEnd);
            }
        #endif
            // DEBUG END

            // Reset reconfiguration flag
            nReconfigure        = 0;
        }

        void Crossover::process(const float *in, size_t samples)
        {
            reconfigure();

            for (size_t sample=0; sample < samples; )
            {
                size_t to_do        = lsp_min(samples - sample, nBufSize);
                band_t *left        = &vBands[0];
                const float *src    = in;

                if (nPlanSize > 0)
                {
                    // Process each band except last
                    for (size_t i=0; i<nPlanSize; ++i)
                    {
                        split_t *sp         = vPlan[i];
                        band_t *right       = &vBands[sp->nBandId];

                        // Perform split first
                        if (left->pFunc != NULL)
                            sp->sLPF.process(vLpfBuf, src, to_do);
                        sp->sHPF.process(vHpfBuf, src, to_do);

                        // Now call handlers
                        if (left->pFunc != NULL)
                            left->pFunc(left->pObject, left->pSubject, left->nId, vLpfBuf, sample, to_do);

                        src                 = vHpfBuf;
                        left                = right;
                    }

                    // Process last band
                    if (left->pFunc != NULL)
                        left->pFunc(left->pObject, left->pSubject, left->nId, vHpfBuf, sample, to_do);
                }
                else if (left->pFunc != NULL)
                {
                    dsp::mul_k3(vLpfBuf, src, vBands[0].fGain, to_do);
                    left->pFunc(left->pObject, left->pSubject, left->nId, vLpfBuf, sample, to_do);
                }

                // Update pointers
                in                 += to_do;
                sample             -= to_do;
            }
        }

        bool Crossover::freq_chart(size_t band, float *re, float *im, const float *f, size_t count)
        {
            // Valid index of the band?
            if (band > nSplits)
                return false;

            // Reconfigure
            reconfigure();

            // Band is enabled ?
            band_t *b       = &vBands[band];
            if (!b->bEnabled)
            {
                dsp::fill_zero(re, count);
                dsp::fill_zero(im, count);
            }
            else if (nPlanSize == 0)
            {
                dsp::fill_one(re, count);
                dsp::fill_zero(im, count);
            }
            else if (b->pEnd == NULL)
                b->pStart->sHPF.freq_chart(re, im, f, count);
            else if (b->pStart == NULL)
                b->pEnd->sLPF.freq_chart(re, im, f, count);
            else
            {
                // Compute frequency chart with chunks of maximum nBufSize size
                while (count > 0)
                {
                    size_t to_do    = lsp_min(count, nBufSize);

                    // Apply frequency chart
                    b->pStart->sHPF.freq_chart(re, im, f, to_do);
                    b->pEnd->sLPF.freq_chart(size_t(0), vLpfBuf, vHpfBuf, f, to_do);
                    dsp::complex_mul2(re, im, vLpfBuf, vHpfBuf, to_do);

                    // Update pointers
                    re             += to_do;
                    im             += to_do;
                    f              += to_do;
                    count          -= to_do;
                }
            }

            return true;
        }

        bool Crossover::freq_chart(size_t band, float *c, const float *f, size_t count)
        {
            // Valid index of the band?
            if (band > nSplits)
                return false;

            // Reconfigure
            reconfigure();

            // Band is enabled ?
            band_t *b       = &vBands[band];
            if (!b->bEnabled)
                dsp::pcomplex_fill_ri(c, 0.0f, 0.0f, count);
            else if (nPlanSize == 0)
                dsp::pcomplex_fill_ri(c, 1.0f, 0.0f, count);
            else if (b->pEnd == NULL)
                b->pStart->sHPF.freq_chart(c, f, count);
            else if (b->pStart == NULL)
                b->pEnd->sLPF.freq_chart(c, f, count);
            else
            {
                // Compute frequency chart with chunks of maximum nBufSize size
                while (count > 0)
                {
                    // We can go out of vLpfBuf because vHpfBuf is there after it
                    size_t to_do    = lsp_min(count, nBufSize);

                    // Apply frequency chart
                    b->pStart->sHPF.freq_chart(c, f, to_do);
                    b->pEnd->sLPF.freq_chart(size_t(0), vLpfBuf, f, to_do);
                    dsp::pcomplex_mul2(c, vLpfBuf, to_do);

                    // Update pointers
                    c              += to_do * 2;
                    f              += to_do;
                    count          -= to_do;
                }
            }

            return true;
        }

        void Crossover::dump(IStateDumper *v) const
        {
            v->write("nReconfigure", nReconfigure);
            v->write("nSplits", nSplits);
            v->write("nBufSize", nBufSize);
            v->write("nSampleRate", nSampleRate);
            v->write("nPlanSize", nPlanSize);

            v->begin_array("vBands", vBands, nSplits+1);
            for (size_t i=0; i<=nSplits; ++i)
            {
                band_t *b   = &vBands[i];
                v->begin_object(b, sizeof(band_t));
                {
                    v->write("fGain", b->fGain);
                    v->write("fStart", b->fStart);
                    v->write("fEnd", b->fEnd);
                    v->write("bEnabled", b->bEnabled);
                    v->write("pStart", b->pStart);
                    v->write("pEnd", b->pEnd);

                    v->write("pFunc", b->pFunc);
                    v->write("pOpbject", b->pObject);
                    v->write("pSubject", b->pSubject);
                    v->write("nId", b->nId);
                }
                v->end_object();
            }
            v->end_array();

            v->begin_array("vSplit", vSplit, nSplits);
            for (size_t i=0; i < nSplits; ++i)
            {
                split_t *s  = &vSplit[i];
                v->begin_object(s, sizeof(split_t));
                {
                    v->write_object("sLPF", &s->sLPF);
                    v->write_object("sHPF", &s->sHPF);

                    v->write("nBandId", s->nBandId);
                    v->write("nSlopw", s->nSlope);
                    v->write("fFreq", s->fFreq);
                    v->write("nMode", s->nMode);
                }
                v->end_object();
            }
            v->end_array();

            v->writev("vPlan", vPlan, nPlanSize);

            v->write("vLpfBuf", vLpfBuf);
            v->write("vHpfBuf", vHpfBuf);
            v->write("pData", pData);
        }
    } /* namespace dspu */
} /* namespace lsp */
