annotate src/OPERATORS/op-range.cc @ 7789:82be108cc558

First attempt at single precision tyeps * * * corrections to qrupdate single precision routines * * * prefer demotion to single over promotion to double * * * Add single precision support to log2 function * * * Trivial PROJECT file update * * * Cache optimized hermitian/transpose methods * * * Add tests for tranpose/hermitian and ChangeLog entry for new transpose code
author David Bateman <dbateman@free.fr>
date Sun, 27 Apr 2008 22:34:17 +0200
parents a1dbe9d80eee
children 283989f2da9b
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
1 /*
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
2
7017
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
3 Copyright (C) 1996, 1997, 1998, 2000, 2002, 2004, 2005, 2007
a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe]
jwe
parents: 7016
diff changeset
4 John W. Eaton
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
5
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
6 This file is part of Octave.
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
7
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
8 Octave is free software; you can redistribute it and/or modify it
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
9 under the terms of the GNU General Public License as published by the
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
10 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
11 option) any later version.
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
12
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
13 Octave is distributed in the hope that it will be useful, but WITHOUT
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
16 for more details.
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
17
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
18 You should have received a copy of the GNU General Public License
7016
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
19 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 5307
diff changeset
20 <http://www.gnu.org/licenses/>.
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
21
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
22 */
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
23
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
24 #ifdef HAVE_CONFIG_H
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
25 #include <config.h>
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
26 #endif
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
27
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
28 #include "gripes.h"
4055
7690958e7726 [project @ 2002-08-18 01:56:58 by jwe]
jwe
parents: 3538
diff changeset
29 #include "oct-obj.h"
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
30 #include "ov.h"
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
31 #include "ov-range.h"
4915
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
32 #include "ov-ch-mat.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
33 #include "ov-scalar.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
34 #include "ov-re-mat.h"
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
35 #include "ov-flt-re-mat.h"
4915
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
36 #include "ov-complex.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
37 #include "ov-cx-mat.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
38 #include "ov-bool.h"
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
39 #include "ov-bool-mat.h"
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
40 #include "ov-typeinfo.h"
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
41 #include "ops.h"
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
42
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
43 // range unary ops.
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
44
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
45 DEFUNOP (not, range)
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
46 {
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
47 CAST_UNOP_ARG (const octave_range&);
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
48
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
49 return octave_value (! v.matrix_value());
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
50 }
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
51
4965
c0d8e8afa82f [project @ 2004-09-06 20:19:57 by jwe]
jwe
parents: 4915
diff changeset
52 DEFUNOP_OP (uplus, range, /* no-op */)
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
53 DEFUNOP_OP (uminus, range, -)
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
54
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
55 DEFUNOP (transpose, range)
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
56 {
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
57 CAST_UNOP_ARG (const octave_range&);
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
58
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
59 return octave_value (v.matrix_value().transpose ());
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
60 }
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
61
4915
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
62 DEFNDCATOP_FN (r_r, range, range, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
63 DEFNDCATOP_FN (r_s, range, scalar, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
64 DEFNDCATOP_FN (r_m, range, matrix, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
65 DEFNDCATOP_FN (r_cs, range, complex, array, complex_array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
66 DEFNDCATOP_FN (r_cm, range, complex_matrix, array, complex_array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
67 DEFNDCATOP_FN (r_b, range, bool, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
68 DEFNDCATOP_FN (r_bm, range, bool_matrix, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
69 DEFNDCATOP_FN (r_chm, range, char_matrix, array, char_array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
70 DEFNDCATOP_FN (s_r, scalar, range, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
71 DEFNDCATOP_FN (m_r, matrix, range, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
72 DEFNDCATOP_FN (cs_r, complex, range, complex_array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
73 DEFNDCATOP_FN (cm_r, complex_matrix, range, complex_array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
74 DEFNDCATOP_FN (b_r, bool, range, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
75 DEFNDCATOP_FN (bm_r, bool_matrix, range, array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
76 DEFNDCATOP_FN (chm_r, char_matrix, range, char_array, array, concat)
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
77
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
78 CONVDECL (range_to_float_matrix)
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
79 {
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
80 CAST_CONV_ARG (const octave_range&);
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
81
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
82 return new octave_float_matrix (FloatNDArray (v.array_value ()));
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
83 }
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
84
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
85 void
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
86 install_range_ops (void)
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
87 {
3538
0ff7323dab8b [project @ 2000-02-02 12:36:25 by jwe]
jwe
parents: 3207
diff changeset
88 INSTALL_UNOP (op_not, octave_range, not);
4965
c0d8e8afa82f [project @ 2004-09-06 20:19:57 by jwe]
jwe
parents: 4915
diff changeset
89 INSTALL_UNOP (op_uplus, octave_range, uplus);
3538
0ff7323dab8b [project @ 2000-02-02 12:36:25 by jwe]
jwe
parents: 3207
diff changeset
90 INSTALL_UNOP (op_uminus, octave_range, uminus);
0ff7323dab8b [project @ 2000-02-02 12:36:25 by jwe]
jwe
parents: 3207
diff changeset
91 INSTALL_UNOP (op_transpose, octave_range, transpose);
0ff7323dab8b [project @ 2000-02-02 12:36:25 by jwe]
jwe
parents: 3207
diff changeset
92 INSTALL_UNOP (op_hermitian, octave_range, transpose);
4915
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
93
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
94 INSTALL_CATOP (octave_range, octave_range, r_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
95 INSTALL_CATOP (octave_range, octave_scalar, r_s);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
96 INSTALL_CATOP (octave_range, octave_matrix, r_m);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
97 INSTALL_CATOP (octave_range, octave_complex, r_cs);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
98 INSTALL_CATOP (octave_range, octave_complex_matrix, r_cm);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
99 INSTALL_CATOP (octave_range, octave_bool, r_b);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
100 INSTALL_CATOP (octave_range, octave_bool_matrix, r_bm);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
101 INSTALL_CATOP (octave_range, octave_char_matrix, r_chm);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
102 INSTALL_CATOP (octave_scalar, octave_range, s_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
103 INSTALL_CATOP (octave_matrix, octave_range, m_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
104 INSTALL_CATOP (octave_complex, octave_range, cs_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
105 INSTALL_CATOP (octave_complex_matrix, octave_range, cm_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
106 INSTALL_CATOP (octave_bool, octave_range, b_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
107 INSTALL_CATOP (octave_bool_matrix, octave_range, bm_r);
c638c144d4da [project @ 2004-07-23 19:01:22 by jwe]
jwe
parents: 4192
diff changeset
108 INSTALL_CATOP (octave_char_matrix, octave_range, chm_r);
7789
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
109
82be108cc558 First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents: 7017
diff changeset
110 INSTALL_CONVOP (octave_range, octave_float_matrix, range_to_float_matrix);
3207
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
111 }
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
112
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
113 /*
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
114 ;;; Local Variables: ***
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
115 ;;; mode: C++ ***
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
116 ;;; End: ***
dddfaa93a99c [project @ 1998-11-03 04:04:37 by jwe]
jwe
parents:
diff changeset
117 */