Hi, Surfing the source code of some functions of cblas implementation I found: 1- In source_scal_r.h, _c_s.h and _c.h I think that the sentence ix = OFFSET(N, incX); is not needed because the functions only works when incX is greater than 0 and, in this case, OFFSET(N, incX) allways returns 0. 2- I can see that, as in the original BLAS, *axpy functions has loop unrolling by hand. Is this nessesary today with the optimizations of the copiler? I can see in other functions of the gsl-cblas that unrolling is not used (though in the original fortran BLAS is still used).
-- ***************************************** José Luis García Pallero [email protected] (o< / / \ V_/_ Use Debian GNU/Linux and enjoy! ***************************************** _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
