## Interior Point Algorithm for Linear Programming. Mathematica 5.0 (2003)

This is an implementation of an infeasible predictor-corrector interior point algorithm which is able to handle very large scale linear programming problems. It is available through the LinearProgramming function. Here is an example of solving a linear programming problem of 16558 constraints and 48763 variables. It took 83 seconds.

In[4]:= pds10=ExampleData[{"LinearProgramming","pds-10"}]; In[5]:= solution=LinearProgramming@@pds10;//Timing Out[5]= {83.1124, Null} In[6]:= pds10[[2]] Out[6]= SparseArray[<106436>, {16558, 48763}]