================================================================================
Date: Mon, 2 Mar 2009 01:58:12 -0700
From: Vasile Sima
Cc: Peter Benner
Subject: [Lapack] Minor bug in DTRSYL

Hello,

I would like to draw your attention to a very minor bug in the file
dtrsyl.f (version 3.2, November 2006).  Specifically, the line 167,

      SCALE = ONE

must be moved before the line 152 (Section "Quick return if possible"),
since this should be an output even if min(M,N) = 0.

Sincerely yours,

Vasile Sima
================================================================================
Date: Tue, 3 Mar 2009 10:16:43 -0700
From: Jason Riedy
To: Vasile Sima
Cc: Peter Benner
Subject: Re: [Lapack] Minor bug in DTRSYL

And Vasile Sima writes:

> I would like to draw your attention to a very minor bug in the file
> dtrsyl.f (version 3.2, November 2006).  Specifically, the line 167,
> SCALE = ONE
> must be moved before the line 152 (Section "Quick return if possible"),
> since this should be an output even if min(M,N) = 0.

Committed that fix, thanks!

I doubt if LAPACK is completely consistent with zero-dimension matrices,
though.  Any other such fixes happily accepted!

Jason
================================================================================
Author: Jason Riedy
Date: Tue Mar  3 12:16:37 2009
New Revision: 644
URL: https://icl.cs.utk.edu/trac/lapack-dev/changeset/644
Log: Return SCALE = ONE in xTRSYL for zero-dimension matrices.

I doubt if LAPACK is consistent about zero-dim matrices, but it's
worth a shot.  Plenty of discussion about their semantics on the
GNU Octave development lists.

Reported by Vasile Sima

Signed-off-by: Jason Riedy

Modified:
    lapack/trunk/SRC/ctrsyl.f
    lapack/trunk/SRC/dtrsyl.f
    lapack/trunk/SRC/strsyl.f
    lapack/trunk/SRC/ztrsyl.f
================================================================================