Three-monotone interpolation

A function $f\colon\mathbb R\to\mathbb R$ is called \emph{$k$-monotone} if it is $(k-2)$-times differentiable and its $(k-2)$nd derivative is convex. A point set $P\subset\mathbb R^2$ is \emph{$k$-monotone interpolable} if it lies on a graph of a $k$-monotone function. These notions have been studied in analysis, approximation theory etc. since the 1940s. We show that 3-monotone interpolability is very non-local: we exhibit an arbitrarily large finite $P$ for which every proper subset is $3$-monotone interpolable but $P$ itself is not. On the other hand, we prove a Ramsey-type result: for every $n$ there exists $N$ such that every $N$-point $P$ with distinct $x$-coordinates contains an $n$-point $Q$ such that $Q$ or its vertical mirror reflection are $3$-monotone interpolable. The analogs for $k$-monotone interpolability with $k=1$ and $k=2$ are classical theorems of Erd\H{o}s and Szekeres, while the cases with $k\ge4$ remain open. We also investigate the computational complexity of deciding $3$-monotone interpolability of a given point set. Using a known characterization, this decision problem can be stated as an instance of polynomial optimization and reformulated as a semidefinite program. We exhibit an example for which this semidefinite program has only doubly exponentially large feasible solutions, and thus known algorithms cannot solve it in polynomial time. While such phenomena have been well known for semidefinite programming in general, ours seems to be the first such example in polynomial optimization, and it involves only univariate quadratic polynomials.