December 3rd, 2014, 8:47 am
QuoteOriginally posted by: outrunIt can be done (I've been explained, never tried it myself) with expression templatesWhat is ideal, is BLAS like this, directly in application codeQuoteGeneric applications can be written in a natural notation, e.g. v += A*q - w;, while the library dispatches to the appropriate algorithms: matrix vector products vs. matrix products vs. vector scalar products etcetera. The goal is to encapsulate performance issues inside the library and provide scientists an intuitive interface. MTL4 is used in different finite element and finite volume packages, e.g.