================================================================================ 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 ================================================================================