IIR.cpp

Go to the documentation of this file.

#include "IIR.hpp"
 
namespace MayaFlux::Nodes::Filters {
 
IIR::IIR(const std::shared_ptr<Node>& input, const std::vector<double>& a_coef, const std::vector<double>& b_coef)
    : Filter(input, a_coef, b_coef)
{
}
 
IIR::IIR(const std::vector<double>& a_coef, const std::vector<double>& b_coef)
    : Filter(a_coef, b_coef)
{
}
 
double IIR::process_sample(double input)
{
    if (is_bypass_enabled()) {
        return input;
    }
 
    double processed_input = input;
    if (m_input_node) {
        atomic_inc_modulator_count(m_input_node->m_modulator_count, 1);
        uint32_t state = m_input_node->m_state.load();
        if (state & NodeState::PROCESSED) {
            processed_input += m_input_node->get_last_output();
        } else {
            processed_input += m_input_node->process_sample(input);
            atomic_add_flag(m_input_node->m_state, NodeState::PROCESSED);
        }
    }
 
    if (m_use_external_input_context) {
        build_input_history(processed_input);
    } else {
        update_inputs(processed_input);
    }
 
    double output = 0.;
    const size_t num_feedforward = std::min(m_coef_b.size(), m_input_history.size());
    for (size_t i = 0; i < num_feedforward; ++i) {
        output += m_coef_b[i] * m_input_history[i];
    }
 
    const size_t num_feedback = std::min(m_coef_a.size(), m_output_history.size());
    for (size_t i = 1; i < num_feedback; ++i) {
        output -= m_coef_a[i] * m_output_history[i];
    }
 
    // output /= m_coef_a[0];
 
    update_outputs(output);
 
    if ((!m_state_saved || (m_state_saved && m_fire_events_during_snapshot))
        && !m_networked_node) {
        notify_tick(output);
    }
 
    if (m_input_node) {
        atomic_dec_modulator_count(m_input_node->m_modulator_count, 1);
        try_reset_processed_state(m_input_node);
    }
 
    return output * get_gain();
}
 
void IIR::save_state()
{
    m_saved_input_history = m_input_history;
    m_saved_output_history = m_output_history;
 
    if (m_input_node)
        m_input_node->save_state();
 
    m_state_saved = true;
}
 
void IIR::restore_state()
{
    m_input_history = m_saved_input_history;
    m_output_history = m_saved_output_history;
 
    if (m_input_node)
        m_input_node->restore_state();
 
    m_state_saved = false;
}
 
}