Mercurial > octave-nkf
annotate src/ov-scalar.cc @ 7961:a5d1e27ee1f4 ss-3-1-51
3.1.51 snapshot
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Tue, 22 Jul 2008 11:40:48 -0400 |
parents | f336dd8e96d0 |
children | c777f3ce02d8 |
rev | line source |
---|---|
2376 | 1 /* |
2 | |
7017 | 3 Copyright (C) 1996, 1997, 1998, 2000, 2002, 2003, 2004, 2005, 2006, |
4 2007 John W. Eaton | |
2376 | 5 |
6 This file is part of Octave. | |
7 | |
8 Octave is free software; you can redistribute it and/or modify it | |
9 under the terms of the GNU General Public License as published by the | |
7016 | 10 Free Software Foundation; either version 3 of the License, or (at your |
11 option) any later version. | |
2376 | 12 |
13 Octave is distributed in the hope that it will be useful, but WITHOUT | |
14 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
15 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
16 for more details. | |
17 | |
18 You should have received a copy of the GNU General Public License | |
7016 | 19 along with Octave; see the file COPYING. If not, see |
20 <http://www.gnu.org/licenses/>. | |
2376 | 21 |
22 */ | |
23 | |
24 #ifdef HAVE_CONFIG_H | |
25 #include <config.h> | |
26 #endif | |
27 | |
3503 | 28 #include <iostream> |
2901 | 29 |
4944 | 30 #include "data-conv.h" |
31 #include "mach-info.h" | |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
32 #include "lo-specfun.h" |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
33 #include "lo-mappers.h" |
4944 | 34 |
2376 | 35 #include "defun.h" |
36 #include "gripes.h" | |
37 #include "oct-obj.h" | |
4944 | 38 #include "oct-stream.h" |
2376 | 39 #include "ov-scalar.h" |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
40 #include "ov-float.h" |
3223 | 41 #include "ov-base.h" |
42 #include "ov-base-scalar.h" | |
43 #include "ov-base-scalar.cc" | |
2423 | 44 #include "ov-re-mat.h" |
2376 | 45 #include "ov-typeinfo.h" |
46 #include "pr-output.h" | |
47 #include "xdiv.h" | |
48 #include "xpow.h" | |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
49 #include "ops.h" |
2376 | 50 |
4687 | 51 #include "ls-oct-ascii.h" |
52 #include "ls-hdf5.h" | |
53 | |
3223 | 54 template class octave_base_scalar<double>; |
55 | |
3219 | 56 DEFINE_OCTAVE_ALLOCATOR (octave_scalar); |
2376 | 57 |
4612 | 58 DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA (octave_scalar, "scalar", "double"); |
2376 | 59 |
7789
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
60 static octave_base_value * |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
61 default_numeric_demotion_function (const octave_base_value& a) |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
62 { |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
63 CAST_CONV_ARG (const octave_scalar&); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
64 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
65 return new octave_float_scalar (v.float_value ()); |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
66 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
67 |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
68 octave_base_value::type_conv_fcn |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
69 octave_scalar::numeric_demotion_function (void) const |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
70 { |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
71 return default_numeric_demotion_function; |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
72 } |
82be108cc558
First attempt at single precision tyeps
David Bateman <dbateman@free.fr>
parents:
7740
diff
changeset
|
73 |
2376 | 74 octave_value |
5885 | 75 octave_scalar::do_index_op (const octave_value_list& idx, bool resize_ok) |
2376 | 76 { |
77 octave_value retval; | |
78 | |
3933 | 79 if (idx.valid_scalar_indices ()) |
2376 | 80 retval = scalar; |
81 else | |
82 { | |
5775 | 83 // FIXME -- this doesn't solve the problem of |
2376 | 84 // |
85 // a = 1; a([1,1], [1,1], [1,1]) | |
86 // | |
87 // and similar constructions. Hmm... | |
88 | |
5775 | 89 // FIXME -- using this constructor avoids narrowing the |
2423 | 90 // 1x1 matrix back to a scalar value. Need a better solution |
91 // to this problem. | |
92 | |
93 octave_value tmp (new octave_matrix (matrix_value ())); | |
2376 | 94 |
3933 | 95 retval = tmp.do_index_op (idx, resize_ok); |
2376 | 96 } |
97 | |
98 return retval; | |
99 } | |
100 | |
5731 | 101 octave_value |
102 octave_scalar::resize (const dim_vector& dv, bool fill) const | |
103 { | |
104 if (fill) | |
105 { | |
106 NDArray retval (dv, NDArray::resize_fill_value()); | |
107 | |
108 if (dv.numel ()) | |
109 retval(0) = scalar; | |
110 | |
111 return retval; | |
112 } | |
113 else | |
114 { | |
115 NDArray retval (dv); | |
116 | |
117 if (dv.numel ()) | |
118 retval(0) = scalar; | |
119 | |
120 return retval; | |
121 } | |
122 } | |
123 | |
2376 | 124 octave_value |
5279 | 125 octave_scalar::convert_to_str_internal (bool, bool, char type) const |
2376 | 126 { |
127 octave_value retval; | |
128 | |
129 if (xisnan (scalar)) | |
130 ::error ("invalid conversion from NaN to character"); | |
131 else | |
132 { | |
4100 | 133 int ival = NINT (scalar); |
134 | |
135 if (ival < 0 || ival > UCHAR_MAX) | |
136 { | |
5775 | 137 // FIXME -- is there something better we could do? |
4100 | 138 |
139 ival = 0; | |
140 | |
141 ::warning ("range error for conversion to character value"); | |
142 } | |
143 | |
5279 | 144 retval = octave_value (std::string (1, static_cast<char> (ival)), type); |
2376 | 145 } |
146 | |
147 return retval; | |
148 } | |
149 | |
4687 | 150 bool |
6974 | 151 octave_scalar::save_ascii (std::ostream& os) |
4687 | 152 { |
153 double d = double_value (); | |
154 | |
5958 | 155 octave_write_double (os, d); |
156 | |
157 os << "\n"; | |
4687 | 158 |
159 return true; | |
160 } | |
161 | |
162 bool | |
163 octave_scalar::load_ascii (std::istream& is) | |
164 { | |
165 scalar = octave_read_double (is); | |
166 if (!is) | |
167 { | |
168 error ("load: failed to load scalar constant"); | |
169 return false; | |
170 } | |
171 | |
172 return true; | |
173 } | |
174 | |
175 bool | |
176 octave_scalar::save_binary (std::ostream& os, bool& /* save_as_floats */) | |
177 { | |
5760 | 178 char tmp = LS_DOUBLE; |
179 os.write (reinterpret_cast<char *> (&tmp), 1); | |
4687 | 180 double dtmp = double_value (); |
5760 | 181 os.write (reinterpret_cast<char *> (&dtmp), 8); |
4687 | 182 |
183 return true; | |
184 } | |
185 | |
186 bool | |
187 octave_scalar::load_binary (std::istream& is, bool swap, | |
4815 | 188 oct_mach_info::float_format fmt) |
4687 | 189 { |
190 char tmp; | |
5760 | 191 if (! is.read (reinterpret_cast<char *> (&tmp), 1)) |
4687 | 192 return false; |
193 | |
194 double dtmp; | |
5760 | 195 read_doubles (is, &dtmp, static_cast<save_type> (tmp), 1, swap, fmt); |
4687 | 196 if (error_state || ! is) |
197 return false; | |
198 | |
199 scalar = dtmp; | |
200 return true; | |
201 } | |
202 | |
203 #if defined (HAVE_HDF5) | |
4944 | 204 |
4687 | 205 bool |
206 octave_scalar::save_hdf5 (hid_t loc_id, const char *name, | |
207 bool /* save_as_floats */) | |
208 { | |
209 hsize_t dimens[3]; | |
210 hid_t space_hid = -1, data_hid = -1; | |
211 bool retval = true; | |
212 | |
4815 | 213 space_hid = H5Screate_simple (0, dimens, 0); |
4687 | 214 if (space_hid < 0) return false; |
215 | |
216 data_hid = H5Dcreate (loc_id, name, H5T_NATIVE_DOUBLE, space_hid, | |
217 H5P_DEFAULT); | |
218 if (data_hid < 0) | |
219 { | |
220 H5Sclose (space_hid); | |
221 return false; | |
222 } | |
223 | |
224 double tmp = double_value (); | |
225 retval = H5Dwrite (data_hid, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, | |
4815 | 226 H5P_DEFAULT, &tmp) >= 0; |
4687 | 227 |
228 H5Dclose (data_hid); | |
229 H5Sclose (space_hid); | |
4837 | 230 |
4687 | 231 return retval; |
232 } | |
233 | |
234 bool | |
235 octave_scalar::load_hdf5 (hid_t loc_id, const char *name, | |
236 bool /* have_h5giterate_bug */) | |
237 { | |
238 hid_t data_hid = H5Dopen (loc_id, name); | |
239 hid_t space_id = H5Dget_space (data_hid); | |
240 | |
241 hsize_t rank = H5Sget_simple_extent_ndims (space_id); | |
242 | |
243 if (rank != 0) | |
244 { | |
245 H5Dclose (data_hid); | |
246 return false; | |
247 } | |
248 | |
249 double dtmp; | |
250 if (H5Dread (data_hid, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, | |
4815 | 251 H5P_DEFAULT, &dtmp) < 0) |
4687 | 252 { |
253 H5Dclose (data_hid); | |
254 return false; | |
255 } | |
256 | |
257 scalar = dtmp; | |
4837 | 258 |
4687 | 259 H5Dclose (data_hid); |
4837 | 260 |
4687 | 261 return true; |
262 } | |
4944 | 263 |
4687 | 264 #endif |
265 | |
5900 | 266 mxArray * |
267 octave_scalar::as_mxArray (void) const | |
268 { | |
269 mxArray *retval = new mxArray (mxDOUBLE_CLASS, 1, 1, mxREAL); | |
270 | |
271 double *pr = static_cast<double *> (retval->get_data ()); | |
272 | |
273 pr[0] = scalar; | |
274 | |
275 return retval; | |
276 } | |
277 | |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
278 #define SCALAR_MAPPER(MAP, FCN) \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
279 octave_value \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
280 octave_scalar::MAP (void) const \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
281 { \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
282 return octave_value (FCN (scalar)); \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
283 } |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
284 |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
285 #define CD_SCALAR_MAPPER(MAP, RFCN, CFCN, L1, L2) \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
286 octave_value \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
287 octave_scalar::MAP (void) const \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
288 { \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
289 return (scalar < L1 || scalar > L2 \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
290 ? octave_value (CFCN (Complex (scalar))) \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
291 : octave_value (RFCN (scalar))); \ |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
292 } |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
293 |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
294 static double |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
295 xconj (double x) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
296 { |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
297 return x; |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
298 } |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
299 |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
300 SCALAR_MAPPER (erf, ::erf) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
301 SCALAR_MAPPER (erfc, ::erfc) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
302 SCALAR_MAPPER (gamma, xgamma) |
7601
8a939b217863
Treat negative values to lgamma and beta correctly
David Bateman <dbateman@free.fr>
parents:
7529
diff
changeset
|
303 CD_SCALAR_MAPPER (lgamma, xlgamma, xlgamma, 0.0, octave_Inf) |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
304 SCALAR_MAPPER (abs, ::fabs) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
305 SCALAR_MAPPER (acos, ::acos) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
306 CD_SCALAR_MAPPER (acosh, ::acosh, ::acosh, 1.0, octave_Inf) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
307 SCALAR_MAPPER (angle, ::arg) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
308 SCALAR_MAPPER (arg, ::arg) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
309 CD_SCALAR_MAPPER (asin, ::asin, ::asin, -1.0, 1.0) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
310 SCALAR_MAPPER (asinh, ::asinh) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
311 SCALAR_MAPPER (atan, ::atan) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
312 CD_SCALAR_MAPPER (atanh, ::atanh, ::atanh, -1.0, 1.0) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
313 SCALAR_MAPPER (ceil, ::ceil) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
314 SCALAR_MAPPER (conj, xconj) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
315 SCALAR_MAPPER (cos, ::cos) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
316 SCALAR_MAPPER (cosh, ::cosh) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
317 SCALAR_MAPPER (exp, ::exp) |
7638
2df457529cfa
implement expm1 and log1p functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7636
diff
changeset
|
318 SCALAR_MAPPER (expm1, ::expm1) |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
319 SCALAR_MAPPER (fix, ::fix) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
320 SCALAR_MAPPER (floor, ::floor) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
321 SCALAR_MAPPER (imag, ::imag) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
322 CD_SCALAR_MAPPER (log, ::log, std::log, 0.0, octave_Inf) |
7740 | 323 CD_SCALAR_MAPPER (log2, xlog2, xlog2, 0.0, octave_Inf) |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
324 CD_SCALAR_MAPPER (log10, ::log10, std::log10, 0.0, octave_Inf) |
7638
2df457529cfa
implement expm1 and log1p functions
Jaroslav Hajek <highegg@gmail.com>
parents:
7636
diff
changeset
|
325 CD_SCALAR_MAPPER (log1p, ::log1p, ::log1p, -1.0, octave_Inf) |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
326 SCALAR_MAPPER (real, ::real) |
7529
7e1b042c5418
octave_scalar::round: use xround instead of ::round (MSVC does not have ::round).
Michael Goffioul <michael.goffioul@gmail.com>
parents:
7503
diff
changeset
|
327 SCALAR_MAPPER (round, xround) |
7636
99c410f7f0b0
implement mapper function for banker's rounding
Jaroslav Hajek <highegg@gmail.com>
parents:
7601
diff
changeset
|
328 SCALAR_MAPPER (roundb, xroundb) |
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
329 SCALAR_MAPPER (signum, ::signum) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
330 SCALAR_MAPPER (sin, ::sin) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
331 SCALAR_MAPPER (sinh, ::sinh) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
332 CD_SCALAR_MAPPER (sqrt, ::sqrt, std::sqrt, 0.0, octave_Inf) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
333 SCALAR_MAPPER (tan, ::tan) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
334 SCALAR_MAPPER (tanh, ::tanh) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
335 SCALAR_MAPPER (finite, xfinite) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
336 SCALAR_MAPPER (isinf, xisinf) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
337 SCALAR_MAPPER (isna, octave_is_NA) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
338 SCALAR_MAPPER (isnan, xisnan) |
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
339 |
2376 | 340 /* |
341 ;;; Local Variables: *** | |
342 ;;; mode: C++ *** | |
343 ;;; End: *** | |
344 */ |