polesnzeroes.cc

Go to the documentation of this file.

#define TF_POLESNZEROES_CC_VERSION \
  "TF_POLESNZEROES_CC   V1.0   "

#include <fourier/polesnzeroes.h>

namespace fourier {

  /*----------------------------------------------------------------------*/

  void PolesNZeroes::clear() 
  {
    Mpoles.clear(); Mzeroes.clear(); 
    Mnumerator=1.; Mdenominator=1.;
    Mnormal=true;
  }

  /*----------------------------------------------------------------------*/

  void PolesNZeroes::setpole(const PolesNZeroes::Tcvalue& pole) 
  {
    if (isnormal()) 
    { Mpoles.push_back(pole); }
    else
    { Mzeroes.push_back(pole); }
  }

  /*----------------------------------------------------------------------*/

  void PolesNZeroes::setzero(const PolesNZeroes::Tcvalue& zero) 
  {
    if (isnormal()) 
    { Mzeroes.push_back(zero); }
    else
    { Mpoles.push_back(zero); }
  }

  /*----------------------------------------------------------------------*/

  void PolesNZeroes::numfactor(const PolesNZeroes::Tcvalue& factor) 
  {
    if (isnormal()) 
    { Mnumerator *= factor; }
    else
    { Mdenominator *= factor; }
  }

  /*----------------------------------------------------------------------*/

  void PolesNZeroes::denfactor(const PolesNZeroes::Tcvalue& factor) 
  {
    if (isnormal()) 
    { Mdenominator *= factor; }
    else
    { Mnumerator *= factor; }
  }

  /*----------------------------------------------------------------------*/

  PolesNZeroes::Tcvalue PolesNZeroes::operator()(const Tcvalue& omega)  const
  {
    const double waterlevel=1.e-30;
    Tcvalue om = std::abs(omega) > waterlevel ? omega : waterlevel;
    Tcvalue numerator=Mnumerator;
    Tcvalue denominator=Mdenominator;
    Tlist::const_iterator in=Mzeroes.begin();
    while (in != Mzeroes.end())
    { numerator *= (om - *in); ++in; }
    Tlist::const_iterator id=Mpoles.begin();
    while (id != Mpoles.end())
    { denominator *= (om - *id); ++id; }
    return(numerator/denominator);
  }

} // namespace fourier

/* ----- END OF polesnzeroes.cc ----- */