--On 7 mai 2013 09:27:07 -0400 Raul Miller <[email protected]> wrote:
My initial impression, here, is BLAS is written in fortran, suggesting that there's something wrong with the C language and/or compilers, for generating efficient numerical code.
BLAS was originally written in Fortran, but today most machines have BLAS libraries optimized beyond what any compiler would do. They are either hand-optimized or auto-tuned (look at Atlas for example). Optimization involves not only code generation, but also the choice of optimal chunk sizes to get the best possible cache performance.
BLAS is very hard to beat, so the practical question is: could one make J use BLAS for matrix multiplication (automatically, from standard J verbs) and other operations where BLAS offers support?
---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
