Mercurial > octave-nkf

annotate src/DLD-FUNCTIONS/str2double.cc @ 10378:eeb6c09ec51a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
avoid g++ extension in str2double
author Jaroslav Hajek <highegg@gmail.com>
date Mon, 01 Mar 2010 09:49:56 +0100
parents 7658cd4bdcf2
children f578e6468d0c
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
1 /*
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
2
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
3 Copyright (C) 2010 Jaroslav Hajek
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
4 Copyright (C) 2010 VZLU Prague
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
5
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
6 This file is part of Octave.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
7
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
8 Octave is free software; you can redistribute it and/or modify it
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
9 under the terms of the GNU General Public License as published by the
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
10 Free Software Foundation; either version 3 of the License, or (at your
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
11 option) any later version.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
12
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
13 Octave is distributed in the hope that it will be useful, but WITHOUT
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
16 for more details.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
17
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
18 You should have received a copy of the GNU General Public License
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
19 along with Octave; see the file COPYING. If not, see
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
20 <http://www.gnu.org/licenses/>.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
21
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
22 */
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
23
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
24 #ifdef HAVE_CONFIG_H
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
25 #include <config.h>
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
26 #endif
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
27
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
28 #include <string>
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
29 #include <cctype>
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
30 #include <sstream>
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
31 #include <algorithm>
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
32
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
33 #include "lo-ieee.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
34
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
35 #include "Cell.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
36 #include "ov.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
37 #include "defun-dld.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
38 #include "gripes.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
39 #include "utils.h"
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
40
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
41 static inline bool
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
42 is_imag_unit (int c)
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
43 { return c == 'i' || c == 'j'; }
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
44
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
45 static std::istringstream&
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
46 single_num (std::istringstream& is, double& num)
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
47 {
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
48 char c = is.peek ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
49 if (c == 'I')
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
50 {
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
51 // It's infinity.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
52 is.get ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
53 char c1 = is.get (), c2 = is.get ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
54 if (c1 == 'n' && c2 == 'f')
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
55 {
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
56 num = octave_Inf;
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
57 is.peek (); // Sets eof bit.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
58 }
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
59 else
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
60 is.setstate (std::ios::failbit); // indicate that read has failed.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
61 }
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
62 else if (c == 'N')
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
63 {
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
64 // It's NaN.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
65 is.get ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
66 char c1 = is.get (), c2 = is.get ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
67 if (c1 == 'a' && c2 == 'N')
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
68 {
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
69 num = octave_NaN;
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
70 is.peek (); // Sets eof bit.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
71 }
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
72 else
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
73 is.setstate (std::ios::failbit); // indicate that read has failed.
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
74 }
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
75 else
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
76 is >> num;
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
77
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
78 return is;
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
79 }
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
80
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
81 static std::istringstream&
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
82 extract_num (std::istringstream& is, double& num, bool& imag, bool& have_sign)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
83 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
84 have_sign = imag = false;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
85
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
86 char c = is.peek ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
87 bool negative = false;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
88
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
89 // Accept leading sign.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
90 if (c == '+' || c == '-')
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
91 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
92 negative = c == '-';
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
93 is.get ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
94 c = is.peek ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
95 have_sign = true;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
96 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
97
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
98 if (is_imag_unit (c))
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
99 {
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
100 c = is.get ();
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
101 // It's i*num or just i.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
102 imag = true;
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
103 char cn = is.peek ();
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
104 if (cn == '*')
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
105 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
106 // Multiplier follows, we extract it as a number.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
107 is.get ();
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
108 single_num (is, num);
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
109 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
110 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
111 num = 1.0;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
112 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
113 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
114 {
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
115 single_num (is, num);
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
116 if (is.good ())
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
117 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
118 c = is.peek ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
119 if (c == '*')
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
120 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
121 is.get ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
122 c = is.get ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
123 if (is_imag_unit (c))
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
124 imag = true;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
125 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
126 is.setstate (std::ios::failbit); // indicate that read has failed.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
127 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
128 else if (is_imag_unit (c))
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
129 imag = true;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
130 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
131 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
132
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
133 if (negative)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
134 num = -num;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
135
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
136 return is;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
137 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
138
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
139 static Complex
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
140 str2double1 (std::string str)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
141 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
142 Complex val (0.0, 0.0);
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
143 std::string::iterator se = str.end ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
144 // Remove commas (thousand separators) and spaces.
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
145 se = std::remove (str.begin (), se, ',');
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
146 se = std::remove (str.begin (), se, ' ');
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
147 str.erase (se, str.end ());
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
148
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
149 std::istringstream is (str);
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
150 double num;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
151 bool i1, i2, s1, s2;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
152
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
153 if (is.eof ())
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
154 val = octave_NaN;
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
155 else if (! extract_num (is, num, i1, s1))
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
156 val = octave_NaN;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
157 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
158 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
159 if (i1)
10378
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
160 #ifdef __GNUC__
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
161 val.imag () = num; // GNU C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
162 #else
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
163 val = Complex (val.real (), num); // ISO C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
164 #endif
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
165 else
10378
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
166 #ifdef __GNUC__
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
167 val.real () = num; // GNU C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
168 #else
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
169 val = Complex (num, val.imag ()); // ISO C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
170 #endif
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
171
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
172 if (! is.eof ())
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
173 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
174 if (! extract_num (is, num, i2, s2) || i1 == i2 || ! s2)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
175 val = octave_NaN;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
176 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
177 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
178 if (i2)
10378
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
179 #ifdef __GNUC__
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
180 val.imag () = num; // GNU C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
181 #else
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
182 val = Complex (val.real (), num); // ISO C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
183 #endif
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
184 else
10378
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
185 #ifdef __GNUC__
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
186 val.real () = num; // GNU C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
187 #else
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
188 val = Complex (num, val.imag ()); // ISO C++
eeb6c09ec51a avoid g++ extension in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10357
diff changeset
189 #endif
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
190 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
191 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
192 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
193
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
194 return val;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
195 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
196
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
197 DEFUN_DLD (str2double, args, ,
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
198 "-*- texinfo -*-\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
199 @deftypefn {Built-in Function} {} str2double (@var{s})\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
200 Converts a string to real or complex number.\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
201 A complex number should be in one of the formats:\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
202 \n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
203 @itemize\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
204 @item a + bi\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
205 @item a + b*i\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
206 @item a + i*b\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
207 @item bi + a\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
208 @item b*i + a\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
209 @item i*b + a\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
210 @end itemize\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
211 \n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
212 It is also possible to use @code{j} instead of @code{i}, or write just\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
213 @code{i} instead of @code{1*i}.\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
214 @code{a} and @code{b} should be real numbers\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
215 in a standard format.\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
216 @var{s} can also be a character matrix, in which case the conversion is repeated\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
217 for each row, or a cell array of strings, in which case each element is converted\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
218 and an array of the same dimensions is returned.\n\
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
219 @end deftypefn")
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
220 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
221 octave_value retval;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
222
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
223 if (args.length () == 1)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
224 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
225 if (args(0).is_string ())
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
226 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
227 if (args(0).rows () == 1 && args(0).ndims () == 2)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
228 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
229 retval = str2double1 (args(0).string_value ());
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
230 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
231 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
232 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
233 const string_vector sv = args(0).all_strings ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
234 if (! error_state)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
235 retval = sv.map<Complex> (str2double1);
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
236 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
237 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
238 else if (args(0).is_cellstr ())
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
239 {
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
240 Array<std::string> sa = args(0).cellstr_value ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
241 retval = sa.map<Complex> (str2double1);
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
242 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
243 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
244 gripe_wrong_type_arg ("str2double", args(0));
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
245 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
246 else
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
247 print_usage ();
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
248
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
249 return retval;
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
250 }
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
251
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
252 /*
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
253
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
254 %!assert (str2double ("1"), 1)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
255 %!assert (str2double ("-.1e-5"), -1e-6)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
256 %!assert (str2double ("1,222.5"), 1222.5)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
257 %!assert (str2double ("i"), i)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
258 %!assert (str2double ("2 + j"), 2+j)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
259 %!assert (str2double ("i*2 + 3"), 3+2i)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
260 %!assert (str2double (".5*i + 3.5"), 3.5+0.5i)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
261 %!assert (str2double ("1e-3 + i*.25"), 1e-3 + i*.25)
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
262 %!assert (str2double (["2 + j";"1.25e-3";"-05"]), [2+i; 1.25e-3; -05])
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
263 %!assert (str2double ({"2 + j","1.25e-3","-05"}), [2+i, 1.25e-3, -05])
10357
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
264 %!assert (str2double ("NaN"), NaN)
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
265 %!assert (str2double ("Inf"), Inf)
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
266 %!assert (str2double ("-Inf"), -Inf)
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
267 %!assert (str2double ("i*Inf"), complex (0, Inf))
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
268 %!assert (str2double ("NaN + Inf*i"), complex (NaN, Inf))
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
269 %!assert (str2double ("Inf - Inf*i"), complex (Inf, -Inf))
7658cd4bdcf2 handle Infs and NaNs in str2double
Jaroslav Hajek <highegg@gmail.com>
parents: 10356
diff changeset
270 %!assert (str2double ("-i*NaN - Inf"), complex (-Inf, -NaN))
10356
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
271
1d7930b77ab9 implement compiled str2double
Jaroslav Hajek <highegg@gmail.com>
parents:
diff changeset
272 */