Mesh2D.h

#ifndef STK_MESH2D_H
#define STK_MESH2D_H
 
#include "Instrmnt.h"
#include "OnePole.h"
 
namespace stk {
 
/***************************************************/
/***************************************************/
 
const unsigned short NXMAX = 12;
const unsigned short NYMAX = 12;
 
class Mesh2D : public Instrmnt
{
 public:
  Mesh2D( unsigned short nX, unsigned short nY );
 
  ~Mesh2D( void );
 
  void clear( void );
 
  void setNX( unsigned short lenX );
 
  void setNY( unsigned short lenY );
 
  void setInputPosition( StkFloat xFactor, StkFloat yFactor );
 
  void setDecay( StkFloat decayFactor );
 
  void noteOn( StkFloat frequency, StkFloat amplitude );
 
  void noteOff( StkFloat amplitude );
 
  StkFloat energy( void );
 
  StkFloat inputTick( StkFloat input );
 
  void controlChange( int number, StkFloat value );
 
  StkFloat tick( unsigned int channel = 0 );
 
 
  StkFrames& tick( StkFrames& frames, unsigned int channel = 0 );
 
 protected:
 
  StkFloat tick0();
  StkFloat tick1();
  void clearMesh();
 
  unsigned short NX_, NY_;
  unsigned short xInput_, yInput_;
  OnePole  filterX_[NXMAX];
  OnePole  filterY_[NYMAX];
  StkFloat v_[NXMAX-1][NYMAX-1]; // junction velocities
  StkFloat vxp_[NXMAX][NYMAX];   // positive-x velocity wave
  StkFloat vxm_[NXMAX][NYMAX];   // negative-x velocity wave
  StkFloat vyp_[NXMAX][NYMAX];   // positive-y velocity wave
  StkFloat vym_[NXMAX][NYMAX];   // negative-y velocity wave
 
  // Alternate buffers
  StkFloat vxp1_[NXMAX][NYMAX];  // positive-x velocity wave
  StkFloat vxm1_[NXMAX][NYMAX];  // negative-x velocity wave
  StkFloat vyp1_[NXMAX][NYMAX];  // positive-y velocity wave
  StkFloat vym1_[NXMAX][NYMAX];  // negative-y velocity wave
 
  int counter_; // time in samples
};
 
inline StkFrames& Mesh2D :: tick( StkFrames& frames, unsigned int channel )
{
  unsigned int nChannels = lastFrame_.channels();
#if defined(_STK_DEBUG_)
  if ( channel > frames.channels() - nChannels ) {
    oStream_ << "Mesh2D::tick(): channel and StkFrames arguments are incompatible!";
    handleError( StkError::FUNCTION_ARGUMENT );
  }
#endif
 
  StkFloat *samples = &frames[channel];
  unsigned int j, hop = frames.channels() - nChannels;
  if ( nChannels == 1 ) {
    for ( unsigned int i=0; i<frames.frames(); i++, samples += hop )
      *samples++ = tick();
  }
  else {
    for ( unsigned int i=0; i<frames.frames(); i++, samples += hop ) {
      *samples++ = tick();
      for ( j=1; j<nChannels; j++ )
        *samples++ = lastFrame_[j];
    }
  }
 
  return frames;
}
 
} // stk namespace
 
#endif