Mereni

test.jl

using MovingWeightedLeastSquares

f(x) = sin(2x)/2 + 10

function t()
  xs = collect(-6:0.2:6)
  fs = [f(x) for x in xs]
  ε = 0.4
  w = (d, eps) -> (exp(-d^2))
  obj = mwls_naive(xs, fs, ε, w, maxDegree = 5)
  for x in xs
    res = obj(x)
    ref = f(x)
    @printf("%lf\t%lf\t%lf\t%lf\n", x, res, ref, abs(res - ref))
  end
end