Systolic L-U decomposition array with a new reciprocal cell
V.K. Jain, David L. Landis, C.E. Alvarez · 2003
A systolic architecture for L-U decomposition using a recently developed reciprocal cell is presented. The arithmetic of this new cell is based on second-order polynomial interpolation, which results in high speed and inherent stability in inversion. In part, the speed of the cell is realized by use of a special empirical mapping. For the 16-b mantissa of a floating point number, the lookup table size required is 64 words, and the RMS value of the error is approximately one-third of the LSB of the 16-b result.>