operator()() [2/2]

PolesNZeroes::Tcvalue fourier::PolesNZeroes::operator()

(

const Tcvalue &

omega

)

const

Definition at line 92 of file polesnzeroes.cc.

References Mdenominator, Mnumerator, Mpoles, and Mzeroes.

  {
    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);
  }