annotate liboctave/util/lo-array-gripes.cc @ 16779:8fce0ed4894a

Specialize is_empty and numel methods for sparse matrices (debian bug #706376) * ov-base.h (virtual bool is_empty (void) const) : Make method virtual * ov-base-sparse.h (bool is_empty (void) const)) : Declare new method (octave_idx_type numel (void) const): New method. * ov-base-sparse.cc (template <class T> bool octave_base_sparse<T>:is_empty (void) const)) : Define new method * ov-bool-sparse.h (idx_vector index_vector (void) const): Start specialization for sparse matrices * lo-array-gripes.cc (void gripe_index_value (void)): Clarify error message * data.cc : Add test to isempty * idx-vector.cc (idx_vector::idx_vector_rep::idx_vector_rep (const Sparse<bool>&) : Fix sparse indexing bug * __sprand_impl__.m : Don't use randparm when there are more than sizemax() elements. * sprand.m, sprandn.m : Add a test for large, low density matrices
author David Bateman <dbateman@free.fr>
date Thu, 20 Jun 2013 02:17:25 +0200
parents 648dabbb4c6b
children d63878346099
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
11135
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
1 /*
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
2
14138
72c96de7a403 maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents: 11590
diff changeset
3 Copyright (C) 2003-2012 John W. Eaton
11135
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
4 Copyright (C) 2009 VZLU Prague
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
5
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
6 This file is part of Octave.
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
7
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
8 Octave is free software; you can redistribute it and/or modify it
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
9 under the terms of the GNU General Public License as published by the
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
10 Free Software Foundation; either version 3 of the License, or (at your
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
11 option) any later version.
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
12
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
13 Octave is distributed in the hope that it will be useful, but WITHOUT
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
16 for more details.
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
17
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
18 You should have received a copy of the GNU General Public License
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
19 along with Octave; see the file COPYING. If not, see
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
20 <http://www.gnu.org/licenses/>.
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
21
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
22 */
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
23
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
24 #ifdef HAVE_CONFIG_H
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
25 #include <config.h>
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
26 #endif
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
27
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
28 #include "lo-array-gripes.h"
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
29 #include "lo-error.h"
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
30
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
31 const char *error_id_nonconformant_args = "Octave:nonconformant-args";
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
32
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
33 const char *error_id_index_out_of_bounds = "Octave:index-out-of-bounds";
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
34
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
35 const char *error_id_invalid_index = "Octave:invalid-index";
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
36
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
37 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
38 gripe_nan_to_logical_conversion (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
39 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
40 (*current_liboctave_error_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
41 ("invalid conversion from NaN to logical");
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
42 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
43
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
44 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
45 gripe_nan_to_character_conversion (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
46 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
47 (*current_liboctave_error_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
48 ("invalid conversion from NaN to character");
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
49 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
50
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
51 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
52 gripe_nonconformant (const char *op, octave_idx_type op1_len,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
53 octave_idx_type op2_len)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
54 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
55 const char *err_id = error_id_nonconformant_args;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
56
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
57 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
58 (err_id, "%s: nonconformant arguments (op1 len: %d, op2 len: %d)",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
59 op, op1_len, op2_len);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
60 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
61
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
62 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
63 gripe_nonconformant (const char *op,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
64 octave_idx_type op1_nr, octave_idx_type op1_nc,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
65 octave_idx_type op2_nr, octave_idx_type op2_nc)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
66 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
67 const char *err_id = error_id_nonconformant_args;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
68
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
69 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
70 (err_id, "%s: nonconformant arguments (op1 is %dx%d, op2 is %dx%d)",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
71 op, op1_nr, op1_nc, op2_nr, op2_nc);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
72 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
73
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
74 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
75 gripe_nonconformant (const char *op, const dim_vector& op1_dims,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
76 const dim_vector& op2_dims)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
77 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
78 const char *err_id = error_id_nonconformant_args;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
79
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
80 std::string op1_dims_str = op1_dims.str ();
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
81 std::string op2_dims_str = op2_dims.str ();
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
82
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
83 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
84 (err_id, "%s: nonconformant arguments (op1 is %s, op2 is %s)",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
85 op, op1_dims_str.c_str (), op2_dims_str.c_str ());
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
86 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
87
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
88 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
89 gripe_index_out_of_range (int nd, int dim, octave_idx_type idx,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
90 octave_idx_type ext)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
91 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
92 const char *err_id = error_id_index_out_of_bounds;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
93
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
94 switch (nd)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
95 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
96 case 1:
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
97 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
98 (err_id, "A(I): index out of bounds; value %d out of bound %d",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
99 idx, ext);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
100 break;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
101
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
102 case 2:
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
103 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
104 (err_id, "A(I,J): %s index out of bounds; value %d out of bound %d",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
105 (dim == 1) ? "row" : "column", idx, ext);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
106 break;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
107
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
108 default:
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
109 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
110 (err_id, "A(I,J,...): index to dimension %d out of bounds; value %d out of bound %d",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
111 dim, idx, ext);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
112 break;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
113 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
114 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
115
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
116 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
117 gripe_del_index_out_of_range (bool is1d, octave_idx_type idx,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
118 octave_idx_type ext)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
119 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
120 const char *err_id = error_id_index_out_of_bounds;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
121
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
122 (*current_liboctave_error_with_id_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
123 (err_id, "A(%s) = []: index out of bounds; value %d out of bound %d",
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
124 is1d ? "I" : "..,I,..", idx, ext);
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
125 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
126
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
127 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
128 gripe_invalid_index (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
129 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
130 const char *err_id = error_id_invalid_index;
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
131
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
132 (*current_liboctave_error_with_id_handler)
16779
8fce0ed4894a Specialize is_empty and numel methods for sparse matrices (debian bug #706376)
David Bateman <dbateman@free.fr>
parents: 15271
diff changeset
133 #ifdef USE_64_BIT_IDX_T
8fce0ed4894a Specialize is_empty and numel methods for sparse matrices (debian bug #706376)
David Bateman <dbateman@free.fr>
parents: 15271
diff changeset
134 (err_id, "subscript indices must be either positive integers less than 2^63 or logicals");
8fce0ed4894a Specialize is_empty and numel methods for sparse matrices (debian bug #706376)
David Bateman <dbateman@free.fr>
parents: 15271
diff changeset
135 #else
8fce0ed4894a Specialize is_empty and numel methods for sparse matrices (debian bug #706376)
David Bateman <dbateman@free.fr>
parents: 15271
diff changeset
136 (err_id, "subscript indices must be either positive integers less than 2^31 or logicals");
8fce0ed4894a Specialize is_empty and numel methods for sparse matrices (debian bug #706376)
David Bateman <dbateman@free.fr>
parents: 15271
diff changeset
137 #endif
11135
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
138 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
139
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
140 // FIXME -- the following is a common error message to resize,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
141 // regardless of whether it's called from assign or elsewhere. It
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
142 // seems OK to me, but eventually the gripe can be specialized.
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
143 // Anyway, propagating various error messages into procedure is, IMHO,
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
144 // a nonsense. If anything, we should change error handling here (and
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
145 // throughout liboctave) to allow custom handling of errors
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
146
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
147 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
148 gripe_invalid_resize (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
149 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
150 (*current_liboctave_error_with_id_handler)
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11523
diff changeset
151 ("Octave:invalid-resize",
11590
4ced6b90fffb style fixes for warning and error messages in source files
John W. Eaton <jwe@octave.org>
parents: 11586
diff changeset
152 "Invalid resizing operation or ambiguous assignment to an out-of-bounds array element");
11135
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
153 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
154
11586
12df7854fa7c strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents: 11523
diff changeset
155 void
11135
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
156 gripe_invalid_assignment_size (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
157 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
158 (*current_liboctave_error_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
159 ("A(I) = X: X must have the same size as I");
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
160 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
161
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
162 void
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
163 gripe_assignment_dimension_mismatch (void)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
164 {
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
165 (*current_liboctave_error_handler)
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
166 ("A(I,J,...) = X: dimensions mismatch");
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
167 }
047b0e877a14 add files for previous change
John W. Eaton <jwe@octave.org>
parents:
diff changeset
168