================================================================================ Date: Sat, 28 Feb 2009 23:07:19 -0700 From: Daniel Waggoner Subject: [Lapack] a possible but in ?tgsen Dear Sirs: I am using LAPACK 3.2. I have found that the routine dtgsen() occasionally crashes. Upon further investigation, on line 704 of dtgsen.f, it appears that the matrix Q is referenced even if WANTQ is false. I suspect that there should be an if statement and that Q should only be changed only if WANTQ is true. Hope you find this of use and, as you know, LAPACK 3.2 is a great collection of programs. Daniel Waggoner ================================================================================ Date: Tue, 3 Mar 2009 10:15:50 -0700 From: Jason Riedy To: Daniel Waggoner Subject: Re: [Lapack] a possible but in ?tgsen And Daniel Waggoner writes: > I am using LAPACK 3.2. I have found that the routine dtgsen() > occasionally crashes. Upon further investigation, on line 704 of > dtgsen.f, it appears that the matrix Q is referenced even if WANTQ is > false. I suspect that there should be an if statement and that Q should > only be changed only if WANTQ is true. Committed that fix, thanks! Jason ================================================================================ [Lapack SVN] r643 - in /lapack/trunk/SRC: dtgsen.f stgsen.f Author: Jason Riedy Date: Tue Mar 3 12:12:04 2009 New Revision: 643 Log: Only touch Q in {d,s}tgsen.f if WANTQ is true. Reported by Daniel Waggoner Signed-off-by: Jason Riedy Modified: lapack/trunk/SRC/dtgsen.f lapack/trunk/SRC/stgsen.f ================================================================================