In this work, a specialized interior point method is developed for planning cancer treatment by radiotherapy. This is a large-scale problem with a specific matrix structure. That structure is explored in an efficient way reducing the dimension of the linear system, which must be solved at each iteration. Moreover, the system obtained is sparse, symmetric and positive definite. Matlab numerical results show that large-scale problems can be solved in few iterations and short computational time by this approach.
interior point methods; linear programming; radiotherapy