Improving numerical software

It turns out both are true, for different groups of users. Traditionally, users have asked for and been provided with “black box” software in the form of mathematical libraries such as LAPACK, LINPACK, NAG, and IMSL. More recently, the high-performance community has discovered that they must write custom software for their problem. Their reasons include inadequate functionality of existing software libraries, data structures that are not natural or convenient for a particular problem and overly general software that sacrifices too much performance when applied to a special case of interest. Can we improve numerical libraries? The answer is yes and we are going to investigate in this paper some methods.