annotate scripts/general/__isequal__.m @ 7016:93c65f2a5668

[project @ 2007-10-12 06:40:56 by jwe]
author jwe
date Fri, 12 Oct 2007 06:41:26 +0000
parents 6bbf56a9718a
children a1dbe9d80eee
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
1 ## Copyright (C) 2000, 2006 Paul Kienzle
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
2 ##
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
3 ## This file is part of Octave.
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
4 ##
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
5 ## Octave is free software; you can redistribute it and/or modify it
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
6 ## under the terms of the GNU General Public License as published by
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
7 ## the Free Software Foundation; either version 3 of the License, or (at
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
8 ## your option) any later version.
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
9 ##
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
10 ## Octave is distributed in the hope that it will be useful, but
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
11 ## WITHOUT ANY WARRANTY; without even the implied warranty of
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
12 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
13 ## General Public License for more details.
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
14 ##
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
15 ## 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: 6945
diff changeset
16 ## along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6945
diff changeset
17 ## <http://www.gnu.org/licenses/>.
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
18
6945
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe]
jwe
parents: 6609
diff changeset
19 ## Undocumented internal function.
6bbf56a9718a [project @ 2007-10-02 20:47:22 by jwe]
jwe
parents: 6609
diff changeset
20
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
21 ## -*- texinfo -*-
5550
815926a781f6 [project @ 2005-11-30 03:22:53 by jwe]
jwe
parents: 5549
diff changeset
22 ## @deftypefn {Function File} {} __isequal__ (@var{nans_compare_equal}, @var{x1}, @var{x2}, @dots{})
815926a781f6 [project @ 2005-11-30 03:22:53 by jwe]
jwe
parents: 5549
diff changeset
23 ## Return true if @var{x1}, @var{x2}, @dots{} are all equal and
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
24 ## @var{nans_compare_equal} evaluates to false.
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
25 ##
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
26 ## If @var{nans_compare_equal} evaluates to true, then assume NaN == NaN.
5642
2618a0750ae6 [project @ 2006-03-06 21:26:48 by jwe]
jwe
parents: 5583
diff changeset
27 ## @seealso{isequal, isequalwithequalnans}
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
28 ## @end deftypefn
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
29
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
30 ## Modified by: William Poetra Yoga Hadisoeseno
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
31
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
32 ## Algorithm:
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
33 ##
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
34 ## 1. Determine the class of x
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
35 ## 2. If x is of the struct, cell, list or char class, for each
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
36 ## argument after x, determine whether it has the same class
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
37 ## and size as x.
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
38 ## Otherwise, for each argument after x, verify that it is not
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
39 ## of the struct, cell, list or char class, and that it has
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
40 ## the same size as x.
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
41 ## 3. For each argument after x, compare it for equality with x:
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
42 ## a. struct compare each member by name, not by order (recursive)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
43 ## b. cell/list compare each member by order (recursive)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
44 ## c. char compare each member with strcmp
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
45 ## d. <other> compare each nonzero member, and assume NaN == NaN
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
46 ## if nans_compare_equal is nonzero.
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
47
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
48 function t = __isequal__ (nans_compare_equal, x, varargin)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
49
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
50 if (nargin < 3)
6046
34f96dd5441b [project @ 2006-10-10 16:10:25 by jwe]
jwe
parents: 6023
diff changeset
51 print_usage ();
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
52 endif
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
53
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
54 l_v = nargin - 2;
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
55
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
56 ## Generic tests.
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
57
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
58 ## Give an error for a list (that will make the code simpler and lists
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
59 ## are deprecated anyway.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
60 if (islist (x))
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
61 error ("__isequal__: list objects are deprecated and cannot be tested for equality here; use cell arrays instead");
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
62 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
63
6609
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
64 ## All arguments must either be of the same class or they must be
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
65 ## numeric values.
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
66 t = (all (strcmp (class(x),
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
67 cellfun (@class, varargin, "UniformOutput", false)))
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
68 || (isnumeric (x)
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
69 && all (cellfun (@isnumeric, varargin, "UniformOutput", true))));
30891d1d0c86 [project @ 2007-05-09 02:12:04 by jwe]
jwe
parents: 6157
diff changeset
70
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
71 if (t)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
72 ## Test that everything has the same number of dimensions.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
73 s_x = size (x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
74 s_v = cellfun (@size, varargin, "UniformOutput", false);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
75 t = all (length (s_x) == cellfun (@length, s_v));
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
76 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
77
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
78 if (t)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
79 ## Test that everything is the same size since it has the same
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
80 ## dimensionality.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
81 l_x = length (s_x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
82 s_v = reshape ([s_v{:}], length (s_x), []);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
83 idx = 0;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
84 while (t && idx < l_x)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
85 idx++;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
86 t = all (s_x(idx) == s_v(idx, :));
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
87 endwhile
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
88 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
89
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
90 if (t)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
91 ## Check individual classes.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
92 if (isstruct (x))
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
93 ## Test the number of fields.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
94 fn_x = fieldnames (x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
95 l_fn_x = length (fn_x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
96 fn_v = cellfun (@fieldnames, varargin, "UniformOutput", false);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
97 t = all (l_fn_x == cellfun (@length, fn_v));
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
98
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
99 ## Test that all the names are equal.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
100 idx = 0;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
101 s_fn_x = sort (fn_x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
102 while (t && idx < l_v)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
103 idx++;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
104 ## We'll allow the fieldnames to be in a different order.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
105 t = all (strcmp (s_fn_x, sort (fn_v{idx})));
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
106 endwhile
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
107
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
108 idx = 0;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
109 while (t && idx < l_fn_x)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
110 ## Test that all field values are equal.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
111 idx++;
6157
045038e0108a [project @ 2006-11-13 22:22:53 by jwe]
jwe
parents: 6130
diff changeset
112 args = {nans_compare_equal, {x.(fn_x{idx})}};
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
113 for argn = 1:l_v
6157
045038e0108a [project @ 2006-11-13 22:22:53 by jwe]
jwe
parents: 6130
diff changeset
114 args{argn+2} = {varargin{argn}.(fn_x{idx})};
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
115 endfor
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
116 ## Minimize function calls by calling for all the arguments at
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
117 ## once.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
118 t = __isequal__ (args{:});
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
119 endwhile
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
120
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
121 elseif (iscell (x))
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
122 ## Check that each element of a cell is equal.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
123 l_x = numel (x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
124 idx = 0;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
125 while (t && idx < l_x)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
126 idx++;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
127 args = {nans_compare_equal, x{idx}};
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
128 for p = 1:l_v
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
129 args{p+2} = varargin{p}{idx};
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
130 endfor
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
131 t = __isequal__ (args{:});
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
132 endwhile
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
133
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
134 elseif (ischar (x))
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
135
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
136 ## Sizes are equal already, so we can just make everything into a
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
137 ## row and test the rows.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
138 for i = 1:l_v
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
139 strings{i} = reshape (varargin{i}, 1, []);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
140 endfor
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
141 t = all (strcmp (reshape (x, 1, []), strings));
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
142
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
143 else
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
144 ## Check the numeric types.
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
145
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
146 if (issparse (x))
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
147 f_x = spfind (x);
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
148 else
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
149 f_x = find (x);
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
150 endif
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
151 l_f_x = length (f_x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
152 x = x(f_x);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
153 for argn = 1:l_v
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
154 y = varargin{argn};
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
155 if (issparse (y))
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
156 f_y = spfind (y);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
157 else
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
158 f_y = find (y);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
159 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
160
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
161 t = (l_f_x == length (f_y)) && all (f_x == f_y);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
162 if (!t)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
163 return;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
164 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
165
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
166 y = y(f_y);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
167 m = (x == y);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
168 t = all (m);
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
169
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
170 if (!t)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
171 if (nans_compare_equal)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
172 t = isnan (x(!m)) && isnan (y(!m));
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
173 else
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
174 return;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
175 endif
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
176 endif
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
177 endfor
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
178
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
179 endif
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
180 endif
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
181
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
182 endfunction
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
183
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
184 ## test size and shape
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
185 %!assert(__isequal__(0,[1,2,3,4],[1,2,3,4]), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
186 %!assert(__isequal__(0,[1;2;3;4],[1;2;3;4]), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
187 %!assert(__isequal__(0,[1,2,3,4],[1;2;3;4]), false)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
188 %!assert(__isequal__(0,[1,2,3,4],[1,2;3,4]), false)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
189 %!assert(__isequal__(0,[1,2,3,4],[1,3;2,4]), false)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
190
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
191 %!test
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
192 %! A = 1:8;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
193 %! B = reshape (A, 2, 2, 2);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
194 %! assert (__isequal__ (0, A, B), false);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
195
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
196 %!test
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
197 %! A = reshape (1:8, 2, 2, 2);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
198 %! B = A;
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
199 %! assert (__isequal__ (0, A, B), true);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
200
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
201 %!test
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
202 %! A = reshape (1:8, 2, 4);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
203 %! B = reshape (A, 2, 2, 2);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
204 %! assert (__isequal__ (0, A, B), false);
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
205
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
206 ## test for equality
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
207 %!assert(__isequal__(0,[1,2,3,4],[1,2,3,4]), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
208 %!assert(__isequal__(1,{1,2,NaN,4},{1,2,NaN,4}), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
209 %!assert(__isequal__(1,[1,2,NaN,4],[1,2,NaN,4]), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
210 %!assert(__isequal__(0,['a','b','c','d'],['a','b','c','d']), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
211 ## Test multi-line strings
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
212 %!assert(__isequal__(0,["test";"strings"],["test";"strings"],["test";"strings"]), true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
213 ## test for inequality
5549
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
214 %!assert(__isequal__(0,[1,2,3,4],[1;2;3;4]),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
215 %!assert(__isequal__(0,{1,2,3,4},[1,2,3,4]),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
216 %!assert(__isequal__(0,[1,2,3,4],{1,2,3,4}),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
217 %!assert(__isequal__(0,[1,2,NaN,4],[1,2,NaN,4]),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
218 %!assert(__isequal__(1,[1,2,NaN,4],[1,NaN,3,4]),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
219 %!assert(__isequal__(1,[1,2,NaN,4],[1,2,3,4]),false)
6db3a5df1eab [project @ 2005-11-30 03:15:19 by jwe]
jwe
parents:
diff changeset
220 %!assert(__isequal__(0,['a','b','c','d'],['a';'b';'c';'d']),false)
6130
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
221 %!assert(__isequal__(0,{'a','b','c','d'},{'a';'b';'c';'d'}),false)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
222 ## test for equality (struct)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
223 %!assert(__isequal__(0,struct('a',1,'b',2),struct('a',1,'b',2)),true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
224 %!assert(__isequal__(0,struct('a',1,'b',2),struct('a',1,'b',2),struct('a',1,'b',2)),true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
225 %!assert(__isequal__(0,struct('a','abc','b',2),struct('a','abc','b',2)),true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
226 %!assert(__isequal__(1,struct('a',NaN,'b',2),struct('a',NaN,'b',2),struct('a',NaN,'b',2)),true)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
227 ## test for inequality (struct)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
228 %!assert(__isequal__(0,struct('a',NaN,'b',2),struct('a',NaN,'b',2),struct('a',NaN,'b',2)),false)
6eba20084f8f [project @ 2006-11-01 16:54:04 by jwe]
jwe
parents: 6046
diff changeset
229