From ripley at stats.ox.ac.uk  Mon Aug 29 10:05:51 2005
From: ripley at stats.ox.ac.uk (ripley@stats.ox.ac.uk)
Date: Mon, 29 Aug 2005 10:05:51 +0200 (CEST)
Subject: [Rd] xerbla called from BLAS routine (PR#8100)
Message-ID: <20050829080551.212E718B8A@slim.kubism.ku.dk>

  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.

--27464147-1060265079-1125299740=:28726
Content-Type: TEXT/PLAIN; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: QUOTED-PRINTABLE

Please check out

http://www.netlib.org/blas/index.html

file=09xerbla.f
for=09error handler for the Level 2 and 3 BLAS routines.
,=09(Auxiliary Routine)

and the refefence code which says

*  XERBLA  is an error handler for the LAPACK routines.
*  It is called by an LAPACK routine if an input parameter has an
*  invalid value.  A message is printed and execution stops.

So, XERBLA is an LAPACK handler, intended to be called from LAPACK=20
routines -- the level 2 and level 3 BLAS routines are considered to be
part of LAPACK in the reference implementation.

On Mon, 29 Aug 2005, Prof Brian Ripley wrote:

> I don't believe your message actually came from R's XERBLA.  Note that th=
e=20
> name is not quoted and that the name of the routine is not 6 characters l=
ong.
>
> Here is a genuine example from R:
>
> Error in La.svd(x, nu, nv) : LAPACK routine 'DGEBRD' gave error code -10
>
> There is no R 2.2.0, BTW.  If referring to R-devel, please do so with an =
SVN=20
> revision.
>
> On Sun, 28 Aug 2005 gb at stat.umu.se wrote:
>
>> Full_Name: G=F6ran Brostr=F6m
>> Version: R-2.1.1, 2.2.0
>> OS: Debian unstable
>> Submission from: (NULL) (213.65.9.59)
>>=20
>>=20
>> Some BLAS routines call xerbla for error messages, which results in a=20
>> message
>> like
>> "LAPACK routine DGER gave error code -9".
>>=20
>> Suggested solution: In
>>=20
>> void F77_NAME(xerbla)(char *srname, int *info)
>> {
>>   /* srname is not null-terminated.  It should be 6 characters. */
>>    char buf[7];
>>    strncpy(buf, srname, 6);
>>    buf[6] =3D '\0';
>>    error(_("LAPACK routine '%6s' gave error code %d"), buf, -(*info));
>> }
>>=20
>> change 'LAPACK' to 'The' (or 'The BLAS/LAPACK').
>>=20
>> ______________________________________________
>> R-devel at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-devel
>>=20
>>=20
>
> --=20
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595

--=20
Brian D. Ripley,                  ripley at stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595
--27464147-1060265079-1125299740=:28726--



[Rd] xerbla called from BLAS routine (PR#8100)
ripley@stats.ox.ac.uk