Abstract
The mathematical description of the algorithm presented in a previous paper and a more efficient software implementation is given.