annotate src/DLD-FUNCTIONS/minmax.cc @ 3800:115cbfcd067e

[project @ 2001-02-22 21:55:15 by jwe]
author jwe
date Thu, 22 Feb 2001 21:55:17 +0000
parents 13905c3a24af
children cbac6756967e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
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
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
3 Copyright (C) 1996, 1997 John W. Eaton
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
4
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
6
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
7 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
8 under the terms of the GNU General Public License as published by the
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
9 Free Software Foundation; either version 2, or (at your option) any
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
10 later version.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
11
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
12 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
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
14 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
15 for more details.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
16
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
18 along with Octave; see the file COPYING. If not, write to the Free
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
19 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
20
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 #ifdef HAVE_CONFIG_H
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
24 #include <config.h>
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
25 #endif
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
26
3145
0d640dc625c7 [project @ 1998-02-05 08:44:59 by jwe]
jwe
parents: 3014
diff changeset
27 #include <cmath>
0d640dc625c7 [project @ 1998-02-05 08:44:59 by jwe]
jwe
parents: 3014
diff changeset
28
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
29 #include "lo-ieee.h"
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
30 #include "lo-mappers.h"
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
31
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
32 #include "defun-dld.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
33 #include "error.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
34 #include "gripes.h"
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
35 #include "oct-obj.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 // XXX FIXME XXX -- it would be nice to share code among the min/max
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
38 // functions below.
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
39
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
40 #define EMPTY_RETURN_CHECK(T) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
41 if (nr == 0 || nc == 0) \
3800
115cbfcd067e [project @ 2001-02-22 21:55:15 by jwe]
jwe
parents: 3775
diff changeset
42 return T (nr, nc);
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
43
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
44 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
45 min (double d, const Matrix& m)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
46 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
47 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
48 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
49
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
50 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
51
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
52 Matrix result (nr, nc);
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 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
55 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
56 result (i, j) = xmin (d, m (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
57
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
58 return result;
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
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
61 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
62 min (const Matrix& m, double d)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
63 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
64 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
65 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
66
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
67 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
68
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
69 Matrix result (nr, nc);
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 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
72 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
73 result (i, j) = xmin (m (i, j), d);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
74
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
75 return result;
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
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
78 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
79 min (const Complex& c, const ComplexMatrix& m)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
80 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
81 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
82 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
83
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
84 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
85
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
86 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
87
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
88 for (int j = 0; j < nc; j++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
89 for (int i = 0; i < nr; i++)
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
90 result (i, j) = xmin (c, m (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
91
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
92 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
93 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
94
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
95 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
96 min (const ComplexMatrix& m, const Complex& c)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
97 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
98 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
99 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
100
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
101 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
102
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
103 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
104
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
105 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
106 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
107 result (i, j) = xmin (m (i, j), c);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
108
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
109 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
110 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
111
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
112 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
113 min (const Matrix& a, const Matrix& b)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
114 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
115 int nr = a.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
116 int nc = a.columns ();
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
117
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
118 if (nr != b.rows () || nc != b.columns ())
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
119 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
120 error ("two-arg min expecting args of same size");
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
121 return Matrix ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
122 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
123
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
124 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
125
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
126 Matrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
127
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
128 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
129 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
130 result (i, j) = xmin (a (i, j), b (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
131
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
132 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
133 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
134
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
135 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
136 min (const ComplexMatrix& a, const ComplexMatrix& b)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
137 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
138 int nr = a.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
139 int nc = a.columns ();
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
140
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
141 if (nr != b.rows () || nc != b.columns ())
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
142 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
143 error ("two-arg min expecting args of same size");
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
144 return ComplexMatrix ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
145 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
146
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
147 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
148
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
149 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
150
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
151 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
152 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
153 int columns_are_real_only = 1;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
154 for (int i = 0; i < nr; i++)
3243
dd00769643ae [project @ 1999-05-28 04:19:00 by jwe]
jwe
parents: 3145
diff changeset
155 if (imag (a (i, j)) != 0.0 || imag (b (i, j)) != 0.0)
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
156 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
157 columns_are_real_only = 0;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
158 break;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
159 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
160
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
161 if (columns_are_real_only)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
162 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
163 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
164 result (i, j) = xmin (real (a (i, j)), real (b (i, j)));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
165 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
166 else
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 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
169 result (i, j) = xmin (a (i, j), b (i, j));
2928
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
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
173 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
174 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
175
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
176 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
177 max (double d, const Matrix& m)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
178 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
179 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
180 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
181
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
182 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
183
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
184 Matrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
185
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
186 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
187 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
188 result (i, j) = xmax (d, m (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
189
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
190 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
191 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
192
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
193 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
194 max (const Matrix& m, double d)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
195 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
196 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
197 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
198
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
199 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
200
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
201 Matrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
202
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
203 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
204 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
205 result (i, j) = xmax (m (i, j), d);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
206
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
207 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
208 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
209
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
210 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
211 max (const Complex& c, const ComplexMatrix& m)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
212 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
213 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
214 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
215
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
216 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
217
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
218 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
219
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
220 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
221 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
222 result (i, j) = xmax (c, m (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
223
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
224 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
225 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
226
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
227 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
228 max (const ComplexMatrix& m, const Complex& c)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
229 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
230 int nr = m.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
231 int nc = m.columns ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
232
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
233 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
234
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
235 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
236
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
237 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
238 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
239 result (i, j) = xmax (m (i, j), c);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
240
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
241 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
242 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
243
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
244 static Matrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
245 max (const Matrix& a, const Matrix& b)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
246 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
247 int nr = a.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
248 int nc = a.columns ();
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
249
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
250 if (nr != b.rows () || nc != b.columns ())
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
251 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
252 error ("two-arg max expecting args of same size");
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
253 return Matrix ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
254 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
255
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
256 EMPTY_RETURN_CHECK (Matrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
257
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
258 Matrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
259
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
260 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
261 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
262 result (i, j) = xmax (a (i, j), b (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
263
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
264 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
265 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
266
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
267 static ComplexMatrix
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
268 max (const ComplexMatrix& a, const ComplexMatrix& b)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
269 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
270 int nr = a.rows ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
271 int nc = a.columns ();
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
272
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
273 if (nr != b.rows () || nc != b.columns ())
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
274 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
275 error ("two-arg max expecting args of same size");
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
276 return ComplexMatrix ();
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
277 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
278
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
279 EMPTY_RETURN_CHECK (ComplexMatrix);
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
280
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
281 ComplexMatrix result (nr, nc);
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
282
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
283 for (int j = 0; j < nc; j++)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
284 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
285 int columns_are_real_only = 1;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
286 for (int i = 0; i < nr; i++)
3243
dd00769643ae [project @ 1999-05-28 04:19:00 by jwe]
jwe
parents: 3145
diff changeset
287 if (imag (a (i, j)) != 0.0 || imag (b (i, j)) != 0.0)
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
288 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
289 columns_are_real_only = 0;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
290 break;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
291 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
292
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
293 if (columns_are_real_only)
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
294 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
295 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
296 result (i, j) = xmax (real (a (i, j)), real (b (i, j)));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
297 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
298 else
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
299 {
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
300 for (int i = 0; i < nr; i++)
3248
68259f410026 [project @ 1999-07-13 03:34:54 by jwe]
jwe
parents: 3243
diff changeset
301 result (i, j) = xmax (a (i, j), b (i, j));
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
302 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
303 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
304
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
305 return result;
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
306 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
307
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
308 #define MINMAX_BODY(FCN) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
309 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
310 octave_value_list retval; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
311 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
312 int nargin = args.length (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
313 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
314 if (nargin < 1 || nargin > 2 || nargout > 2) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
315 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
316 print_usage (#FCN); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
317 return retval; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
318 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
319 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
320 octave_value arg1; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
321 octave_value arg2; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
322 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
323 switch (nargin) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
324 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
325 case 2: \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
326 arg2 = args(1); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
327 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
328 case 1: \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
329 arg1 = args(0); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
330 break; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
331 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
332 default: \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
333 panic_impossible (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
334 break; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
335 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
336 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
337 if (nargin == 1 && (nargout == 1 || nargout == 0)) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
338 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
339 if (arg1.is_real_type ()) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
340 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
341 Matrix m = arg1.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
342 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
343 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
344 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
345 if (m.rows () == 1) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
346 retval(0) = m.row_ ## FCN (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
347 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
348 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
349 if (m.rows () == 0 || m.columns () == 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
350 retval(0) = Matrix (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
351 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
352 retval(0) = m.column_ ## FCN (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
353 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
354 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
355 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
356 else if (arg1.is_complex_type ()) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
357 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
358 ComplexMatrix m = arg1.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
359 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
360 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
361 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
362 if (m.rows () == 1) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
363 retval(0) = m.row_ ## FCN (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
364 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
365 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
366 if (m.rows () == 0 || m.columns () == 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
367 retval(0) = Matrix (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
368 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
369 retval(0) = m.column_ ## FCN (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
370 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
371 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
372 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
373 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
374 gripe_wrong_type_arg (#FCN, arg1); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
375 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
376 else if (nargin == 1 && nargout == 2) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
377 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
378 Array<int> index; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
379 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
380 if (arg1.is_real_type ()) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
381 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
382 Matrix m = arg1.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
383 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
384 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
385 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
386 retval.resize (2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
387 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
388 if (m.rows () == 1) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
389 retval(0) = m.row_ ## FCN (index); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
390 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
391 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
392 if (m.rows () == 0 || m.columns () == 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
393 retval(0) = Matrix (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
394 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
395 retval(0) = m.column_ ## FCN (index); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
396 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
397 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
398 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
399 else if (arg1.is_complex_type ()) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
400 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
401 ComplexMatrix m = arg1.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
402 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
403 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
404 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
405 retval.resize (2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
406 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
407 if (m.rows () == 1) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
408 retval(0) = m.row_ ## FCN (index); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
409 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
410 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
411 if (m.rows () == 0 || m.columns () == 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
412 retval(0) = Matrix (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
413 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
414 retval(0) = m.column_ ## FCN (index); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
415 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
416 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
417 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
418 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
419 gripe_wrong_type_arg (#FCN, arg1); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
420 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
421 int len = index.length (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
422 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
423 if (len > 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
424 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
425 RowVector idx (len); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
426 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
427 for (int i = 0; i < len; i++) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
428 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
429 int tmp = index.elem (i) + 1; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
430 idx.elem (i) = (tmp <= 0) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
431 ? octave_NaN : static_cast<double> (tmp); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
432 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
433 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
434 retval(1) = idx; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
435 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
436 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
437 retval(1) = Matrix (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
438 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
439 else if (nargin == 2) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
440 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
441 int arg1_is_scalar = arg1.is_scalar_type (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
442 int arg2_is_scalar = arg2.is_scalar_type (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
443 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
444 int arg1_is_complex = arg1.is_complex_type (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
445 int arg2_is_complex = arg2.is_complex_type (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
446 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
447 if (arg1_is_scalar) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
448 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
449 if (arg1_is_complex || arg2_is_complex) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
450 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
451 Complex c1 = arg1.complex_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
452 ComplexMatrix m2 = arg2.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
453 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
454 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
455 ComplexMatrix result = FCN (c1, m2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
456 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
457 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
458 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
459 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
460 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
461 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
462 double d1 = arg1.double_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
463 Matrix m2 = arg2.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
464 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
465 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
466 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
467 Matrix result = FCN (d1, m2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
468 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
469 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
470 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
471 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
472 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
473 else if (arg2_is_scalar) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
474 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
475 if (arg1_is_complex || arg2_is_complex) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
476 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
477 ComplexMatrix m1 = arg1.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
478 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
479 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
480 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
481 Complex c2 = arg2.complex_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
482 ComplexMatrix result = FCN (m1, c2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
483 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
484 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
485 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
486 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
487 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
488 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
489 Matrix m1 = arg1.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
490 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
491 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
492 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
493 double d2 = arg2.double_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
494 Matrix result = FCN (m1, d2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
495 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
496 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
497 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
498 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
499 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
500 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
501 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
502 if (arg1_is_complex || arg2_is_complex) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
503 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
504 ComplexMatrix m1 = arg1.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
505 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
506 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
507 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
508 ComplexMatrix m2 = arg2.complex_matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
509 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
510 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
511 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
512 ComplexMatrix result = FCN (m1, m2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
513 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
514 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
515 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
516 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
517 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
518 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
519 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
520 Matrix m1 = arg1.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
521 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
522 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
523 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
524 Matrix m2 = arg2.matrix_value (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
525 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
526 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
527 { \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
528 Matrix result = FCN (m1, m2); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
529 if (! error_state) \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
530 retval(0) = result; \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
531 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
532 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
533 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
534 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
535 } \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
536 else \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
537 panic_impossible (); \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
538 \
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
539 return retval
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
540
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
541 DEFUN_DLD (min, args, nargout,
3443
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
542 "-*- texinfo -*-\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
543 For a vector argument, return the minimum value. For a matrix\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
544 argument, return the minimum value from each column, as a row\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
545 vector. Thus,\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
546 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
547 @example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
548 min (min (@var{x}))\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
549 @end example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
550 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
551 @noindent\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
552 returns the smallest element of @var{x}.\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
553 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
554 For complex arguments, the magnitude of the elements are used for\n\
3657
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
555 comparison.\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
556 \n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
557 If called with two output arguments, also returns the index of the\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
558 minimum value(s). Thus,\n\
3775
13905c3a24af [project @ 2001-02-06 01:56:59 by jwe]
jwe
parents: 3747
diff changeset
559 @example\n\
3657
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
560 [x, ix] = min ([1, 3, 0, 2, 5])\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
561 @end example\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
562 \n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
563 @noindent\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
564 returns @var{x} = 0 and @var{ix} = 3.")
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
565 {
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
566 MINMAX_BODY (min);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
567 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
568
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
569 DEFUN_DLD (max, args, nargout,
3443
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
570 "-*- texinfo -*-\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
571 For a vector argument, return the maximum value. For a matrix\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
572 argument, return the maximum value from each column, as a row\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
573 vector. Thus,\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
574 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
575 @example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
576 max (max (@var{x}))\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
577 @end example\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
578 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
579 @noindent\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
580 returns the largest element of @var{x}.\n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
581 \n\
bf22deaf04ed [project @ 2000-01-17 08:28:27 by jwe]
jwe
parents: 3418
diff changeset
582 For complex arguments, the magnitude of the elements are used for\n\
3775
13905c3a24af [project @ 2001-02-06 01:56:59 by jwe]
jwe
parents: 3747
diff changeset
583 comparison.\n\
3657
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
584 \n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
585 If called with two output arguments, also returns the index of the\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
586 maximum value(s). Thus,\n\
3775
13905c3a24af [project @ 2001-02-06 01:56:59 by jwe]
jwe
parents: 3747
diff changeset
587 @example\n\
3657
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
588 [x, ix] = max([1, 3, 5, 2, 5])\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
589 @end example\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
590 \n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
591 @noindent\n\
a908150a3a32 [project @ 2000-04-11 19:02:03 by jwe]
jwe
parents: 3443
diff changeset
592 returns @var{x} = 5 and @var{ix} = 3.")
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
593 {
3747
50f30e40abca [project @ 2000-12-07 05:47:06 by jwe]
jwe
parents: 3657
diff changeset
594 MINMAX_BODY (max);
2928
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
595 }
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
596
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
597 /*
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
598 ;;; Local Variables: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
599 ;;; mode: C++ ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
600 ;;; End: ***
295f037b4b3e [project @ 1997-05-05 05:32:33 by jwe]
jwe
parents:
diff changeset
601 */