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
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions