Quadratic Programming Interior Point Method