From bug-octave-request at bevo dot che dot wisc dot edu  Wed Jan 26 00:15:48 2000
Subject: AIX compilation problem
From: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
To: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
Cc: "G. B|rger" <buerger at pik-potsdam dot de>, bug-octave@bevo.che.wisc.edu
Date: Wed, 26 Jan 2000 00:16:12 -0600 (CST)

On 14-May-1999, John W. Eaton <jwe at bevo dot che dot wisc dot edu> wrote:

| On 14-May-1999, G. B|rger <buerger at pik-potsdam dot de> wrote:
| 
| | Another question:
| | The function `pinv` behaves strange. For example, for a matrix
| | 
| | a=rand(10000,3)
| | 
| | both b=pinv(a) and b=pinv(a') generates a memory failure. On the other
| | hand, b=inv(a'*a)*a' is fine.
| 
| The pinv function uses the SVD to compute the solution.
| Unfortunately, it doesn't try to be smart about memory, so it attempts
| to generate the full SVD.  In your case, that means attempting to
| generate a 10000x10000 matrix (in either case).
| 
| Better implementations of the Matrix::psuedo_inverse and
| ComplexMatrix::psuedo_inverse functions are needed.

This problem is fixed in the development sources now.

Thanks,

jwe



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

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