#ifndef STK_FILTER_H#define STK_FILTER_H #include "Stk.h"#include <vector>#include <cmath> namespace stk { /***************************************************//*! \class Filter    \brief STK abstract filter class.     This class provides limited common functionality for STK digital    filter subclasses.  It is general enough to support both    monophonic and polyphonic input/output classes.     by Perry R. Cook and Gary P. Scavone, 1995-2012.*//***************************************************/ class Filter : public Stk{public:  //! Class constructor.  Filter( void ) { gain_ = 1.0; channelsIn_ = 1; lastFrame_.resize( 1, 1, 0.0 ); };   //! Return the number of input channels for the class.  unsigned int channelsIn( void ) const { return channelsIn_; };   //! Return the number of output channels for the class.  unsigned int channelsOut( void ) const { return lastFrame_.channels(); };   //! Clears all internal states of the filter.  virtual void clear( void );   //! Set the filter gain.  /*!    The gain is applied at the filter input and does not affect the    coefficient values.  The default gain value is 1.0.   */  void setGain( StkFloat gain ) { gain_ = gain; };   //! Return the current filter gain.  StkFloat getGain( void ) const { return gain_; };   //! Return the filter phase delay at the specified frequency.  /*!    Note that the phase delay calculation accounts for the filter    gain.  The frequency value should be greater than 0.0 and less    than or equal to one-half the sample rate.   */  StkFloat phaseDelay( StkFloat frequency );   //! Return an StkFrames reference to the last output sample frame.  const StkFrames& lastFrame( void ) const { return lastFrame_; };   //! Take a channel of the StkFrames object as inputs to the filter and replace with corresponding outputs.  /*!    The StkFrames argument reference is returned.  The \c channel    argument must be less than the number of channels in the    StkFrames argument (the first channel is specified by 0).    However, range checking is only performed if _STK_DEBUG_ is    defined during compilation, in which case an out-of-range value    will trigger an StkError exception.  */  virtual StkFrames& tick( StkFrames& frames, unsigned int channel = 0 ) = 0; protected:   StkFloat gain_;  unsigned int channelsIn_;  StkFrames lastFrame_;   std::vector<StkFloat> b_;  std::vector<StkFloat> a_;  StkFrames outputs_;  StkFrames inputs_; }; inline void Filter :: clear( void ){  unsigned int i;  for ( i=0; i<inputs_.size(); i++ )    inputs_[i] = 0.0;  for ( i=0; i<outputs_.size(); i++ )    outputs_[i] = 0.0;  for ( i=0; i<lastFrame_.size(); i++ )    lastFrame_[i] = 0.0;  } inline StkFloat Filter :: phaseDelay( StkFloat frequency ){  if ( frequency <= 0.0 || frequency > 0.5 * Stk::sampleRate() ) {    oStream_ << "Filter::phaseDelay: argument (" << frequency << ") is out of range!";    handleError( StkError::WARNING ); return 0.0;  }   StkFloat omegaT = 2 * PI * frequency / Stk::sampleRate();  StkFloat real = 0.0, imag = 0.0;  for ( unsigned int i=0; i<b_.size(); i++ ) {    real += b_[i] * std::cos( i * omegaT );    imag -= b_[i] * std::sin( i * omegaT );  }  real *= gain_;  imag *= gain_;   StkFloat phase = atan2( imag, real );   real = 0.0, imag = 0.0;  for ( unsigned int i=0; i<a_.size(); i++ ) {    real += a_[i] * std::cos( i * omegaT );    imag -= a_[i] * std::sin( i * omegaT );  }   phase -= std::atan2( imag, real );  phase = std::fmod( -phase, 2 * PI );  return phase / omegaT;} } // stk namespace #endif