Why not use a convex hull as the projection function?

[zalo]

The fit function might require a library, but it will be tighter and the projection function is the same (just a series of plane constraints).

Also, just wanted to say that this is a super neat trick, and I appreciate your writing style of including snippets in the posts; I try do something similar in my interactive posts.

( @orangeduck in-case the notifications aren't set up)