[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Help-gsl] About invalid parameters in cblas implementation

From: Brian Gough
Subject: Re: [Help-gsl] About invalid parameters in cblas implementation
Date: Tue, 30 Jun 2009 21:15:05 +0100
User-agent: Wanderlust/2.14.0 (Africa) Emacs/22.2 Mule/5.0 (SAKAKI)

At Sun, 28 Jun 2009 17:26:39 +0200,
José Luis García Pallero wrote:
> 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?

Good point, I *think* I considered this issue at the time but I don't
remember and may have overlooked these.  I will have to dig out the
original BLAS paper and check whether the use of XERBLA was optional
or part of the specification.

Brian Gough
(GSL Maintainer)

Support freedom by joining the FSF

reply via email to

[Prev in Thread] Current Thread [Next in Thread]