/Users/kfitz/Projects/Loris/Loris development/src/Filter.h

#ifndef INCLUDE_FILTER_H
#define INCLUDE_FILTER_H
/*
 * This is the Loris C++ Class Library, implementing analysis, 
 * manipulation, and synthesis of digitized sounds using the Reassigned 
 * Bandwidth-Enhanced Additive Sound Model.
 *
 * Loris is Copyright (c) 1999-2010 by Kelly Fitz and Lippold Haken
 *
 * This program 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 2 of the License, or
 * (at your option) any later version.
 *
 * This program 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 this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 *
 * Filter.h
 *
 * Definition of class Loris::Filter, a generic digital filter of 
 * arbitrary order having both feed-forward and feedback coefficients. 
 *
 * Kelly Fitz, 1 Sept 1999
 * revised 9 Oct 2009
 * loris@cerlsoundgroup.org
 *
 * http://www.cerlsoundgroup.org/Loris/
 *
 */

#include "LorisExceptions.h"
#include "Notifier.h"

#include <algorithm>
#include <deque>
#include <vector>

//  begin namespace
namespace Loris {

// ---------------------------------------------------------------------------
//  class Filter
//
//
class Filter
{
public:

//  --- lifecycle ---

    //  default construction
    Filter( void );
    
    //  initialized construction
    //
#if !defined(NO_TEMPLATE_MEMBERS)
    template<typename IterT1, typename IterT2>
    Filter( IterT1 ffwdbegin, IterT1 ffwdend,   //  feed-forward coeffs
            IterT2 fbackbegin, IterT2 fbackend, //  feedback coeffs
            double gain = 1. );
#else
    Filter( const double * ffwdbegin, const double * ffwdend, //    feed-forward coeffs
            const double * fbackbegin, const double * fbackend, //  feedback coeffs
            double gain = 1. );
#endif

    
    // copy constructor 
    Filter( const Filter & other );
    
    //  assignment operator
    Filter & operator=( const Filter & rhs );
    
    ~Filter( void );
    

//  --- filtering ---

    double apply( double input );

    double operator() ( double input ) { return apply(input); }    

//  --- access/mutation ---

        
        std::vector< double > numerator( void );

        
        const std::vector< double > numerator( void ) const;
    
        
        std::vector< double > denominator( void );

        
        const std::vector< double > denominator( void ) const;
        
    
    void clear( void );

    
private:    
    
//  --- implementation ---

    std::deque< double > m_delayline;
        
    std::vector< double > m_ffwdcoefs;  

    std::vector< double > m_fbackcoefs; 
    
    double m_gain;      

};  //  end of class Filter



// ---------------------------------------------------------------------------
//  constructor
// ---------------------------------------------------------------------------
//
#if !defined(NO_TEMPLATE_MEMBERS)
template<typename IterT1, typename IterT2>
Filter::Filter( IterT1 ffwdbegin, IterT1 ffwdend,   //  feed-forward coeffs
                IterT2 fbackbegin, IterT2 fbackend, //  feedback coeffs
                double gain ) :
#else
inline 
Filter::Filter( const double * ffwdbegin, const double * ffwdend, //    feed-forward coeffs
                const double * fbackbegin, const double * fbackend, //  feedback coeffs
                double gain ) :
#endif
    m_ffwdcoefs( ffwdbegin, ffwdend ),
    m_fbackcoefs( fbackbegin, fbackend ),
    m_delayline( std::max( ffwdend-ffwdbegin, fbackend-fbackbegin ) - 1, 0. ),
    m_gain( gain )
{
    if ( *fbackbegin == 0. )
    {
        Throw( InvalidObject, 
               "Tried to create a Filter with feeback coefficient at zero delay equal to 0.0" );
    }

    //  normalize the coefficients by 1/a[0], if a[0] is not equal to 1.0
    //  (already checked for a[0] == 0 above)
    if ( *fbackbegin != 1. )
    {
        //  scale all filter coefficients by a[0]:
        std::transform( m_ffwdcoefs.begin(), m_ffwdcoefs.end(), m_ffwdcoefs.begin(),
                        std::bind2nd( std::divides<double>(), *fbackbegin ) );
        std::transform( m_fbackcoefs.begin(), m_fbackcoefs.end(), m_fbackcoefs.begin(), 
                        std::bind2nd( std::divides<double>(), *fbackbegin ) );
        m_fbackcoefs[0] = 1.;
    }
}


}   //  end of namespace Loris

#endif /* ndef INCLUDE_FILTER_H */

---