|
7017
|
1 ## Copyright (C) 2006, 2007 John W. Eaton |
|
7016
|
2 ## |
|
|
3 ## This file is part of Octave. |
|
|
4 ## |
|
|
5 ## Octave is free software; you can redistribute it and/or modify it |
|
|
6 ## under the terms of the GNU General Public License as published by |
|
|
7 ## the Free Software Foundation; either version 3 of the License, or (at |
|
|
8 ## your option) any later version. |
|
|
9 ## |
|
|
10 ## Octave is distributed in the hope that it will be useful, but |
|
|
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
|
|
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
|
|
13 ## General Public License for more details. |
|
|
14 ## |
|
|
15 ## You should have received a copy of the GNU General Public License |
|
|
16 ## along with Octave; see the file COPYING. If not, see |
|
|
17 ## <http://www.gnu.org/licenses/>. |
|
|
18 |
|
5590
|
19 %% Automatically generated from DejaGNU files |
|
|
20 |
|
|
21 %% test/octave.test/matrix/all-1.m |
|
|
22 %!test |
|
|
23 %! x = ones (3); |
|
|
24 %! x(1,1) = 0; |
|
|
25 %! assert((all (all (rand (3) + 1) == [1, 1, 1]) == 1 |
|
|
26 %! && all (all (x) == [0, 1, 1]) == 1 |
|
|
27 %! && all (x, 1) == [0, 1, 1] |
|
|
28 %! && all (x, 2) == [0; 1; 1])); |
|
|
29 |
|
|
30 %% test/octave.test/matrix/all-2.m |
|
5751
|
31 %!error <Invalid call to all.*> all (); |
|
5590
|
32 |
|
|
33 %% test/octave.test/matrix/all-3.m |
|
5751
|
34 %!error <Invalid call to all.*> all (1, 2, 3); |
|
5590
|
35 |
|
|
36 %% test/octave.test/matrix/any-1.m |
|
|
37 %!test |
|
|
38 %! x = zeros (3); |
|
|
39 %! x(3,3) = 1; |
|
|
40 %! assert((all (any (x) == [0, 0, 1]) == 1 |
|
|
41 %! && all (any (ones (3)) == [1, 1, 1]) == 1 |
|
|
42 %! && any (x, 1) == [0, 0, 1] |
|
|
43 %! && any (x, 2) == [0; 0; 1])); |
|
|
44 |
|
|
45 %% test/octave.test/matrix/any-2.m |
|
5751
|
46 %!error <Invalid call to any.*> any (); |
|
5590
|
47 |
|
|
48 %% test/octave.test/matrix/any-3.m |
|
5751
|
49 %!error <Invalid call to any.*> any (1, 2, 3); |
|
5590
|
50 |
|
|
51 %% test/octave.test/matrix/diff-1.m |
|
|
52 %!assert((diff ([1, 2, 3, 4]) == [1, 1, 1] |
|
|
53 %! && diff ([1, 3, 7, 19], 2) == [2, 8] |
|
|
54 %! && diff ([1, 2; 5, 4; 8, 7; 9, 6; 3, 1]) == [4, 2; 3, 3; 1, -1; -6, -5] |
|
|
55 %! && diff ([1, 2; 5, 4; 8, 7; 9, 6; 3, 1], 3) == [-1, -5; -5, 0] |
|
|
56 %! && isempty (diff (1)))); |
|
|
57 |
|
|
58 %% test/octave.test/matrix/diff-2.m |
|
|
59 %!error diff ([1, 2; 3, 4], -1); |
|
|
60 |
|
|
61 %% test/octave.test/matrix/diff-3.m |
|
|
62 %!error diff ("foo"); |
|
|
63 |
|
|
64 %% test/octave.test/matrix/diff-4.m |
|
|
65 %!error diff (); |
|
|
66 |
|
|
67 %% test/octave.test/matrix/diff-5.m |
|
|
68 %!error diff (1, 2, 3, 4); |
|
|
69 |
|
|
70 %% test/octave.test/matrix/find-1.m |
|
|
71 %!assert((find ([1, 0, 1, 0, 1]) == [1, 3, 5] |
|
|
72 %! && find ([1; 0; 3; 0; 1]) == [1; 3; 5] |
|
|
73 %! && find ([0, 0, 2; 0, 3, 0; -1, 0, 0]) == [3; 5; 7])); |
|
|
74 |
|
|
75 %% test/octave.test/matrix/find-2.m |
|
|
76 %!test |
|
|
77 %! [i, j, v] = find ([0, 0, 2; 0, 3, 0; -1, 0, 0]); |
|
|
78 %! |
|
|
79 %! assert(i == [3; 2; 1] && j == [1; 2; 3] && v == [-1; 3; 2]); |
|
|
80 |
|
|
81 %% test/octave.test/matrix/find-3.m |
|
5751
|
82 %!error <Invalid call to find.*> find (); |
|
5590
|
83 |
|
|
84 %% test/octave.test/matrix/fliplr-1.m |
|
|
85 %!assert((fliplr ([1, 2; 3, 4]) == [2, 1; 4, 3] |
|
|
86 %! && fliplr ([1, 2; 3, 4; 5, 6]) == [2, 1; 4, 3; 6, 5] |
|
|
87 %! && fliplr ([1, 2, 3; 4, 5, 6]) == [3, 2, 1; 6, 5, 4])); |
|
|
88 |
|
|
89 %% test/octave.test/matrix/fliplr-2.m |
|
|
90 %!error <usage.*fliplr> fliplr(); |
|
|
91 |
|
|
92 %% test/octave.test/matrix/fliplr-3.m |
|
|
93 %!error fliplr (1, 2); |
|
|
94 |
|
|
95 %% test/octave.test/matrix/flipud-1.m |
|
|
96 %!assert((flipud ([1, 2; 3, 4]) == [3, 4; 1, 2] |
|
|
97 %! && flipud ([1, 2; 3, 4; 5, 6]) == [5, 6; 3, 4; 1, 2] |
|
|
98 %! && flipud ([1, 2, 3; 4, 5, 6]) == [4, 5, 6; 1, 2, 3])); |
|
|
99 |
|
|
100 %% test/octave.test/matrix/flipud-2.m |
|
|
101 %!error flipud (); |
|
|
102 |
|
|
103 %% test/octave.test/matrix/flipud-3.m |
|
|
104 %!error flipud (1, 2); |
|
|
105 |
|
|
106 %% test/octave.test/matrix/rot90-1.m |
|
|
107 %!test |
|
|
108 %! x1 = [1, 2; |
|
|
109 %! 3, 4]; |
|
|
110 %! x2 = [2, 4; |
|
|
111 %! 1, 3]; |
|
|
112 %! x3 = [4, 3; |
|
|
113 %! 2, 1]; |
|
|
114 %! x4 = [3, 1; |
|
|
115 %! 4, 2]; |
|
|
116 %! |
|
|
117 %! assert((rot90 (x1)== x2 && rot90 (x1, 2) == x3 && rot90 (x1, 3) == x4 |
|
|
118 %! && rot90 (x1, 4) == x1 && rot90 (x1, 5) == x2 && rot90 (x1, -1) == x4)); |
|
|
119 |
|
|
120 %% test/octave.test/matrix/rot90-2.m |
|
|
121 %!error rot90 (); |
|
|
122 |
|
|
123 %% test/octave.test/matrix/rot90-3.m |
|
|
124 %!error rot90 (1, 2, 3); |
|
|
125 |
|
|
126 %% test/octave.test/matrix/reshape-1.m |
|
|
127 %!assert((size (reshape (rand (4, 4), 2, 8)) == [2, 8] |
|
|
128 %! && size (reshape (rand (4, 4), 8, 2)) == [8, 2] |
|
|
129 %! && size (reshape (rand (15, 4), 1, 60)) == [1, 60] |
|
|
130 %! && size (reshape (rand (15, 4), 60, 1)) == [60, 1])); |
|
|
131 |
|
|
132 %% test/octave.test/matrix/reshape-2.m |
|
|
133 %!test |
|
|
134 %! s.a = 1; |
|
|
135 %! fail("reshape (s, 2, 3)"); |
|
|
136 |
|
|
137 %% test/octave.test/matrix/reshape-3.m |
|
5751
|
138 %!error <Invalid call to reshape.*> reshape (); |
|
5590
|
139 |
|
|
140 %% test/octave.test/matrix/reshape-4.m |
|
|
141 %!error reshape (1, 2, 3, 4); |
|
|
142 |
|
|
143 %% test/octave.test/matrix/shift-1.m |
|
|
144 %!test |
|
|
145 %! a = [1, 2, 3]; |
|
|
146 %! b = [4, 5, 6]; |
|
|
147 %! c = [7, 8, 9]; |
|
|
148 %! |
|
|
149 %! r = [a, b, c]; |
|
|
150 %! m = [a; b; c]; |
|
|
151 %! |
|
|
152 %! assert((shift (r, 3) == [c, a, b] |
|
|
153 %! && shift (r, -6) == [c, a, b] |
|
|
154 %! && shift (r, -3) == [b, c, a] |
|
|
155 %! && shift (m, 1) == [c; a; b] |
|
|
156 %! && shift (m, -2) == [c; a; b])); |
|
|
157 |
|
|
158 %% test/octave.test/matrix/shift-2.m |
|
|
159 %!error shift (); |
|
|
160 |
|
|
161 %% test/octave.test/matrix/shift-3.m |
|
|
162 %!error shift (1, 2, 3, 4); |
|
|
163 |
|
|
164 %% test/octave.test/matrix/sort-1.m |
|
|
165 %!test |
|
|
166 %! a = [1, 2; 2, 3; 3, 1]; |
|
|
167 %! s = [1, 1; 2, 2; 3, 3]; |
|
|
168 %! i = [1, 3; 2, 1; 3, 2]; |
|
|
169 %! [xs, xi] = sort (a); |
|
|
170 %! assert(sort (a) == s && xs == s && xi == i); |
|
|
171 |
|
|
172 %% test/octave.test/matrix/sort-2.m |
|
5751
|
173 %!error <Invalid call to sort.*> sort (); |
|
5590
|
174 |
|
|
175 %% test/octave.test/matrix/sort-3.m |
|
5751
|
176 %!error <Invalid call to sort.*> sort (1, 2, 3, 4); |
|
5590
|
177 |
|
|
178 %% test/octave.test/matrix/tril-1.m |
|
|
179 %!test |
|
|
180 %! a = [1, 2, 3; 4, 5, 6; 7, 8, 9; 10, 11, 12]; |
|
|
181 %! |
|
|
182 %! l0 = [1, 0, 0; 4, 5, 0; 7, 8, 9; 10, 11, 12]; |
|
|
183 %! l1 = [1, 2, 0; 4, 5, 6; 7, 8, 9; 10, 11, 12]; |
|
|
184 %! l2 = [1, 2, 3; 4, 5, 6; 7, 8, 9; 10, 11, 12]; |
|
|
185 %! lm1 = [0, 0, 0; 4, 0, 0; 7, 8, 0; 10, 11, 12]; |
|
|
186 %! lm2 = [0, 0, 0; 0, 0, 0; 7, 0, 0; 10, 11, 0]; |
|
|
187 %! lm3 = [0, 0, 0; 0, 0, 0; 0, 0, 0; 10, 0, 0]; |
|
|
188 %! lm4 = [0, 0, 0; 0, 0, 0; 0, 0, 0; 0, 0, 0]; |
|
|
189 %! |
|
|
190 %! assert((tril (a, -4) == lm4 && tril (a, -3) == lm3 |
|
|
191 %! && tril (a, -2) == lm2 && tril (a, -1) == lm1 |
|
|
192 %! && tril (a) == l0 && tril (a, 1) == l1 && tril (a, 2) == l2)); |
|
|
193 |
|
|
194 %% test/octave.test/matrix/tril-2.m |
|
|
195 %!error tril (); |
|
|
196 |
|
|
197 %% test/octave.test/matrix/tril-3.m |
|
|
198 %!error tril (1, 2, 3); |
|
|
199 |
|
|
200 %% test/octave.test/matrix/triu-1.m |
|
|
201 %!test |
|
|
202 %! a = [1, 2, 3; 4, 5, 6; 7, 8, 9; 10, 11, 12]; |
|
|
203 %! |
|
|
204 %! u0 = [1, 2, 3; 0, 5, 6; 0, 0, 9; 0, 0, 0]; |
|
|
205 %! u1 = [0, 2, 3; 0, 0, 6; 0, 0, 0; 0, 0, 0]; |
|
|
206 %! u2 = [0, 0, 3; 0, 0, 0; 0, 0, 0; 0, 0, 0]; |
|
|
207 %! u3 = [0, 0, 0; 0, 0, 0; 0, 0, 0; 0, 0, 0]; |
|
|
208 %! um1 = [1, 2, 3; 4, 5, 6; 0, 8, 9; 0, 0, 12]; |
|
|
209 %! um2 = [1, 2, 3; 4, 5, 6; 7, 8, 9; 0, 11, 12]; |
|
|
210 %! um3 = [1, 2, 3; 4, 5, 6; 7, 8, 9; 10, 11, 12]; |
|
|
211 %! |
|
|
212 %! assert((triu (a, -3) == um3 && triu (a, -2) == um2 |
|
|
213 %! && triu (a, -1) == um1 && triu (a) == u0 && triu (a, 1) == u1 |
|
|
214 %! && triu (a, 2) == u2 && triu (a, 3) == u3)); |
|
|
215 |
|
|
216 %% test/octave.test/matrix/triu-2.m |
|
|
217 %!error triu (); |
|
|
218 |
|
|
219 %% test/octave.test/matrix/triu-3.m |
|
|
220 %!error triu (1, 2, 3); |
|
|
221 |
|
|
222 %% test/octave.test/matrix/vec-1.m |
|
|
223 %!assert(vec ([1, 2; 3, 4]) == [1; 3; 2; 4] && vec ([1, 3, 2, 4]) == [1; 3; 2; 4]); |
|
|
224 |
|
|
225 %% test/octave.test/matrix/vec-2.m |
|
|
226 %!error vec (); |
|
|
227 |
|
|
228 %% test/octave.test/matrix/vec-3.m |
|
|
229 %!error vec (1, 2); |
|
|
230 |
|
|
231 %% test/octave.test/matrix/vech-1.m |
|
|
232 %!assert(all (vech ([1, 2, 3; 4, 5, 6; 7, 8, 9]) == [1; 4; 7; 5; 8; 9])); |
|
|
233 |
|
|
234 %% test/octave.test/matrix/vech-2.m |
|
|
235 %!error vech (); |
|
|
236 |
|
|
237 %% test/octave.test/matrix/vech-3.m |
|
|
238 %!error vech (1, 2); |
|
|
239 |
|
|
240 %% test/octave.test/matrix/eye-1.m |
|
|
241 %!test |
|
|
242 %! i33 = [1, 0, 0; 0, 1, 0; 0, 0, 1]; |
|
|
243 %! i23 = [1, 0, 0; 0, 1, 0]; |
|
|
244 %! assert((eye (3) == i33 && eye (size (i33)) == i33 && eye (3, 3) == i33 |
|
|
245 %! && eye (2, 3) == i23 && eye (3, 2) == i23')); |
|
|
246 |
|
|
247 %% test/octave.test/matrix/eye-2.m |
|
5751
|
248 %!error <Invalid call to eye.*> eye (1, 2, 3); |
|
5590
|
249 |
|
|
250 %% test/octave.test/matrix/ones-1.m |
|
|
251 %!test |
|
|
252 %! x33 = [1, 1, 1; 1, 1, 1; 1, 1, 1]; |
|
|
253 %! x23 = [1, 1, 1; 1, 1, 1]; |
|
|
254 %! assert((ones (3) == x33 && ones (size (x33)) == x33 && ones (3, 3) == x33 |
|
|
255 %! && ones (2, 3) == x23 && ones (3, 2) == x23')); |
|
|
256 |
|
|
257 %% test/octave.test/matrix/ones-2.m |
|
|
258 %!assert(all (size (ones (3, 4, 5)) == [3, 4, 5])); |
|
|
259 |
|
|
260 %% test/octave.test/matrix/zeros-1.m |
|
|
261 %!test |
|
|
262 %! x33 = [0, 0, 0; 0, 0, 0; 0, 0, 0]; |
|
|
263 %! x23 = [0, 0, 0; 0, 0, 0]; |
|
|
264 %! assert((zeros (3) == x33 && zeros (size (x33)) == x33 && zeros (3, 3) == x33 |
|
|
265 %! && zeros (2, 3) == x23 && zeros (3, 2) == x23')); |
|
|
266 |
|
|
267 %% test/octave.test/matrix/zeros-2.m |
|
|
268 %!assert(all (size (zeros (3, 4, 5)) == [3, 4, 5])); |
|
|
269 |
|
|
270 %% test/octave.test/matrix/rand-1.m |
|
|
271 %!test |
|
|
272 %! rand ("seed", 0.5); |
|
|
273 %! r1 = rand (100); |
|
|
274 %! rand ("seed", 0.5); |
|
|
275 %! r2 = rand (100); |
|
|
276 %! assert(rand (100) < 1 && rand (100) > 0 && r1 == r2); |
|
|
277 |
|
|
278 %% test/octave.test/matrix/rand-2.m |
|
|
279 %!assert(all (size (rand (1, 2, 3)) == [1, 2, 3])); |
|
|
280 |
|
|
281 %% test/octave.test/matrix/randn-1.m |
|
|
282 %!test |
|
|
283 %! randn ("seed", 0.5); |
|
|
284 %! r1 = randn (100); |
|
|
285 %! randn ("seed", 0.5); |
|
|
286 %! r2 = randn (100); |
|
|
287 %! assert(all (all (r1 == r2))); |
|
|
288 |
|
|
289 %% test/octave.test/matrix/randn-2.m |
|
|
290 %!assert(all (size (randn (1, 2, 3)) == [1, 2, 3])); |
|
|
291 |
|
|
292 %% test/octave.test/matrix/diag-1.m |
|
|
293 %!test |
|
|
294 %! d = [1; 2; 3]; |
|
|
295 %! |
|
|
296 %! d0 = [1, 0, 0; |
|
|
297 %! 0, 2, 0; |
|
|
298 %! 0, 0, 3]; |
|
|
299 %! |
|
|
300 %! d1 = [0, 1, 0, 0; |
|
|
301 %! 0, 0, 2, 0; |
|
|
302 %! 0, 0, 0, 3; |
|
|
303 %! 0, 0, 0, 0]; |
|
|
304 %! |
|
|
305 %! d2 = [0, 0, 1, 0, 0; |
|
|
306 %! 0, 0, 0, 2, 0; |
|
|
307 %! 0, 0, 0, 0, 3; |
|
|
308 %! 0, 0, 0, 0, 0; |
|
|
309 %! 0, 0, 0, 0, 0]; |
|
|
310 %! |
|
|
311 %! dm1 = [0, 0, 0, 0; |
|
|
312 %! 1, 0, 0, 0; |
|
|
313 %! 0, 2, 0, 0; |
|
|
314 %! 0, 0, 3, 0]; |
|
|
315 %! |
|
|
316 %! dm2 = [0, 0, 0, 0, 0; |
|
|
317 %! 0, 0, 0, 0, 0; |
|
|
318 %! 1, 0, 0, 0, 0; |
|
|
319 %! 0, 2, 0, 0, 0; |
|
|
320 %! 0, 0, 3, 0, 0]; |
|
|
321 %! |
|
|
322 %! assert((diag (d) == d0 && diag (d, 1) == d1 && diag (d, 2) == d2 |
|
|
323 %! && diag (d, -1) == dm1 && diag (d, -2) == dm2 |
|
|
324 %! && diag (d0) == d && diag (d1, 1) == d && diag (dm1, -1) == d)); |
|
|
325 |
|
|
326 %% test/octave.test/matrix/diag-2.m |
|
5751
|
327 %!error <Invalid call to diag.*> diag (); |
|
5590
|
328 |
|
|
329 %% test/octave.test/matrix/diag-3.m |
|
5751
|
330 %!error <Invalid call to diag.*> diag (1, 2, 3); |
|
5590
|
331 |
|
|
332 %% test/octave.test/matrix/linspace-1.m |
|
|
333 %!test |
|
|
334 %! x1 = linspace (1, 2); |
|
|
335 %! x2 = linspace (1, 2, 10); |
|
|
336 %! x3 = linspace (1, -2, 10); |
|
|
337 %! assert((size (x1) == [1, 100] && x1(1) == 1 && x1(100) == 2 |
|
|
338 %! && size (x2) == [1, 10] && x2(1) == 1 && x2(10) == 2 |
|
|
339 %! && size (x3) == [1, 10] && x3(1) == 1 && x3(10) == -2)); |
|
|
340 |
|
|
341 %% test/octave.test/matrix/linspace-2.m |
|
|
342 %!test |
|
|
343 %! warn_fortran_indexing = 0; |
|
|
344 %! assert(all (linspace ([1, 2; 3, 4], 5, 6) == linspace (1, 5, 6))); |
|
|
345 |
|
|
346 %% test/octave.test/matrix/linspace-3.m |
|
5751
|
347 %!error <Invalid call to linspace.*> linspace (); |
|
5590
|
348 |
|
|
349 %% test/octave.test/matrix/linspace-4.m |
|
5751
|
350 %!error <Invalid call to linspace.*> linspace (1, 2, 3, 4); |
|
5590
|
351 |
|
|
352 %% test/octave.test/matrix/linspace-5.m |
|
|
353 %!test |
|
|
354 %! warn_fortran_indexing = 1; |
|
|
355 %! fail("linspace ([1, 2; 3, 4], 5, 6)","warning"); |
|
|
356 |
|
|
357 %% test/octave.test/matrix/logspace-1.m |
|
|
358 %!test |
|
|
359 %! x1 = logspace (1, 2); |
|
|
360 %! x2 = logspace (1, 2, 10); |
|
|
361 %! x3 = logspace (1, -2, 10); |
|
|
362 %! x4 = logspace (1, pi, 10); |
|
|
363 %! assert((size (x1) == [1, 50] && x1(1) == 10 && x1(50) == 100 |
|
|
364 %! && size (x2) == [1, 10] && x2(1) == 10 && x2(10) == 100 |
|
|
365 %! && size (x3) == [1, 10] && x3(1) == 10 && x3(10) == 0.01 |
|
|
366 %! && size (x4) == [1, 10] && x4(1) == 10 && abs (x4(10) - pi) < sqrt (eps))); |
|
|
367 |
|
|
368 %% test/octave.test/matrix/logspace-2.m |
|
|
369 %!error logspace ([1, 2; 3, 4], 5, 6); |
|
|
370 |
|
|
371 %% test/octave.test/matrix/logspace-3.m |
|
|
372 %!error logspace (); |
|
|
373 |
|
|
374 %% test/octave.test/matrix/logspace-4.m |
|
|
375 %!error logspace (1, 2, 3, 4); |
|
|
376 |
|
|
377 %% test/octave.test/matrix/sylvester_matrix-1.m |
|
|
378 %!assert((sylvester_matrix (1) == [1, 1; 1, -1] |
|
|
379 %! && (sylvester_matrix (2) |
|
|
380 %! == [1, 1, 1, 1; 1, -1, 1, -1; 1, 1, -1, -1; 1, -1, -1, 1]))); |
|
|
381 |
|
|
382 %% test/octave.test/matrix/sylvester_matrix-2.m |
|
|
383 %!error sylvester_matrix ([1, 2; 3, 4]); |
|
|
384 |
|
|
385 %% test/octave.test/matrix/sylvester_matrix-3.m |
|
|
386 %!error sylvester_matrix (); |
|
|
387 |
|
|
388 %% test/octave.test/matrix/sylvester_matrix-4.m |
|
|
389 %!error sylvester_matrix (1, 2); |
|
|
390 |
|
|
391 %% test/octave.test/matrix/hankel-1.m |
|
|
392 %!assert((hankel (1) == 1 && hankel ([1, 2]) == [1, 2; 2, 0] |
|
|
393 %! && hankel ([1, 2], [2; -1; -3]) == [1, 2, -1; 2, -1, -3])); |
|
|
394 |
|
|
395 %% test/octave.test/matrix/hankel-2.m |
|
|
396 %!error hankel ([1, 2; 3, 4], [1, 2; 3, 4]); |
|
|
397 |
|
|
398 %% test/octave.test/matrix/hankel-3.m |
|
|
399 %!error hankel (); |
|
|
400 |
|
|
401 %% test/octave.test/matrix/hankel-4.m |
|
|
402 %!error hankel (1, 2, 3); |
|
|
403 |
|
|
404 %% test/octave.test/matrix/hilb-1.m |
|
|
405 %!assert((hilb (2) == [1, 1/2; 1/2, 1/3] |
|
|
406 %! && hilb (3) == [1, 1/2, 1/3; 1/2, 1/3, 1/4; 1/3, 1/4, 1/5])); |
|
|
407 |
|
|
408 %% test/octave.test/matrix/hilb-2.m |
|
|
409 %!error hilb (); |
|
|
410 |
|
|
411 %% test/octave.test/matrix/hilb-3.m |
|
|
412 %!error hilb (1, 2); |
|
|
413 |
|
|
414 %% test/octave.test/matrix/invhilb-1.m |
|
|
415 %!test |
|
|
416 %! result4 = [16, -120, 240, -140; |
|
|
417 %! -120, 1200, -2700, 1680; |
|
|
418 %! 240, -2700, 6480, -4200; |
|
|
419 %! -140, 1680, -4200, 2800]; |
|
|
420 %! |
|
|
421 %! assert((invhilb (1) == 1 && invhilb (2) == [4, -6; -6, 12] |
|
|
422 %! && invhilb (4) == result4 |
|
|
423 %! && abs (invhilb (7) * hilb (7) - eye (7)) < sqrt (eps))); |
|
|
424 |
|
|
425 %% test/octave.test/matrix/invhilb-2.m |
|
|
426 %!error invhilb ([1, 2]); |
|
|
427 |
|
|
428 %% test/octave.test/matrix/invhilb-3.m |
|
|
429 %!error invhilb (); |
|
|
430 |
|
|
431 %% test/octave.test/matrix/invhilb-4.m |
|
|
432 %!error invhilb (1, 2); |
|
|
433 |
|
|
434 %% test/octave.test/matrix/toeplitz-1.m |
|
|
435 %!assert((toeplitz (1) == 1 |
|
|
436 %! && toeplitz ([1, 2, 3], [1; -3; -5]) == [1, -3, -5; 2, 1, -3; 3, 2, 1] |
|
|
437 %! && toeplitz ([1, 2, 3], [1; -3i; -5i]) == [1, -3i, -5i; 2, 1, -3i; 3, 2, 1])); |
|
|
438 |
|
|
439 %% test/octave.test/matrix/toeplitz-2.m |
|
|
440 %!error toeplitz ([1, 2; 3, 4], 1); |
|
|
441 |
|
|
442 %% test/octave.test/matrix/toeplitz-3.m |
|
|
443 %!error toeplitz (); |
|
|
444 |
|
|
445 %% test/octave.test/matrix/toeplitz-4.m |
|
|
446 %!error toeplitz (1, 2, 3); |
|
|
447 |
|
|
448 %% test/octave.test/matrix/vander-1.m |
|
|
449 %!assert((vander (1) == 1 && vander ([1, 2, 3]) == vander ([1; 2; 3]) |
|
|
450 %! && vander ([1, 2, 3]) == [1, 1, 1; 4, 2, 1; 9, 3, 1] |
|
|
451 %! && vander ([1, 2, 3]*i) == [-1, i, 1; -4, 2i, 1; -9, 3i, 1])); |
|
|
452 |
|
|
453 %% test/octave.test/matrix/vander-2.m |
|
|
454 %!error vander ([1, 2; 3, 4]); |
|
|
455 |
|
|
456 %% test/octave.test/matrix/vander-3.m |
|
|
457 %!error vander (); |
|
|
458 |
|
|
459 %% test/octave.test/matrix/vander-4.m |
|
|
460 %!error vander (1, 2); |
|
|
461 |
