Make gram_matrix and lagrangian_multipliers before optimize is called

[blegat]

Then the user has to call value(::GramMatrix) in order to get the actual values. See https://discourse.julialang.org/t/when-adding-a-sumofsquares-constraint-is-it-possible-to-get-the-jump-variable-corresponding-to-the-gram-matrix/127966