From bug-octave-request at bevo dot che dot wisc dot edu Thu Nov 14 09:14:54 2002
Subject: problems with svd, pinv an other  functions
From: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
To: "Eugenio SEGUIN" <eseguin at idecom dot unsj dot edu dot ar>
Cc: bug-octave at bevo dot che dot wisc dot edu
Date: Thu, 14 Nov 2002 09:14:43 -0600

On 14-Nov-2002, Eugenio SEGUIN <eseguin at idecom dot unsj dot edu dot ar> wrote:

| To: bug-octave at bevo dot che dot wisc dot edu
| Cc: eseguin
| Subject: [50 character or so descriptive subject here (for reference)]
| --------
| Bug report for Octave 2.1.36 configured for alphaev67-dec-osf5.1
| 
| Description:
| -----------
| The svd routine not work, the inv or inverse, pinv and x=A\b routine 
| not work 
| with complex matrix. The message when I run octave-bug is:
| Warning: Cannot convert string "-dt-interface user-medium-r-normal-
| s*-*-*-*-*-*-*-*-*" to type FontStruct
| Warning: Cannot convert string "-*-Menu-Medium-R-Normal--*-120-*-
| *-P-*-ISO8859-1" to type FontStruct
| 
|  
| Repeat-By:
| ---------
| [u s v]= svd(hilb(3))
|  ** On entry to DGESVD parameter number 13 had an illegal value
| error: exception encountered in Fortran subroutine dgesvd_
| error: unrecoverable error in dgesvd
| 
| A=ceil(10*randn(5,3))

This happens because you are linked to an old version of lapack that
does not support all the features that Octave expects (in this case,
being able to pass a negative value in for the length of the work
array and have dgesvd return the required value).

Either install ATLAS and rebuild Octave, or rebuild Octave using the
--with-lapack=no option when running configure.

jwe



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------