// Copyright (C) 2025 - zsliu98
// This file is part of ZLSplitter
//
// ZLSplitter is free software: you can redistribute it and/or modify it under the terms of the GNU Affero General Public License Version 3 as published by the Free Software Foundation.
//
// ZLSplitter 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 Affero General Public License for more details.
//
// You should have received a copy of the GNU Affero General Public License along with ZLSplitter. If not, see <https://www.gnu.org/licenses/>.

#pragma once
#pragma clang diagnostic push
#pragma clang diagnostic ignored "-Wall"
#pragma clang diagnostic ignored "-Weverything"
#include <kfr/kfr.h>
#include <kfr/dsp.hpp>
#include <kfr/dft.hpp>
#pragma clang diagnostic pop

namespace zldsp::fft {
    template<typename FloatType>
    void fillCycleHanningWindow(kfr::univector<FloatType> &window, const size_t size) {
        kfr::univector<float> temp_window;
        temp_window.resize(size + 1);
        temp_window = kfr::window_hann<FloatType>(size + 1);
        auto actual_window = kfr::make_univector(temp_window.data(), size);
        window = actual_window;
    }

    template<typename FloatType>
    class KFREngine {
    public:
        KFREngine() = default;

        void setOrder(const size_t order) {
            fft_size_ = static_cast<size_t>(1) << order;
            fft_plan_ = std::make_unique<kfr::dft_plan_real<FloatType> >(fft_size_);
            temp_buffer_.resize(fft_plan_->temp_size);
        }

        void forward(FloatType *in_buffer, std::complex<FloatType> *out_buffer) {
            fft_plan_->execute(out_buffer, in_buffer, temp_buffer_.data());
        }

        void forward(FloatType *in_buffer, FloatType *float_out_buffer) {
            auto out_buffer = reinterpret_cast<std::complex<FloatType> *>(float_out_buffer);
            forward(in_buffer, out_buffer);
        }

        void backward(std::complex<FloatType> *out_buffer, FloatType *in_buffer) {
            fft_plan_->execute(in_buffer, out_buffer, temp_buffer_.data());
        }

        void backward(FloatType *float_out_buffer, FloatType *in_buffer) {
            auto out_buffer = reinterpret_cast<std::complex<FloatType> *>(float_out_buffer);
            backward(out_buffer, in_buffer);
        }

        void forwardMagnitudeOnly(FloatType *buffer) {
            forward(buffer, buffer);
            auto *out = reinterpret_cast<std::complex<FloatType> *>(buffer);
            for (size_t i = 0; i < (fft_size_ / 2) + 1; ++i) {
                buffer[i] = std::abs(out[i]);
            }
        }

        [[nodiscard]] size_t getSize() const { return fft_size_; }

    private:
        size_t fft_size_{0};
        std::unique_ptr<kfr::dft_plan_real<FloatType> > fft_plan_;
        kfr::univector<kfr::u8> temp_buffer_;
    };
}
