annotate liboctave/array/Range.h @ 20603:f61c67865d9f

Don't return A for inv (A) when A is a singular Diagonal matrix (bug #46103). * dDiagMatrix.cc (inverse), fDiagMatrix.cc (inverse): return octave_Inf for the inverse of a diagonal element whose value is 0. * test/diag-perm.tst: Add tests for new behavior.
author Rik <rik@octave.org>
date Wed, 07 Oct 2015 07:41:44 -0700
parents 96153b16febe
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
1 /*
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
2
19731
4197fc428c7d maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents: 19322
diff changeset
3 Copyright (C) 1993-2015 John W. Eaton
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
4
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
5 This file is part of Octave.
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
6
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
7 Octave is free software; you can redistribute it and/or modify it
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
8 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: 6108
diff changeset
9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6108
diff changeset
10 option) any later version.
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
11
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
12 Octave is distributed in the hope that it will be useful, but WITHOUT
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
15 for more details.
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
16
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
17 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: 6108
diff changeset
18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe]
jwe
parents: 6108
diff changeset
19 <http://www.gnu.org/licenses/>.
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
20
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
21 */
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
22
382
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
23 #if !defined (octave_Range_h)
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
24 #define octave_Range_h 1
e02d6b664394 [project @ 1994-03-09 21:08:03 by jwe]
jwe
parents: 238
diff changeset
25
8950
d865363208d6 include <iosfwd> instead of <iostream> in header files
John W. Eaton <jwe@octave.org>
parents: 8920
diff changeset
26 #include <iosfwd>
7458
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
27
4810
72a6d410a14a [project @ 2004-03-03 18:49:39 by jwe]
jwe
parents: 4192
diff changeset
28 #include "dMatrix.h"
7458
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
29 #include "oct-sort.h"
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
30
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
31 class
6108
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe]
jwe
parents: 5347
diff changeset
32 OCTAVE_API
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
33 Range
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
34 {
17769
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
35 public:
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
36
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
37 Range (void)
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
38 : rng_base (0), rng_limit (0), rng_inc (0), rng_numel (0), cache (1, 0) { }
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
39
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
40 Range (const Range& r)
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
41 : rng_base (r.rng_base), rng_limit (r.rng_limit), rng_inc (r.rng_inc),
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
42 rng_numel (r.rng_numel), cache (r.cache) { }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
43
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
44 Range (double b, double l)
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
45 : rng_base (b), rng_limit (l), rng_inc (1),
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
46 rng_numel (numel_internal ()), cache ()
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
47 {
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
48 double tmplimit = rng_limit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
49
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
50 if (rng_inc > 0)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
51 tmplimit = max ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
52 else
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
53 tmplimit = min ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
54
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
55 if (tmplimit != rng_limit)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
56 rng_limit = tmplimit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
57 }
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
58
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
59 Range (double b, double l, double i)
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
60 : rng_base (b), rng_limit (l), rng_inc (i),
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
61 rng_numel (numel_internal ()), cache ()
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
62 {
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
63 double tmplimit = rng_limit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
64
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
65 if (rng_inc > 0)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
66 tmplimit = max ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
67 else
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
68 tmplimit = min ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
69
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
70 if (tmplimit != rng_limit)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
71 rng_limit = tmplimit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
72 }
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
73
8589
0131fa223dbc make length invariant in range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 8553
diff changeset
74 // For operators' usage (to preserve element count).
16169
0303fda3e929 Fix range behavior with -0 endpoints (bug #38423)
Rik <rik@octave.org>
parents: 15271
diff changeset
75 Range (double b, double i, octave_idx_type n)
0303fda3e929 Fix range behavior with -0 endpoints (bug #38423)
Rik <rik@octave.org>
parents: 15271
diff changeset
76 : rng_base (b), rng_limit (b + (n-1) * i), rng_inc (i),
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
77 rng_numel (n), cache ()
17769
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
78 {
19322
d60b44acfff7 Change a bitwise OR to a logical OR in Range.h.
Daniel J Sebald <daniel.sebald@ieee.org>
parents: 18886
diff changeset
79 if (! xfinite (b) || ! xfinite (i) || ! xfinite (rng_limit))
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
80 rng_numel = -2;
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
81 else
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
82 {
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
83 // Code below is only needed if the resulting range must be 100%
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
84 // correctly constructed. If the Range object created is only
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
85 // a temporary one used by operators this may be unnecessary.
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
86 double tmplimit = rng_limit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
87
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
88 if (rng_inc > 0)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
89 tmplimit = max ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
90 else
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
91 tmplimit = min ();
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
92
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
93 if (tmplimit != rng_limit)
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
94 rng_limit = tmplimit;
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
95 }
17769
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
96 }
8589
0131fa223dbc make length invariant in range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 8553
diff changeset
97
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
98 double base (void) const { return rng_base; }
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
99 double limit (void) const { return rng_limit; }
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
100 double inc (void) const { return rng_inc; }
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
101 GCC_ATTR_DEPRECATED octave_idx_type nelem (void) const { return numel (); }
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
102 octave_idx_type numel (void) const { return rng_numel; }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
103
2383
096529066838 [project @ 1996-10-12 17:50:47 by jwe]
jwe
parents: 1993
diff changeset
104 bool all_elements_are_ints (void) const;
096529066838 [project @ 1996-10-12 17:50:47 by jwe]
jwe
parents: 1993
diff changeset
105
645
21cb3c7dc9e3 [project @ 1994-08-24 06:19:54 by jwe]
jwe
parents: 461
diff changeset
106 Matrix matrix_value (void) const;
21cb3c7dc9e3 [project @ 1994-08-24 06:19:54 by jwe]
jwe
parents: 461
diff changeset
107
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
108 double min (void) const;
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
109 double max (void) const;
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
110
7458
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
111 void sort_internal (bool ascending = true);
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
112 void sort_internal (Array<octave_idx_type>& sidx, bool ascending = true);
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
113
7620
36594d5bbe13 Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents: 7463
diff changeset
114 Matrix diag (octave_idx_type k = 0) const;
36594d5bbe13 Move diag function into the octave_value class
David Bateman <dbateman@free.fr>
parents: 7463
diff changeset
115
7463
2467639bd8c0 eliminate UNDEFINED sort mode
John W. Eaton <jwe@octave.org>
parents: 7458
diff changeset
116 Range sort (octave_idx_type dim = 0, sortmode mode = ASCENDING) const;
7458
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
117
1032e24f199f make sort work for ranges
John W. Eaton <jwe@octave.org>
parents: 7017
diff changeset
118 Range sort (Array<octave_idx_type>& sidx, octave_idx_type dim = 0,
10312
cbc402e64d83 untabify liboctave header files
John W. Eaton <jwe@octave.org>
parents: 10158
diff changeset
119 sortmode mode = ASCENDING) const;
208
380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe]
jwe
parents: 3
diff changeset
120
8721
e9cb742df9eb imported patch sort3.diff
Jaroslav Hajek <highegg@gmail.com>
parents: 8589
diff changeset
121 sortmode is_sorted (sortmode mode = ASCENDING) const;
e9cb742df9eb imported patch sort3.diff
Jaroslav Hajek <highegg@gmail.com>
parents: 8589
diff changeset
122
9986
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
123 // Support for single-index subscripting, without generating matrix cache.
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
124
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
125 double checkelem (octave_idx_type i) const;
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
126
16169
0303fda3e929 Fix range behavior with -0 endpoints (bug #38423)
Rik <rik@octave.org>
parents: 15271
diff changeset
127 double elem (octave_idx_type i) const;
9986
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
128
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
129 Array<double> index (const idx_vector& i) const;
672e1b49e01e optimize indexing of ranges by single subscripts
Jaroslav Hajek <highegg@gmail.com>
parents: 8971
diff changeset
130
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
131 void set_base (double b);
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
132
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
133 void set_limit (double l);
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
134
20548
96153b16febe Overhaul Range object in liboctave.
Rik <rik@octave.org>
parents: 20263
diff changeset
135 void set_inc (double i);
1528
dc527156c38c [project @ 1995-10-05 01:44:18 by jwe]
jwe
parents: 1315
diff changeset
136
17769
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
137 friend OCTAVE_API std::ostream& operator << (std::ostream& os,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
138 const Range& r);
6108
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe]
jwe
parents: 5347
diff changeset
139 friend OCTAVE_API std::istream& operator >> (std::istream& is, Range& r);
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
140
8971
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
141 friend OCTAVE_API Range operator - (const Range& r);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
142 friend OCTAVE_API Range operator + (double x, const Range& r);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
143 friend OCTAVE_API Range operator + (const Range& r, double x);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
144 friend OCTAVE_API Range operator - (double x, const Range& r);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
145 friend OCTAVE_API Range operator - (const Range& r, double x);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
146 friend OCTAVE_API Range operator * (double x, const Range& r);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
147 friend OCTAVE_API Range operator * (const Range& r, double x);
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
148
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
149 void print_range (void);
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
150
17769
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents: 17744
diff changeset
151 private:
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
152
208
380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe]
jwe
parents: 3
diff changeset
153 double rng_base;
380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe]
jwe
parents: 3
diff changeset
154 double rng_limit;
380cb9ccc889 [project @ 1993-11-12 10:35:05 by jwe]
jwe
parents: 3
diff changeset
155 double rng_inc;
1860
821870c30840 [project @ 1996-02-04 10:04:58 by jwe]
jwe
parents: 1528
diff changeset
156
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
157 octave_idx_type rng_numel;
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
158
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
159 mutable Matrix cache;
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
160
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
161 octave_idx_type numel_internal (void) const;
4811
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
162
198f3712c692 [project @ 2004-03-03 19:24:33 by jwe]
jwe
parents: 4810
diff changeset
163 void clear_cache (void) const { cache.resize (0, 0); }
8971
967a692ddfe2 fix range arithmetics
Jaroslav Hajek <highegg@gmail.com>
parents: 8950
diff changeset
164
18886
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
165 protected:
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
166
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
167 // For operators' usage (to allow all values to be set directly).
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
168 Range (double b, double l, double i, octave_idx_type n)
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
169 : rng_base (b), rng_limit (l), rng_inc (i),
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
170 rng_numel (n), cache ()
18886
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
171 {
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
172 if (! xfinite (b) || ! xfinite (i) || ! xfinite (l))
20263
00cf2847355d Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents: 19731
diff changeset
173 rng_numel = -2;
18886
47d4b680d0e0 improve accuracy of range/scalar arithmetic (bug #42589)
John W. Eaton <jwe@octave.org>
parents: 17769
diff changeset
174 }
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
175 };
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
176
6108
143b556ce725 [project @ 2006-10-27 01:45:54 by jwe]
jwe
parents: 5347
diff changeset
177 extern OCTAVE_API Range operator - (const Range& r);
2599
9705ea40f0a2 [project @ 1996-12-19 04:49:09 by jwe]
jwe
parents: 2383
diff changeset
178
8553
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
179 extern OCTAVE_API Range operator + (double x, const Range& r);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
180
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
181 extern OCTAVE_API Range operator + (const Range& r, double x);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
182
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
183 extern OCTAVE_API Range operator - (double x, const Range& r);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
184
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
185 extern OCTAVE_API Range operator - (const Range& r, double x);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
186
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
187 extern OCTAVE_API Range operator * (double x, const Range& r);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
188
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
189 extern OCTAVE_API Range operator * (const Range& r, double x);
c7ff200e45f5 optimize range-scalar ops
Jaroslav Hajek <highegg@gmail.com>
parents: 7620
diff changeset
190
3
9a4c07481e61 [project @ 1993-08-08 01:20:23 by jwe]
jwe
parents:
diff changeset
191 #endif