search

JuliaMath / QuadGK.jl

adaptive 1d numerical Gauss–Kronrod integration in Julia
MIT License
268 stars 37 forks source link

order=1 seems børked #66

Closed stevengj closed 1 year ago

stevengj commented 1 year ago

This seems wrong:

julia> quadgk_count(x -> 1, 0, 1, order=1)
(2.4444444444019076, 0.5555555555356156, 23333345)
stevengj commented 1 year ago

Okay, the problem is that the unrolling in evalrule assumed a higher order.