ApproxFun.jl icon indicating copy to clipboard operation
ApproxFun.jl copied to clipboard

Multiple Integrals give different results

Open Veenty opened this issue 10 months ago • 1 comments

So I'm doing iterated integrals and there is something funky when I composite them So for example Evaluation(Chebyshev(), -1)*Derivative(Chebyshev(),1)*Integral(Chebyshev(), 1) - Evaluation(Chebyshev(), -1)

Gives me the 0 vector, however

(Evaluation(Chebyshev(), -1)*Derivative(Chebyshev(),1)*Integral(Chebyshev(), 1))*Integral(Chebyshev(), 1)*Integral(Chebyshev(), 1) - Evaluation(Chebyshev(), -1)*Integral(Chebyshev(), 2)

Does not give me the 0 vector, I get

PlusOperator : Chebyshev() → ConstantSpace(Point(-1)) 0.0 0.125 2.77556e-17 -0.125 3.46945e-18 -1.73472e-17 2.1684e-17 -4.33681e-19 6.72205e-18 4.77049e-18 ⋯

Veenty avatar Feb 10 '25 16:02 Veenty

So it turns out that applying iterated integrals is not the same as applying Integrate(Chebyshev(), 2). However Integrate(Chebyshev(), 2) = Integrate(Chebyshev(), 1)*Integrate(Chebyshev(), 1), which honestly I think is quite confusing

Veenty avatar Feb 11 '25 20:02 Veenty