operator()()

LinearInterpolator::Tvalue ts::ipo::LinearInterpolator::operator() ( const libtime::TAbsoluteTime &  t ) const

virtual

calculate interpolated value

Parameters

t	time for value to return

Returns

new sample value at time t

Implements ts::ipo::Interpolator.

Definition at line 79 of file ipolin.cc.

References ts::ipo::Interpolator::debug(), ts::TimeSeries< S, H >::header, Madjustedfirst, Mdt, Minputwindow, ts::ipo::Interpolator::Mts, TSXX_assert, TSXX_debug, and TSXX_value.

    {
      using libtime::TAbsoluteTime;
      using libtime::TRelativeTime;
      const Tconst_series& s=this->Mts;
      const Theader& hd=this->Mts.header;
      //aff::Tsubscript ileft=s.f()+int(libtime::time2double(t-first)/hd.dt);
      /*
       * Notice that the division in libtime is defined to return the number
       * of samples that will come closest to the time interval. That means
       * that the operator round to the nearest integer and not to the lower
       * integer. For this reason we use the adjusted time of first sample.
       */
      aff::Tsubscript ileft=s.first()+(t-Madjustedfirst)/Mdt;
      TSXX_assert((Minputwindow.includes(t) && (ileft<=s.last())),
                  "requested sample lies outside time series");
      aff::Tsubscript iright= (ileft < s.last()) ? ileft+1 : ileft;
      double f=libtime::time2double(t-sff::wid2isample(hd, ileft-s.f()))/hd.dt;
      TSXX_debug(this->debug(), "LinearInterpolator::operator():",
           TSXX_value(ileft) << " " <<
           TSXX_value(iright) << " " <<
           TSXX_value(s.f()) << " " <<
           TSXX_value(f) << " " <<
           TSXX_value(hd.dt) << std::endl <<
           "                                " << TSXX_value(Mdt.timestring())
           << std::endl <<
           "                     " << TSXX_value(Madjustedfirst.timestring())
           << std::endl <<
           "  " << TSXX_value(sff::wid2isample(hd, ileft-s.f()).timestring())
           << std::endl <<
           "                                  " << TSXX_value(t.timestring())
           << std::endl <<
           " " << TSXX_value(sff::wid2isample(hd, iright-s.f()).timestring())
           << std::endl <<
           "                 " << TSXX_value((t-Madjustedfirst).timestring())
           << std::endl <<
           "                 " << TSXX_value((t-Madjustedfirst)/Mdt)
                );
      Tvalue retval=((1.-f)*s(ileft)+f*s(iright));
      return(retval);
    }

Here is the call graph for this function: