octave: src/DLD-FUNCTIONS/eig.cc annotate

annotate src/DLD-FUNCTIONS/eig.cc @ 7789:82be108cc558

First attempt at single precision tyeps * * * corrections to qrupdate single precision routines * * * prefer demotion to single over promotion to double * * * Add single precision support to log2 function * * * Trivial PROJECT file update * * * Cache optimized hermitian/transpose methods * * * Add tests for tranpose/hermitian and ChangeLog entry for new transpose code

author	David Bateman <dbateman@free.fr>
date	Sun, 27 Apr 2008 22:34:17 +0200
parents	91f8446ce4ae
children	87865ed7405f

rev	line source
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	1 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	2
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	3 Copyright (C) 1996, 1997, 1999, 2000, 2003, 2004, 2005, 2006, 2007
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	4 John W. Eaton
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	5
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	6 This file is part of Octave.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	7
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	8 Octave is free software; you can redistribute it and/or modify it
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	9 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5829 diff changeset	10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5829 diff changeset	11 option) any later version.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	12
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	13 Octave is distributed in the hope that it will be useful, but WITHOUT
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	16 for more details.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	17
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	18 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5829 diff changeset	19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 5829 diff changeset	20 <http://www.gnu.org/licenses/>.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	21
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	22 */
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	23
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	24 #ifdef HAVE_CONFIG_H
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	25 #include <config.h>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	26 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	27
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	28 #include "EIG.h"
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	29 #include "fEIG.h"
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	30
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	31 #include "defun-dld.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	32 #include "error.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	33 #include "gripes.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	34 #include "oct-obj.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	35 #include "utils.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	36
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	37 DEFUN_DLD (eig, args, nargout,
3548 ab7fa5a8f23f [project @ 2000-02-03 01:17:15 by jwe] jwe parents: 3372 diff changeset	38 "-- texinfo --\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	39 @deftypefn {Loadable Function} {@var{lambda} =} eig (@var{a})\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	40 @deftypefnx {Loadable Function} {[@var{v}, @var{lambda}] =} eig (@var{a})\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	41 The eigenvalues (and eigenvectors) of a matrix are computed in a several\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	42 step process which begins with a Hessenberg decomposition, followed by a\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	43 Schur decomposition, from which the eigenvalues are apparent. The\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	44 eigenvectors, when desired, are computed by further manipulations of the\n\
f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	45 Schur decomposition.\n\
5829 93785a1b0f97 [project @ 2006-05-23 18:34:10 by jwe] jwe parents: 5823 diff changeset	46 \n\
93785a1b0f97 [project @ 2006-05-23 18:34:10 by jwe] jwe parents: 5823 diff changeset	47 The eigenvalues returned by @code{eig} are not ordered.\n\
3372 f16c2ce14886 [project @ 1999-11-23 19:07:09 by jwe] jwe parents: 3014 diff changeset	48 @end deftypefn")
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	49 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	50 octave_value_list retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	51
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	52 int nargin = args.length ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	53
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	54 if (nargin != 1 \|\| nargout > 2)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	55 {
5823 080c08b192d8 [project @ 2006-05-19 05:32:17 by jwe] jwe parents: 5307 diff changeset	56 print_usage ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	57 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	58 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	59
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	60 octave_value arg = args(0);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	61
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4725 diff changeset	62 octave_idx_type nr = arg.rows ();
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4725 diff changeset	63 octave_idx_type nc = arg.columns ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	64
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	65 int arg_is_empty = empty_arg ("eig", nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	66 if (arg_is_empty < 0)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	67 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	68 else if (arg_is_empty > 0)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	69 return octave_value_list (2, Matrix ());
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	70
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	71 if (nr != nc)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	72 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	73 gripe_square_matrix_required ("eig");
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	74 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	75 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	76
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	77 Matrix tmp;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	78 ComplexMatrix ctmp;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	79 FloatMatrix ftmp;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	80 FloatComplexMatrix fctmp;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	81
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	82 if (arg.is_single_type ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	83 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	84 FloatEIG result;
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	85
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	86 if (arg.is_real_type ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	87 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	88 ftmp = arg.float_matrix_value ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	89
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	90 if (error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	91 return retval;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	92 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	93 result = FloatEIG (ftmp, nargout > 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	94 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	95 else if (arg.is_complex_type ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	96 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	97 fctmp = arg.float_complex_matrix_value ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	98
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	99 if (error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	100 return retval;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	101 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	102 result = FloatEIG (fctmp, nargout > 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	103 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	104
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	105 if (! error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	106 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	107 if (nargout == 0 \|\| nargout == 1)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	108 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	109 retval(0) = result.eigenvalues ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	110 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	111 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	112 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	113 // Blame it on Matlab.
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	114
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	115 FloatComplexDiagMatrix d (result.eigenvalues ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	116
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	117 retval(1) = d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	118 retval(0) = result.eigenvectors ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	119 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	120 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	121 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	122 else
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	123 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	124 EIG result;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	125
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	126 if (arg.is_real_type ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	127 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	128 tmp = arg.matrix_value ();
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	129
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	130 if (error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	131 return retval;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	132 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	133 result = EIG (tmp, nargout > 1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	134 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	135 else if (arg.is_complex_type ())
7578 91f8446ce4ae handle possible error from EIG John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	136 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	137 ctmp = arg.complex_matrix_value ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	138
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	139 if (error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	140 return retval;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	141 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	142 result = EIG (ctmp, nargout > 1);
7578 91f8446ce4ae handle possible error from EIG John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	143 }
91f8446ce4ae handle possible error from EIG John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	144 else
91f8446ce4ae handle possible error from EIG John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	145 {
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	146 gripe_wrong_type_arg ("eig", tmp);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	147 return retval;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	148 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	149
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	150 if (! error_state)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	151 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	152 if (nargout == 0 \|\| nargout == 1)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	153 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	154 retval(0) = result.eigenvalues ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	155 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	156 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	157 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	158 // Blame it on Matlab.
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	159
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	160 ComplexDiagMatrix d (result.eigenvalues ());
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	161
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	162 retval(1) = d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	163 retval(0) = result.eigenvectors ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7578 diff changeset	164 }
7578 91f8446ce4ae handle possible error from EIG John W. Eaton <jwe@octave.org> parents: 7017 diff changeset	165 }
2928 295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	166 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	167
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	168 return retval;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	169 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	170
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	171 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	172 ;;; Local Variables: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	173 ;;; mode: C++ ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	174 ;;; End: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe] jwe parents: diff changeset	175 */

Mercurial > octave

annotate src/DLD-FUNCTIONS/eig.cc @ 7789:82be108cc558