Interior point method