1
|
1 /* |
|
2 |
1827
|
3 Copyright (C) 1996 John W. Eaton |
1
|
4 |
|
5 This file is part of Octave. |
|
6 |
|
7 Octave is free software; you can redistribute it and/or modify it |
|
8 under the terms of the GNU General Public License as published by the |
|
9 Free Software Foundation; either version 2, or (at your option) any |
|
10 later version. |
|
11 |
|
12 Octave is distributed in the hope that it will be useful, but WITHOUT |
|
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
|
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License |
|
15 for more details. |
|
16 |
|
17 You should have received a copy of the GNU General Public License |
|
18 along with Octave; see the file COPYING. If not, write to the Free |
1315
|
19 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
1
|
20 |
|
21 */ |
|
22 |
1299
|
23 #if defined (__GNUG__) |
|
24 #pragma implementation |
|
25 #endif |
|
26 |
240
|
27 #ifdef HAVE_CONFIG_H |
1192
|
28 #include <config.h> |
1
|
29 #endif |
|
30 |
1558
|
31 #include <cctype> |
1355
|
32 #include <cstring> |
|
33 |
1728
|
34 #include <string> |
|
35 |
1558
|
36 #include <fstream.h> |
581
|
37 #include <iostream.h> |
|
38 |
1742
|
39 #include <SLList.h> |
|
40 |
2181
|
41 #include "Array-flags.h" |
|
42 |
1558
|
43 #include "mx-base.h" |
|
44 #include "Range.h" |
1755
|
45 #include "str-vec.h" |
1558
|
46 |
|
47 #include "arith-ops.h" |
2181
|
48 #include "defun.h" |
1
|
49 #include "error.h" |
|
50 #include "gripes.h" |
1558
|
51 #include "idx-vector.h" |
1742
|
52 #include "mappers.h" |
747
|
53 #include "oct-map.h" |
1742
|
54 #include "oct-obj.h" |
1355
|
55 #include "pager.h" |
1558
|
56 #include "pr-output.h" |
|
57 #include "sysdep.h" |
1742
|
58 #include "pt-const.h" |
2124
|
59 #include "pt-walk.h" |
1558
|
60 #include "unwind-prot.h" |
1430
|
61 #include "utils.h" |
1558
|
62 #include "variables.h" |
|
63 |
2086
|
64 #ifndef OCT_VAL_REP |
|
65 #define OCT_VAL_REP octave_value::octave_value_rep |
1558
|
66 #endif |
|
67 |
|
68 #ifndef MAX |
|
69 #define MAX(a,b) ((a) > (b) ? (a) : (b)) |
|
70 #endif |
|
71 |
2086
|
72 #ifndef OCT_VAL_REP |
|
73 #define OCT_VAL_REP octave_value::octave_value_rep |
1558
|
74 #endif |
|
75 |
|
76 #ifndef MAX |
|
77 #define MAX(a,b) ((a) > (b) ? (a) : (b)) |
|
78 #endif |
|
79 |
|
80 // The following three variables could be made static members of the |
2086
|
81 // OCT_VAL_REP class. |
1558
|
82 |
|
83 // Pointer to the blocks of memory we manage. |
2086
|
84 static OCT_VAL_REP *tc_rep_newlist = 0; |
1558
|
85 |
|
86 // Multiplier for allocating new blocks. |
|
87 static const int tc_rep_newlist_grow_size = 128; |
|
88 |
2181
|
89 // If TRUE, allow assignments like |
|
90 // |
|
91 // octave> A(1) = 3; A(2) = 5 |
|
92 // |
|
93 // for A already defined and a matrix type. |
|
94 static bool Vdo_fortran_indexing; |
|
95 |
|
96 // Should we allow things like: |
|
97 // |
|
98 // octave> 'abc' + 0 |
|
99 // 97 98 99 |
|
100 // |
|
101 // to happen? A positive value means yes. A negative value means |
|
102 // yes, but print a warning message. Zero means it should be |
|
103 // considered an error. |
|
104 int Vimplicit_str_to_num_ok; |
|
105 |
|
106 // Should we allow silent conversion of complex to real when a real |
|
107 // type is what we're really looking for? A positive value means yes. |
|
108 // A negative value means yes, but print a warning message. Zero |
|
109 // means it should be considered an error. |
|
110 static int Vok_to_lose_imaginary_part; |
|
111 |
|
112 // If TRUE, create column vectors when doing assignments like: |
|
113 // |
|
114 // octave> A(1) = 3; A(2) = 5 |
|
115 // |
|
116 // (for A undefined). Only matters when resize_on_range_error is also |
|
117 // TRUE. |
|
118 static bool Vprefer_column_vectors; |
|
119 |
|
120 // If TRUE, prefer logical (zore-one) indexing over normal indexing |
|
121 // when there is a conflice. For example, given a = [2, 3], the |
|
122 // expression a ([1, 1]) would return [2 3] (instead of [2 2], which |
|
123 // would be returned if prefer_zero_one_indxing were FALSE). |
|
124 static bool Vprefer_zero_one_indexing; |
|
125 |
2187
|
126 // If TRUE, print the name along with the value. |
|
127 static bool Vprint_answer_id_name; |
|
128 |
2181
|
129 // Should operations on empty matrices return empty matrices or an |
|
130 // error? A positive value means yes. A negative value means yes, |
|
131 // but print a warning message. Zero means it should be considered an |
|
132 // error. |
|
133 int Vpropagate_empty_matrices; |
|
134 |
|
135 // If TRUE, resize matrices when performing and indexed assignment and |
|
136 // the indices are outside the current bounds. |
|
137 bool Vresize_on_range_error; |
|
138 |
|
139 // How many levels of structure elements should we print? |
|
140 static int Vstruct_levels_to_print; |
|
141 |
1558
|
142 // Indentation level for structures. |
1973
|
143 static int struct_indent = 0; |
1558
|
144 |
|
145 static void |
1973
|
146 increment_struct_indent (void) |
1558
|
147 { |
1973
|
148 struct_indent += 2; |
1558
|
149 } |
|
150 |
|
151 static void |
1973
|
152 decrement_struct_indent (void) |
1558
|
153 { |
1973
|
154 struct_indent -= 2; |
1558
|
155 } |
|
156 |
2338
|
157 // XXX FIXME XXX -- these should be member functions. |
|
158 |
1827
|
159 static bool |
1558
|
160 any_element_is_complex (const ComplexMatrix& a) |
|
161 { |
|
162 int nr = a.rows (); |
|
163 int nc = a.columns (); |
1827
|
164 |
1558
|
165 for (int j = 0; j < nc; j++) |
|
166 for (int i = 0; i < nr; i++) |
2305
|
167 if (imag (a (i, j)) != 0.0) |
1827
|
168 return true; |
|
169 |
|
170 return false; |
1558
|
171 } |
747
|
172 |
2338
|
173 static bool |
|
174 any_element_is_complex (const ComplexDiagMatrix& d) |
|
175 { |
|
176 int len = d.length (); |
|
177 |
|
178 for (int i = 0; i < len; i++) |
|
179 if (imag (d (i, i)) != 0.0) |
|
180 return true; |
|
181 |
|
182 return false; |
|
183 } |
|
184 |
|
185 static bool |
|
186 any_element_is_complex (const ComplexRowVector& v) |
|
187 { |
|
188 int len = v.length (); |
|
189 |
|
190 for (int i = 0; i < len; i++) |
|
191 if (imag (v (i)) != 0.0) |
|
192 return true; |
|
193 |
|
194 return false; |
|
195 } |
|
196 |
|
197 static bool |
|
198 any_element_is_complex (const ComplexColumnVector& v) |
|
199 { |
|
200 int len = v.length (); |
|
201 |
|
202 for (int i = 0; i < len; i++) |
|
203 if (imag (v (i)) != 0.0) |
|
204 return true; |
|
205 |
|
206 return false; |
|
207 } |
|
208 |
1168
|
209 // The following three variables could be made static members of the |
2086
|
210 // octave_value class. |
1168
|
211 |
|
212 // Pointer to the blocks of memory we manage. |
2086
|
213 static octave_value *tc_newlist = 0; |
1168
|
214 |
|
215 // Multiplier for allocating new blocks. |
1299
|
216 static const int tc_newlist_grow_size = 128; |
1168
|
217 |
747
|
218 Octave_map |
2086
|
219 octave_value::map_value (void) const |
747
|
220 { |
|
221 return rep->map_value (); |
|
222 } |
1
|
223 |
2086
|
224 octave_value::~octave_value (void) |
1
|
225 { |
|
226 #if defined (MDEBUG) |
2086
|
227 cerr << "~octave_value: rep: " << rep |
1
|
228 << " rep->count: " << rep->count << "\n"; |
|
229 #endif |
|
230 |
|
231 if (--rep->count <= 0) |
|
232 { |
|
233 delete rep; |
531
|
234 rep = 0; |
1
|
235 } |
|
236 } |
|
237 |
|
238 void * |
2086
|
239 octave_value::operator new (size_t size) |
1
|
240 { |
2086
|
241 assert (size == sizeof (octave_value)); |
1168
|
242 |
1299
|
243 if (! tc_newlist) |
|
244 { |
2086
|
245 int block_size = tc_newlist_grow_size * sizeof (octave_value); |
|
246 tc_newlist = (octave_value *) new char [block_size]; |
1299
|
247 |
1321
|
248 int i = 0; |
|
249 |
|
250 for (i = 0; i < tc_newlist_grow_size - 1; i++) |
1299
|
251 tc_newlist[i].freeptr = &tc_newlist[i+1]; |
|
252 |
|
253 tc_newlist[i].freeptr = 0; |
|
254 } |
|
255 |
2086
|
256 octave_value *tmp = tc_newlist; |
1299
|
257 tc_newlist = tc_newlist->freeptr; |
1168
|
258 return tmp; |
1
|
259 } |
|
260 |
|
261 void |
2086
|
262 octave_value::operator delete (void *p, size_t /* size */) |
1
|
263 { |
2086
|
264 octave_value *tmp = (octave_value *) p; |
1299
|
265 tmp->freeptr = tc_newlist; |
|
266 tc_newlist = tmp; |
1
|
267 } |
|
268 |
747
|
269 // Simple assignment. |
|
270 |
2086
|
271 octave_value |
|
272 octave_value::operator = (const octave_value& a) |
747
|
273 { |
|
274 if (rep != a.rep) |
|
275 { |
|
276 if (--rep->count <= 0) |
|
277 delete rep; |
|
278 rep = a.rep; |
|
279 rep->count++; |
|
280 } |
|
281 return *this; |
|
282 } |
|
283 |
2086
|
284 octave_value |
|
285 octave_value::lookup_map_element (const string& ref, bool insert, |
1827
|
286 bool silent) |
1277
|
287 { |
2086
|
288 octave_value retval; |
1277
|
289 |
1755
|
290 if (! ref.empty ()) |
1277
|
291 { |
1755
|
292 SLList<string> list; |
|
293 |
|
294 size_t beg = 0; |
|
295 size_t end; |
|
296 |
1277
|
297 do |
|
298 { |
1755
|
299 end = ref.find ('.', beg); |
|
300 |
|
301 string tmp = (end == NPOS) |
2246
|
302 ? ref.substr (beg) : ref.substr (beg, end - beg); |
1755
|
303 |
|
304 list.append (tmp); |
1277
|
305 } |
1755
|
306 while (end != NPOS && (beg = end + 1)); |
1277
|
307 |
|
308 retval = lookup_map_element (list, insert, silent); |
|
309 } |
|
310 |
|
311 return retval; |
|
312 } |
|
313 |
2086
|
314 octave_value |
|
315 octave_value::lookup_map_element (SLList<string>& list, bool insert, |
1827
|
316 bool silent) |
747
|
317 { |
2086
|
318 octave_value retval; |
|
319 |
|
320 octave_value_rep *tmp_rep = rep; |
747
|
321 |
|
322 Pix p = list.first (); |
|
323 while (p) |
|
324 { |
1755
|
325 string elt = list (p); |
747
|
326 |
|
327 list.next (p); |
|
328 |
2086
|
329 octave_value tmp; |
1277
|
330 |
|
331 tmp = tmp_rep->lookup_map_element (elt, insert, silent); |
747
|
332 |
|
333 if (error_state) |
|
334 break; |
|
335 |
|
336 tmp_rep = tmp.rep; |
|
337 |
|
338 if (! p) |
|
339 retval = tmp; |
|
340 } |
|
341 |
|
342 return retval; |
|
343 } |
|
344 |
1199
|
345 void |
2086
|
346 octave_value::print (void) |
1199
|
347 { |
2095
|
348 print (octave_stdout); |
1199
|
349 } |
|
350 |
1755
|
351 void |
2086
|
352 octave_value::print_with_name (const string& name, bool print_padding) |
1755
|
353 { |
2095
|
354 print_with_name (octave_stdout, name, print_padding); |
1755
|
355 } |
|
356 |
|
357 void |
2086
|
358 octave_value::print_with_name (ostream& output_buf, const string& name, |
2095
|
359 bool print_padding) |
1755
|
360 { |
1827
|
361 bool pad_after = false; |
1755
|
362 |
2187
|
363 if (Vprint_answer_id_name) |
1755
|
364 { |
1973
|
365 if (print_as_scalar ()) |
1755
|
366 output_buf << name << " = "; |
1973
|
367 else if (print_as_structure ()) |
|
368 { |
|
369 pad_after = true; |
|
370 output_buf << name << " ="; |
|
371 } |
1755
|
372 else |
|
373 { |
1827
|
374 pad_after = true; |
1755
|
375 output_buf << name << " =\n\n"; |
|
376 } |
|
377 } |
|
378 |
|
379 print (output_buf); |
|
380 |
|
381 if (print_padding && pad_after) |
|
382 output_buf << "\n"; |
|
383 } |
|
384 |
747
|
385 // Simple structure assignment. |
|
386 |
|
387 void |
2086
|
388 octave_value::make_unique (void) |
747
|
389 { |
|
390 if (rep->count > 1) |
|
391 { |
|
392 --rep->count; |
2086
|
393 rep = new octave_value_rep (*rep); |
747
|
394 rep->count = 1; |
|
395 } |
|
396 |
|
397 if (rep->is_map ()) |
|
398 { |
|
399 for (Pix p = rep->a_map->first (); p != 0; rep->a_map->next (p)) |
|
400 { |
|
401 rep->a_map->contents (p) . make_unique (); |
|
402 } |
|
403 } |
|
404 } |
|
405 |
2086
|
406 octave_value::octave_value_rep * |
|
407 octave_value::make_unique_map (void) |
747
|
408 { |
|
409 if (! rep->is_map ()) |
|
410 { |
|
411 if (--rep->count <= 0) |
|
412 delete rep; |
|
413 |
|
414 Octave_map m; |
2086
|
415 rep = new octave_value_rep (m); |
747
|
416 rep->count = 1; |
|
417 } |
|
418 |
|
419 make_unique (); |
|
420 |
|
421 return rep; |
|
422 } |
|
423 |
2086
|
424 octave_value |
|
425 octave_value::assign_map_element (SLList<string>& list, |
|
426 octave_value& rhs) |
747
|
427 { |
2086
|
428 octave_value_rep *tmp_rep = make_unique_map (); |
747
|
429 |
|
430 if (rhs.is_map ()) |
|
431 rhs.make_unique (); |
|
432 |
|
433 Pix p = list.first (); |
|
434 while (p) |
|
435 { |
1755
|
436 string elt = list (p); |
747
|
437 |
|
438 list.next (p); |
|
439 |
2086
|
440 octave_value& tmp = tmp_rep->lookup_map_element (elt, 1); |
747
|
441 |
|
442 if (! p) |
|
443 { |
|
444 tmp = rhs; |
|
445 return tmp; |
|
446 } |
|
447 |
|
448 tmp_rep = tmp.make_unique_map (); |
|
449 } |
|
450 |
2086
|
451 return octave_value (); |
747
|
452 } |
|
453 |
|
454 // Indexed structure assignment. |
|
455 |
2086
|
456 octave_value |
|
457 octave_value::assign_map_element (SLList<string>& list, |
|
458 octave_value& rhs, |
|
459 const octave_value_list& args) |
747
|
460 { |
2086
|
461 octave_value_rep *tmp_rep = make_unique_map (); |
747
|
462 |
|
463 if (rhs.is_map ()) |
|
464 rhs.make_unique (); |
|
465 |
|
466 Pix p = list.first (); |
|
467 while (p) |
|
468 { |
1755
|
469 string elt = list (p); |
747
|
470 |
|
471 list.next (p); |
|
472 |
2086
|
473 octave_value& tmp = tmp_rep->lookup_map_element (elt, 1); |
747
|
474 |
|
475 if (! p) |
|
476 { |
|
477 tmp.assign (rhs, args); |
|
478 return tmp; |
|
479 } |
|
480 |
|
481 tmp_rep = tmp.make_unique_map (); |
|
482 } |
|
483 |
2086
|
484 return octave_value (); |
747
|
485 } |
|
486 |
2086
|
487 octave_value_list |
|
488 octave_value::eval (bool print, int, const octave_value_list& args) |
1742
|
489 { |
2086
|
490 octave_value_list retval; |
1742
|
491 |
|
492 if (args.length () > 0) |
|
493 retval(0) = rep->do_index (args); |
|
494 else |
|
495 retval(0) = *this; |
|
496 |
|
497 if (retval(0).is_defined ()) |
|
498 retval(0).eval (print); |
|
499 |
|
500 return retval; |
|
501 } |
|
502 |
747
|
503 void |
2124
|
504 octave_value::accept (tree_walker& tw) |
747
|
505 { |
2124
|
506 tw.visit_octave_value (*this); |
747
|
507 } |
|
508 |
1299
|
509 // The real representation of constants. |
|
510 |
2086
|
511 OCT_VAL_REP::octave_value_rep (void) |
1299
|
512 { |
|
513 type_tag = unknown_constant; |
|
514 } |
|
515 |
2086
|
516 OCT_VAL_REP::octave_value_rep (double d) |
1299
|
517 { |
|
518 scalar = d; |
|
519 type_tag = scalar_constant; |
|
520 } |
|
521 |
2338
|
522 // XXX FIXME XXX -- perhaps these constructors should just do the |
|
523 // obvious thing and then call maybe_mutate() instead of duplicating |
|
524 // most of that logic several times here... |
|
525 |
2086
|
526 OCT_VAL_REP::octave_value_rep (const Matrix& m) |
1299
|
527 { |
|
528 if (m.rows () == 1 && m.columns () == 1) |
|
529 { |
2305
|
530 scalar = m (0, 0); |
1299
|
531 type_tag = scalar_constant; |
|
532 } |
|
533 else |
|
534 { |
|
535 matrix = new Matrix (m); |
|
536 type_tag = matrix_constant; |
|
537 } |
|
538 } |
|
539 |
2086
|
540 OCT_VAL_REP::octave_value_rep (const DiagMatrix& d) |
1299
|
541 { |
|
542 if (d.rows () == 1 && d.columns () == 1) |
|
543 { |
2305
|
544 scalar = d (0, 0); |
1299
|
545 type_tag = scalar_constant; |
|
546 } |
|
547 else |
|
548 { |
|
549 matrix = new Matrix (d); |
|
550 type_tag = matrix_constant; |
|
551 } |
|
552 } |
|
553 |
2086
|
554 OCT_VAL_REP::octave_value_rep (const RowVector& v, int prefer_column_vector) |
1299
|
555 { |
|
556 int len = v.capacity (); |
|
557 if (len == 1) |
|
558 { |
2305
|
559 scalar = v (0); |
1299
|
560 type_tag = scalar_constant; |
|
561 } |
|
562 else |
|
563 { |
|
564 int pcv = (prefer_column_vector < 0) |
2181
|
565 ? Vprefer_column_vectors |
1299
|
566 : prefer_column_vector; |
|
567 |
|
568 if (pcv) |
|
569 { |
|
570 Matrix m (len, 1); |
|
571 for (int i = 0; i < len; i++) |
2305
|
572 m (i, 0) = v (i); |
1299
|
573 matrix = new Matrix (m); |
|
574 type_tag = matrix_constant; |
|
575 } |
|
576 else |
|
577 { |
|
578 Matrix m (1, len); |
|
579 for (int i = 0; i < len; i++) |
2305
|
580 m (0, i) = v (i); |
1299
|
581 matrix = new Matrix (m); |
|
582 type_tag = matrix_constant; |
|
583 } |
|
584 } |
|
585 } |
|
586 |
2086
|
587 OCT_VAL_REP::octave_value_rep (const ColumnVector& v, int prefer_column_vector) |
1299
|
588 { |
|
589 int len = v.capacity (); |
|
590 if (len == 1) |
|
591 { |
2305
|
592 scalar = v (0); |
1299
|
593 type_tag = scalar_constant; |
|
594 } |
|
595 else |
|
596 { |
|
597 int pcv = (prefer_column_vector < 0) |
2181
|
598 ? Vprefer_column_vectors |
1299
|
599 : prefer_column_vector; |
|
600 |
|
601 if (pcv) |
|
602 { |
|
603 Matrix m (len, 1); |
|
604 for (int i = 0; i < len; i++) |
2305
|
605 m (i, 0) = v (i); |
1299
|
606 matrix = new Matrix (m); |
|
607 type_tag = matrix_constant; |
|
608 } |
|
609 else |
|
610 { |
|
611 Matrix m (1, len); |
|
612 for (int i = 0; i < len; i++) |
2305
|
613 m (0, i) = v (i); |
1299
|
614 matrix = new Matrix (m); |
|
615 type_tag = matrix_constant; |
|
616 } |
|
617 } |
|
618 } |
|
619 |
2086
|
620 OCT_VAL_REP::octave_value_rep (const Complex& c) |
1299
|
621 { |
1615
|
622 if (::imag (c) == 0.0) |
|
623 { |
|
624 scalar = ::real (c); |
|
625 type_tag = scalar_constant; |
|
626 } |
|
627 else |
|
628 { |
|
629 complex_scalar = new Complex (c); |
|
630 type_tag = complex_scalar_constant; |
|
631 } |
1299
|
632 } |
|
633 |
2086
|
634 OCT_VAL_REP::octave_value_rep (const ComplexMatrix& m) |
1299
|
635 { |
|
636 if (m.rows () == 1 && m.columns () == 1) |
|
637 { |
2305
|
638 Complex c = m (0, 0); |
1615
|
639 |
|
640 if (::imag (c) == 0.0) |
|
641 { |
|
642 scalar = ::real (c); |
|
643 type_tag = scalar_constant; |
|
644 } |
|
645 else |
|
646 { |
|
647 complex_scalar = new Complex (c); |
|
648 type_tag = complex_scalar_constant; |
|
649 } |
1299
|
650 } |
2338
|
651 else if (! any_element_is_complex (*complex_matrix)) |
|
652 { |
|
653 matrix = new Matrix (::real (m)); |
|
654 type_tag = matrix_constant; |
|
655 } |
1299
|
656 else |
|
657 { |
|
658 complex_matrix = new ComplexMatrix (m); |
|
659 type_tag = complex_matrix_constant; |
|
660 } |
|
661 } |
|
662 |
2086
|
663 OCT_VAL_REP::octave_value_rep (const ComplexDiagMatrix& d) |
1299
|
664 { |
|
665 if (d.rows () == 1 && d.columns () == 1) |
|
666 { |
2305
|
667 Complex c = d (0, 0); |
1615
|
668 |
|
669 if (::imag (c) == 0.0) |
|
670 { |
|
671 scalar = ::real (c); |
|
672 type_tag = scalar_constant; |
|
673 } |
|
674 else |
|
675 { |
|
676 complex_scalar = new Complex (c); |
|
677 type_tag = complex_scalar_constant; |
|
678 } |
1299
|
679 } |
2338
|
680 else if (! any_element_is_complex (d)) |
|
681 { |
|
682 matrix = new Matrix (::real (d)); |
|
683 type_tag = matrix_constant; |
|
684 } |
1299
|
685 else |
|
686 { |
|
687 complex_matrix = new ComplexMatrix (d); |
|
688 type_tag = complex_matrix_constant; |
|
689 } |
|
690 } |
|
691 |
2086
|
692 OCT_VAL_REP::octave_value_rep (const ComplexRowVector& v, |
1299
|
693 int prefer_column_vector) |
|
694 { |
|
695 int len = v.capacity (); |
|
696 if (len == 1) |
|
697 { |
2305
|
698 Complex c = v (0); |
1615
|
699 |
|
700 if (::imag (c) == 0.0) |
|
701 { |
|
702 scalar = ::real (c); |
|
703 type_tag = scalar_constant; |
|
704 } |
|
705 else |
|
706 { |
|
707 complex_scalar = new Complex (c); |
|
708 type_tag = complex_scalar_constant; |
|
709 } |
1299
|
710 } |
|
711 else |
|
712 { |
|
713 int pcv = (prefer_column_vector < 0) |
2181
|
714 ? Vprefer_column_vectors |
1299
|
715 : prefer_column_vector; |
|
716 |
|
717 if (pcv) |
|
718 { |
2338
|
719 if (! any_element_is_complex (v)) |
|
720 { |
|
721 Matrix m (len, 1); |
|
722 for (int i = 0; i < len; i++) |
|
723 m (i, 0) = ::real (v (i)); |
|
724 matrix = new Matrix (m); |
|
725 type_tag = matrix_constant; |
|
726 } |
|
727 else |
|
728 { |
|
729 ComplexMatrix m (len, 1); |
|
730 for (int i = 0; i < len; i++) |
|
731 m (i, 0) = v (i); |
|
732 complex_matrix = new ComplexMatrix (m); |
|
733 type_tag = complex_matrix_constant; |
|
734 } |
1299
|
735 } |
|
736 else |
|
737 { |
2338
|
738 if (! any_element_is_complex (v)) |
|
739 { |
|
740 Matrix m (len, 1); |
|
741 for (int i = 0; i < len; i++) |
|
742 m (0, i) = ::real (v (i)); |
|
743 matrix = new Matrix (m); |
|
744 type_tag = matrix_constant; |
|
745 } |
|
746 else |
|
747 { |
|
748 ComplexMatrix m (1, len); |
|
749 for (int i = 0; i < len; i++) |
|
750 m (0, i) = v (i); |
|
751 complex_matrix = new ComplexMatrix (m); |
|
752 type_tag = complex_matrix_constant; |
|
753 } |
1299
|
754 } |
|
755 } |
|
756 } |
|
757 |
2086
|
758 OCT_VAL_REP::octave_value_rep (const ComplexColumnVector& v, int |
1299
|
759 prefer_column_vector) |
|
760 { |
|
761 int len = v.capacity (); |
|
762 if (len == 1) |
|
763 { |
2305
|
764 Complex c = v (0); |
1615
|
765 |
|
766 if (::imag (c) == 0.0) |
|
767 { |
|
768 scalar = ::real (c); |
|
769 type_tag = scalar_constant; |
|
770 } |
|
771 else |
|
772 { |
|
773 complex_scalar = new Complex (c); |
|
774 type_tag = complex_scalar_constant; |
|
775 } |
1299
|
776 } |
|
777 else |
|
778 { |
|
779 int pcv = (prefer_column_vector < 0) |
2181
|
780 ? Vprefer_column_vectors |
1299
|
781 : prefer_column_vector; |
|
782 |
|
783 if (pcv) |
|
784 { |
2338
|
785 if (! any_element_is_complex (v)) |
|
786 { |
|
787 Matrix m (len, 1); |
|
788 for (int i = 0; i < len; i++) |
|
789 m (i, 0) = ::real (v (i)); |
|
790 matrix = new Matrix (m); |
|
791 type_tag = matrix_constant; |
|
792 } |
|
793 else |
|
794 { |
|
795 ComplexMatrix m (len, 1); |
|
796 for (int i = 0; i < len; i++) |
|
797 m (i, 0) = v (i); |
|
798 complex_matrix = new ComplexMatrix (m); |
|
799 type_tag = complex_matrix_constant; |
|
800 } |
1299
|
801 } |
|
802 else |
|
803 { |
2338
|
804 if (! any_element_is_complex (v)) |
|
805 { |
|
806 Matrix m (len, 1); |
|
807 for (int i = 0; i < len; i++) |
|
808 m (0, i) = ::real (v (i)); |
|
809 matrix = new Matrix (m); |
|
810 type_tag = matrix_constant; |
|
811 } |
|
812 else |
|
813 { |
|
814 ComplexMatrix m (1, len); |
|
815 for (int i = 0; i < len; i++) |
|
816 m (0, i) = v (i); |
|
817 complex_matrix = new ComplexMatrix (m); |
|
818 type_tag = complex_matrix_constant; |
|
819 } |
1299
|
820 } |
|
821 } |
|
822 } |
|
823 |
2086
|
824 OCT_VAL_REP::octave_value_rep (const char *s) |
1299
|
825 { |
1572
|
826 char_matrix = new charMatrix (s); |
|
827 type_tag = char_matrix_constant_str; |
1355
|
828 } |
|
829 |
2086
|
830 OCT_VAL_REP::octave_value_rep (const string& s) |
1731
|
831 { |
|
832 char_matrix = new charMatrix (s); |
|
833 type_tag = char_matrix_constant_str; |
1755
|
834 } |
|
835 |
2086
|
836 OCT_VAL_REP::octave_value_rep (const string_vector& s) |
1755
|
837 { |
|
838 int nr = s.length (); |
|
839 int nc = s.max_length (); |
|
840 char_matrix = new charMatrix (nr, nc, 0); |
|
841 for (int i = 0; i < nr; i++) |
|
842 { |
|
843 nc = s[i].length (); |
|
844 for (int j = 0; j < nc; j++) |
2305
|
845 (*char_matrix) (i, j) = s[i][j]; |
1755
|
846 } |
|
847 type_tag = char_matrix_constant_str; |
1731
|
848 } |
|
849 |
2086
|
850 OCT_VAL_REP::octave_value_rep (const charMatrix& chm, bool is_str) |
1355
|
851 { |
1572
|
852 char_matrix = new charMatrix (chm); |
|
853 type_tag = is_str ? char_matrix_constant_str : char_matrix_constant; |
1299
|
854 } |
|
855 |
2086
|
856 OCT_VAL_REP::octave_value_rep (double b, double l, double i) |
1299
|
857 { |
|
858 range = new Range (b, l, i); |
|
859 int nel = range->nelem (); |
|
860 if (nel > 1) |
|
861 type_tag = range_constant; |
|
862 else |
|
863 { |
|
864 delete range; |
|
865 if (nel == 1) |
|
866 { |
|
867 scalar = b; |
|
868 type_tag = scalar_constant; |
|
869 } |
|
870 else if (nel == 0) |
|
871 { |
|
872 matrix = new Matrix (); |
|
873 type_tag = matrix_constant; |
|
874 } |
|
875 else |
|
876 { |
|
877 type_tag = unknown_constant; |
|
878 if (nel == -1) |
|
879 ::error ("number of elements in range exceeds INT_MAX"); |
|
880 else |
|
881 ::error ("invalid range"); |
|
882 } |
|
883 } |
|
884 } |
|
885 |
2086
|
886 OCT_VAL_REP::octave_value_rep (const Range& r) |
1299
|
887 { |
|
888 int nel = r.nelem (); |
|
889 if (nel > 1) |
|
890 { |
|
891 range = new Range (r); |
|
892 type_tag = range_constant; |
|
893 } |
|
894 else if (nel == 1) |
|
895 { |
|
896 scalar = r.base (); |
|
897 type_tag = scalar_constant; |
|
898 } |
|
899 else if (nel == 0) |
|
900 { |
|
901 matrix = new Matrix (); |
|
902 type_tag = matrix_constant; |
|
903 } |
|
904 else |
|
905 { |
|
906 type_tag = unknown_constant; |
|
907 if (nel == -1) |
|
908 ::error ("number of elements in range exceeds INT_MAX"); |
|
909 else |
|
910 ::error ("invalid range"); |
|
911 } |
|
912 } |
|
913 |
2086
|
914 OCT_VAL_REP::octave_value_rep (const Octave_map& m) |
1299
|
915 { |
|
916 a_map = new Octave_map (m); |
|
917 type_tag = map_constant; |
|
918 } |
|
919 |
2086
|
920 OCT_VAL_REP::octave_value_rep (OCT_VAL_REP::constant_type t) |
1299
|
921 { |
|
922 assert (t == magic_colon || t == all_va_args); |
|
923 type_tag = t; |
|
924 } |
|
925 |
2086
|
926 OCT_VAL_REP::octave_value_rep (const octave_value_rep& t) |
1299
|
927 { |
|
928 type_tag = t.type_tag; |
|
929 |
|
930 switch (t.type_tag) |
|
931 { |
|
932 case unknown_constant: |
|
933 break; |
|
934 |
|
935 case scalar_constant: |
|
936 scalar = t.scalar; |
|
937 break; |
|
938 |
|
939 case matrix_constant: |
|
940 matrix = new Matrix (*(t.matrix)); |
|
941 break; |
|
942 |
1572
|
943 case char_matrix_constant: |
|
944 char_matrix = new charMatrix (*(t.char_matrix)); |
|
945 break; |
|
946 |
|
947 case char_matrix_constant_str: |
|
948 char_matrix = new charMatrix (*(t.char_matrix)); |
1299
|
949 break; |
|
950 |
|
951 case complex_matrix_constant: |
|
952 complex_matrix = new ComplexMatrix (*(t.complex_matrix)); |
|
953 break; |
|
954 |
|
955 case complex_scalar_constant: |
|
956 complex_scalar = new Complex (*(t.complex_scalar)); |
|
957 break; |
|
958 |
|
959 case range_constant: |
|
960 range = new Range (*(t.range)); |
|
961 break; |
|
962 |
|
963 case map_constant: |
|
964 a_map = new Octave_map (*(t.a_map)); |
|
965 break; |
|
966 |
|
967 case magic_colon: |
|
968 case all_va_args: |
|
969 break; |
|
970 } |
|
971 |
1755
|
972 orig_text = t.orig_text; |
1299
|
973 } |
|
974 |
2086
|
975 OCT_VAL_REP::~octave_value_rep (void) |
1299
|
976 { |
|
977 switch (type_tag) |
|
978 { |
|
979 case matrix_constant: |
|
980 delete matrix; |
|
981 break; |
|
982 |
|
983 case complex_scalar_constant: |
|
984 delete complex_scalar; |
|
985 break; |
|
986 |
|
987 case complex_matrix_constant: |
|
988 delete complex_matrix; |
|
989 break; |
|
990 |
1572
|
991 case char_matrix_constant: |
|
992 case char_matrix_constant_str: |
|
993 delete char_matrix; |
1299
|
994 break; |
|
995 |
|
996 case range_constant: |
|
997 delete range; |
|
998 break; |
|
999 |
|
1000 case map_constant: |
|
1001 delete a_map; |
|
1002 break; |
|
1003 |
|
1004 case unknown_constant: |
|
1005 case scalar_constant: |
|
1006 case magic_colon: |
|
1007 case all_va_args: |
|
1008 break; |
|
1009 } |
|
1010 } |
|
1011 |
|
1012 void * |
2086
|
1013 OCT_VAL_REP::operator new (size_t size) |
1299
|
1014 { |
2086
|
1015 assert (size == sizeof (OCT_VAL_REP)); |
1299
|
1016 |
|
1017 if (! tc_rep_newlist) |
|
1018 { |
2086
|
1019 int block_size = tc_rep_newlist_grow_size * sizeof (OCT_VAL_REP); |
|
1020 tc_rep_newlist = (OCT_VAL_REP *) new char [block_size]; |
1299
|
1021 |
1321
|
1022 int i = 0; |
|
1023 |
|
1024 for (i = 0; i < tc_rep_newlist_grow_size - 1; i++) |
1299
|
1025 tc_rep_newlist[i].freeptr = &tc_rep_newlist[i+1]; |
|
1026 |
|
1027 tc_rep_newlist[i].freeptr = 0; |
|
1028 } |
|
1029 |
2086
|
1030 OCT_VAL_REP *tmp = tc_rep_newlist; |
1299
|
1031 tc_rep_newlist = tc_rep_newlist->freeptr; |
|
1032 return tmp; |
|
1033 } |
|
1034 |
|
1035 void |
2086
|
1036 OCT_VAL_REP::operator delete (void *p, size_t /* size */) |
1299
|
1037 { |
2086
|
1038 OCT_VAL_REP *tmp = (OCT_VAL_REP *) p; |
1299
|
1039 tmp->freeptr = tc_rep_newlist; |
|
1040 tc_rep_newlist = tmp; |
|
1041 } |
|
1042 |
|
1043 int |
2086
|
1044 OCT_VAL_REP::rows (void) const |
1299
|
1045 { |
|
1046 int retval = -1; |
|
1047 |
|
1048 switch (type_tag) |
|
1049 { |
|
1050 case scalar_constant: |
|
1051 case complex_scalar_constant: |
|
1052 retval = 1; |
|
1053 break; |
|
1054 |
1572
|
1055 case char_matrix_constant: |
|
1056 case char_matrix_constant_str: |
|
1057 retval = char_matrix->rows (); |
1355
|
1058 break; |
|
1059 |
1299
|
1060 case range_constant: |
|
1061 retval = (columns () > 0); |
|
1062 break; |
|
1063 |
|
1064 case matrix_constant: |
|
1065 retval = matrix->rows (); |
|
1066 break; |
|
1067 |
|
1068 case complex_matrix_constant: |
|
1069 retval = complex_matrix->rows (); |
|
1070 break; |
|
1071 |
|
1072 default: |
|
1073 break; |
|
1074 } |
|
1075 |
|
1076 return retval; |
|
1077 } |
|
1078 |
|
1079 int |
2086
|
1080 OCT_VAL_REP::columns (void) const |
1299
|
1081 { |
|
1082 int retval = -1; |
|
1083 |
|
1084 switch (type_tag) |
|
1085 { |
|
1086 case scalar_constant: |
|
1087 case complex_scalar_constant: |
|
1088 retval = 1; |
|
1089 break; |
|
1090 |
|
1091 case matrix_constant: |
|
1092 retval = matrix->columns (); |
|
1093 break; |
|
1094 |
|
1095 case complex_matrix_constant: |
|
1096 retval = complex_matrix->columns (); |
|
1097 break; |
|
1098 |
1572
|
1099 case char_matrix_constant: |
|
1100 case char_matrix_constant_str: |
|
1101 retval = char_matrix->columns (); |
1299
|
1102 break; |
|
1103 |
|
1104 case range_constant: |
|
1105 retval = range->nelem (); |
|
1106 break; |
|
1107 |
|
1108 default: |
|
1109 break; |
|
1110 } |
|
1111 |
|
1112 return retval; |
|
1113 } |
|
1114 |
2086
|
1115 octave_value |
|
1116 OCT_VAL_REP::all (void) const |
1299
|
1117 { |
2086
|
1118 octave_value retval; |
1299
|
1119 |
|
1120 if (error_state) |
|
1121 return retval; |
|
1122 |
|
1123 if (! is_numeric_type ()) |
|
1124 { |
2086
|
1125 octave_value tmp = make_numeric (); |
1299
|
1126 |
|
1127 if (error_state) |
|
1128 return retval; |
|
1129 |
|
1130 return tmp.all (); |
|
1131 } |
|
1132 |
|
1133 switch (type_tag) |
|
1134 { |
|
1135 case scalar_constant: |
1355
|
1136 retval = (double) (scalar != 0.0); |
1299
|
1137 break; |
|
1138 |
|
1139 case matrix_constant: |
1355
|
1140 retval = matrix->all (); |
1299
|
1141 break; |
|
1142 |
|
1143 case complex_scalar_constant: |
1355
|
1144 retval = (double) (*complex_scalar != 0.0); |
1299
|
1145 break; |
|
1146 |
|
1147 case complex_matrix_constant: |
1355
|
1148 retval = complex_matrix->all (); |
1299
|
1149 break; |
|
1150 |
|
1151 default: |
|
1152 gripe_wrong_type_arg ("all", *this); |
|
1153 break; |
|
1154 } |
|
1155 |
|
1156 return retval; |
|
1157 } |
|
1158 |
2086
|
1159 octave_value |
|
1160 OCT_VAL_REP::any (void) const |
1299
|
1161 { |
2086
|
1162 octave_value retval; |
1299
|
1163 |
|
1164 if (error_state) |
|
1165 return retval; |
|
1166 |
|
1167 if (! is_numeric_type ()) |
|
1168 { |
2086
|
1169 octave_value tmp = make_numeric (); |
1299
|
1170 |
|
1171 if (error_state) |
|
1172 return retval; |
|
1173 |
|
1174 return tmp.any (); |
|
1175 } |
|
1176 |
|
1177 switch (type_tag) |
|
1178 { |
|
1179 case scalar_constant: |
1355
|
1180 retval = (double) (scalar != 0.0); |
1299
|
1181 break; |
|
1182 |
|
1183 case matrix_constant: |
1355
|
1184 retval = matrix->any (); |
1299
|
1185 break; |
|
1186 |
|
1187 case complex_scalar_constant: |
1355
|
1188 retval = (double) (*complex_scalar != 0.0); |
1299
|
1189 break; |
|
1190 |
|
1191 case complex_matrix_constant: |
1355
|
1192 retval = complex_matrix->any (); |
1299
|
1193 break; |
|
1194 |
|
1195 default: |
|
1196 gripe_wrong_type_arg ("any", *this); |
|
1197 break; |
|
1198 } |
|
1199 |
|
1200 return retval; |
|
1201 } |
|
1202 |
1827
|
1203 bool |
2086
|
1204 OCT_VAL_REP::valid_as_scalar_index (void) const |
1299
|
1205 { |
|
1206 return (type_tag == magic_colon |
|
1207 || (type_tag == scalar_constant |
|
1208 && ! xisnan (scalar) |
|
1209 && NINT (scalar) == 1) |
|
1210 || (type_tag == range_constant |
|
1211 && range->nelem () == 1 |
|
1212 && ! xisnan (range->base ()) |
|
1213 && NINT (range->base ()) == 1)); |
|
1214 } |
|
1215 |
1827
|
1216 bool |
2086
|
1217 OCT_VAL_REP::valid_as_zero_index (void) const |
1299
|
1218 { |
|
1219 return ((type_tag == scalar_constant |
|
1220 && ! xisnan (scalar) |
|
1221 && NINT (scalar) == 0) |
|
1222 || (type_tag == matrix_constant |
|
1223 && matrix->rows () == 0 |
|
1224 && matrix->columns () == 0) |
|
1225 || (type_tag == range_constant |
|
1226 && range->nelem () == 1 |
|
1227 && ! xisnan (range->base ()) |
|
1228 && NINT (range->base ()) == 0)); |
|
1229 } |
|
1230 |
1827
|
1231 bool |
2086
|
1232 OCT_VAL_REP::is_true (void) const |
1299
|
1233 { |
1827
|
1234 int retval = false; |
1299
|
1235 |
|
1236 if (error_state) |
|
1237 return retval; |
|
1238 |
|
1239 if (! is_numeric_type ()) |
|
1240 { |
2086
|
1241 octave_value tmp = make_numeric (); |
1299
|
1242 |
|
1243 if (error_state) |
|
1244 return retval; |
|
1245 |
|
1246 return tmp.is_true (); |
|
1247 } |
|
1248 |
|
1249 switch (type_tag) |
|
1250 { |
|
1251 case scalar_constant: |
|
1252 retval = (scalar != 0.0); |
|
1253 break; |
|
1254 |
|
1255 case matrix_constant: |
|
1256 { |
|
1257 Matrix m = (matrix->all ()) . all (); |
|
1258 retval = (m.rows () == 1 |
|
1259 && m.columns () == 1 |
2305
|
1260 && m (0, 0) != 0.0); |
1299
|
1261 } |
|
1262 break; |
|
1263 |
|
1264 case complex_scalar_constant: |
|
1265 retval = (*complex_scalar != 0.0); |
|
1266 break; |
|
1267 |
|
1268 case complex_matrix_constant: |
|
1269 { |
|
1270 Matrix m = (complex_matrix->all ()) . all (); |
|
1271 retval = (m.rows () == 1 |
|
1272 && m.columns () == 1 |
2305
|
1273 && m (0, 0) != 0.0); |
1299
|
1274 } |
|
1275 break; |
|
1276 |
|
1277 default: |
|
1278 gripe_wrong_type_arg (0, *this); |
|
1279 break; |
|
1280 } |
|
1281 |
|
1282 return retval; |
|
1283 } |
|
1284 |
|
1285 static void |
|
1286 warn_implicit_conversion (const char *from, const char *to) |
|
1287 { |
|
1288 warning ("implicit conversion from %s to %s", from, to); |
|
1289 } |
|
1290 |
2067
|
1291 // XXX FIXME XXX -- we need a better way of handling conversions. |
|
1292 |
1299
|
1293 double |
2086
|
1294 OCT_VAL_REP::double_value (bool force_string_conv) const |
1299
|
1295 { |
|
1296 double retval = octave_NaN; |
|
1297 |
|
1298 switch (type_tag) |
|
1299 { |
|
1300 case scalar_constant: |
|
1301 retval = scalar; |
|
1302 break; |
|
1303 |
|
1304 case matrix_constant: |
|
1305 { |
2181
|
1306 if (Vdo_fortran_indexing && rows () > 0 && columns () > 0) |
2305
|
1307 retval = (*matrix) (0, 0); |
1299
|
1308 else |
|
1309 gripe_invalid_conversion ("real matrix", "real scalar"); |
|
1310 } |
|
1311 break; |
|
1312 |
|
1313 case complex_matrix_constant: |
|
1314 case complex_scalar_constant: |
|
1315 { |
2181
|
1316 int flag = Vok_to_lose_imaginary_part; |
1299
|
1317 |
|
1318 if (flag < 0) |
|
1319 warn_implicit_conversion ("complex scalar", "real scalar"); |
|
1320 |
|
1321 if (flag) |
|
1322 { |
|
1323 if (type_tag == complex_scalar_constant) |
|
1324 retval = ::real (*complex_scalar); |
|
1325 else if (type_tag == complex_matrix_constant) |
|
1326 { |
2181
|
1327 if (Vdo_fortran_indexing |
1299
|
1328 && rows () > 0 && columns () > 0) |
2305
|
1329 retval = ::real ((*complex_matrix) (0, 0)); |
1299
|
1330 else |
|
1331 gripe_invalid_conversion ("complex matrix", "real scalar"); |
|
1332 } |
|
1333 else |
|
1334 panic_impossible (); |
|
1335 } |
|
1336 else |
|
1337 gripe_invalid_conversion ("complex scalar", "real scalar"); |
|
1338 } |
|
1339 break; |
|
1340 |
1572
|
1341 case char_matrix_constant: |
1299
|
1342 { |
1572
|
1343 int len = char_matrix->rows (); |
|
1344 if ((char_matrix->rows () == 1 && len == 1) |
2181
|
1345 || (len > 1 && Vdo_fortran_indexing)) |
2305
|
1346 retval = toascii ((int) (*char_matrix) (0, 0)); |
1572
|
1347 else |
|
1348 gripe_invalid_conversion ("char matrix", "real scalar"); |
|
1349 } |
|
1350 break; |
|
1351 |
|
1352 case char_matrix_constant_str: |
|
1353 { |
1827
|
1354 int flag = force_string_conv; |
1299
|
1355 if (! flag) |
2181
|
1356 flag = Vimplicit_str_to_num_ok; |
1299
|
1357 |
|
1358 if (flag < 0) |
|
1359 warn_implicit_conversion ("string", "real scalar"); |
|
1360 |
1572
|
1361 int len = char_matrix->rows (); |
1355
|
1362 if (flag |
1572
|
1363 && ((char_matrix->rows () == 1 && len == 1) |
2181
|
1364 || (len > 1 && Vdo_fortran_indexing))) |
2305
|
1365 retval = toascii ((int) (*char_matrix) (0, 0)); |
1299
|
1366 else |
|
1367 gripe_invalid_conversion ("string", "real scalar"); |
|
1368 } |
|
1369 break; |
|
1370 |
|
1371 case range_constant: |
|
1372 { |
|
1373 int nel = range->nelem (); |
2181
|
1374 if (nel == 1 || (nel > 1 && Vdo_fortran_indexing)) |
1299
|
1375 retval = range->base (); |
|
1376 else |
|
1377 gripe_invalid_conversion ("range", "real scalar"); |
|
1378 } |
|
1379 break; |
|
1380 |
|
1381 default: |
|
1382 gripe_invalid_conversion (type_as_string (), "real scalar"); |
|
1383 break; |
|
1384 } |
|
1385 |
|
1386 return retval; |
|
1387 } |
|
1388 |
|
1389 Matrix |
2086
|
1390 OCT_VAL_REP::matrix_value (bool force_string_conv) const |
1299
|
1391 { |
|
1392 Matrix retval; |
|
1393 |
|
1394 switch (type_tag) |
|
1395 { |
|
1396 case scalar_constant: |
|
1397 retval = Matrix (1, 1, scalar); |
|
1398 break; |
|
1399 |
|
1400 case matrix_constant: |
|
1401 retval = *matrix; |
|
1402 break; |
|
1403 |
|
1404 case complex_scalar_constant: |
|
1405 case complex_matrix_constant: |
|
1406 { |
2181
|
1407 int flag = Vok_to_lose_imaginary_part; |
1299
|
1408 if (flag < 0) |
|
1409 warn_implicit_conversion ("complex matrix", "real matrix"); |
|
1410 |
|
1411 if (flag) |
|
1412 { |
|
1413 if (type_tag == complex_scalar_constant) |
|
1414 retval = Matrix (1, 1, ::real (*complex_scalar)); |
|
1415 else if (type_tag == complex_matrix_constant) |
|
1416 retval = ::real (*complex_matrix); |
|
1417 else |
|
1418 panic_impossible (); |
|
1419 } |
|
1420 else |
|
1421 gripe_invalid_conversion ("complex matrix", "real matrix"); |
|
1422 } |
|
1423 break; |
|
1424 |
1572
|
1425 case char_matrix_constant: |
|
1426 retval = Matrix (*char_matrix); |
|
1427 break; |
|
1428 |
|
1429 case char_matrix_constant_str: |
1299
|
1430 { |
1827
|
1431 int flag = force_string_conv; |
1299
|
1432 if (! flag) |
2181
|
1433 flag = Vimplicit_str_to_num_ok; |
1299
|
1434 |
|
1435 if (flag < 0) |
|
1436 warn_implicit_conversion ("string", "real matrix"); |
|
1437 |
|
1438 if (flag) |
1572
|
1439 retval = Matrix (*char_matrix); |
1299
|
1440 else |
|
1441 gripe_invalid_conversion ("string", "real matrix"); |
|
1442 } |
|
1443 break; |
|
1444 |
|
1445 case range_constant: |
|
1446 retval = range->matrix_value (); |
|
1447 break; |
|
1448 |
|
1449 default: |
|
1450 gripe_invalid_conversion (type_as_string (), "real matrix"); |
|
1451 break; |
|
1452 } |
|
1453 |
|
1454 return retval; |
|
1455 } |
|
1456 |
|
1457 Complex |
2086
|
1458 OCT_VAL_REP::complex_value (bool force_string_conv) const |
1299
|
1459 { |
|
1460 Complex retval (octave_NaN, octave_NaN); |
|
1461 |
|
1462 switch (type_tag) |
|
1463 { |
|
1464 case complex_scalar_constant: |
|
1465 retval = *complex_scalar; |
|
1466 break; |
|
1467 |
|
1468 case scalar_constant: |
|
1469 retval = scalar; |
|
1470 break; |
|
1471 |
|
1472 case complex_matrix_constant: |
|
1473 case matrix_constant: |
|
1474 { |
2181
|
1475 if (Vdo_fortran_indexing && rows () > 0 && columns () > 0) |
1299
|
1476 { |
|
1477 if (type_tag == complex_matrix_constant) |
2305
|
1478 retval = (*complex_matrix) (0, 0); |
1299
|
1479 else |
2305
|
1480 retval = (*matrix) (0, 0); |
1299
|
1481 } |
|
1482 else |
|
1483 gripe_invalid_conversion ("real matrix", "real scalar"); |
|
1484 } |
|
1485 break; |
|
1486 |
1572
|
1487 case char_matrix_constant: |
1299
|
1488 { |
1572
|
1489 int len = char_matrix->cols (); |
|
1490 if ((char_matrix->rows () == 1 && len == 1) |
2181
|
1491 || (len > 1 && Vdo_fortran_indexing)) |
2305
|
1492 retval = toascii ((int) (*char_matrix) (0, 0)); |
1572
|
1493 else |
|
1494 gripe_invalid_conversion ("char matrix", "complex scalar"); |
|
1495 } |
|
1496 break; |
|
1497 |
|
1498 case char_matrix_constant_str: |
|
1499 { |
1827
|
1500 int flag = force_string_conv; |
1299
|
1501 if (! flag) |
2181
|
1502 flag = Vimplicit_str_to_num_ok; |
1299
|
1503 |
|
1504 if (flag < 0) |
|
1505 warn_implicit_conversion ("string", "complex scalar"); |
|
1506 |
1572
|
1507 int len = char_matrix->cols (); |
1355
|
1508 if (flag |
1572
|
1509 && ((char_matrix->rows () == 1 && len == 1) |
2181
|
1510 || (len > 1 && Vdo_fortran_indexing))) |
2305
|
1511 retval = toascii ((int) (*char_matrix) (0, 0)); |
1299
|
1512 else |
|
1513 gripe_invalid_conversion ("string", "complex scalar"); |
|
1514 } |
|
1515 break; |
|
1516 |
|
1517 case range_constant: |
|
1518 { |
|
1519 int nel = range->nelem (); |
2181
|
1520 if (nel == 1 || (nel > 1 && Vdo_fortran_indexing)) |
1299
|
1521 retval = range->base (); |
|
1522 else |
|
1523 gripe_invalid_conversion ("range", "complex scalar"); |
|
1524 } |
|
1525 break; |
|
1526 |
|
1527 default: |
|
1528 gripe_invalid_conversion (type_as_string (), "complex scalar"); |
|
1529 break; |
|
1530 } |
|
1531 |
|
1532 return retval; |
|
1533 } |
|
1534 |
|
1535 ComplexMatrix |
2086
|
1536 OCT_VAL_REP::complex_matrix_value (bool force_string_conv) const |
1299
|
1537 { |
|
1538 ComplexMatrix retval; |
|
1539 |
|
1540 switch (type_tag) |
|
1541 { |
|
1542 case scalar_constant: |
|
1543 retval = ComplexMatrix (1, 1, Complex (scalar)); |
|
1544 break; |
|
1545 |
|
1546 case complex_scalar_constant: |
|
1547 retval = ComplexMatrix (1, 1, *complex_scalar); |
|
1548 break; |
|
1549 |
|
1550 case matrix_constant: |
|
1551 retval = ComplexMatrix (*matrix); |
|
1552 break; |
|
1553 |
|
1554 case complex_matrix_constant: |
|
1555 retval = *complex_matrix; |
|
1556 break; |
|
1557 |
1572
|
1558 case char_matrix_constant: |
|
1559 retval = ComplexMatrix (*char_matrix); |
|
1560 break; |
|
1561 |
|
1562 case char_matrix_constant_str: |
1299
|
1563 { |
1827
|
1564 int flag = force_string_conv; |
1299
|
1565 if (! flag) |
2181
|
1566 flag = Vimplicit_str_to_num_ok; |
1299
|
1567 |
|
1568 if (flag < 0) |
|
1569 warn_implicit_conversion ("string", "complex matrix"); |
|
1570 |
|
1571 if (flag) |
1572
|
1572 retval = ComplexMatrix (*char_matrix); |
1299
|
1573 else |
1572
|
1574 gripe_invalid_conversion ("complex", "real matrix"); |
1299
|
1575 } |
|
1576 break; |
|
1577 |
|
1578 case range_constant: |
|
1579 retval = range->matrix_value (); |
|
1580 break; |
|
1581 |
|
1582 default: |
|
1583 gripe_invalid_conversion (type_as_string (), "complex matrix"); |
|
1584 break; |
|
1585 } |
|
1586 |
|
1587 return retval; |
|
1588 } |
|
1589 |
1572
|
1590 // XXX FIXME XXX -- this needs to try to do some conversions... |
|
1591 |
|
1592 charMatrix |
2086
|
1593 OCT_VAL_REP::char_matrix_value (bool force_string_conv) const |
1572
|
1594 { |
|
1595 charMatrix retval; |
|
1596 |
1827
|
1597 int flag = force_string_conv; |
1572
|
1598 if (! flag) |
2181
|
1599 flag = Vimplicit_str_to_num_ok; |
1572
|
1600 |
|
1601 switch (type_tag) |
|
1602 { |
|
1603 case char_matrix_constant: |
|
1604 case char_matrix_constant_str: |
|
1605 retval = *char_matrix; |
|
1606 break; |
|
1607 |
|
1608 default: |
1652
|
1609 if (! (rows () == 0 && columns () == 0)) |
|
1610 gripe_invalid_conversion (type_as_string (), "string"); |
1572
|
1611 break; |
|
1612 } |
|
1613 |
|
1614 return retval; |
|
1615 } |
|
1616 |
|
1617 charMatrix |
2086
|
1618 OCT_VAL_REP::all_strings (void) const |
1355
|
1619 { |
1572
|
1620 if (type_tag == char_matrix_constant_str) |
|
1621 return *char_matrix; |
1355
|
1622 else |
|
1623 { |
|
1624 gripe_invalid_conversion (type_as_string (), "string"); |
|
1625 return 0; |
|
1626 } |
|
1627 } |
|
1628 |
1728
|
1629 string |
2086
|
1630 OCT_VAL_REP::string_value (void) const |
1299
|
1631 { |
1728
|
1632 string retval; |
|
1633 |
1572
|
1634 if (type_tag == char_matrix_constant_str) |
1728
|
1635 retval = char_matrix->row_as_string (0); // XXX FIXME??? XXX |
1299
|
1636 else |
1728
|
1637 gripe_invalid_conversion (type_as_string (), "string"); |
|
1638 |
|
1639 return retval; |
1299
|
1640 } |
|
1641 |
|
1642 Range |
2086
|
1643 OCT_VAL_REP::range_value (void) const |
1299
|
1644 { |
|
1645 assert (type_tag == range_constant); |
|
1646 return *range; |
|
1647 } |
|
1648 |
|
1649 Octave_map |
2086
|
1650 OCT_VAL_REP::map_value (void) const |
1299
|
1651 { |
|
1652 assert (type_tag == map_constant); |
|
1653 return *a_map; |
|
1654 } |
|
1655 |
2086
|
1656 octave_value& |
|
1657 OCT_VAL_REP::lookup_map_element (const string& name, bool insert, bool silent) |
1299
|
1658 { |
2086
|
1659 static octave_value retval; |
1299
|
1660 |
|
1661 if (type_tag == map_constant) |
|
1662 { |
|
1663 Pix idx = a_map->seek (name); |
|
1664 |
|
1665 if (idx) |
|
1666 return a_map->contents (idx); |
|
1667 else if (insert) |
|
1668 return (*a_map) [name]; |
|
1669 else if (! silent) |
1755
|
1670 error ("structure has no member `%s'", name.c_str ()); |
1299
|
1671 } |
|
1672 else if (! silent) |
|
1673 error ("invalid structure access attempted"); |
|
1674 |
|
1675 return retval; |
|
1676 } |
|
1677 |
|
1678 // This could be made more efficient by doing all the work here rather |
|
1679 // than relying on matrix_value() to do any possible type conversions. |
|
1680 |
|
1681 ColumnVector |
2086
|
1682 OCT_VAL_REP::vector_value (bool force_string_conv, |
1827
|
1683 bool force_vector_conversion) const |
1299
|
1684 { |
|
1685 ColumnVector retval; |
|
1686 |
1827
|
1687 Matrix m = matrix_value (force_string_conv); |
1299
|
1688 |
|
1689 if (error_state) |
|
1690 return retval; |
|
1691 |
|
1692 int nr = m.rows (); |
|
1693 int nc = m.columns (); |
|
1694 if (nr == 1) |
|
1695 { |
|
1696 retval.resize (nc); |
|
1697 for (int i = 0; i < nc; i++) |
2305
|
1698 retval (i) = m (0, i); |
1299
|
1699 } |
|
1700 else if (nc == 1) |
|
1701 { |
|
1702 retval.resize (nr); |
|
1703 for (int i = 0; i < nr; i++) |
2305
|
1704 retval (i) = m (i, 0); |
1299
|
1705 } |
|
1706 else if (nr > 0 && nc > 0 |
2181
|
1707 && (Vdo_fortran_indexing || force_vector_conversion)) |
1299
|
1708 { |
|
1709 retval.resize (nr * nc); |
|
1710 int k = 0; |
|
1711 for (int j = 0; j < nc; j++) |
|
1712 for (int i = 0; i < nr; i++) |
2305
|
1713 retval (k++) = m (i, j); |
1299
|
1714 } |
|
1715 else |
|
1716 gripe_invalid_conversion ("real matrix", "real vector"); |
|
1717 |
|
1718 return retval; |
|
1719 } |
|
1720 |
|
1721 // This could be made more efficient by doing all the work here rather |
|
1722 // than relying on complex_matrix_value() to do any possible type |
|
1723 // conversions. |
|
1724 |
|
1725 ComplexColumnVector |
2086
|
1726 OCT_VAL_REP::complex_vector_value (bool force_string_conv, |
1827
|
1727 bool force_vector_conversion) const |
1299
|
1728 { |
|
1729 ComplexColumnVector retval; |
|
1730 |
1827
|
1731 ComplexMatrix m = complex_matrix_value (force_string_conv); |
1299
|
1732 |
|
1733 if (error_state) |
|
1734 return retval; |
|
1735 |
|
1736 int nr = m.rows (); |
|
1737 int nc = m.columns (); |
|
1738 if (nr == 1) |
|
1739 { |
|
1740 retval.resize (nc); |
|
1741 for (int i = 0; i < nc; i++) |
2305
|
1742 retval (i) = m (0, i); |
1299
|
1743 } |
|
1744 else if (nc == 1) |
|
1745 { |
|
1746 retval.resize (nr); |
|
1747 for (int i = 0; i < nr; i++) |
2305
|
1748 retval (i) = m (i, 0); |
1299
|
1749 } |
|
1750 else if (nr > 0 && nc > 0 |
2181
|
1751 && (Vdo_fortran_indexing || force_vector_conversion)) |
1299
|
1752 { |
|
1753 retval.resize (nr * nc); |
|
1754 int k = 0; |
|
1755 for (int j = 0; j < nc; j++) |
|
1756 for (int i = 0; i < nr; i++) |
2305
|
1757 retval (k++) = m (i, j); |
1299
|
1758 } |
|
1759 else |
|
1760 gripe_invalid_conversion ("complex matrix", "complex vector"); |
|
1761 |
|
1762 return retval; |
|
1763 } |
|
1764 |
2086
|
1765 octave_value |
|
1766 OCT_VAL_REP::convert_to_str (void) const |
1299
|
1767 { |
2086
|
1768 octave_value retval; |
1299
|
1769 |
|
1770 switch (type_tag) |
|
1771 { |
|
1772 case complex_scalar_constant: |
|
1773 case scalar_constant: |
|
1774 { |
|
1775 double d = double_value (); |
|
1776 |
|
1777 if (xisnan (d)) |
|
1778 { |
|
1779 ::error ("invalid conversion from NaN to character"); |
|
1780 return retval; |
|
1781 } |
|
1782 else |
|
1783 { |
1358
|
1784 // XXX FIXME XXX -- warn about out of range conversions? |
|
1785 |
1299
|
1786 int i = NINT (d); |
|
1787 char s[2]; |
|
1788 s[0] = (char) i; |
|
1789 s[1] = '\0'; |
2086
|
1790 retval = octave_value (s, 1); |
1299
|
1791 } |
|
1792 } |
|
1793 break; |
|
1794 |
|
1795 case complex_matrix_constant: |
|
1796 case matrix_constant: |
|
1797 { |
|
1798 if (rows () == 0 && columns () == 0) |
|
1799 { |
|
1800 char s = '\0'; |
2086
|
1801 retval = octave_value (&s, 1); |
1299
|
1802 } |
|
1803 else |
|
1804 { |
1355
|
1805 Matrix m = matrix_value (); |
|
1806 |
|
1807 int nr = m.rows (); |
|
1808 int nc = m.columns (); |
|
1809 |
|
1810 if (nr == 0 || nc == 0) |
1299
|
1811 { |
|
1812 char s = '\0'; |
2086
|
1813 retval = octave_value (&s, 1); |
1299
|
1814 } |
|
1815 else |
|
1816 { |
1572
|
1817 charMatrix chm (nr, nc); |
|
1818 |
|
1819 for (int j = 0; j < nc; j++) |
1299
|
1820 { |
1572
|
1821 for (int i = 0; i < nr; i++) |
1299
|
1822 { |
2305
|
1823 double d = m (i, j); |
1355
|
1824 |
|
1825 if (xisnan (d)) |
|
1826 { |
|
1827 ::error ("invalid conversion from NaN to character"); |
|
1828 return retval; |
|
1829 } |
|
1830 else |
|
1831 { |
1358
|
1832 // XXX FIXME XXX -- warn about out of |
|
1833 // range conversions? |
1355
|
1834 |
|
1835 int ival = NINT (d); |
2305
|
1836 chm (i, j) = (char) ival; |
1355
|
1837 } |
1299
|
1838 } |
|
1839 } |
1355
|
1840 |
2086
|
1841 retval = octave_value (chm, 1); |
1299
|
1842 } |
|
1843 } |
|
1844 } |
|
1845 break; |
|
1846 |
|
1847 case range_constant: |
|
1848 { |
|
1849 Range r = range_value (); |
|
1850 double b = r.base (); |
|
1851 double incr = r.inc (); |
|
1852 int nel = r.nelem (); |
|
1853 char *s = new char [nel+1]; |
|
1854 s[nel] = '\0'; |
|
1855 for (int i = 0; i < nel; i++) |
|
1856 { |
|
1857 double d = b + i * incr; |
|
1858 |
|
1859 if (xisnan (d)) |
|
1860 { |
|
1861 ::error ("invalid conversion from NaN to character"); |
|
1862 delete [] s; |
|
1863 return retval; |
|
1864 } |
|
1865 else |
|
1866 { |
1358
|
1867 // XXX FIXME XXX -- warn about out of range |
|
1868 // conversions? |
|
1869 |
1299
|
1870 int ival = NINT (d); |
|
1871 s[i] = (char) ival; |
|
1872 } |
|
1873 } |
2086
|
1874 retval = octave_value (s, 1); |
1299
|
1875 delete [] s; |
|
1876 } |
|
1877 break; |
|
1878 |
1572
|
1879 case char_matrix_constant: |
2086
|
1880 retval = octave_value (*char_matrix, 1); |
1572
|
1881 break; |
|
1882 |
|
1883 case char_matrix_constant_str: |
2086
|
1884 retval = octave_value (*char_matrix, 1); |
1299
|
1885 break; |
|
1886 |
|
1887 default: |
|
1888 gripe_invalid_conversion (type_as_string (), "string"); |
|
1889 break; |
|
1890 } |
|
1891 |
|
1892 return retval; |
|
1893 } |
|
1894 |
|
1895 void |
2086
|
1896 OCT_VAL_REP::convert_to_row_or_column_vector (void) |
1299
|
1897 { |
|
1898 assert (type_tag == matrix_constant || type_tag == complex_matrix_constant); |
|
1899 |
|
1900 int nr = rows (); |
|
1901 int nc = columns (); |
|
1902 |
|
1903 if (nr == 1 || nc == 1) |
|
1904 return; |
|
1905 |
|
1906 int len = nr * nc; |
|
1907 |
|
1908 assert (len > 0); |
|
1909 |
|
1910 int new_nr = 1; |
|
1911 int new_nc = 1; |
|
1912 |
2181
|
1913 if (Vprefer_column_vectors) |
1299
|
1914 new_nr = len; |
|
1915 else |
|
1916 new_nc = len; |
|
1917 |
|
1918 if (type_tag == matrix_constant) |
|
1919 { |
|
1920 Matrix *m = new Matrix (new_nr, new_nc); |
|
1921 |
|
1922 double *cop_out = matrix->fortran_vec (); |
|
1923 |
|
1924 for (int i = 0; i < len; i++) |
|
1925 { |
|
1926 if (new_nr == 1) |
2305
|
1927 (*m) (0, i) = *cop_out++; |
1299
|
1928 else |
2305
|
1929 (*m) (i, 0) = *cop_out++; |
1299
|
1930 } |
|
1931 |
|
1932 delete matrix; |
|
1933 matrix = m; |
|
1934 } |
|
1935 else |
|
1936 { |
|
1937 ComplexMatrix *cm = new ComplexMatrix (new_nr, new_nc); |
|
1938 |
|
1939 Complex *cop_out = complex_matrix->fortran_vec (); |
|
1940 |
|
1941 for (int i = 0; i < len; i++) |
|
1942 { |
|
1943 if (new_nr == 1) |
2305
|
1944 (*cm) (0, i) = *cop_out++; |
1299
|
1945 else |
2305
|
1946 (*cm) (i, 0) = *cop_out++; |
1299
|
1947 } |
|
1948 |
|
1949 delete complex_matrix; |
|
1950 complex_matrix = cm; |
|
1951 } |
|
1952 } |
|
1953 |
|
1954 void |
2086
|
1955 OCT_VAL_REP::convert_to_matrix_type (bool make_complex) |
1558
|
1956 { |
|
1957 switch (type_tag) |
|
1958 { |
|
1959 case complex_scalar_constant: |
|
1960 { |
|
1961 Complex *old_complex = complex_scalar; |
|
1962 complex_matrix = new ComplexMatrix (1, 1, *complex_scalar); |
|
1963 type_tag = complex_matrix_constant; |
|
1964 delete old_complex; |
|
1965 } |
|
1966 break; |
|
1967 |
|
1968 case scalar_constant: |
|
1969 { |
1581
|
1970 if (make_complex) |
|
1971 { |
|
1972 complex_matrix = new ComplexMatrix (1, 1, scalar); |
|
1973 type_tag = complex_matrix_constant; |
|
1974 } |
|
1975 else |
|
1976 { |
|
1977 matrix = new Matrix (1, 1, scalar); |
|
1978 type_tag = matrix_constant; |
|
1979 } |
1558
|
1980 } |
|
1981 break; |
|
1982 |
|
1983 case unknown_constant: |
|
1984 { |
1581
|
1985 if (make_complex) |
|
1986 { |
|
1987 complex_matrix = new ComplexMatrix (); |
|
1988 type_tag = complex_matrix_constant; |
|
1989 } |
|
1990 else |
|
1991 { |
|
1992 matrix = new Matrix (); |
|
1993 type_tag = matrix_constant; |
|
1994 } |
1558
|
1995 } |
|
1996 break; |
|
1997 |
1572
|
1998 case range_constant: |
|
1999 { |
1581
|
2000 if (make_complex) |
|
2001 { |
|
2002 ComplexMatrix *tmp = new ComplexMatrix (range->matrix_value ()); |
|
2003 delete range; |
|
2004 complex_matrix = tmp; |
|
2005 type_tag = complex_matrix_constant; |
|
2006 } |
|
2007 else |
|
2008 { |
|
2009 Matrix *tmp = new Matrix (range->matrix_value ()); |
|
2010 delete range; |
|
2011 matrix = tmp; |
|
2012 type_tag = matrix_constant; |
|
2013 } |
1572
|
2014 } |
|
2015 break; |
|
2016 |
1558
|
2017 default: |
|
2018 panic_impossible (); |
|
2019 break; |
|
2020 } |
|
2021 } |
|
2022 |
|
2023 void |
2086
|
2024 OCT_VAL_REP::force_numeric (bool force_string_conv) |
1299
|
2025 { |
|
2026 switch (type_tag) |
|
2027 { |
|
2028 case scalar_constant: |
|
2029 case matrix_constant: |
|
2030 case complex_scalar_constant: |
|
2031 case complex_matrix_constant: |
1572
|
2032 case char_matrix_constant: |
1299
|
2033 break; |
|
2034 |
1572
|
2035 case char_matrix_constant_str: |
1299
|
2036 { |
2181
|
2037 if (! force_string_conv && ! Vimplicit_str_to_num_ok) |
1299
|
2038 { |
1490
|
2039 ::error ("string to numeric conversion failed --\ |
|
2040 default conversion turned off"); |
1299
|
2041 return; |
|
2042 } |
|
2043 |
1572
|
2044 int nr = char_matrix->rows (); |
|
2045 int nc = char_matrix->cols (); |
1355
|
2046 |
|
2047 if (nr == 1 && nc == 1) |
1299
|
2048 { |
|
2049 type_tag = scalar_constant; |
2305
|
2050 double tmp = toascii ((int) (*char_matrix) (0, 0)); |
1572
|
2051 delete char_matrix; |
|
2052 scalar = tmp; |
1299
|
2053 } |
1355
|
2054 else if (nr == 0 || nc == 0) |
1299
|
2055 { |
1572
|
2056 delete char_matrix; |
1299
|
2057 type_tag = matrix_constant; |
|
2058 matrix = new Matrix (0, 0); |
|
2059 } |
1355
|
2060 else if (nr > 0 && nc > 0) |
|
2061 { |
|
2062 type_tag = matrix_constant; |
|
2063 |
|
2064 Matrix *tm = new Matrix (nr, nc); |
|
2065 |
|
2066 for (int i = 0; i < nr; i++) |
|
2067 { |
|
2068 for (int j = 0; j < nc; j++) |
|
2069 { |
2305
|
2070 int c = (int) (*char_matrix) (i, j); |
|
2071 (*tm) (i, j) = toascii (c); |
1355
|
2072 } |
|
2073 } |
1572
|
2074 delete char_matrix; |
1355
|
2075 matrix = tm; |
|
2076 } |
1299
|
2077 else |
|
2078 panic_impossible (); |
|
2079 } |
|
2080 break; |
|
2081 |
|
2082 case range_constant: |
|
2083 { |
|
2084 int len = range->nelem (); |
|
2085 if (len > 1) |
|
2086 { |
|
2087 type_tag = matrix_constant; |
|
2088 Matrix *tm = new Matrix (1, len); |
|
2089 double b = range->base (); |
|
2090 double increment = range->inc (); |
|
2091 for (int i = 0; i < len; i++) |
2305
|
2092 (*tm) (0, i) = b + i * increment; |
1572
|
2093 delete range; |
1299
|
2094 matrix = tm; |
|
2095 } |
|
2096 else if (len == 1) |
|
2097 { |
|
2098 type_tag = scalar_constant; |
|
2099 scalar = range->base (); |
|
2100 } |
|
2101 } |
|
2102 break; |
|
2103 |
|
2104 default: |
|
2105 gripe_invalid_conversion (type_as_string (), "numeric type"); |
|
2106 break; |
|
2107 } |
|
2108 } |
|
2109 |
2086
|
2110 octave_value |
|
2111 OCT_VAL_REP::make_numeric (bool force_string_conv) const |
1299
|
2112 { |
2086
|
2113 octave_value retval; |
1299
|
2114 |
|
2115 switch (type_tag) |
|
2116 { |
|
2117 case scalar_constant: |
1355
|
2118 retval = scalar; |
1299
|
2119 break; |
|
2120 |
|
2121 case matrix_constant: |
1355
|
2122 retval = *matrix; |
1299
|
2123 break; |
|
2124 |
|
2125 case complex_scalar_constant: |
1355
|
2126 retval = *complex_scalar; |
1299
|
2127 break; |
|
2128 |
|
2129 case complex_matrix_constant: |
1355
|
2130 retval = *complex_matrix; |
1299
|
2131 break; |
|
2132 |
1572
|
2133 case char_matrix_constant: |
|
2134 retval = *char_matrix; |
|
2135 break; |
|
2136 |
|
2137 case char_matrix_constant_str: |
|
2138 { |
1827
|
2139 int flag = force_string_conv; |
1572
|
2140 if (! flag) |
2181
|
2141 flag = Vimplicit_str_to_num_ok; |
1572
|
2142 |
|
2143 if (flag < 0) |
|
2144 warn_implicit_conversion ("string", "char matrix"); |
|
2145 |
|
2146 if (flag) |
|
2147 { |
2260
|
2148 retval = octave_value (*char_matrix, true); |
1827
|
2149 retval.force_numeric (force_string_conv); |
1572
|
2150 } |
|
2151 else |
|
2152 gripe_invalid_conversion ("string", "char matrix"); |
|
2153 } |
1299
|
2154 break; |
|
2155 |
|
2156 case range_constant: |
1355
|
2157 retval = *range; |
1827
|
2158 retval.force_numeric (force_string_conv); |
1299
|
2159 break; |
|
2160 |
|
2161 default: |
|
2162 gripe_invalid_conversion (type_as_string (), "numeric value"); |
|
2163 break; |
|
2164 } |
|
2165 |
|
2166 return retval; |
|
2167 } |
|
2168 |
|
2169 void |
2086
|
2170 OCT_VAL_REP::bump_value (tree_expression::type etype) |
1299
|
2171 { |
|
2172 switch (etype) |
|
2173 { |
|
2174 case tree_expression::increment: |
|
2175 switch (type_tag) |
|
2176 { |
|
2177 case scalar_constant: |
|
2178 scalar++; |
|
2179 break; |
|
2180 |
|
2181 case matrix_constant: |
|
2182 *matrix = *matrix + 1.0; |
|
2183 break; |
|
2184 |
|
2185 case complex_scalar_constant: |
|
2186 *complex_scalar = *complex_scalar + 1.0; |
|
2187 break; |
|
2188 |
|
2189 case complex_matrix_constant: |
|
2190 *complex_matrix = *complex_matrix + 1.0; |
|
2191 break; |
|
2192 |
|
2193 case range_constant: |
|
2194 range->set_base (range->base () + 1.0); |
|
2195 range->set_limit (range->limit () + 1.0); |
|
2196 break; |
|
2197 |
|
2198 default: |
|
2199 gripe_wrong_type_arg ("operator ++", type_as_string ()); |
|
2200 break; |
|
2201 } |
|
2202 break; |
|
2203 |
|
2204 case tree_expression::decrement: |
|
2205 switch (type_tag) |
|
2206 { |
|
2207 case scalar_constant: |
|
2208 scalar--; |
|
2209 break; |
|
2210 |
|
2211 case matrix_constant: |
|
2212 *matrix = *matrix - 1.0; |
|
2213 break; |
|
2214 |
|
2215 case range_constant: |
|
2216 range->set_base (range->base () - 1.0); |
|
2217 range->set_limit (range->limit () - 1.0); |
|
2218 break; |
|
2219 |
|
2220 default: |
|
2221 gripe_wrong_type_arg ("operator --", type_as_string ()); |
|
2222 break; |
|
2223 } |
|
2224 break; |
|
2225 |
|
2226 default: |
|
2227 panic_impossible (); |
|
2228 break; |
|
2229 } |
|
2230 } |
|
2231 |
|
2232 void |
2086
|
2233 OCT_VAL_REP::resize (int i, int j) |
1299
|
2234 { |
|
2235 switch (type_tag) |
|
2236 { |
|
2237 case matrix_constant: |
|
2238 matrix->resize (i, j); |
|
2239 break; |
|
2240 |
|
2241 case complex_matrix_constant: |
|
2242 complex_matrix->resize (i, j); |
|
2243 break; |
|
2244 |
|
2245 default: |
|
2246 gripe_wrong_type_arg ("resize", type_as_string ()); |
|
2247 break; |
|
2248 } |
|
2249 } |
|
2250 |
|
2251 void |
2086
|
2252 OCT_VAL_REP::resize (int i, int j, double val) |
1299
|
2253 { |
|
2254 switch (type_tag) |
|
2255 { |
|
2256 case matrix_constant: |
|
2257 matrix->resize (i, j, val); |
|
2258 break; |
|
2259 |
|
2260 case complex_matrix_constant: |
|
2261 complex_matrix->resize (i, j, val); |
|
2262 break; |
|
2263 |
|
2264 default: |
|
2265 gripe_wrong_type_arg ("resize", type_as_string ()); |
|
2266 break; |
|
2267 } |
|
2268 } |
|
2269 |
|
2270 void |
2086
|
2271 OCT_VAL_REP::stash_original_text (const string &s) |
1299
|
2272 { |
1755
|
2273 orig_text = s; |
1299
|
2274 } |
|
2275 |
|
2276 void |
2086
|
2277 OCT_VAL_REP::maybe_mutate (void) |
1299
|
2278 { |
|
2279 switch (type_tag) |
|
2280 { |
|
2281 case complex_scalar_constant: |
|
2282 if (::imag (*complex_scalar) == 0.0) |
|
2283 { |
|
2284 double d = ::real (*complex_scalar); |
|
2285 delete complex_scalar; |
|
2286 scalar = d; |
|
2287 type_tag = scalar_constant; |
|
2288 } |
|
2289 break; |
|
2290 |
|
2291 case complex_matrix_constant: |
|
2292 if (! any_element_is_complex (*complex_matrix)) |
|
2293 { |
|
2294 Matrix *m = new Matrix (::real (*complex_matrix)); |
|
2295 delete complex_matrix; |
|
2296 matrix = m; |
|
2297 type_tag = matrix_constant; |
|
2298 } |
|
2299 break; |
|
2300 |
|
2301 default: |
|
2302 break; |
|
2303 } |
|
2304 |
1358
|
2305 // Avoid calling rows() and columns() for things like magic_colon. |
1299
|
2306 |
|
2307 int nr = 1; |
|
2308 int nc = 1; |
|
2309 if (type_tag == matrix_constant |
|
2310 || type_tag == complex_matrix_constant |
|
2311 || type_tag == range_constant) |
|
2312 { |
|
2313 nr = rows (); |
|
2314 nc = columns (); |
|
2315 } |
|
2316 |
|
2317 switch (type_tag) |
|
2318 { |
|
2319 case matrix_constant: |
|
2320 if (nr == 1 && nc == 1) |
|
2321 { |
2305
|
2322 double d = (*matrix) (0, 0); |
1299
|
2323 delete matrix; |
|
2324 scalar = d; |
|
2325 type_tag = scalar_constant; |
|
2326 } |
|
2327 break; |
|
2328 |
|
2329 case complex_matrix_constant: |
|
2330 if (nr == 1 && nc == 1) |
|
2331 { |
2305
|
2332 Complex c = (*complex_matrix) (0, 0); |
1299
|
2333 delete complex_matrix; |
|
2334 complex_scalar = new Complex (c); |
|
2335 type_tag = complex_scalar_constant; |
|
2336 } |
|
2337 break; |
|
2338 |
|
2339 case range_constant: |
|
2340 if (nr == 1 && nc == 1) |
|
2341 { |
|
2342 double d = range->base (); |
|
2343 delete range; |
|
2344 scalar = d; |
|
2345 type_tag = scalar_constant; |
|
2346 } |
|
2347 break; |
|
2348 |
|
2349 default: |
|
2350 break; |
|
2351 } |
|
2352 } |
|
2353 |
|
2354 void |
2086
|
2355 OCT_VAL_REP::print (ostream& output_buf) |
1299
|
2356 { |
|
2357 if (error_state) |
|
2358 return; |
|
2359 |
|
2360 switch (type_tag) |
|
2361 { |
|
2362 case scalar_constant: |
1973
|
2363 octave_print_internal (output_buf, scalar, false); |
1299
|
2364 break; |
|
2365 |
|
2366 case matrix_constant: |
1973
|
2367 octave_print_internal (output_buf, *matrix, false, |
|
2368 struct_indent); |
1299
|
2369 break; |
|
2370 |
|
2371 case complex_scalar_constant: |
1973
|
2372 octave_print_internal (output_buf, *complex_scalar, false); |
1299
|
2373 break; |
|
2374 |
|
2375 case complex_matrix_constant: |
1973
|
2376 octave_print_internal (output_buf, *complex_matrix, false, |
|
2377 struct_indent); |
1299
|
2378 break; |
|
2379 |
1572
|
2380 case char_matrix_constant: |
1973
|
2381 octave_print_internal (output_buf, *char_matrix, false, |
|
2382 struct_indent); |
1572
|
2383 break; |
|
2384 |
|
2385 case char_matrix_constant_str: |
1973
|
2386 octave_print_internal (output_buf, *char_matrix, false, true, |
|
2387 struct_indent); |
1299
|
2388 break; |
|
2389 |
|
2390 case range_constant: |
1973
|
2391 octave_print_internal (output_buf, *range, false, struct_indent); |
1299
|
2392 break; |
|
2393 |
|
2394 case map_constant: |
|
2395 { |
1358
|
2396 // XXX FIXME XXX -- would be nice to print the output in some |
|
2397 // standard order. Maybe all substructures first, maybe |
|
2398 // alphabetize entries, etc. |
|
2399 |
2086
|
2400 begin_unwind_frame ("OCT_VAL_REP_print"); |
1299
|
2401 |
1973
|
2402 unwind_protect_int (struct_indent); |
2181
|
2403 unwind_protect_int (Vstruct_levels_to_print); |
|
2404 |
|
2405 if (Vstruct_levels_to_print-- > 0) |
1299
|
2406 { |
1973
|
2407 output_buf.form ("\n%*s{\n", struct_indent, ""); |
|
2408 |
|
2409 increment_struct_indent (); |
|
2410 |
|
2411 Pix p = a_map->first (); |
|
2412 |
|
2413 while (p) |
1299
|
2414 { |
1973
|
2415 bool pad_after = false; |
|
2416 |
1755
|
2417 string key = a_map->key (p); |
2086
|
2418 octave_value val = a_map->contents (p); |
1299
|
2419 |
1973
|
2420 a_map->next (p); |
|
2421 |
|
2422 output_buf.form ("%*s%s =", struct_indent, |
1755
|
2423 "", key.c_str ()); |
|
2424 |
1973
|
2425 if (val.print_as_scalar ()) |
|
2426 output_buf << " "; |
|
2427 else if (val.print_as_structure ()) |
|
2428 { |
|
2429 if (p) |
|
2430 pad_after = true; |
|
2431 } |
|
2432 else |
|
2433 { |
|
2434 if (p) |
|
2435 pad_after = true; |
|
2436 |
|
2437 output_buf << "\n\n"; |
|
2438 } |
1299
|
2439 |
|
2440 val.print (output_buf); |
1973
|
2441 |
|
2442 if (pad_after) |
|
2443 output_buf << "\n"; |
1299
|
2444 } |
|
2445 |
1973
|
2446 decrement_struct_indent (); |
|
2447 |
|
2448 output_buf.form ("%*s%s", struct_indent, "", "}\n"); |
1299
|
2449 } |
|
2450 else |
1973
|
2451 output_buf << " <structure>\n"; |
1299
|
2452 |
2086
|
2453 run_unwind_frame ("OCT_VAL_REP_print"); |
1299
|
2454 } |
|
2455 break; |
|
2456 |
|
2457 case unknown_constant: |
|
2458 case magic_colon: |
|
2459 case all_va_args: |
|
2460 panic_impossible (); |
|
2461 break; |
|
2462 } |
|
2463 } |
|
2464 |
|
2465 void |
2086
|
2466 OCT_VAL_REP::gripe_wrong_type_arg (const char *name, |
|
2467 const octave_value_rep& tcr) const |
1299
|
2468 { |
|
2469 if (name) |
|
2470 ::error ("%s: wrong type argument `%s'", name, tcr.type_as_string ()); |
|
2471 else |
|
2472 ::error ("wrong type argument `%s'", name, tcr.type_as_string ()); |
|
2473 } |
|
2474 |
|
2475 char * |
2086
|
2476 OCT_VAL_REP::type_as_string (void) const |
1299
|
2477 { |
|
2478 switch (type_tag) |
|
2479 { |
|
2480 case scalar_constant: |
|
2481 return "real scalar"; |
|
2482 |
|
2483 case matrix_constant: |
|
2484 return "real matrix"; |
|
2485 |
|
2486 case complex_scalar_constant: |
|
2487 return "complex scalar"; |
|
2488 |
|
2489 case complex_matrix_constant: |
|
2490 return "complex matrix"; |
|
2491 |
1572
|
2492 case char_matrix_constant: |
|
2493 return "char matrix"; |
|
2494 |
|
2495 case char_matrix_constant_str: |
1299
|
2496 return "string"; |
|
2497 |
|
2498 case range_constant: |
|
2499 return "range"; |
|
2500 |
|
2501 case map_constant: |
|
2502 return "structure"; |
|
2503 |
|
2504 default: |
|
2505 return "<unknown type>"; |
|
2506 } |
|
2507 } |
|
2508 |
2086
|
2509 octave_value |
|
2510 do_binary_op (octave_value& a, octave_value& b, tree_expression::type t) |
1299
|
2511 { |
2086
|
2512 octave_value retval; |
1299
|
2513 |
1827
|
2514 bool first_empty = (a.rows () == 0 || a.columns () == 0); |
|
2515 bool second_empty = (b.rows () == 0 || b.columns () == 0); |
1299
|
2516 |
|
2517 if (first_empty || second_empty) |
|
2518 { |
2181
|
2519 int flag = Vpropagate_empty_matrices; |
1299
|
2520 if (flag < 0) |
|
2521 warning ("binary operation on empty matrix"); |
|
2522 else if (flag == 0) |
|
2523 { |
|
2524 ::error ("invalid binary operation on empty matrix"); |
|
2525 return retval; |
|
2526 } |
|
2527 } |
|
2528 |
2260
|
2529 int force = (a.is_string () && b.is_string () |
|
2530 && (t == tree_expression::cmp_lt |
|
2531 || t == tree_expression::cmp_le |
|
2532 || t == tree_expression::cmp_eq |
|
2533 || t == tree_expression::cmp_ge |
|
2534 || t == tree_expression::cmp_gt |
|
2535 || t == tree_expression::cmp_ne)); |
|
2536 |
|
2537 octave_value tmp_a = a.make_numeric (force); |
1299
|
2538 |
|
2539 if (error_state) |
|
2540 return retval; |
|
2541 |
2260
|
2542 octave_value tmp_b = b.make_numeric (force); |
1299
|
2543 |
|
2544 if (error_state) |
|
2545 return retval; |
|
2546 |
2086
|
2547 OCT_VAL_REP::constant_type a_type = tmp_a.const_type (); |
|
2548 OCT_VAL_REP::constant_type b_type = tmp_b.const_type (); |
1299
|
2549 |
|
2550 double d1, d2; |
|
2551 Matrix m1, m2; |
|
2552 Complex c1, c2; |
|
2553 ComplexMatrix cm1, cm2; |
|
2554 |
|
2555 switch (a_type) |
|
2556 { |
2086
|
2557 case OCT_VAL_REP::scalar_constant: |
1299
|
2558 |
|
2559 d1 = tmp_a.double_value (); |
|
2560 |
|
2561 switch (b_type) |
|
2562 { |
2086
|
2563 case OCT_VAL_REP::scalar_constant: |
1299
|
2564 d2 = tmp_b.double_value (); |
|
2565 retval = do_binary_op (d1, d2, t); |
|
2566 break; |
|
2567 |
2086
|
2568 case OCT_VAL_REP::matrix_constant: |
|
2569 case OCT_VAL_REP::char_matrix_constant: |
1299
|
2570 m2 = tmp_b.matrix_value (); |
|
2571 retval = do_binary_op (d1, m2, t); |
|
2572 break; |
|
2573 |
2086
|
2574 case OCT_VAL_REP::complex_scalar_constant: |
1299
|
2575 c2 = tmp_b.complex_value (); |
|
2576 retval = do_binary_op (d1, c2, t); |
|
2577 break; |
|
2578 |
2086
|
2579 case OCT_VAL_REP::complex_matrix_constant: |
1299
|
2580 cm2 = tmp_b.complex_matrix_value (); |
|
2581 retval = do_binary_op (d1, cm2, t); |
|
2582 break; |
|
2583 |
|
2584 default: |
|
2585 gripe_wrong_type_arg_for_binary_op (tmp_b); |
|
2586 break; |
|
2587 } |
|
2588 break; |
|
2589 |
2086
|
2590 case OCT_VAL_REP::matrix_constant: |
|
2591 case OCT_VAL_REP::char_matrix_constant: |
1299
|
2592 |
|
2593 m1 = tmp_a.matrix_value (); |
|
2594 |
|
2595 switch (b_type) |
|
2596 { |
2086
|
2597 case OCT_VAL_REP::scalar_constant: |
1299
|
2598 d2 = tmp_b.double_value (); |
|
2599 retval = do_binary_op (m1, d2, t); |
|
2600 break; |
|
2601 |
2086
|
2602 case OCT_VAL_REP::matrix_constant: |
|
2603 case OCT_VAL_REP::char_matrix_constant: |
1299
|
2604 m2 = tmp_b.matrix_value (); |
|
2605 retval = do_binary_op (m1, m2, t); |
|
2606 break; |
|
2607 |
2086
|
2608 case OCT_VAL_REP::complex_scalar_constant: |
1299
|
2609 c2 = tmp_b.complex_value (); |
|
2610 retval = do_binary_op (m1, c2, t); |
|
2611 break; |
|
2612 |
2086
|
2613 case OCT_VAL_REP::complex_matrix_constant: |
1299
|
2614 cm2 = tmp_b.complex_matrix_value (); |
|
2615 retval = do_binary_op (m1, cm2, t); |
|
2616 break; |
|
2617 |
|
2618 default: |
|
2619 gripe_wrong_type_arg_for_binary_op (tmp_b); |
|
2620 break; |
|
2621 } |
|
2622 break; |
|
2623 |
2086
|
2624 case OCT_VAL_REP::complex_scalar_constant: |
1299
|
2625 |
|
2626 c1 = tmp_a.complex_value (); |
|
2627 |
|
2628 switch (b_type) |
|
2629 { |
2086
|
2630 case OCT_VAL_REP::scalar_constant: |
1299
|
2631 d2 = tmp_b.double_value (); |
|
2632 retval = do_binary_op (c1, d2, t); |
|
2633 break; |
|
2634 |
2086
|
2635 case OCT_VAL_REP::matrix_constant: |
|
2636 case OCT_VAL_REP::char_matrix_constant: |
1299
|
2637 m2 = tmp_b.matrix_value (); |
|
2638 retval = do_binary_op (c1, m2, t); |
|
2639 break; |
|
2640 |
2086
|
2641 case OCT_VAL_REP::complex_scalar_constant: |
1299
|
2642 c2 = tmp_b.complex_value (); |
|
2643 retval = do_binary_op (c1, c2, t); |
|
2644 break; |
|
2645 |
2086
|
2646 case OCT_VAL_REP::complex_matrix_constant: |
1299
|
2647 cm2 = tmp_b.complex_matrix_value (); |
|
2648 retval = do_binary_op (c1, cm2, t); |
|
2649 break; |
|
2650 |
|
2651 default: |
|
2652 gripe_wrong_type_arg_for_binary_op (tmp_b); |
|
2653 break; |
|
2654 } |
|
2655 break; |
|
2656 |
2086
|
2657 case OCT_VAL_REP::complex_matrix_constant: |
1299
|
2658 |
|
2659 cm1 = tmp_a.complex_matrix_value (); |
|
2660 |
|
2661 switch (b_type) |
|
2662 { |
2086
|
2663 case OCT_VAL_REP::scalar_constant: |
1299
|
2664 d2 = tmp_b.double_value (); |
|
2665 retval = do_binary_op (cm1, d2, t); |
|
2666 break; |
|
2667 |
2086
|
2668 case OCT_VAL_REP::matrix_constant: |
|
2669 case OCT_VAL_REP::char_matrix_constant: |
1299
|
2670 m2 = tmp_b.matrix_value (); |
|
2671 retval = do_binary_op (cm1, m2, t); |
|
2672 break; |
|
2673 |
2086
|
2674 case OCT_VAL_REP::complex_scalar_constant: |
1299
|
2675 c2 = tmp_b.complex_value (); |
|
2676 retval = do_binary_op (cm1, c2, t); |
|
2677 break; |
|
2678 |
2086
|
2679 case OCT_VAL_REP::complex_matrix_constant: |
1299
|
2680 cm2 = tmp_b.complex_matrix_value (); |
|
2681 retval = do_binary_op (cm1, cm2, t); |
|
2682 break; |
|
2683 |
|
2684 default: |
|
2685 gripe_wrong_type_arg_for_binary_op (tmp_b); |
|
2686 break; |
|
2687 } |
|
2688 break; |
|
2689 |
|
2690 default: |
|
2691 gripe_wrong_type_arg_for_binary_op (tmp_a); |
|
2692 break; |
|
2693 } |
|
2694 |
|
2695 return retval; |
|
2696 } |
|
2697 |
2086
|
2698 octave_value |
|
2699 do_unary_op (octave_value& a, tree_expression::type t) |
1299
|
2700 { |
2086
|
2701 octave_value retval; |
1299
|
2702 |
|
2703 if (a.rows () == 0 || a.columns () == 0) |
|
2704 { |
2181
|
2705 int flag = Vpropagate_empty_matrices; |
1299
|
2706 if (flag < 0) |
|
2707 warning ("unary operation on empty matrix"); |
|
2708 else if (flag == 0) |
|
2709 { |
|
2710 ::error ("invalid unary operation on empty matrix"); |
|
2711 return retval; |
|
2712 } |
|
2713 } |
|
2714 |
2258
|
2715 // XXX FIXME XXX -- it is very unlikely that this is the correct |
|
2716 // place for this special case... |
|
2717 |
|
2718 if (a.const_type () == OCT_VAL_REP::char_matrix_constant_str |
2266
|
2719 && (t == tree_expression::transpose |
|
2720 || t == tree_expression::hermitian)) |
2258
|
2721 { |
|
2722 charMatrix chm = a.all_strings (); |
|
2723 |
|
2724 if (! error_state) |
|
2725 retval = octave_value (chm.transpose (), true); |
|
2726 } |
|
2727 else |
1299
|
2728 { |
2258
|
2729 octave_value tmp_a = a.make_numeric (); |
|
2730 |
|
2731 if (error_state) |
|
2732 return retval; |
|
2733 |
|
2734 switch (tmp_a.const_type ()) |
|
2735 { |
|
2736 case OCT_VAL_REP::scalar_constant: |
|
2737 retval = do_unary_op (tmp_a.double_value (), t); |
|
2738 break; |
|
2739 |
|
2740 case OCT_VAL_REP::matrix_constant: |
|
2741 { |
|
2742 Matrix m = tmp_a.matrix_value (); |
|
2743 retval = do_unary_op (m, t); |
|
2744 } |
|
2745 break; |
|
2746 |
|
2747 case OCT_VAL_REP::complex_scalar_constant: |
|
2748 retval = do_unary_op (tmp_a.complex_value (), t); |
|
2749 break; |
|
2750 |
|
2751 case OCT_VAL_REP::complex_matrix_constant: |
|
2752 { |
|
2753 ComplexMatrix m = tmp_a.complex_matrix_value (); |
|
2754 retval = do_unary_op (m, t); |
|
2755 } |
|
2756 break; |
|
2757 |
|
2758 default: |
|
2759 gripe_wrong_type_arg_for_unary_op (tmp_a); |
|
2760 break; |
|
2761 } |
1299
|
2762 } |
|
2763 |
|
2764 return retval; |
|
2765 } |
|
2766 |
|
2767 // Indexing operations for the tree-constant representation class. |
1558
|
2768 |
|
2769 void |
2086
|
2770 OCT_VAL_REP::clear_index (void) |
1558
|
2771 { |
|
2772 switch (type_tag) |
|
2773 { |
|
2774 case matrix_constant: |
|
2775 matrix->clear_index (); |
|
2776 break; |
|
2777 |
2086
|
2778 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2779 complex_matrix->clear_index (); |
|
2780 break; |
|
2781 |
1572
|
2782 case char_matrix_constant: |
|
2783 case char_matrix_constant_str: |
|
2784 char_matrix->clear_index (); |
|
2785 break; |
|
2786 |
1558
|
2787 default: |
|
2788 panic_impossible (); |
|
2789 break; |
|
2790 } |
|
2791 } |
|
2792 |
|
2793 #if 0 |
|
2794 void |
2086
|
2795 OCT_VAL_REP::set_index (double d) |
1558
|
2796 { |
|
2797 switch (type_tag) |
|
2798 { |
|
2799 case matrix_constant: |
|
2800 matrix->set_index (d); |
|
2801 break; |
|
2802 |
2086
|
2803 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2804 complex_matrix->set_index (d); |
|
2805 break; |
|
2806 |
2086
|
2807 case OCT_VAL_REP::char_matrix_constant: |
|
2808 case OCT_VAL_REP::char_matrix_constant_str: |
1572
|
2809 char_matrix->set_index (d); |
|
2810 break; |
|
2811 |
1558
|
2812 default: |
|
2813 panic_impossible (); |
|
2814 break; |
|
2815 } |
|
2816 } |
|
2817 #endif |
|
2818 |
|
2819 void |
2086
|
2820 OCT_VAL_REP::set_index (const Range& r) |
1558
|
2821 { |
|
2822 switch (type_tag) |
|
2823 { |
|
2824 case matrix_constant: |
|
2825 matrix->set_index (r); |
|
2826 break; |
|
2827 |
2086
|
2828 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2829 complex_matrix->set_index (r); |
|
2830 break; |
|
2831 |
2086
|
2832 case OCT_VAL_REP::char_matrix_constant: |
|
2833 case OCT_VAL_REP::char_matrix_constant_str: |
1572
|
2834 char_matrix->set_index (r); |
|
2835 break; |
|
2836 |
1558
|
2837 default: |
|
2838 panic_impossible (); |
|
2839 break; |
|
2840 } |
|
2841 } |
|
2842 |
|
2843 void |
2086
|
2844 OCT_VAL_REP::set_index (const ColumnVector& v) |
1558
|
2845 { |
|
2846 switch (type_tag) |
|
2847 { |
|
2848 case matrix_constant: |
|
2849 matrix->set_index (v); |
|
2850 break; |
|
2851 |
2086
|
2852 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2853 complex_matrix->set_index (v); |
|
2854 break; |
|
2855 |
2086
|
2856 case OCT_VAL_REP::char_matrix_constant: |
|
2857 case OCT_VAL_REP::char_matrix_constant_str: |
1572
|
2858 char_matrix->set_index (v); |
|
2859 break; |
|
2860 |
1558
|
2861 default: |
|
2862 panic_impossible (); |
|
2863 break; |
|
2864 } |
|
2865 } |
|
2866 |
|
2867 void |
2086
|
2868 OCT_VAL_REP::set_index (const Matrix& m) |
1558
|
2869 { |
|
2870 int nr = m.rows (); |
|
2871 int nc = m.cols (); |
|
2872 |
|
2873 if (nr <= 1 || nc <= 1 |
2181
|
2874 || Vdo_fortran_indexing) |
1558
|
2875 { |
|
2876 switch (type_tag) |
|
2877 { |
|
2878 case matrix_constant: |
|
2879 matrix->set_index (m); |
|
2880 break; |
|
2881 |
2086
|
2882 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2883 complex_matrix->set_index (m); |
|
2884 break; |
|
2885 |
2086
|
2886 case OCT_VAL_REP::char_matrix_constant: |
|
2887 case OCT_VAL_REP::char_matrix_constant_str: |
1572
|
2888 char_matrix->set_index (m); |
|
2889 break; |
|
2890 |
1558
|
2891 default: |
|
2892 panic_impossible (); |
|
2893 break; |
|
2894 } |
|
2895 } |
|
2896 else |
|
2897 ::error ("invalid matrix used as index"); |
|
2898 } |
|
2899 |
|
2900 // XXX FIXME XXX -- this should probably be handled some other way... |
|
2901 // The arg here is expected to be ':'. |
|
2902 void |
2086
|
2903 OCT_VAL_REP::set_index (char c) |
1558
|
2904 { |
|
2905 switch (type_tag) |
|
2906 { |
|
2907 case matrix_constant: |
|
2908 matrix->set_index (c); |
|
2909 break; |
|
2910 |
2086
|
2911 case OCT_VAL_REP::complex_matrix_constant: |
1558
|
2912 complex_matrix->set_index (c); |
|
2913 break; |
|
2914 |
2086
|
2915 case OCT_VAL_REP::char_matrix_constant: |
|
2916 case OCT_VAL_REP::char_matrix_constant_str: |
1572
|
2917 char_matrix->set_index (c); |
|
2918 break; |
|
2919 |
1558
|
2920 default: |
|
2921 panic_impossible (); |
|
2922 break; |
|
2923 } |
|
2924 } |
|
2925 |
|
2926 void |
2086
|
2927 OCT_VAL_REP::set_index (const octave_value_list& args, bool rhs_is_complex) |
1558
|
2928 { |
2115
|
2929 // XXX FIXME XXX -- it's not good that we have to list all the types |
|
2930 // that can be indexed here. |
|
2931 |
1558
|
2932 switch (type_tag) |
|
2933 { |
|
2934 case unknown_constant: |
|
2935 case scalar_constant: |
|
2936 case complex_scalar_constant: |
|
2937 case range_constant: |
1581
|
2938 convert_to_matrix_type (rhs_is_complex); |
1558
|
2939 break; |
|
2940 |
2115
|
2941 case matrix_constant: |
|
2942 case complex_matrix_constant: |
|
2943 case char_matrix_constant: |
|
2944 case char_matrix_constant_str: |
|
2945 break; |
|
2946 |
1558
|
2947 default: |
2115
|
2948 ::error ("indexing %s type not implemented", type_as_string ()); |
1558
|
2949 break; |
|
2950 } |
|
2951 |
2115
|
2952 if (! error_state) |
1558
|
2953 { |
2115
|
2954 int n = args.length (); |
|
2955 |
|
2956 for (int i = 0; i < n; i++) |
1558
|
2957 { |
2115
|
2958 octave_value arg = args (i); |
|
2959 |
|
2960 switch (arg.const_type ()) |
|
2961 { |
|
2962 case range_constant: |
|
2963 set_index (arg.range_value ()); |
|
2964 break; |
|
2965 |
|
2966 case magic_colon: |
|
2967 set_index (':'); |
|
2968 break; |
|
2969 |
|
2970 default: |
|
2971 set_index (arg.matrix_value ()); |
|
2972 break; |
|
2973 } |
|
2974 |
|
2975 if (error_state) |
|
2976 { |
|
2977 clear_index (); |
|
2978 break; |
|
2979 } |
1558
|
2980 } |
|
2981 } |
|
2982 } |
|
2983 |
1827
|
2984 static inline bool |
2086
|
2985 valid_scalar_indices (const octave_value_list& args) |
1558
|
2986 { |
|
2987 int nargin = args.length (); |
|
2988 |
|
2989 for (int i = 0; i < nargin; i++) |
|
2990 if (! args(i).valid_as_scalar_index ()) |
1827
|
2991 return false; |
|
2992 |
|
2993 return true; |
1558
|
2994 } |
1299
|
2995 |
2086
|
2996 octave_value |
|
2997 OCT_VAL_REP::do_index (const octave_value_list& args) |
1299
|
2998 { |
2086
|
2999 octave_value retval; |
1299
|
3000 |
|
3001 if (error_state) |
|
3002 return retval; |
|
3003 |
1827
|
3004 bool originally_scalar_type = is_scalar_type (); |
1558
|
3005 |
|
3006 if (originally_scalar_type && valid_scalar_indices (args)) |
1299
|
3007 { |
1558
|
3008 switch (type_tag) |
1299
|
3009 { |
1558
|
3010 case scalar_constant: |
|
3011 retval = scalar; |
|
3012 break; |
|
3013 |
|
3014 case complex_scalar_constant: |
|
3015 retval = *complex_scalar; |
1299
|
3016 break; |
|
3017 |
|
3018 default: |
1558
|
3019 panic_impossible (); |
1299
|
3020 break; |
|
3021 } |
|
3022 } |
|
3023 else |
|
3024 { |
1558
|
3025 set_index (args); |
|
3026 |
|
3027 if (! error_state) |
|
3028 { |
|
3029 switch (type_tag) |
|
3030 { |
|
3031 case range_constant: |
|
3032 force_numeric (); |
|
3033 // Fall through... |
|
3034 |
|
3035 case matrix_constant: |
|
3036 retval = Matrix (matrix->value ()); |
|
3037 break; |
|
3038 |
|
3039 case complex_matrix_constant: |
|
3040 retval = ComplexMatrix (complex_matrix->value ()); |
|
3041 break; |
|
3042 |
1572
|
3043 case char_matrix_constant: |
|
3044 retval = charMatrix (char_matrix->value ()); |
|
3045 break; |
|
3046 |
|
3047 case char_matrix_constant_str: |
2275
|
3048 { |
|
3049 // Kluge to prevent s([]) from turning into a string |
|
3050 // with no rows... |
|
3051 charMatrix tmp (char_matrix->value ()); |
|
3052 |
|
3053 if (tmp.rows () == 0 && tmp.columns () == 0) |
|
3054 tmp.resize (1, 0); |
|
3055 |
|
3056 retval = octave_value (tmp, 1); |
|
3057 } |
1572
|
3058 break; |
|
3059 |
1558
|
3060 default: |
|
3061 error ("can't index %s variables", type_as_string ()); |
|
3062 break; |
|
3063 } |
2090
|
3064 |
|
3065 // We may have converted this value from a scalar to a |
|
3066 // matrix to allow indexing to work. |
|
3067 |
|
3068 if (! error_state) |
|
3069 maybe_mutate (); |
1558
|
3070 } |
1299
|
3071 } |
|
3072 |
|
3073 return retval; |
|
3074 } |
|
3075 |
1558
|
3076 void |
2086
|
3077 OCT_VAL_REP::maybe_widen (OCT_VAL_REP::constant_type rhs_type) |
1299
|
3078 { |
1558
|
3079 switch (type_tag) |
1299
|
3080 { |
|
3081 case matrix_constant: |
1558
|
3082 switch (rhs_type) |
|
3083 { |
|
3084 case complex_scalar_constant: |
|
3085 case complex_matrix_constant: |
1299
|
3086 { |
1558
|
3087 ComplexMatrix *cm = new ComplexMatrix (*matrix); |
|
3088 delete matrix; |
|
3089 complex_matrix = cm; |
|
3090 type_tag = complex_matrix_constant; |
1299
|
3091 } |
1558
|
3092 break; |
|
3093 |
|
3094 default: |
|
3095 break; |
|
3096 } |
1299
|
3097 break; |
|
3098 |
1572
|
3099 case char_matrix_constant: |
|
3100 switch (rhs_type) |
|
3101 { |
|
3102 case scalar_constant: |
|
3103 case matrix_constant: |
|
3104 { |
|
3105 Matrix *m = new Matrix (*char_matrix); |
|
3106 delete matrix; |
|
3107 matrix = m; |
|
3108 type_tag = matrix_constant; |
|
3109 } |
|
3110 break; |
|
3111 |
|
3112 case complex_scalar_constant: |
|
3113 case complex_matrix_constant: |
|
3114 { |
|
3115 ComplexMatrix *cm = new ComplexMatrix (*char_matrix); |
|
3116 delete matrix; |
|
3117 complex_matrix = cm; |
|
3118 type_tag = complex_matrix_constant; |
|
3119 } |
|
3120 break; |
|
3121 |
|
3122 default: |
|
3123 break; |
|
3124 } |
|
3125 break; |
|
3126 |
1299
|
3127 default: |
|
3128 break; |
|
3129 } |
|
3130 } |
|
3131 |
|
3132 // Assignment operations for the tree-constant representation class. |
|
3133 |
|
3134 // Top-level tree-constant function that handles assignments. Only |
|
3135 // decide if the left-hand side is currently a scalar or a matrix and |
|
3136 // hand off to other functions to do the real work. |
|
3137 |
1572
|
3138 // XXX FIXME XXX -- need some other way to make these functions |
|
3139 // visible here (they should be in some header file...) |
|
3140 |
1558
|
3141 extern void assign (Array2<Complex>&, const Array2<Complex>&); |
|
3142 extern void assign (Array2<Complex>&, const Array2<double>&); |
1572
|
3143 extern void assign (Array2<Complex>&, const Array2<char>&); |
|
3144 |
1558
|
3145 extern void assign (Array2<double>&, const Array2<double>&); |
1572
|
3146 extern void assign (Array2<double>&, const Array2<char>&); |
|
3147 |
|
3148 extern void assign (Array2<char>&, const Array2<char>&); |
1558
|
3149 |
1299
|
3150 void |
2086
|
3151 OCT_VAL_REP::assign (octave_value& rhs, const octave_value_list& args) |
1299
|
3152 { |
1593
|
3153 // XXX FIXME XXX -- we should probably have special cases for rhs |
|
3154 // being a range type, since converting to a matrix can waste a lot |
|
3155 // of memory. |
|
3156 |
2086
|
3157 octave_value rhs_tmp = rhs; |
1593
|
3158 |
1652
|
3159 if (! (is_string () |
|
3160 && (rhs_tmp.is_string () |
|
3161 || rhs_tmp.is_zero_by_zero ()))) |
|
3162 { |
2269
|
3163 rhs_tmp = rhs_tmp.make_numeric (); |
1652
|
3164 |
|
3165 if (error_state) |
|
3166 return; |
|
3167 } |
|
3168 |
|
3169 if (rhs_tmp.is_string |
|
3170 && rhs_tmp.rows () == 1 |
|
3171 && rhs_tmp.columns () == 0) |
|
3172 { |
2269
|
3173 rhs_tmp = rhs_tmp.make_numeric (1); |
1652
|
3174 |
|
3175 if (error_state) |
|
3176 return; |
|
3177 } |
1299
|
3178 |
1572
|
3179 // An assignment to a range will normally require a conversion to a |
|
3180 // vector in the end anyway, since it will normally destroy the |
|
3181 // equally-spaced property of the range elements. This is not as |
|
3182 // memory efficient as possible, but it is much simpler than writing |
|
3183 // additional indexing and assignment functions especially for |
|
3184 // Ranges. |
|
3185 |
|
3186 if (is_defined () && ! (is_numeric_type () || is_string ())) |
1652
|
3187 { |
|
3188 force_numeric (); |
|
3189 |
|
3190 if (error_state) |
|
3191 return; |
|
3192 } |
|
3193 |
|
3194 if (! rhs_tmp.is_zero_by_zero ()) |
|
3195 { |
|
3196 maybe_widen (rhs_tmp.const_type ()); |
|
3197 |
|
3198 if (error_state) |
|
3199 return; |
|
3200 } |
|
3201 |
|
3202 set_index (args, rhs_tmp.is_complex_type ()); |
1299
|
3203 |
|
3204 if (error_state) |
|
3205 return; |
|
3206 |
1652
|
3207 switch (type_tag) |
1299
|
3208 { |
1652
|
3209 case complex_matrix_constant: |
|
3210 { |
|
3211 switch (rhs_tmp.const_type ()) |
1299
|
3212 { |
1652
|
3213 case complex_scalar_constant: |
|
3214 case complex_matrix_constant: |
|
3215 ::assign (*complex_matrix, rhs_tmp.complex_matrix_value ()); |
|
3216 break; |
|
3217 |
|
3218 case scalar_constant: |
|
3219 case matrix_constant: |
|
3220 ::assign (*complex_matrix, rhs_tmp.matrix_value ()); |
|
3221 break; |
|
3222 |
|
3223 default: |
|
3224 panic_impossible ();; |
|
3225 break; |
1299
|
3226 } |
1652
|
3227 } |
|
3228 break; |
|
3229 |
|
3230 case scalar_constant: |
|
3231 case matrix_constant: |
|
3232 { |
|
3233 switch (rhs_tmp.const_type ()) |
1572
|
3234 { |
1652
|
3235 case scalar_constant: |
|
3236 case matrix_constant: |
|
3237 ::assign (*matrix, rhs_tmp.matrix_value ()); |
|
3238 break; |
|
3239 |
|
3240 case char_matrix_constant: |
|
3241 ::assign (*matrix, rhs_tmp.char_matrix_value ()); |
|
3242 break; |
|
3243 |
|
3244 default: |
|
3245 panic_impossible (); |
|
3246 break; |
1572
|
3247 } |
1652
|
3248 } |
|
3249 break; |
|
3250 |
|
3251 case char_matrix_constant: |
|
3252 ::assign (*char_matrix, rhs_tmp.char_matrix_value ()); |
|
3253 break; |
|
3254 |
|
3255 case char_matrix_constant_str: |
|
3256 ::assign (*char_matrix, rhs_tmp.char_matrix_value ()); |
|
3257 if (char_matrix->rows () == 0 && char_matrix->columns () == 0) |
|
3258 char_matrix->resize (1, 0); |
|
3259 break; |
|
3260 |
|
3261 default: |
|
3262 panic_impossible (); |
|
3263 break; |
1299
|
3264 } |
1688
|
3265 |
|
3266 // Do the right thing for assignments like `x(1) = pi' when x is |
|
3267 // undefined before the assignment. |
|
3268 |
|
3269 if (is_matrix_type () || is_range ()) |
|
3270 maybe_mutate (); |
1299
|
3271 } |
|
3272 |
1827
|
3273 bool |
2086
|
3274 OCT_VAL_REP::print_as_scalar (void) |
1755
|
3275 { |
|
3276 int nr = rows (); |
|
3277 int nc = columns (); |
1827
|
3278 |
1755
|
3279 return (is_scalar_type () |
|
3280 || (is_string () && nr <= 1) |
|
3281 || (is_matrix_type () |
|
3282 && ((nr == 1 && nc == 1) |
|
3283 || nr == 0 |
|
3284 || nc == 0))); |
|
3285 } |
|
3286 |
1827
|
3287 bool |
2086
|
3288 OCT_VAL_REP::print_as_structure (void) |
1755
|
3289 { |
|
3290 return is_map (); |
|
3291 } |
|
3292 |
2181
|
3293 static int |
|
3294 do_fortran_indexing (void) |
|
3295 { |
|
3296 Vdo_fortran_indexing = check_preference ("do_fortran_indexing"); |
|
3297 |
|
3298 liboctave_dfi_flag = Vdo_fortran_indexing; |
|
3299 |
|
3300 return 0; |
|
3301 } |
|
3302 |
|
3303 static int |
|
3304 implicit_str_to_num_ok (void) |
|
3305 { |
|
3306 Vimplicit_str_to_num_ok = check_preference ("implicit_str_to_num_ok"); |
|
3307 |
|
3308 return 0; |
|
3309 } |
|
3310 |
|
3311 static int |
|
3312 ok_to_lose_imaginary_part (void) |
|
3313 { |
|
3314 Vok_to_lose_imaginary_part = check_preference ("ok_to_lose_imaginary_part"); |
|
3315 |
|
3316 return 0; |
|
3317 } |
|
3318 |
|
3319 static int |
|
3320 prefer_column_vectors (void) |
|
3321 { |
|
3322 Vprefer_column_vectors |
|
3323 = check_preference ("prefer_column_vectors"); |
|
3324 |
|
3325 liboctave_pcv_flag = Vprefer_column_vectors; |
|
3326 |
|
3327 return 0; |
|
3328 } |
|
3329 |
|
3330 static int |
|
3331 prefer_zero_one_indexing (void) |
|
3332 { |
|
3333 Vprefer_zero_one_indexing = check_preference ("prefer_zero_one_indexing"); |
|
3334 |
|
3335 liboctave_pzo_flag = Vprefer_zero_one_indexing; |
|
3336 |
|
3337 return 0; |
|
3338 } |
|
3339 |
|
3340 static int |
2187
|
3341 print_answer_id_name (void) |
|
3342 { |
|
3343 Vprint_answer_id_name = check_preference ("print_answer_id_name"); |
|
3344 |
|
3345 return 0; |
|
3346 } |
|
3347 |
|
3348 static int |
2181
|
3349 propagate_empty_matrices (void) |
|
3350 { |
|
3351 Vpropagate_empty_matrices = check_preference ("propagate_empty_matrices"); |
|
3352 |
|
3353 return 0; |
|
3354 } |
|
3355 |
|
3356 static int |
|
3357 resize_on_range_error (void) |
|
3358 { |
|
3359 Vresize_on_range_error = check_preference ("resize_on_range_error"); |
|
3360 |
|
3361 liboctave_rre_flag = Vresize_on_range_error; |
|
3362 |
|
3363 return 0; |
|
3364 } |
|
3365 |
|
3366 static int |
|
3367 struct_levels_to_print (void) |
|
3368 { |
|
3369 double val; |
|
3370 if (builtin_real_scalar_variable ("struct_levels_to_print", val) |
|
3371 && ! xisnan (val)) |
|
3372 { |
|
3373 int ival = NINT (val); |
|
3374 if (ival >= 0 && (double) ival == val) |
|
3375 { |
|
3376 Vstruct_levels_to_print = ival; |
|
3377 return 0; |
|
3378 } |
|
3379 } |
|
3380 gripe_invalid_value_specified ("struct_levels_to_print"); |
|
3381 return -1; |
|
3382 } |
|
3383 |
|
3384 void |
|
3385 symbols_of_pt_const (void) |
|
3386 { |
|
3387 DEFVAR (do_fortran_indexing, 0.0, 0, do_fortran_indexing, |
|
3388 "allow single indices for matrices"); |
|
3389 |
|
3390 DEFVAR (implicit_str_to_num_ok, 0.0, 0, implicit_str_to_num_ok, |
|
3391 "allow implicit string to number conversion"); |
|
3392 |
|
3393 DEFVAR (ok_to_lose_imaginary_part, "warn", 0, ok_to_lose_imaginary_part, |
|
3394 "silently convert from complex to real by dropping imaginary part"); |
|
3395 |
|
3396 DEFVAR (prefer_column_vectors, 1.0, 0, prefer_column_vectors, |
|
3397 "prefer column/row vectors"); |
|
3398 |
|
3399 DEFVAR (prefer_zero_one_indexing, 0.0, 0, prefer_zero_one_indexing, |
|
3400 "when there is a conflict, prefer zero-one style indexing"); |
|
3401 |
2187
|
3402 DEFVAR (print_answer_id_name, 1.0, 0, print_answer_id_name, |
|
3403 "set output style to print `var_name = ...'"); |
|
3404 |
2181
|
3405 DEFVAR (propagate_empty_matrices, 1.0, 0, propagate_empty_matrices, |
|
3406 "operations on empty matrices return an empty matrix, not an error"); |
|
3407 |
|
3408 DEFVAR (resize_on_range_error, 1.0, 0, resize_on_range_error, |
|
3409 "enlarge matrices on assignment"); |
|
3410 |
|
3411 DEFVAR (struct_levels_to_print, 2.0, 0, struct_levels_to_print, |
|
3412 "number of levels of structure elements to print"); |
|
3413 } |
|
3414 |
96
|
3415 /* |
1
|
3416 ;;; Local Variables: *** |
|
3417 ;;; mode: C++ *** |
|
3418 ;;; End: *** |
|
3419 */ |