julia> D = Interval(-2,-1), Interval(1,2);
S = JacobiWeight(-0.5,-0.5, Chebyshev(D[1])) ∪ JacobiWeight(-0.5,-0.5, Chebyshev(D[2]))
μ = [
DefiniteIntegral(component(S,1));
DefiniteIntegral(component(S,2));
real(Hilbert(S))
] \ [-0.5;1;0]
@show sum(components(μ)[1])
@show sum(components(μ)[2])
;
sum((components(μ))[1]) = -5.505172187825942e15
sum((components(μ))[2]) = 5.505172187825942e15
The output on the last two lines should have been
sum((components(μ))[1]) = -0.5
sum((components(μ))[2]) = 1
I believe this used to work at some point, though I am not sure I remember correctly.
The output on the last two lines should have been
I believe this used to work at some point, though I am not sure I remember correctly.