octave: liboctave/lo-utils.cc annotate

annotate liboctave/lo-utils.cc @ 10447:12d25a1d84bf

Use gnulib headers instead of C++ std headers where needed.

author	David Grundberg <davidg@cs.umu.se>
date	Wed, 24 Mar 2010 17:27:58 +0100
parents	479cc8a0a846
children	bbe99b2a5ba7

rev	line source
1993 1b57120c997b [project @ 1996-03-03 01:16:15 by jwe] jwe parents: 1967 diff changeset	1 // utils.cc
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	2 /*
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	3
8920 eb63fbe60fab update copyright notices John W. Eaton <jwe@octave.org> parents: 7789 diff changeset	4 Copyright (C) 1996, 1997, 2000, 2002, 2003, 2004, 2005, 2006, 2007, 2008
7017 a1dbe9d80eee [project @ 2007-10-12 21:27:11 by jwe] jwe parents: 7016 diff changeset	5 John W. Eaton
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	6
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	7 This file is part of Octave.
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 Octave is free software; you can redistribute it and/or modify it
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	10 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6907 diff changeset	11 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6907 diff changeset	12 option) any later version.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	13
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	14 Octave is distributed in the hope that it will be useful, but WITHOUT
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	15 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	16 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	17 for more details.
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	18
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	19 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	20 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6907 diff changeset	21 <http://www.gnu.org/licenses/>.
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	22
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	23 */
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 #ifdef HAVE_CONFIG_H
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	26 #include <config.h>
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	27 #endif
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	28
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	29 #include <cctype>
10447 12d25a1d84bf Use gnulib headers instead of C++ std headers where needed. David Grundberg <davidg@cs.umu.se> parents: 10411 diff changeset	30 #include "stdlib.h"
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	31 #include <cstdio>
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
6490 0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	34 #include <limits>
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	35 #include <string>
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	36
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	37 #include <sys/types.h>
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	38 #include <unistd.h>
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
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	47 // Convert X to the nearest integer value. Should not pass NaN to
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	48 // this function.
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	49
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	50 // Sometimes you need a large integer, but not always.
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	51
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	52 octave_idx_type
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	53 NINTbig (double x)
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	54 {
6490 0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	55 if (x > std::numeric_limits<octave_idx_type>::max ())
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	56 return std::numeric_limits<octave_idx_type>::max ();
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	57 else if (x < std::numeric_limits<octave_idx_type>::min ())
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	58 return std::numeric_limits<octave_idx_type>::min ();
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	59 else
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	60 return static_cast<octave_idx_type> ((x > 0) ? (x + 0.5) : (x - 0.5));
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	61 }
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 4784 diff changeset	62
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	63 octave_idx_type
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	64 NINTbig (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	65 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	66 if (x > std::numeric_limits<octave_idx_type>::max ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	67 return std::numeric_limits<octave_idx_type>::max ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	68 else if (x < std::numeric_limits<octave_idx_type>::min ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	69 return std::numeric_limits<octave_idx_type>::min ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	70 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	71 return static_cast<octave_idx_type> ((x > 0) ? (x + 0.5) : (x - 0.5));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	72 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	73
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	74 int
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	75 NINT (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	76 {
6490 0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	77 if (x > std::numeric_limits<int>::max ())
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	78 return std::numeric_limits<int>::max ();
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	79 else if (x < std::numeric_limits<int>::min ())
0ad7655cf2bc [project @ 2007-04-05 02:44:34 by jwe] jwe parents: 6202 diff changeset	80 return std::numeric_limits<int>::min ();
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	81 else
4784 743ef6154f8a [project @ 2004-02-17 19:59:23 by jwe] jwe parents: 4527 diff changeset	82 return static_cast<int> ((x > 0) ? (x + 0.5) : (x - 0.5));
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	83 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	84
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	85 int
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	86 NINT (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	87 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	88 if (x > std::numeric_limits<int>::max ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	89 return std::numeric_limits<int>::max ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	90 else if (x < std::numeric_limits<int>::min ())
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	91 return std::numeric_limits<int>::min ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	92 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	93 return static_cast<int> ((x > 0) ? (x + 0.5) : (x - 0.5));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	94 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	95
1967 4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	96 double
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	97 D_NINT (double x)
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	98 {
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	99 if (xisinf (x) \|\| xisnan (x))
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	100 return x;
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	101 else
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	102 return floor (x + 0.5);
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	103 }
4cb6aaa0f2c7 [project @ 1996-02-17 02:54:24 by jwe] jwe parents: diff changeset	104
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	105 float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	106 F_NINT (float x)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	107 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	108 if (xisinf (x) \|\| xisnan (x))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	109 return x;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	110 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	111 return floor (x + 0.5);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	112 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	113
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	114 // Save a string.
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	115
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	116 char *
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	117 strsave (const char *s)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	118 {
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	119 if (! s)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	120 return 0;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	121
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	122 int len = strlen (s);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	123 char *tmp = new char [len+1];
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	124 tmp = strcpy (tmp, s);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	125 return tmp;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	126 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	127
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	128 // This function was adapted from xputenv from Karl Berry's kpathsearch
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	129 // library.
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	130
5775 ace8d8d26933 [project @ 2006-04-24 19:13:06 by jwe] jwe parents: 5307 diff changeset	131 // FIXME -- make this do the right thing if we don't have a
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	132 // SMART_PUTENV.
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	133
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	134 void
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2926 diff changeset	135 octave_putenv (const std::string& name, const std::string& value)
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	136 {
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	137 int new_len = name.length () + value.length () + 2;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	138
10411 479cc8a0a846 use gnulib namespace John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	139 char new_item = static_cast<char> (gnulib::malloc (new_len));
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	140
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	141 sprintf (new_item, "%s=%s", name.c_str (), value.c_str ());
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	142
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	143 // As far as I can see there's no way to distinguish between the
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	144 // various errors; putenv doesn't have errno values.
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	145
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	146 if (putenv (new_item) < 0)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	147 (*current_liboctave_error_handler) ("putenv (%s) failed", new_item);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	148 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	149
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2926 diff changeset	150 std::string
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	151 octave_fgets (FILE *f)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	152 {
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	153 bool eof;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	154 return octave_fgets (f, eof);
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	155 }
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	156
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	157 std::string
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	158 octave_fgets (FILE *f, bool& eof)
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	159 {
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	160 eof = false;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	161
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2926 diff changeset	162 std::string retval;
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 int grow_size = 1024;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	165 int max_size = grow_size;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	166
10411 479cc8a0a846 use gnulib namespace John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	167 char buf = static_cast<char > (gnulib::malloc (max_size));
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	168 char *bufptr = buf;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	169 int len = 0;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	170
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	171 do
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	172 {
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	173 if (fgets (bufptr, grow_size, f))
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	174 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	175 len = strlen (bufptr);
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	176
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	177 if (len == grow_size - 1)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	178 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	179 int tmp = bufptr - buf + grow_size - 1;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	180 grow_size *= 2;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	181 max_size += grow_size;
10411 479cc8a0a846 use gnulib namespace John W. Eaton <jwe@octave.org> parents: 10314 diff changeset	182 buf = static_cast<char *> (gnulib::realloc (buf, max_size));
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	183 bufptr = buf + tmp;
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	184
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	185 if (*(bufptr-1) == '\n')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	186 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	187 *bufptr = '\0';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	188 retval = buf;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	189 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	190 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	191 else if (bufptr[len-1] != '\n')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	192 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	193 bufptr[len++] = '\n';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	194 bufptr[len] = '\0';
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	195 retval = buf;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	196 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	197 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	198 retval = buf;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	199 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	200 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	201 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	202 if (len == 0)
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	203 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	204 eof = true;
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	205
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	206 free (buf);
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	207
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	208 buf = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	209 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	210
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	211 break;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	212 }
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	213 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	214 while (retval.empty ());
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	215
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	216 if (buf)
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	217 free (buf);
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	218
10142 829e69ec3110 make OCTAVE_QUIT a function Jaroslav Hajek <highegg@gmail.com> parents: 10068 diff changeset	219 octave_quit ();
10068 ca93f583573d handle interrupts octave_fgets Jaroslav Hajek <highegg@gmail.com> parents: 9469 diff changeset	220
2926 66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	221 return retval;
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	222 }
66ef74ee5d9f [project @ 1997-05-05 03:20:52 by jwe] jwe parents: 2847 diff changeset	223
3970 4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	224 std::string
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	225 octave_fgetl (FILE *f)
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	226 {
4527 c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	227 bool eof;
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	228 return octave_fgetl (f, eof);
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	229 }
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	230
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	231 std::string
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	232 octave_fgetl (FILE *f, bool& eof)
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	233 {
c0a23a13eea2 [project @ 2003-10-03 02:52:46 by jwe] jwe parents: 4130 diff changeset	234 std::string retval = octave_fgets (f, eof);
3970 4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	235
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	236 size_t len = retval.length ();
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	237
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	238 if (retval[len-1] == '\n')
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	239 retval.resize (len-1);
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	240
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	241 return retval;
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	242 }
4f884e25aab9 [project @ 2002-07-02 04:14:35 by jwe] jwe parents: 3504 diff changeset	243
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	244 static inline double
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	245 read_inf_nan_na (std::istream& is, char c, char sign = '+')
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	246 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	247 double d = 0.0;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	248
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	249 switch (c)
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	250 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	251 case 'i': case 'I':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	252 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	253 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	254 if (c == 'n' \|\| c == 'N')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	255 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	256 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	257 if (c == 'f' \|\| c == 'F')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	258 d = sign == '-' ? -octave_Inf : octave_Inf;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	259 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	260 is.putback (c);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	261 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	262 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	263 is.putback (c);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	264 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	265 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	266
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	267 case 'n': case 'N':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	268 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	269 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	270 if (c == 'a' \|\| c == 'A')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	271 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	272 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	273 if (c == 'n' \|\| c == 'N')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	274 d = octave_NaN;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	275 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	276 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	277 is.putback (c);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	278 d = octave_NA;
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 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	281 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	282 is.putback (c);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	283 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	284 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	285
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	286 default:
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	287 abort ();
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	288 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	289
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	290 return d;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	291 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	292
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	293 template <>
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	294 double
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	295 octave_read_value (std::istream& is)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	296 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	297 double d = 0.0;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	298
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	299 char c1 = ' ';
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	300
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	301 while (isspace (c1))
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	302 c1 = is.get ();
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	303
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	304 switch (c1)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	305 {
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	306 case '-':
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	307 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	308 char c2 = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	309 c2 = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	310 if (c2 == 'i' \|\| c2 == 'I')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	311 d = read_inf_nan_na (is, c2, c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	312 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	313 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	314 is.putback (c2);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	315 is.putback (c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	316 is >> d;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	317 }
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	318 }
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	319 break;
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	320
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	321 case '+':
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	322 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	323 char c2 = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	324 c2 = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	325 if (c2 == 'i' \|\| c2 == 'I')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	326 d = read_inf_nan_na (is, c2, c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	327 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	328 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	329 is.putback (c2);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	330 is.putback (c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	331 is >> d;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	332 }
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	333 }
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	334 break;
0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	335
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	336 case 'i': case 'I':
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	337 case 'n': case 'N':
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	338 d = read_inf_nan_na (is, c1);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	339 break;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	340
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	341 default:
6194 0d2ff0dfb159 [project @ 2006-12-01 05:25:38 by jwe] jwe parents: 5775 diff changeset	342 is.putback (c1);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	343 is >> d;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	344 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	345
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	346 return d;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	347 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	348
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	349 template <>
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	350 Complex
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	351 octave_read_value (std::istream& is)
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	352 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	353 double re = 0.0, im = 0.0;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	354
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	355 Complex cx = 0.0;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	356
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	357 char ch = ' ';
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	358
6907 c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	359 while (isspace (ch))
c8081f4c094f [project @ 2007-09-17 17:31:10 by jwe] jwe parents: 6897 diff changeset	360 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	361
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	362 if (ch == '(')
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	363 {
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	364 re = octave_read_value<double> (is);
6897 935d23e16951 [project @ 2007-09-13 19:13:56 by jwe] jwe parents: 6490 diff changeset	365 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	366
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	367 if (ch == ',')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	368 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	369 im = octave_read_value<double> (is);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	370 ch = is.get ();
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	371
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	372 if (ch == ')')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	373 cx = Complex (re, im);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	374 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	375 is.setstate (std::ios::failbit);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	376 }
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	377 else if (ch == ')')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	378 cx = re;
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	379 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	380 is.setstate (std::ios::failbit);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	381 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	382 else
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	383 {
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	384 is.putback (ch);
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	385 cx = octave_read_value<double> (is);
4130 7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	386 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	387
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	388 return cx;
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	389
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	390 }
7d9bda865012 [project @ 2002-10-28 21:05:30 by jwe] jwe parents: 3970 diff changeset	391
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	392 static inline float
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	393 read_float_inf_nan_na (std::istream& is, char c, char sign = '+')
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	394 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	395 float d = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	396
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	397 switch (c)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	398 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	399 case 'i': case 'I':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	400 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	401 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	402 if (c == 'n' \|\| c == 'N')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	403 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	404 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	405 if (c == 'f' \|\| c == 'F')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	406 d = sign == '-' ? -octave_Inf : octave_Inf;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	407 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	408 is.putback (c);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	409 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	410 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	411 is.putback (c);
7789 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 break;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	414
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	415 case 'n': case 'N':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	416 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	417 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	418 if (c == 'a' \|\| c == 'A')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	419 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	420 c = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	421 if (c == 'n' \|\| c == 'N')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	422 d = octave_NaN;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	423 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	424 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	425 is.putback (c);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	426 d = octave_NA;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	427 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	428 }
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	429 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	430 is.putback (c);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	431 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	432 break;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	433
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	434 default:
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	435 abort ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	436 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	437
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	438 return d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	439 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	440
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	441 template <>
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	442 float
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	443 octave_read_value (std::istream& is)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	444 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	445 float d = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	446
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	447 char c1 = ' ';
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	448
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	449 while (isspace (c1))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	450 c1 = is.get ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	451
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	452 switch (c1)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	453 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	454 case '-':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	455 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	456 char c2 = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	457 c2 = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	458 if (c2 == 'i' \|\| c2 == 'I')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	459 d = read_float_inf_nan_na (is, c2, c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	460 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	461 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	462 is.putback (c2);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	463 is.putback (c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	464 is >> d;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	465 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	466 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	467 break;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	468
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	469 case '+':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	470 {
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	471 char c2 = 0;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	472 c2 = is.get ();
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	473 if (c2 == 'i' \|\| c2 == 'I')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	474 d = read_float_inf_nan_na (is, c2, c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	475 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	476 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	477 is.putback (c2);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	478 is.putback (c1);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	479 is >> d;
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	480 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	481 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	482 break;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	483
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	484 case 'i': case 'I':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	485 case 'n': case 'N':
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	486 d = read_float_inf_nan_na (is, c1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	487 break;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	488
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	489 default:
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	490 is.putback (c1);
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	491 is >> d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	492 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	493
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	494 return d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	495 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	496
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	497 template <>
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	498 FloatComplex
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	499 octave_read_value (std::istream& is)
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	500 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	501 float re = 0.0, im = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	502
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	503 FloatComplex cx = 0.0;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	504
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	505 char ch = ' ';
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	506
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	507 while (isspace (ch))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	508 ch = is.get ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	509
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	510 if (ch == '(')
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	511 {
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	512 re = octave_read_value<float> (is);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	513 ch = is.get ();
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	514
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	515 if (ch == ',')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	516 {
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	517 im = octave_read_value<float> (is);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	518 ch = is.get ();
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	519
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	520 if (ch == ')')
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	521 cx = FloatComplex (re, im);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	522 else
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	523 is.setstate (std::ios::failbit);
07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	524 }
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	525 else if (ch == ')')
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	526 cx = re;
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	527 else
10314 07ebe522dac2 untabify liboctave C++ sources John W. Eaton <jwe@octave.org> parents: 10182 diff changeset	528 is.setstate (std::ios::failbit);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	529 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	530 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	531 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	532 is.putback (ch);
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	533 cx = octave_read_value<float> (is);
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	534 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	535
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	536 return cx;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	537
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	538 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	539
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	540 void
9469 c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	541 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	542 {
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	543 if (lo_ieee_is_NA (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	544 os << "NA";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	545 else if (lo_ieee_isnan (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	546 os << "NaN";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	547 else if (lo_ieee_isinf (d))
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	548 os << (d < 0 ? "-Inf" : "Inf");
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	549 else
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	550 os << d;
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	551 }
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	552
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	553 void
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	554 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	555 {
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	556 os << "(";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	557 octave_write_double (os, real (c));
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	558 os << ",";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	559 octave_write_double (os, imag (c));
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	560 os << ")";
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	561 }
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	562
c6edba80dfae sanity checks for loading sparse matrices John W. Eaton <jwe@octave.org> parents: 8920 diff changeset	563 void
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	564 octave_write_float (std::ostream& os, float d)
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	565 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	566 if (lo_ieee_is_NA (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	567 os << "NA";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	568 else if (lo_ieee_isnan (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	569 os << "NaN";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	570 else if (lo_ieee_isinf (d))
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	571 os << (d < 0 ? "-Inf" : "Inf");
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	572 else
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	573 os << d;
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	574 }
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	575
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	576 void
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	577 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	578 {
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	579 os << "(";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	580 octave_write_float (os, real (c));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	581 os << ",";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	582 octave_write_float (os, imag (c));
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	583 os << ")";
82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7048 diff changeset	584 }

Mercurial > octave

annotate liboctave/lo-utils.cc @ 10447:12d25a1d84bf