示例#1
0
  def test___init__(self):
    ''' Just run the linear and quadratic interpolations '''

    fseries, signalseries, errorseries = self.ServeInput()
    interpolator = Interpolator(fseries, signalseries, running_regr_type = 'linear', windowlength=0)
    y_interpolated = interpolator.y(fseries)
    e_interpolated = interpolator.e(fseries)
    self.assertAlmostEqual(numpy.sum(signalseries-y_interpolated), 0.0, places=10)
    self.assertAlmostEqual(numpy.sum(e_interpolated), 0.0, places=10)

    fseries, signalseries, errorseries = self.ServeInput()
    interpolator = Interpolator(fseries, signalseries, running_regr_type = 'linear', windowlength=5)
    y_interpolated = interpolator.y(fseries)
    e_interpolated = interpolator.e(fseries)

    fseries, signalseries, errorseries = self.ServeInput(create_errorseries=True)
    interpolator = Interpolator(fseries, signalseries, running_regr_type = 'linear', windowlength=5)
    y_interpolated = interpolator.y(fseries)
    e_interpolated = interpolator.e(fseries)

    fseries, signalseries, errorseries = self.ServeInput()
    interpolator = Interpolator(fseries, signalseries, running_regr_type = 'quadratic', windowlength=5)
    y_interpolated = interpolator.y(fseries)
    e_interpolated = interpolator.e(fseries)

    fseries, signalseries, errorseries = self.ServeInput(create_errorseries=True)
    interpolator = Interpolator(fseries, signalseries, running_regr_type = 'quadratic', windowlength=5)
    y_interpolated = interpolator.y(fseries)
    e_interpolated = interpolator.e(fseries)