  Numerical software in computational science and engineering often relies on
highly-optimized building blocks from libraries such as BLAS and LAPACK, and
while such libraries provide portable performance for a wide range of computing
architectures, they still present limitations in terms of flexibility. We
advocate a domain-specific program generator capable of producing library
routines tailored to the specific needs of the application in terms of sizes,
interface, and target architecture.
