Im new to mathematical optimisation, so I’m looking for a simple example I can learn from …

My (initial) problem looks like

1000 cartesian points in Float16

I want to find a polynomial p(x) (to be evaluated in Float16) to fit these points.

Conceptually

my constraints are |p(x) - y| <= \text{ulp}(y) and

the objective is \text{Min} of \text{norm(|p(x) - y|, Inf)}

But Im not seeing how to declare it within JuMP and I haven’t found a suitable example to learn from.

Any pointers ?

thanks in advance,

JT

Edit: I corrected the constraints refer to \text{ulp}(y) instead of x.