Due to the fact that the LP infrastructure
imposes a positivity constraint on the solution, we are forced to
solve for a model space twice as large as we would choose with both
negativity and positivity constraints, and then combine the two. We
will use the hyperbolic radon transform (HRT) as the analysis
operator, in equation (2) or *A* in equation
(3). This operator is normally approximately full and,
therefore, a tractable operator to use for this method. As
such, the programming is implemented with a sparse matrix approach rather
than operators as it is traditional at SEP.

10/14/2003