if you are interested, I did something similar in C++ using QuadProg++, a library that implements the algorithm of Goldfarb and Idnani for the solution of a (convex) Quadratic Programming problem by means of an active-set dual method.
I didn’t yet know Julia at that time… I plan to port it in Julia, but I’m too busy at the moment 