Mercurial > octave-antonio
annotate liboctave/util/lo-utils.cc @ 17744:d63878346099
maint: Update copyright notices for release.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Wed, 23 Oct 2013 22:09:27 -0400 |
parents | 7b305b36b87e |
children | 49a5a4be04a1 |
rev | line source |
---|---|
1993 | 1 // utils.cc |
1967 | 2 /* |
3 | |
17744
d63878346099
maint: Update copyright notices for release.
John W. Eaton <jwe@octave.org>
parents:
17623
diff
changeset
|
4 Copyright (C) 1996-2013 John W. Eaton |
1967 | 5 |
6 This file is part of Octave. | |
7 | |
8 Octave is free software; you can redistribute it and/or modify it | |
9 under the terms of the GNU General Public License as published by the | |
7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
11 option) any later version. | |
1967 | 12 |
13 Octave is distributed in the hope that it will be useful, but WITHOUT | |
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
16 for more details. | |
17 | |
18 You should have received a copy of the GNU General Public License | |
7016 | 19 along with Octave; see the file COPYING. If not, see |
20 <http://www.gnu.org/licenses/>. | |
1967 | 21 |
22 */ | |
23 | |
24 #ifdef HAVE_CONFIG_H | |
25 #include <config.h> | |
26 #endif | |
27 | |
6907 | 28 #include <cctype> |
17623
7b305b36b87e
use setenv to implement putenv and avoid memory leak (bug #40230)
John W. Eaton <jwe@octave.org>
parents:
16971
diff
changeset
|
29 #include <cerrno> |
10463
bbe99b2a5ba7
undo recent gnulib-related changes
John W. Eaton <jwe@octave.org>
parents:
10447
diff
changeset
|
30 #include <cstdlib> |
2926 | 31 #include <cstdio> |
7048 | 32 #include <cstring> |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
33 #include <cfloat> |
1967 | 34 |
6490 | 35 #include <limits> |
2926 | 36 #include <string> |
37 | |
38 #include <sys/types.h> | |
39 #include <unistd.h> | |
40 | |
10068
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
41 #include "quit.h" |
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
42 |
2926 | 43 #include "lo-error.h" |
4130 | 44 #include "lo-ieee.h" |
1967 | 45 #include "lo-mappers.h" |
46 #include "lo-utils.h" | |
47 | |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
48 bool xis_int_or_inf_or_nan (double x) |
11013
63f79f798a14
fix small typos in new tests
Jaroslav Hajek <highegg@gmail.com>
parents:
11010
diff
changeset
|
49 { return xisnan (x) || D_NINT (x) == x; } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
50 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
51 bool xis_one_or_zero (double x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
52 { return x == 0 || x == 1; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
53 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
54 bool xis_zero (double x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
55 { return x == 0; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
56 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
57 bool xtoo_large_for_float (double x) |
15213
336f42406671
use numeric_limits functions instead of DBL_MIN, DBL_MAX, etc.
John W. Eaton <jwe@octave.org>
parents:
15212
diff
changeset
|
58 { |
16971
259c1f295a1e
Use xfinite to replace some (isinf || isnan) instances in C++ code.
Rik <rik@octave.org>
parents:
15454
diff
changeset
|
59 return (xfinite (x) && fabs (x) > std::numeric_limits<float>::max ()); |
15213
336f42406671
use numeric_limits functions instead of DBL_MIN, DBL_MAX, etc.
John W. Eaton <jwe@octave.org>
parents:
15212
diff
changeset
|
60 } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
61 |
15212
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
62 bool xtoo_large_for_float (const Complex& x) |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
63 { |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
64 return (xtoo_large_for_float (x.real ()) |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
65 || xtoo_large_for_float (x.imag ())); |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
66 } |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
67 |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
68 bool xis_int_or_inf_or_nan (float x) |
11013
63f79f798a14
fix small typos in new tests
Jaroslav Hajek <highegg@gmail.com>
parents:
11010
diff
changeset
|
69 { return xisnan (x) || D_NINT (x) == x; } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
70 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
71 bool xis_one_or_zero (float x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
72 { return x == 0 || x == 1; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
73 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
74 bool xis_zero (float x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
75 { return x == 0; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
76 |
2926 | 77 // Save a string. |
78 | |
79 char * | |
80 strsave (const char *s) | |
81 { | |
82 if (! s) | |
83 return 0; | |
84 | |
85 int len = strlen (s); | |
86 char *tmp = new char [len+1]; | |
87 tmp = strcpy (tmp, s); | |
88 return tmp; | |
89 } | |
90 | |
91 void | |
3504 | 92 octave_putenv (const std::string& name, const std::string& value) |
2926 | 93 { |
17623
7b305b36b87e
use setenv to implement putenv and avoid memory leak (bug #40230)
John W. Eaton <jwe@octave.org>
parents:
16971
diff
changeset
|
94 if (gnulib::setenv (name.c_str (), value.c_str (), true) < 0) |
7b305b36b87e
use setenv to implement putenv and avoid memory leak (bug #40230)
John W. Eaton <jwe@octave.org>
parents:
16971
diff
changeset
|
95 (*current_liboctave_error_handler) ("putenv: %s", |
7b305b36b87e
use setenv to implement putenv and avoid memory leak (bug #40230)
John W. Eaton <jwe@octave.org>
parents:
16971
diff
changeset
|
96 gnulib::strerror (errno)); |
2926 | 97 } |
98 | |
3504 | 99 std::string |
2926 | 100 octave_fgets (FILE *f) |
101 { | |
4527 | 102 bool eof; |
103 return octave_fgets (f, eof); | |
104 } | |
105 | |
106 std::string | |
107 octave_fgets (FILE *f, bool& eof) | |
108 { | |
109 eof = false; | |
110 | |
3504 | 111 std::string retval; |
2926 | 112 |
113 int grow_size = 1024; | |
114 int max_size = grow_size; | |
115 | |
10411 | 116 char *buf = static_cast<char *> (gnulib::malloc (max_size)); |
2926 | 117 char *bufptr = buf; |
118 int len = 0; | |
119 | |
120 do | |
121 { | |
12918
f3a8d1efe2c1
use gnulib:: qualifiers for more stdio functions
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
122 if (gnulib::fgets (bufptr, grow_size, f)) |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
123 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
124 len = strlen (bufptr); |
2926 | 125 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
126 if (len == grow_size - 1) |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
127 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
128 int tmp = bufptr - buf + grow_size - 1; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
129 grow_size *= 2; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
130 max_size += grow_size; |
10411 | 131 buf = static_cast<char *> (gnulib::realloc (buf, max_size)); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
132 bufptr = buf + tmp; |
2926 | 133 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
134 if (*(bufptr-1) == '\n') |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
135 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
136 *bufptr = '\0'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
137 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
138 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
139 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
140 else if (bufptr[len-1] != '\n') |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
141 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
142 bufptr[len++] = '\n'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
143 bufptr[len] = '\0'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
144 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
145 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
146 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
147 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
148 } |
2926 | 149 else |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
150 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
151 if (len == 0) |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
152 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
153 eof = true; |
4527 | 154 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
155 free (buf); |
2926 | 156 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
157 buf = 0; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
158 } |
2926 | 159 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
160 break; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
161 } |
2926 | 162 } |
163 while (retval.empty ()); | |
164 | |
165 if (buf) | |
166 free (buf); | |
167 | |
10142
829e69ec3110
make OCTAVE_QUIT a function
Jaroslav Hajek <highegg@gmail.com>
parents:
10068
diff
changeset
|
168 octave_quit (); |
10068
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
169 |
2926 | 170 return retval; |
171 } | |
172 | |
3970 | 173 std::string |
174 octave_fgetl (FILE *f) | |
175 { | |
4527 | 176 bool eof; |
177 return octave_fgetl (f, eof); | |
178 } | |
179 | |
180 std::string | |
181 octave_fgetl (FILE *f, bool& eof) | |
182 { | |
183 std::string retval = octave_fgets (f, eof); | |
3970 | 184 |
185 size_t len = retval.length (); | |
186 | |
187 if (retval[len-1] == '\n') | |
188 retval.resize (len-1); | |
189 | |
190 return retval; | |
191 } | |
192 | |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
193 // Note that the caller is responsible for repositioning the stream on |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
194 // failure. |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
195 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
196 template <typename T> |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
197 T |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
198 read_inf_nan_na (std::istream& is, char c0) |
4130 | 199 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
200 T val = 0.0; |
4130 | 201 |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
202 switch (c0) |
4130 | 203 { |
204 case 'i': case 'I': | |
205 { | |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
206 char c1 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
207 if (c1 == 'n' || c1 == 'N') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
208 { |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
209 char c2 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
210 if (c2 == 'f' || c2 == 'F') |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
211 val = std::numeric_limits<T>::infinity (); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
212 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
213 is.setstate (std::ios::failbit); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
214 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
215 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
216 is.setstate (std::ios::failbit); |
4130 | 217 } |
218 break; | |
219 | |
220 case 'n': case 'N': | |
221 { | |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
222 char c1 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
223 if (c1 == 'a' || c1 == 'A') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
224 { |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
225 char c2 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
226 if (c2 == 'n' || c2 == 'N') |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
227 val = std::numeric_limits<T>::quiet_NaN (); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
228 else |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
229 val = octave_numeric_limits<T>::NA (); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
230 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
231 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
232 is.setstate (std::ios::failbit); |
4130 | 233 } |
234 break; | |
235 | |
236 default: | |
237 abort (); | |
238 } | |
239 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
240 return val; |
4130 | 241 } |
242 | |
12936
b74cb659e757
accept but discard sign when reading NA and NaN values
John W. Eaton <jwe@octave.org>
parents:
12918
diff
changeset
|
243 // Read a double value. Discard any sign on NaN and NA. |
b74cb659e757
accept but discard sign when reading NA and NaN values
John W. Eaton <jwe@octave.org>
parents:
12918
diff
changeset
|
244 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
245 template <typename T> |
4130 | 246 double |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
247 octave_read_fp_value (std::istream& is) |
4130 | 248 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
249 T val = 0.0; |
4130 | 250 |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
251 // FIXME -- resetting stream position is likely to fail unless we are |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
252 // reading from a file. |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
253 std::ios::streampos pos = is.tellg (); |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
254 |
6907 | 255 char c1 = ' '; |
4130 | 256 |
6907 | 257 while (isspace (c1)) |
258 c1 = is.get (); | |
259 | |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
260 bool neg = false; |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
261 |
6194 | 262 switch (c1) |
4130 | 263 { |
6194 | 264 case '-': |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
265 neg = true; |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
266 // fall through... |
6194 | 267 |
268 case '+': | |
269 { | |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
270 char c2 = 0; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
271 c2 = is.get (); |
12936
b74cb659e757
accept but discard sign when reading NA and NaN values
John W. Eaton <jwe@octave.org>
parents:
12918
diff
changeset
|
272 if (c2 == 'i' || c2 == 'I' || c2 == 'n' || c2 == 'N') |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
273 val = read_inf_nan_na<T> (is, c2); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
274 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
275 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
276 is.putback (c2); |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
277 is >> val; |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
278 } |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
279 |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
280 if (neg && ! is.fail ()) |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
281 val = -val; |
6194 | 282 } |
283 break; | |
284 | |
4130 | 285 case 'i': case 'I': |
286 case 'n': case 'N': | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
287 val = read_inf_nan_na<T> (is, c1); |
4130 | 288 break; |
289 | |
290 default: | |
6194 | 291 is.putback (c1); |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
292 is >> val; |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
293 break; |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
294 } |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
295 |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
296 std::ios::iostate status = is.rdstate (); |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
297 if (status & std::ios::failbit) |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
298 { |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
299 is.clear (); |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
300 is.seekg (pos); |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
301 is.setstate (status); |
4130 | 302 } |
303 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
304 return val; |
4130 | 305 } |
306 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
307 template <typename T> |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
308 std::complex<T> |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
309 octave_read_cx_fp_value (std::istream& is) |
4130 | 310 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
311 T re = 0.0, im = 0.0; |
4130 | 312 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
313 std::complex<T> cx = 0.0; |
4130 | 314 |
6907 | 315 char ch = ' '; |
4130 | 316 |
6907 | 317 while (isspace (ch)) |
318 ch = is.get (); | |
4130 | 319 |
320 if (ch == '(') | |
321 { | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
322 re = octave_read_value<T> (is); |
6897 | 323 ch = is.get (); |
4130 | 324 |
325 if (ch == ',') | |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
326 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
327 im = octave_read_value<T> (is); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
328 ch = is.get (); |
4130 | 329 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
330 if (ch == ')') |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
331 cx = std::complex<T> (re, im); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
332 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
333 is.setstate (std::ios::failbit); |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
334 } |
4130 | 335 else if (ch == ')') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
336 cx = re; |
4130 | 337 else |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
338 is.setstate (std::ios::failbit); |
4130 | 339 } |
340 else | |
341 { | |
342 is.putback (ch); | |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
343 cx = octave_read_value<double> (is); |
4130 | 344 } |
345 | |
346 return cx; | |
347 } | |
348 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
349 template <> OCTAVE_API double octave_read_value (std::istream& is) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
350 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
351 return octave_read_fp_value<double> (is); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
352 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
353 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
354 template <> OCTAVE_API Complex octave_read_value (std::istream& is) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
355 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
356 return octave_read_cx_fp_value<double> (is); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
357 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
358 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
359 template <> OCTAVE_API float octave_read_value (std::istream& is) |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
360 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
361 return octave_read_fp_value<float> (is); |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
362 } |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
363 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
364 template <> OCTAVE_API FloatComplex octave_read_value (std::istream& is) |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
365 { |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
366 return octave_read_cx_fp_value<float> (is); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
367 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
368 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
369 void |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
370 octave_write_double (std::ostream& os, double d) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
371 { |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
372 if (lo_ieee_is_NA (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
373 os << "NA"; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
374 else if (lo_ieee_isnan (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
375 os << "NaN"; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
376 else if (lo_ieee_isinf (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
377 os << (d < 0 ? "-Inf" : "Inf"); |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
378 else |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
379 os << d; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
380 } |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
381 |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
382 void |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
383 octave_write_complex (std::ostream& os, const Complex& c) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
384 { |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
385 os << "("; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
386 octave_write_double (os, real (c)); |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
387 os << ","; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
388 octave_write_double (os, imag (c)); |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
389 os << ")"; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
390 } |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
391 |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
392 void |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
393 octave_write_float (std::ostream& os, float d) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
394 { |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
395 if (lo_ieee_is_NA (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
396 os << "NA"; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
397 else if (lo_ieee_isnan (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
398 os << "NaN"; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
399 else if (lo_ieee_isinf (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
400 os << (d < 0 ? "-Inf" : "Inf"); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
401 else |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
402 os << d; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
403 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
404 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
405 void |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
406 octave_write_float_complex (std::ostream& os, const FloatComplex& c) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
407 { |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
408 os << "("; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
409 octave_write_float (os, real (c)); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
410 os << ","; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
411 octave_write_float (os, imag (c)); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
412 os << ")"; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
413 } |