pcs.h

#pragma once

#include "core/filters_common.h"

#include "utilities/aligned/vector.h"

namespace wayverb {
namespace waveguide {


template <typename T>
constexpr auto factdbl(T t) {
    double out{1};
    for (auto i = t; i >= 1; i -= 2) {
        out *= i;
    }
    return out;
}

struct offset_signal final {
    util::aligned::vector<double> signal;
    size_t offset;
};

offset_signal maxflat(double f0, uint32_t N, double A, uint32_t h_len);

double compute_g0(double acoustic_impedance,
                  double speed_of_sound,
                  double sample_rate,
                  double radius);

core::filter::biquad::coefficients mech_sphere(double M,
                                               double f0,
                                               double Q,
                                               double T);

offset_signal design_pcs_source(size_t length,
                                double acoustic_impedance,
                                double speed_of_sound,
                                double sample_rate,
                                double radius,
                                double sphere_mass,
                                double low_cutoff_hz,
                                double low_q);

}  // namespace waveguide
}  // namespace wayverb