octave: liboctave/util/lo-utils.cc annotate

annotate liboctave/util/lo-utils.cc @ 27919:1891570abac8

update Octave Project Developers copyright for the new year In files that have the "Octave Project Developers" copyright notice, update for 2020.

author	John W. Eaton <jwe@octave.org>
date	Mon, 06 Jan 2020 22:29:51 -0500
parents	b442ec6dda5c
children	bd51beb6205e

rev	line source
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	1 /*
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	2
27919 1891570abac8 update Octave Project Developers copyright for the new year John W. Eaton <jwe@octave.org> parents: 27918 diff changeset	3 Copyright (C) 1996-2020 The Octave Project Developers
27918 b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27021 diff changeset	4
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27021 diff changeset	5 See the file COPYRIGHT.md in the top-level directory of this distribution
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27021 diff changeset	6 or <https://octave.org/COPYRIGHT.html/>.
b442ec6dda5c use centralized file for copyright info for individual contributors John W. Eaton <jwe@octave.org> parents: 27021 diff changeset	7
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	8
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	9 This file is part of Octave.
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	10
24534 194eb4bd202b maint: Update punctuation for GPL v3 license text. Rik <rik@octave.org> parents: 24026 diff changeset	11 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	12 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	13 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	14 (at your option) any later version.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	15
22755 3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	16 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	17 WITHOUT ANY WARRANTY; without even the implied warranty of
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	18 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
3a2b891d0b33 maint: Standardize Copyright formatting. Rik <rik@octave.org> parents: 22402 diff changeset	19 GNU General Public License for more details.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	20
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	21 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: 6907 diff changeset	22 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	23 <https://www.gnu.org/licenses/>.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	24
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	25 */
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	26
21724 aba2e6293dd8 use "#if ..." consistently instead of "#ifdef" and "#ifndef" John W. Eaton <jwe@octave.org> parents: 21723 diff changeset	27 #if defined (HAVE_CONFIG_H)
21301 40de9f8f23a6 Use '#include "config.h"' rather than <config.h>. Rik <rik@octave.org> parents: 21202 diff changeset	28 # include "config.h"
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	29 #endif
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	30
10463 bbe99b2a5ba7 undo recent gnulib-related changes John W. Eaton <jwe@octave.org> parents: 10447 diff changeset	31 #include <cstdlib>
7048 845ca0affec0 [project @ 2007-10-22 16:55:41 by jwe] jwe parents: 7017 diff changeset	32 #include <cstring>
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	33
23443 3f1bf237908b maint: Eliminate <cfloat.h> header from liboctave files. Rik <rik@octave.org> parents: 23220 diff changeset	34 #include <complex>
25438 cb1606f78f6b prefer <istream>, <ostream>, or <iosfwd> to <iostream> where possible John W. Eaton <jwe@octave.org> parents: 25054 diff changeset	35 #include <istream>
6490 0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	36 #include <limits>
25438 cb1606f78f6b prefer <istream>, <ostream>, or <iosfwd> to <iostream> where possible John W. Eaton <jwe@octave.org> parents: 25054 diff changeset	37 #include <ostream>
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	38 #include <string>
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	39
10068 ca93f583573d handle interrupts octave_fgets Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	40 #include "quit.h"
ca93f583573d handle interrupts octave_fgets Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	41
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	42 #include "lo-error.h"
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	43 #include "lo-ieee.h"
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	44 #include "lo-mappers.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	45 #include "lo-utils.h"
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	46
11010 9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	47 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	48 { 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	49
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	50 bool xis_one_or_zero (double x)
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	51 { return x == 0 \|\| x == 1; }
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	52
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	53 bool xis_zero (double x)
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	54 { return x == 0; }
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	55
9478b216752e simplify more array tests Jaroslav Hajek <highegg@gmail.com> parents: 10463 diff changeset	56 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	57 {
23564 7049da1648c0 Deprecate octave::math::finite in favor of octave::math::isfinite. Rik <rik@octave.org> parents: 23448 diff changeset	58 return (octave::math::isfinite (x)
22402 4caa7b28d183 maint: Style check C++ code in liboctave/ Rik <rik@octave.org> parents: 22323 diff changeset	59 && 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)
21782 2aef506f3fec use namespace for lo-mappers.h functions John W. Eaton <jwe@octave.org> parents: 21724 diff changeset	69 { 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	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 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	77 // Save a string.
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	78
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	79 char *
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	80 strsave (const char *s)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	81 {
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	82 if (! s)
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23564 diff changeset	83 return nullptr;
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	84
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	85 int len = strlen (s);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	86 char *tmp = new char [len+1];
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	87 tmp = strcpy (tmp, s);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	88 return tmp;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	89 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	90
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2926 diff changeset	91 std::string
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	92 octave_fgets (FILE *f)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	93 {
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	94 bool eof;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	95 return octave_fgets (f, eof);
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	96 }
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	97
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	98 std::string
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	99 octave_fgets (FILE *f, bool& eof)
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	100 {
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	101 eof = false;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	102
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2926 diff changeset	103 std::string retval;
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	104
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	105 int grow_size = 1024;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	106 int max_size = grow_size;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	107
21942 aab79a1885cc limit gnulib headers to liboctave/wrappers directory John W. Eaton <jwe@octave.org> parents: 21888 diff changeset	108 char buf = static_cast<char > (std::malloc (max_size));
26439 4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	109 if (! buf)
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	110 (*current_liboctave_error_handler) ("octave_fgets: unable to malloc %d bytes", max_size);
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	111
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	112 char *bufptr = buf;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	113 int len = 0;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	114
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	115 do
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	116 {
21942 aab79a1885cc limit gnulib headers to liboctave/wrappers directory John W. Eaton <jwe@octave.org> parents: 21888 diff changeset	117 if (std::fgets (bufptr, grow_size, f))
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	118 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	119 len = strlen (bufptr);
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	120
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	121 if (len == grow_size - 1)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	122 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	123 int tmp = bufptr - buf + grow_size - 1;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	124 grow_size *= 2;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	125 max_size += grow_size;
26439 4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	126 auto tmpbuf = static_cast<char *> (std::realloc (buf, max_size));
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	127 if (! tmpbuf)
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	128 {
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	129 free (buf);
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	130 (*current_liboctave_error_handler) ("octave_fgets: unable to realloc %d bytes", max_size);
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	131 }
4d5bf84ae249 lo-utils.cc: Fix static analyzer detected issues (bug #55347). Rik <rik@octave.org> parents: 26376 diff changeset	132 buf = tmpbuf;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	133 bufptr = buf + tmp;
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	134
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	135 if (*(bufptr-1) == '\n')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	136 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	137 *bufptr = '\0';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	138 retval = buf;
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 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	141 else if (bufptr[len-1] != '\n')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	142 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	143 bufptr[len++] = '\n';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	144 bufptr[len] = '\0';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	145 retval = buf;
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 else
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 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	150 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	151 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	152 if (len == 0)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	153 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	154 eof = true;
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	155
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	156 free (buf);
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	157
23795 980f39c3ab90 Use C++11 nullptr rather than 0 in code (bug #51565). Rik <rik@octave.org> parents: 23564 diff changeset	158 buf = nullptr;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	159 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	160
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	161 break;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	162 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	163 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	164 while (retval.empty ());
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	165
23878 889b8262e661 maint: eliminate useless if (p) before free (p) Mike Miller <mtmiller@octave.org> parents: 23826 diff changeset	166 free (buf);
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	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 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	170 return retval;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	171 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	172
3970 4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	173 std::string
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	174 octave_fgetl (FILE *f)
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	175 {
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	176 bool eof;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	177 return octave_fgetl (f, eof);
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	178 }
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	179
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	180 std::string
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	181 octave_fgetl (FILE *f, bool& eof)
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	182 {
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	183 std::string retval = octave_fgets (f, eof);
3970 4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	184
23812 057a894914df Use C++11 string fcns back() and pop_back() to simplify code. Rik <rik@octave.org> parents: 23807 diff changeset	185 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	186 retval.pop_back ();
3970 4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	187
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	188 return retval;
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	189 }
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	190
17928 079da800451b Fix loading complex values with 'NA' saved in text format (bug #40575). Rik <rik@octave.org> parents: 17861 diff changeset	191 // 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	192
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	193 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	194 T
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	195 read_inf_nan_na (std::istream& is, char c0)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	196 {
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	197 T val = 0.0;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	198
12995 7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	199 switch (c0)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	200 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	201 case 'i': case 'I':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	202 {
12995 7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	203 char c1 = is.get ();
7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	204 if (c1 == 'n' \|\| c1 == 'N')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	205 {
12995 7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	206 char c2 = is.get ();
7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	207 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	208 val = std::numeric_limits<T>::infinity ();
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	209 else
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	210 is.setstate (std::ios::failbit);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	211 }
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);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	214 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	215 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	216
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	217 case 'n': case 'N':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	218 {
12995 7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	219 char c1 = is.get ();
7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	220 if (c1 == 'a' \|\| c1 == 'A')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	221 {
12995 7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	222 char c2 = is.get ();
7872afb42fab fix scanf problem with reading I (bug #33722) John W. Eaton <jwe@octave.org> parents: 12936 diff changeset	223 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	224 val = std::numeric_limits<T>::quiet_NaN ();
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	225 else
17928 079da800451b Fix loading complex values with 'NA' saved in text format (bug #40575). Rik <rik@octave.org> parents: 17861 diff changeset	226 {
21723 bae585228161 use namespace for numeric_limits John W. Eaton <jwe@octave.org> parents: 21444 diff changeset	227 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	228 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	229 is.putback (c2);
25672 31c0bb386bcd Correctly read NA values in dlmread. Rik <rik@octave.org> parents: 25438 diff changeset	230 else
31c0bb386bcd Correctly read NA values in dlmread. Rik <rik@octave.org> parents: 25438 diff changeset	231 is.clear (is.rdstate () & ~std::ios::failbit);
17928 079da800451b Fix loading complex values with 'NA' saved in text format (bug #40575). Rik <rik@octave.org> parents: 17861 diff changeset	232 }
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	233 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	234 else
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	235 is.setstate (std::ios::failbit);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	236 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	237 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	238
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	239 default:
26569 0e77df67b522 Add static compile-time checking of printf functions in liboctave. Markus Mützel <markus.muetzel@gmx.de> parents: 26439 diff changeset	240 (*current_liboctave_error_handler)
0e77df67b522 Add static compile-time checking of printf functions in liboctave. Markus Mützel <markus.muetzel@gmx.de> parents: 26439 diff changeset	241 ("read_inf_nan_na: invalid character '%c'", c0);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	242 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	243
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	244 return val;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	245 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	246
12936 b74cb659e757 accept but discard sign when reading NA and NaN values John W. Eaton <jwe@octave.org> parents: 12918 diff changeset	247 // 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	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 template <typename T>
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	250 double
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	251 octave_read_fp_value (std::istream& is)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	252 {
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	253 T val = 0.0;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	254
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	255 // 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	256 // reading from a file.
27021 c61301840704 use std::streampos instead of std::ios::streampos (patch #9782) John W. Eaton <jwe@octave.org> parents: 26569 diff changeset	257 std::streampos pos = is.tellg ();
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	258
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	259 char c1 = ' ';
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	260
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	261 while (isspace (c1))
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	262 c1 = is.get ();
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	263
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	264 bool neg = false;
980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	265
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	266 switch (c1)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	267 {
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	268 case '-':
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	269 neg = true;
23826 d69021d58a61 avoid fallthrough warnings John W. Eaton <jwe@octave.org> parents: 23812 diff changeset	270 OCTAVE_FALLTHROUGH;
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	271
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	272 case '+':
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	273 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	274 char c2 = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	275 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	276 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	277 val = read_inf_nan_na<T> (is, c2);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	278 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	279 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	280 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	281 is >> val;
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	282 }
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	283
980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	284 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	285 val = -val;
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	286 }
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	287 break;
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	288
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	289 case 'i': case 'I':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	290 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	291 val = read_inf_nan_na<T> (is, c1);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	292 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	293
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	294 default:
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	295 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	296 is >> val;
14806 980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	297 break;
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
980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	300 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	301 if (status & std::ios::failbit)
980e2d5c83f7 avoid calling putback more than once between reads John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	302 {
24026 1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	303 // 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	304 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	305 {
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	306 if (neg)
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	307 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	308 else
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	309 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	310 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	311 }
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	312 else
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	313 {
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	314 // 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	315 is.clear ();
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	316 is.seekg (pos);
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	317 is.setstate (status);
1cb94f46466f Return Inf/-Inf when reading large numbers from files for Matlab compatibility. Rik <rik@octave.org> parents: 23878 diff changeset	318 }
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	319 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	320
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	321 return val;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	322 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	323
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	324 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	325 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	326 octave_read_cx_fp_value (std::istream& is)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	327 {
18084 8e056300994b Follow coding convention of defining and initializing only 1 variable per line in liboctave. Rik <rik@octave.org> parents: 17928 diff changeset	328 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	329 T im = 0.0;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	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 std::complex<T> cx = 0.0;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	332
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	333 char ch = ' ';
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	334
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	335 while (isspace (ch))
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	336 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	337
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	338 if (ch == '(')
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	339 {
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	340 re = octave_read_value<T> (is);
6897 935d23e16951 [project @ 2007-09-13 19:13:56 by jwe] jwe parents: 6490 diff changeset	341 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	342
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	343 if (ch == ',')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	344 {
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	345 im = octave_read_value<T> (is);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	346 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	347
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	348 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	349 cx = std::complex<T> (re, im);
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	350 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	351 is.setstate (std::ios::failbit);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	352 }
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	353 else if (ch == ')')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	354 cx = re;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	355 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	356 is.setstate (std::ios::failbit);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	357 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	358 else
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	359 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	360 is.putback (ch);
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	361 cx = octave_read_value<double> (is);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	362 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	363
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	364 return cx;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	365 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	366
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	367 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	368 {
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	369 return octave_read_fp_value<double> (is);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	370 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	371
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	372 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	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 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	375 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	376
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	377 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	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 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	380 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	381
14811 52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	382 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	383 {
52cb71787cd1 use templates to avoid code duplication in octave_read_value functions John W. Eaton <jwe@octave.org> parents: 14806 diff changeset	384 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	385 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	386
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	387 void
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	388 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	389 {
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	390 if (lo_ieee_is_NA (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	391 os << "NA";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	392 else if (lo_ieee_isnan (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	393 os << "NaN";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	394 else if (lo_ieee_isinf (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	395 os << (d < 0 ? "-Inf" : "Inf");
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	396 else
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	397 os << d;
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	398 }
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	399
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	400 void
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	401 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	402 {
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23795 diff changeset	403 os << '(';
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	404 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	405 os << ',';
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	406 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	407 os << ')';
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	408 }
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	409
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	410 void
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	411 octave_write_float (std::ostream& os, float d)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	412 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	413 if (lo_ieee_is_NA (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	414 os << "NA";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	415 else if (lo_ieee_isnan (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	416 os << "NaN";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	417 else if (lo_ieee_isinf (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	418 os << (d < 0 ? "-Inf" : "Inf");
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	419 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	420 os << d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	421 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	422
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	423 void
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	424 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	425 {
23807 336f89b6208b Use character literals 'c' rather than string literals "c" when possible. Rik <rik@octave.org> parents: 23795 diff changeset	426 os << '(';
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	427 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	428 os << ',';
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	429 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	430 os << ')';
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	431 }

Mercurial > octave

annotate liboctave/util/lo-utils.cc @ 27919:1891570abac8