Skip to content

DefiniteIntegral not respected when solving Hilbert problem #124

@ettersi

Description

@ettersi
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.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions