Primal-Dual Interior Point Method