At Wed, 17 Jun 2009 13:31:10 +0200, José Luis García Pallero wrote: > 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.
You're right, I will remove that unnecessary line. Thanks. > 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). As you say, it is probably not needed with a recent compiler. I don't know whether it does any harm to leave it though (has anyone tried a comparison?). I am tempted to #ifdef it out and replace it with a plain loop. -- Brian Gough (GSL Maintainer) Support freedom by joining the FSF http://www.fsf.org/associate/support_freedom/join_fsf?referrer=37 _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
