> Completely devectorisation can sometimes be useful, but I have come to > prefer writing linear algebra code close to the style of LAPACK where you > organise the calculations in BLAS calls. Not only for speed, but also for > making the code easier to read and infer the structure of the algorithm. > > Med venlig hilsen > > Andreas Noack > > Organizing linear algebra code into BLAS calls make sense to me too. It is like writing vectorized (easy to read) code but with excellent speed and efficient memory allocation.
Jan
