Home   Information   Classes   Download   Usage   Mail List   Requirements   Links   FAQ   Tutorial


Fir.h
1 #ifndef STK_FIR_H
2 #define STK_FIR_H
3 
4 #include "Filter.h"
5 
6 namespace stk {
7 
8 /***************************************************/
28 /***************************************************/
29 
30 class Fir : public Filter
31 {
32 public:
34  Fir( void );
35 
37 
41  Fir( std::vector<StkFloat> &coefficients );
42 
44  ~Fir( void );
45 
47 
52  void setCoefficients( std::vector<StkFloat> &coefficients, bool clearState = false );
53 
55  StkFloat lastOut( void ) const { return lastFrame_[0]; };
56 
58  StkFloat tick( StkFloat input );
59 
61 
69  StkFrames& tick( StkFrames& frames, unsigned int channel = 0 );
70 
72 
80  StkFrames& tick( StkFrames& iFrames, StkFrames &oFrames, unsigned int iChannel = 0, unsigned int oChannel = 0 );
81 
82 protected:
83 
84 };
85 
86 inline StkFloat Fir :: tick( StkFloat input )
87 {
88  lastFrame_[0] = 0.0;
89  inputs_[0] = gain_ * input;
90 
91  for ( unsigned int i=(unsigned int)(b_.size())-1; i>0; i-- ) {
92  lastFrame_[0] += b_[i] * inputs_[i];
93  inputs_[i] = inputs_[i-1];
94  }
95  lastFrame_[0] += b_[0] * inputs_[0];
96 
97  return lastFrame_[0];
98 }
99 
100 inline StkFrames& Fir :: tick( StkFrames& frames, unsigned int channel )
101 {
102 #if defined(_STK_DEBUG_)
103  if ( channel >= frames.channels() ) {
104  oStream_ << "Fir::tick(): channel and StkFrames arguments are incompatible!";
105  handleError( StkError::FUNCTION_ARGUMENT );
106  }
107 #endif
108 
109  StkFloat *samples = &frames[channel];
110  unsigned int i, hop = frames.channels();
111  for ( unsigned int j=0; j<frames.frames(); j++, samples += hop ) {
112  inputs_[0] = gain_ * *samples;
113  *samples = 0.0;
114 
115  for ( i=(unsigned int)b_.size()-1; i>0; i-- ) {
116  *samples += b_[i] * inputs_[i];
117  inputs_[i] = inputs_[i-1];
118  }
119  *samples += b_[0] * inputs_[0];
120  }
121 
122  lastFrame_[0] = *(samples-hop);
123  return frames;
124 }
125 
126 inline StkFrames& Fir :: tick( StkFrames& iFrames, StkFrames& oFrames, unsigned int iChannel, unsigned int oChannel )
127 {
128 #if defined(_STK_DEBUG_)
129  if ( iChannel >= iFrames.channels() || oChannel >= oFrames.channels() ) {
130  oStream_ << "Fir::tick(): channel and StkFrames arguments are incompatible!";
131  handleError( StkError::FUNCTION_ARGUMENT );
132  }
133 #endif
134 
135  StkFloat *iSamples = &iFrames[iChannel];
136  StkFloat *oSamples = &oFrames[oChannel];
137  unsigned int i, iHop = iFrames.channels(), oHop = oFrames.channels();
138  for ( unsigned int j=0; j<iFrames.frames(); j++, iSamples += iHop, oSamples += oHop ) {
139  inputs_[0] = gain_ * *iSamples;
140  *oSamples = 0.0;
141 
142  for ( i=(unsigned int)b_.size()-1; i>0; i-- ) {
143  *oSamples += b_[i] * inputs_[i];
144  inputs_[i] = inputs_[i-1];
145  }
146  *oSamples += b_[0] * inputs_[0];
147  }
148 
149  lastFrame_[0] = *(oSamples-oHop);
150  return iFrames;
151 }
152 
153 } // stk namespace
154 
155 #endif

The Synthesis ToolKit in C++ (STK)
©1995--2014 Perry R. Cook and Gary P. Scavone. All Rights Reserved.