From owner-bug-octave at bevo dot che dot wisc dot edu  Thu Dec  5 13:11:45 1996
Subject: bug report for octave 1.93
From: "John W. Eaton" <jwe at bevo dot che dot wisc dot edu>
To: pfost at aghl dot ruhr-uni-bochum dot de (Martin Pfost)
cc: bug-octave at bevo dot che dot wisc dot edu
Date: Thu, 5 Dec 1996 13:09:07 -0600

On  5-Dec-1996, Martin Pfost <pfost at aghl dot ruhr-uni-bochum dot de> wrote:

: Octave 1.93 returns something weird for a non-square system 
: of linear equations with complex numbers:
: 
:     octave:1> [i;2i] \ [3i;6i]
:     ans = [](0x0)
: 
: However, everything is fine with reals only:
: 
:     octave:2> [1;2] \ [3;6]
:     ans = 3.0000

Please try the following patch.

Thanks,

jwe


Index: CMatrix.cc
===================================================================
RCS file: /home/jwe/src/master/octave/liboctave/CMatrix.cc,v
retrieving revision 1.42
retrieving revision 1.43
diff -c -r1.42 -r1.43
*** CMatrix.cc	1996/10/30 18:51:26	1.42
--- CMatrix.cc	1996/12/05 19:06:58	1.43
***************
*** 1416,1422 ****
--- 1416,1424 ----
        int len_s = m < n ? m : n;
        Array<double> s (len_s);
        double *ps = s.fortran_vec ();
+ 
        double rcond = -1.0;
+ 
        int lwork;
        if (m < n)
  	lwork = 2*m + (nrhs > n ? nrhs : n);
***************
*** 1439,1445 ****
  	(*current_liboctave_error_handler) ("unrecoverable error in zgelss");
        else
  	{
! 	  ComplexMatrix retval (n, nrhs);
  	  for (int j = 0; j < nrhs; j++)
  	    for (int i = 0; i < n; i++)
  	      retval.elem (i, j) = result.elem (i, j);
--- 1441,1447 ----
  	(*current_liboctave_error_handler) ("unrecoverable error in zgelss");
        else
  	{
! 	  retval.resize (n, nrhs);
  	  for (int j = 0; j < nrhs; j++)
  	    for (int i = 0; i < n; i++)
  	      retval.elem (i, j) = result.elem (i, j);
***************
*** 1519,1525 ****
  	(*current_liboctave_error_handler) ("unrecoverable error in zgelss");
        else
  	{
! 	  ComplexColumnVector retval (n);
  	  for (int i = 0; i < n; i++)
  	    retval.elem (i) = result.elem (i);
  	}
--- 1521,1527 ----
  	(*current_liboctave_error_handler) ("unrecoverable error in zgelss");
        else
  	{
! 	  retval.resize (n);
  	  for (int i = 0; i < n; i++)
  	    retval.elem (i) = result.elem (i);
  	}