search

JuliaIntervals / TaylorModels.jl

Rigorous function approximation using Taylor models in Julia
Other
63 stars 14 forks source link

Error when changing polynomial degree for `TaylorModel1{TaylorN{T} ,T}` #135

Closed lbenet closed 2 years ago

lbenet commented 2 years ago

The following yields an error:

julia> using TaylorModels;

julia> xN, yN = set_variables("x y", order=4, numvars=2);

julia> xtm1_order5 = exp( TaylorModel1( Taylor1(xN, 5), 0..0, 0..0, -0.5 .. 0.5) )
  1.0 + 1.0 x + 0.5 x² + 0.16666666666666666 x³ + 0.041666666666666664 x⁴ + [-4.63936, 3.17106]

julia> xtm1_order4 = TaylorModel1( Taylor1(xN, 4), 0..0, 0..0, -0.5 .. 0.5)
  1.0 x + [0, 0]

julia> new_xtm1_order4, _ = TaylorSeries.fixorder(xtm1_order5, xtm1_order4)
ERROR: MethodError: no method matching TaylorModel1(::Taylor1{TaylorN{Float64}}, ::TaylorN{Interval{Float64}}, ::Interval{Float64}, ::Interval{Float64})
...

cc: @AnderGray