annotate src/qzval.cc @ 2599:9705ea40f0a2

[project @ 1996-12-19 04:49:09 by jwe]
author jwe
date Thu, 19 Dec 1996 04:49:28 +0000
parents 8b0911d576bf
children 8b262e771614
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
1 /*
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
2
1884
e62277bf5fe0 [project @ 1996-02-05 18:17:59 by jwe]
jwe
parents: 1847
diff changeset
3 Copyright (C) 1996 John W. Eaton
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
4
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
6
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
7 Octave is free software; you can redistribute it and/or modify it
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
8 under the terms of the GNU General Public License as published by the
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
9 Free Software Foundation; either version 2, or (at your option) any
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
10 later version.
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
11
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
12 Octave is distributed in the hope that it will be useful, but WITHOUT
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
15 for more details.
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
16
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
18 along with Octave; see the file COPYING. If not, write to the Free
1315
611d403c7f3d [project @ 1995-06-25 19:56:32 by jwe]
jwe
parents: 1255
diff changeset
19 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
20
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
21 */
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
22
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
23 // Written by A. S. Hodel <scotte@eng.auburn.edu>
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
24
240
a99f28f5e351 [project @ 1993-11-30 20:24:36 by jwe]
jwe
parents: 234
diff changeset
25 #ifdef HAVE_CONFIG_H
1192
b6360f2d4fa6 [project @ 1995-03-30 21:38:35 by jwe]
jwe
parents: 1009
diff changeset
26 #include <config.h>
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
27 #endif
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
28
1342
61bb2bdee11e [project @ 1995-09-04 00:19:22 by jwe]
jwe
parents: 1315
diff changeset
29 #include <cfloat>
112
3438ee77fea2 [project @ 1993-09-17 06:11:35 by jwe]
jwe
parents: 99
diff changeset
30
1352
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe]
jwe
parents: 1342
diff changeset
31 #include "defun-dld.h"
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe]
jwe
parents: 1342
diff changeset
32 #include "error.h"
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe]
jwe
parents: 1342
diff changeset
33 #include "gripes.h"
19c10b8657d5 [project @ 1995-09-05 08:11:57 by jwe]
jwe
parents: 1342
diff changeset
34 #include "help.h"
1740
fe9d3b2ded26 [project @ 1996-01-12 11:03:26 by jwe]
jwe
parents: 1684
diff changeset
35 #include "oct-obj.h"
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
36
2465
7ee42ff6536a [project @ 1996-11-03 08:54:14 by jwe]
jwe
parents: 2086
diff changeset
37 DEFUN_DLD (qzval, args, ,
519
b9284136189a [project @ 1994-07-19 14:40:20 by jwe]
jwe
parents: 516
diff changeset
38 "X = qzval (A, B)\n\
b9284136189a [project @ 1994-07-19 14:40:20 by jwe]
jwe
parents: 516
diff changeset
39 \n\
b9284136189a [project @ 1994-07-19 14:40:20 by jwe]
jwe
parents: 516
diff changeset
40 compute generalized eigenvalues of the matrix pencil (A - lambda B).\n\
b9284136189a [project @ 1994-07-19 14:40:20 by jwe]
jwe
parents: 516
diff changeset
41 A and B must be real matrices.")
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
42 {
2086
bfb775fb6fe8 [project @ 1996-04-25 05:55:19 by jwe]
jwe
parents: 1996
diff changeset
43 octave_value retval;
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
44
712
36ba0576bd1b [project @ 1994-09-19 14:18:15 by jwe]
jwe
parents: 701
diff changeset
45 int nargin = args.length ();
36ba0576bd1b [project @ 1994-09-19 14:18:15 by jwe]
jwe
parents: 701
diff changeset
46
1954
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
47 if (nargin == 2)
636
fae2bd91c027 [project @ 1994-08-23 18:39:50 by jwe]
jwe
parents: 628
diff changeset
48 {
2086
bfb775fb6fe8 [project @ 1996-04-25 05:55:19 by jwe]
jwe
parents: 1996
diff changeset
49 octave_value arg_a = args(0);
bfb775fb6fe8 [project @ 1996-04-25 05:55:19 by jwe]
jwe
parents: 1996
diff changeset
50 octave_value arg_b = args(1);
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
51
1954
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
52 Matrix a = arg_a.matrix_value ();
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
53 Matrix b = arg_b.matrix_value ();
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
54
1954
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
55 if (! error_state)
2563
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
56 {
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
57 ComplexColumnVector tmp = Qzval (a, b);
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
58
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
59 if (! error_state)
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
60 retval = tmp;
8b0911d576bf [project @ 1996-12-05 19:06:58 by jwe]
jwe
parents: 2465
diff changeset
61 }
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
62 }
1954
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
63 else
5328a7bc9b9d [project @ 1996-02-14 20:23:31 by jwe]
jwe
parents: 1884
diff changeset
64 print_usage ("qzval");
636
fae2bd91c027 [project @ 1994-08-23 18:39:50 by jwe]
jwe
parents: 628
diff changeset
65
44
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
66 return retval;
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
67 }
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
68
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
69 /*
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
70 ;;; Local Variables: ***
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
71 ;;; mode: C++ ***
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
72 ;;; End: ***
f3215b07c171 [project @ 1993-08-10 22:33:20 by jwe]
jwe
parents:
diff changeset
73 */