From owner-bug-octave at bevo dot che dot wisc dot edu  Mon Nov 13 07:24:18 1995
Subject: ComplexDiagMatrix::inverse(info) 
From: "Petr Mikulik" <mikulik at sci dot muni dot cz>
To: bug-octave at che dot utexas dot edu
Date: Mon, 13 Nov 95 14:23:30 CST


I report the following bug:

ComplexDiagMatrix F(N,N);
F.fill(0);
for (int i=0; i<N;i++) F(i,i)=cos(i)+sin(i);
cout <<"DIAG MATRIX =\n"<<F<<endl;
int info;
ComplexDiagMatrix G( F.inverse(info) );
cout <<"INVERSE DIAG MATRIX =\n"<<G<<endl;
cout <<"info="<<info<<endl;
cout <<"THIS MUST BE UNITY =\n"<<F*G<<endl;
cout <<"BUT IT IS NOT!\n";


I think that in CDiagMatrix.cc, routine ComplexDiagMatrix::inverse(info) 
the last statement should correctly be
  return retval
instead of 
  return *this



***  Petr Mikulik                                                    ***
***    Department of Solid State Physics, Masaryk University         ***
***    Kotlarska 2,  611 37 Brno, Czech Republic                     ***
***  tel: (42-5-)41129378          fax: (42-5-)41211214              ***
***  e-mail: mikulik at sci dot muni dot cz   http://www.sci.muni.cz/~mikulik/  ***
------------------------------------------------------------------------
| ... and a little bit of the theory of fun ...                        |
| Laws of Computerdom according to Golub: 2. A carelessly planned      |
| project takes three times longer to complete than expected; a        |
| carefully planned project takes only twice as long.                  |
------------------------------------------------------------------------