Hello, Working with cblas_dgemv() function I can see various diferences between refblas and gsl cblas implementations about invalid parameters. While in gsl cblas the function only abort the execution because incorrect values in "order" and "Trans" input parameters, in the refblas implementation the function calls xerbla when M<0, N<0, LDA<max(1,M), incX==0 or incY==0 (in this cases gsl cblas do not make operations and return). In the gsl_blas_* wrapper this cases are not catched. This produces different executions if you link your programs with gslcblas or refblas, lapack, atlas, etc. Are these differences a known fact?
Thanks -- ***************************************** 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
