Mercurial > octave
annotate liboctave/util/lo-utils.cc @ 25054:6652d3823428 stable
maint: Update copyright dates in all source files.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 30 Mar 2018 09:19:05 -0400 |
parents | 194eb4bd202b |
children | cb1606f78f6b |
rev | line source |
---|---|
1967 | 1 /* |
2 | |
25054
6652d3823428
maint: Update copyright dates in all source files.
John W. Eaton <jwe@octave.org>
parents:
24534
diff
changeset
|
3 Copyright (C) 1996-2018 John W. Eaton |
1967 | 4 |
5 This file is part of Octave. | |
6 | |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
24026
diff
changeset
|
7 Octave is free software: you can redistribute it and/or modify it |
22755
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
8 under the terms of the GNU General Public License as published by |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
24026
diff
changeset
|
9 the Free Software Foundation, either version 3 of the License, or |
22755
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
10 (at your option) any later version. |
1967 | 11 |
22755
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
12 Octave is distributed in the hope that it will be useful, but |
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
13 WITHOUT ANY WARRANTY; without even the implied warranty of |
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
3a2b891d0b33
maint: Standardize Copyright formatting.
Rik <rik@octave.org>
parents:
22402
diff
changeset
|
15 GNU General Public License for more details. |
1967 | 16 |
17 You should have received a copy of the GNU General Public License | |
7016 | 18 along with Octave; see the file COPYING. If not, see |
24534
194eb4bd202b
maint: Update punctuation for GPL v3 license text.
Rik <rik@octave.org>
parents:
24026
diff
changeset
|
19 <https://www.gnu.org/licenses/>. |
1967 | 20 |
21 */ | |
22 | |
21724
aba2e6293dd8
use "#if ..." consistently instead of "#ifdef" and "#ifndef"
John W. Eaton <jwe@octave.org>
parents:
21723
diff
changeset
|
23 #if defined (HAVE_CONFIG_H) |
21301
40de9f8f23a6
Use '#include "config.h"' rather than <config.h>.
Rik <rik@octave.org>
parents:
21202
diff
changeset
|
24 # include "config.h" |
1967 | 25 #endif |
26 | |
10463
bbe99b2a5ba7
undo recent gnulib-related changes
John W. Eaton <jwe@octave.org>
parents:
10447
diff
changeset
|
27 #include <cstdlib> |
7048 | 28 #include <cstring> |
1967 | 29 |
23443
3f1bf237908b
maint: Eliminate <cfloat.h> header from liboctave files.
Rik <rik@octave.org>
parents:
23220
diff
changeset
|
30 #include <complex> |
6490 | 31 #include <limits> |
2926 | 32 #include <string> |
33 | |
10068
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
34 #include "quit.h" |
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
35 |
2926 | 36 #include "lo-error.h" |
4130 | 37 #include "lo-ieee.h" |
1967 | 38 #include "lo-mappers.h" |
39 #include "lo-utils.h" | |
21888
00f1249f2483
create wrappers for some gnulib functions
John W. Eaton <jwe@octave.org>
parents:
21852
diff
changeset
|
40 #include "putenv-wrapper.h" |
1967 | 41 |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
42 bool xis_int_or_inf_or_nan (double x) |
21782
2aef506f3fec
use namespace for lo-mappers.h functions
John W. Eaton <jwe@octave.org>
parents:
21724
diff
changeset
|
43 { return octave::math::isnan (x) || octave::math::x_nint (x) == x; } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
44 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
45 bool xis_one_or_zero (double x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
46 { return x == 0 || x == 1; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
47 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
48 bool xis_zero (double x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
49 { return x == 0; } |
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 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
|
52 { |
23564
7049da1648c0
Deprecate octave::math::finite in favor of octave::math::isfinite.
Rik <rik@octave.org>
parents:
23448
diff
changeset
|
53 return (octave::math::isfinite (x) |
22402
4caa7b28d183
maint: Style check C++ code in liboctave/
Rik <rik@octave.org>
parents:
22323
diff
changeset
|
54 && 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
|
55 } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
56 |
15212
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
57 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
|
58 { |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
59 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
|
60 || 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
|
61 } |
4bbd3bbb8912
reduce code duplication in too_large_for_float array functions
John W. Eaton <jwe@octave.org>
parents:
14811
diff
changeset
|
62 |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
63 bool xis_int_or_inf_or_nan (float x) |
21782
2aef506f3fec
use namespace for lo-mappers.h functions
John W. Eaton <jwe@octave.org>
parents:
21724
diff
changeset
|
64 { return octave::math::isnan (x) || octave::math::x_nint (x) == x; } |
11010
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
65 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
66 bool xis_one_or_zero (float x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
67 { return x == 0 || x == 1; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
68 |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
69 bool xis_zero (float x) |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
70 { return x == 0; } |
9478b216752e
simplify more array tests
Jaroslav Hajek <highegg@gmail.com>
parents:
10463
diff
changeset
|
71 |
2926 | 72 // Save a string. |
73 | |
74 char * | |
75 strsave (const char *s) | |
76 { | |
77 if (! s) | |
23795
980f39c3ab90
Use C++11 nullptr rather than 0 in code (bug #51565).
Rik <rik@octave.org>
parents:
23564
diff
changeset
|
78 return nullptr; |
2926 | 79 |
80 int len = strlen (s); | |
81 char *tmp = new char [len+1]; | |
82 tmp = strcpy (tmp, s); | |
83 return tmp; | |
84 } | |
85 | |
17788
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
86 // This function was adapted from xputenv from Karl Berry's kpathsearch |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
87 // library. |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
88 |
17861
870f3e12e163
maint: Use phrase "FIXME:" for problem areas in code.
Rik <rik@octave.org>
parents:
17788
diff
changeset
|
89 // FIXME: make this do the right thing if we don't have a SMART_PUTENV. |
17788
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
90 |
2926 | 91 void |
3504 | 92 octave_putenv (const std::string& name, const std::string& value) |
2926 | 93 { |
17788
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
94 int new_len = name.length () + value.length () + 2; |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
95 |
21444
74b07936fd05
lo-utils.cc: Add note about inevitable memory leak from using putenv.
Rik <rik@octave.org>
parents:
21301
diff
changeset
|
96 // FIXME: This leaks memory, but so would a call to setenv. |
74b07936fd05
lo-utils.cc: Add note about inevitable memory leak from using putenv.
Rik <rik@octave.org>
parents:
21301
diff
changeset
|
97 // Short of extreme measures to track memory, altering the environment |
74b07936fd05
lo-utils.cc: Add note about inevitable memory leak from using putenv.
Rik <rik@octave.org>
parents:
21301
diff
changeset
|
98 // always leaks memory, but the saving grace is that the leaks are small. |
23447
db1fdf4384dd
maint: Use convention "static_cast<void *>" for casting of pointers.
Rik <rik@octave.org>
parents:
23443
diff
changeset
|
99 char *new_item = static_cast<char *> (std::malloc (new_len)); |
17788
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
100 |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
101 sprintf (new_item, "%s=%s", name.c_str (), value.c_str ()); |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
102 |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
103 // As far as I can see there's no way to distinguish between the |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
104 // various errors; putenv doesn't have errno values. |
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
105 |
21888
00f1249f2483
create wrappers for some gnulib functions
John W. Eaton <jwe@octave.org>
parents:
21852
diff
changeset
|
106 if (octave_putenv_wrapper (new_item) < 0) |
17788
415583856971
undo 7b305b36b87e to avoid Windows putenv/setenv failure (bug #40381)
John W. Eaton <jwe@octave.org>
parents:
17769
diff
changeset
|
107 (*current_liboctave_error_handler) ("putenv (%s) failed", new_item); |
2926 | 108 } |
109 | |
3504 | 110 std::string |
2926 | 111 octave_fgets (FILE *f) |
112 { | |
4527 | 113 bool eof; |
114 return octave_fgets (f, eof); | |
115 } | |
116 | |
117 std::string | |
118 octave_fgets (FILE *f, bool& eof) | |
119 { | |
120 eof = false; | |
121 | |
3504 | 122 std::string retval; |
2926 | 123 |
124 int grow_size = 1024; | |
125 int max_size = grow_size; | |
126 | |
21942
aab79a1885cc
limit gnulib headers to liboctave/wrappers directory
John W. Eaton <jwe@octave.org>
parents:
21888
diff
changeset
|
127 char *buf = static_cast<char *> (std::malloc (max_size)); |
2926 | 128 char *bufptr = buf; |
129 int len = 0; | |
130 | |
131 do | |
132 { | |
21942
aab79a1885cc
limit gnulib headers to liboctave/wrappers directory
John W. Eaton <jwe@octave.org>
parents:
21888
diff
changeset
|
133 if (std::fgets (bufptr, grow_size, f)) |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
134 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
135 len = strlen (bufptr); |
2926 | 136 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
137 if (len == grow_size - 1) |
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 int tmp = bufptr - buf + grow_size - 1; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
140 grow_size *= 2; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
141 max_size += grow_size; |
21942
aab79a1885cc
limit gnulib headers to liboctave/wrappers directory
John W. Eaton <jwe@octave.org>
parents:
21888
diff
changeset
|
142 buf = static_cast<char *> (std::realloc (buf, max_size)); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
143 bufptr = buf + tmp; |
2926 | 144 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
145 if (*(bufptr-1) == '\n') |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
146 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
147 *bufptr = '\0'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
148 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
149 } |
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 else if (bufptr[len-1] != '\n') |
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 bufptr[len++] = '\n'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
154 bufptr[len] = '\0'; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
155 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
156 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
157 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
158 retval = buf; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
159 } |
2926 | 160 else |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
161 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
162 if (len == 0) |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
163 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
164 eof = true; |
4527 | 165 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
166 free (buf); |
2926 | 167 |
23795
980f39c3ab90
Use C++11 nullptr rather than 0 in code (bug #51565).
Rik <rik@octave.org>
parents:
23564
diff
changeset
|
168 buf = nullptr; |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
169 } |
2926 | 170 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
171 break; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
172 } |
2926 | 173 } |
174 while (retval.empty ()); | |
175 | |
23878
889b8262e661
maint: eliminate useless if (p) before free (p)
Mike Miller <mtmiller@octave.org>
parents:
23826
diff
changeset
|
176 free (buf); |
2926 | 177 |
10142
829e69ec3110
make OCTAVE_QUIT a function
Jaroslav Hajek <highegg@gmail.com>
parents:
10068
diff
changeset
|
178 octave_quit (); |
10068
ca93f583573d
handle interrupts octave_fgets
Jaroslav Hajek <highegg@gmail.com>
parents:
9469
diff
changeset
|
179 |
2926 | 180 return retval; |
181 } | |
182 | |
3970 | 183 std::string |
184 octave_fgetl (FILE *f) | |
185 { | |
4527 | 186 bool eof; |
187 return octave_fgetl (f, eof); | |
188 } | |
189 | |
190 std::string | |
191 octave_fgetl (FILE *f, bool& eof) | |
192 { | |
193 std::string retval = octave_fgets (f, eof); | |
3970 | 194 |
23812
057a894914df
Use C++11 string fcns back() and pop_back() to simplify code.
Rik <rik@octave.org>
parents:
23807
diff
changeset
|
195 if (! retval.empty () && retval.back () == '\n') |
057a894914df
Use C++11 string fcns back() and pop_back() to simplify code.
Rik <rik@octave.org>
parents:
23807
diff
changeset
|
196 retval.pop_back (); |
3970 | 197 |
198 return retval; | |
199 } | |
200 | |
17928
079da800451b
Fix loading complex values with 'NA' saved in text format (bug #40575).
Rik <rik@octave.org>
parents:
17861
diff
changeset
|
201 // Note that the caller is responsible for repositioning the stream on failure. |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
202 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
203 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
|
204 T |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
205 read_inf_nan_na (std::istream& is, char c0) |
4130 | 206 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
207 T val = 0.0; |
4130 | 208 |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
209 switch (c0) |
4130 | 210 { |
211 case 'i': case 'I': | |
212 { | |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
213 char c1 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
214 if (c1 == 'n' || c1 == 'N') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
215 { |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
216 char c2 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
217 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
|
218 val = std::numeric_limits<T>::infinity (); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
219 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
220 is.setstate (std::ios::failbit); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
221 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
222 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
223 is.setstate (std::ios::failbit); |
4130 | 224 } |
225 break; | |
226 | |
227 case 'n': case 'N': | |
228 { | |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
229 char c1 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
230 if (c1 == 'a' || c1 == 'A') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
231 { |
12995
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
232 char c2 = is.get (); |
7872afb42fab
fix scanf problem with reading I (bug #33722)
John W. Eaton <jwe@octave.org>
parents:
12936
diff
changeset
|
233 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
|
234 val = std::numeric_limits<T>::quiet_NaN (); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
235 else |
17928
079da800451b
Fix loading complex values with 'NA' saved in text format (bug #40575).
Rik <rik@octave.org>
parents:
17861
diff
changeset
|
236 { |
21723
bae585228161
use namespace for numeric_limits
John W. Eaton <jwe@octave.org>
parents:
21444
diff
changeset
|
237 val = octave::numeric_limits<T>::NA (); |
21184
3db899f1d54c
use istream::traits_type::eof () instead of EOF
John W. Eaton <jwe@octave.org>
parents:
20148
diff
changeset
|
238 if (c2 != std::istream::traits_type::eof ()) |
20148
4a41b7d426b6
Fix loading of ASCII files with NA values in final column (bug #44967).
Rik <rik@octave.org>
parents:
19697
diff
changeset
|
239 is.putback (c2); |
17928
079da800451b
Fix loading complex values with 'NA' saved in text format (bug #40575).
Rik <rik@octave.org>
parents:
17861
diff
changeset
|
240 } |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
241 } |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
242 else |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
243 is.setstate (std::ios::failbit); |
4130 | 244 } |
245 break; | |
246 | |
247 default: | |
23448
e1c02e7126a9
avoid calling abort if possible; include cstdlib for declaration if needed
John W. Eaton <jwe@octave.org>
parents:
23447
diff
changeset
|
248 (*current_liboctave_error_handler) ("read_inf_nan_na: invalid character '%c'"); |
4130 | 249 } |
250 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
251 return val; |
4130 | 252 } |
253 | |
12936
b74cb659e757
accept but discard sign when reading NA and NaN values
John W. Eaton <jwe@octave.org>
parents:
12918
diff
changeset
|
254 // 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
|
255 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
256 template <typename T> |
4130 | 257 double |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
258 octave_read_fp_value (std::istream& is) |
4130 | 259 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
260 T val = 0.0; |
4130 | 261 |
17769
49a5a4be04a1
maint: Use GNU style coding conventions for code in liboctave/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
262 // FIXME: resetting stream position is likely to fail unless we are |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
263 // reading from a file. |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
264 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
|
265 |
6907 | 266 char c1 = ' '; |
4130 | 267 |
6907 | 268 while (isspace (c1)) |
269 c1 = is.get (); | |
270 | |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
271 bool neg = false; |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
272 |
6194 | 273 switch (c1) |
4130 | 274 { |
6194 | 275 case '-': |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
276 neg = true; |
23826
d69021d58a61
avoid fallthrough warnings
John W. Eaton <jwe@octave.org>
parents:
23812
diff
changeset
|
277 OCTAVE_FALLTHROUGH; |
6194 | 278 |
279 case '+': | |
280 { | |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
281 char c2 = 0; |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
282 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
|
283 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
|
284 val = read_inf_nan_na<T> (is, c2); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
285 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
286 { |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
287 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
|
288 is >> val; |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
289 } |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
290 |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
291 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
|
292 val = -val; |
6194 | 293 } |
294 break; | |
295 | |
4130 | 296 case 'i': case 'I': |
297 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
|
298 val = read_inf_nan_na<T> (is, c1); |
4130 | 299 break; |
300 | |
301 default: | |
6194 | 302 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
|
303 is >> val; |
14806
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
304 break; |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
305 } |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
306 |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
307 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
|
308 if (status & std::ios::failbit) |
980e2d5c83f7
avoid calling putback more than once between reads
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
309 { |
24026
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
310 // Convert MAX_VAL returned by C++ streams for very large numbers to Inf |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
311 if (val == std::numeric_limits<T>::max ()) |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
312 { |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
313 if (neg) |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
314 val = -std::numeric_limits<T>::infinity (); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
315 else |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
316 val = std::numeric_limits<T>::infinity (); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
317 is.clear (status & ~std::ios::failbit); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
318 } |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
319 else |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
320 { |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
321 // True error. Reset stream to original position and pass status on. |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
322 is.clear (); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
323 is.seekg (pos); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
324 is.setstate (status); |
1cb94f46466f
Return Inf/-Inf when reading large numbers from files for Matlab compatibility.
Rik <rik@octave.org>
parents:
23878
diff
changeset
|
325 } |
4130 | 326 } |
327 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
328 return val; |
4130 | 329 } |
330 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
331 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
|
332 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
|
333 octave_read_cx_fp_value (std::istream& is) |
4130 | 334 { |
18084
8e056300994b
Follow coding convention of defining and initializing only 1 variable per line in liboctave.
Rik <rik@octave.org>
parents:
17928
diff
changeset
|
335 T re = 0.0; |
8e056300994b
Follow coding convention of defining and initializing only 1 variable per line in liboctave.
Rik <rik@octave.org>
parents:
17928
diff
changeset
|
336 T im = 0.0; |
4130 | 337 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
338 std::complex<T> cx = 0.0; |
4130 | 339 |
6907 | 340 char ch = ' '; |
4130 | 341 |
6907 | 342 while (isspace (ch)) |
343 ch = is.get (); | |
4130 | 344 |
345 if (ch == '(') | |
346 { | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
347 re = octave_read_value<T> (is); |
6897 | 348 ch = is.get (); |
4130 | 349 |
350 if (ch == ',') | |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
351 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
352 im = octave_read_value<T> (is); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
353 ch = is.get (); |
4130 | 354 |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
355 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
|
356 cx = std::complex<T> (re, im); |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
357 else |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
358 is.setstate (std::ios::failbit); |
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
359 } |
4130 | 360 else if (ch == ')') |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
361 cx = re; |
4130 | 362 else |
10314
07ebe522dac2
untabify liboctave C++ sources
John W. Eaton <jwe@octave.org>
parents:
10182
diff
changeset
|
363 is.setstate (std::ios::failbit); |
4130 | 364 } |
365 else | |
366 { | |
367 is.putback (ch); | |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
368 cx = octave_read_value<double> (is); |
4130 | 369 } |
370 | |
371 return cx; | |
372 } | |
373 | |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
374 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
|
375 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
376 return octave_read_fp_value<double> (is); |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
377 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
378 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
379 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
|
380 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
381 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
|
382 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
383 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
384 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
|
385 { |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
386 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
|
387 } |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
388 |
14811
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
389 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
|
390 { |
52cb71787cd1
use templates to avoid code duplication in octave_read_value functions
John W. Eaton <jwe@octave.org>
parents:
14806
diff
changeset
|
391 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
|
392 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
393 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
394 void |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
395 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
|
396 { |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
397 if (lo_ieee_is_NA (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
398 os << "NA"; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
399 else if (lo_ieee_isnan (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
400 os << "NaN"; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
401 else if (lo_ieee_isinf (d)) |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
402 os << (d < 0 ? "-Inf" : "Inf"); |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
403 else |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
404 os << d; |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
405 } |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
406 |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
407 void |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
408 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
|
409 { |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
410 os << '('; |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
411 octave_write_double (os, real (c)); |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
412 os << ','; |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
413 octave_write_double (os, imag (c)); |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
414 os << ')'; |
9469
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
415 } |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
416 |
c6edba80dfae
sanity checks for loading sparse matrices
John W. Eaton <jwe@octave.org>
parents:
8920
diff
changeset
|
417 void |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
418 octave_write_float (std::ostream& os, float d) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
419 { |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
420 if (lo_ieee_is_NA (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
421 os << "NA"; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
422 else if (lo_ieee_isnan (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
423 os << "NaN"; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
424 else if (lo_ieee_isinf (d)) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
425 os << (d < 0 ? "-Inf" : "Inf"); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
426 else |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
427 os << d; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
428 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
429 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
430 void |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
431 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
|
432 { |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
433 os << '('; |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
434 octave_write_float (os, real (c)); |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
435 os << ','; |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
436 octave_write_float (os, imag (c)); |
23807
336f89b6208b
Use character literals 'c' rather than string literals "c" when possible.
Rik <rik@octave.org>
parents:
23795
diff
changeset
|
437 os << ')'; |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7048
diff
changeset
|
438 } |