Mercurial > octave
annotate src/ov-scalar.cc @ 7503:8c32f95c2639
convert mapper functions to new format
| author | David Bateman <dbateman@free.fr> |
|---|---|
| date | Wed, 20 Feb 2008 04:22:50 -0500 |
| parents | a1dbe9d80eee |
| children | 7e1b042c5418 |
| 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" |
| 3223 | 40 #include "ov-base.h" |
| 41 #include "ov-base-scalar.h" | |
| 42 #include "ov-base-scalar.cc" | |
| 2423 | 43 #include "ov-re-mat.h" |
| 2376 | 44 #include "ov-typeinfo.h" |
| 45 #include "pr-output.h" | |
| 46 #include "xdiv.h" | |
| 47 #include "xpow.h" | |
| 48 | |
| 4687 | 49 #include "ls-oct-ascii.h" |
| 50 #include "ls-hdf5.h" | |
| 51 | |
| 3223 | 52 template class octave_base_scalar<double>; |
| 53 | |
| 3219 | 54 DEFINE_OCTAVE_ALLOCATOR (octave_scalar); |
| 2376 | 55 |
| 4612 | 56 DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA (octave_scalar, "scalar", "double"); |
| 2376 | 57 |
| 58 octave_value | |
| 5885 | 59 octave_scalar::do_index_op (const octave_value_list& idx, bool resize_ok) |
| 2376 | 60 { |
| 61 octave_value retval; | |
| 62 | |
| 3933 | 63 if (idx.valid_scalar_indices ()) |
| 2376 | 64 retval = scalar; |
| 65 else | |
| 66 { | |
| 5775 | 67 // FIXME -- this doesn't solve the problem of |
| 2376 | 68 // |
| 69 // a = 1; a([1,1], [1,1], [1,1]) | |
| 70 // | |
| 71 // and similar constructions. Hmm... | |
| 72 | |
| 5775 | 73 // FIXME -- using this constructor avoids narrowing the |
| 2423 | 74 // 1x1 matrix back to a scalar value. Need a better solution |
| 75 // to this problem. | |
| 76 | |
| 77 octave_value tmp (new octave_matrix (matrix_value ())); | |
| 2376 | 78 |
| 3933 | 79 retval = tmp.do_index_op (idx, resize_ok); |
| 2376 | 80 } |
| 81 | |
| 82 return retval; | |
| 83 } | |
| 84 | |
| 4645 | 85 std::streamoff |
| 86 octave_scalar::streamoff_value (void) const | |
| 87 { | |
| 88 std::streamoff retval (-1); | |
| 89 | |
| 90 if (D_NINT (scalar) == scalar) | |
| 91 retval = std::streamoff (static_cast<long> (scalar)); | |
| 92 else | |
| 93 error ("conversion to streamoff value failed"); | |
| 94 | |
| 95 return retval; | |
| 96 } | |
| 97 | |
| 4701 | 98 streamoff_array |
| 99 octave_scalar::streamoff_array_value (void) const | |
| 100 { | |
| 101 streamoff_array retval; | |
| 102 | |
| 103 std::streamoff soff = streamoff_value (); | |
| 104 | |
| 105 if (! error_state) | |
| 106 retval = streamoff_array (dim_vector (1, 1), soff); | |
| 107 | |
| 108 return retval; | |
| 109 } | |
| 110 | |
| 5731 | 111 octave_value |
| 112 octave_scalar::resize (const dim_vector& dv, bool fill) const | |
| 113 { | |
| 114 if (fill) | |
| 115 { | |
| 116 NDArray retval (dv, NDArray::resize_fill_value()); | |
| 117 | |
| 118 if (dv.numel ()) | |
| 119 retval(0) = scalar; | |
| 120 | |
| 121 return retval; | |
| 122 } | |
| 123 else | |
| 124 { | |
| 125 NDArray retval (dv); | |
| 126 | |
| 127 if (dv.numel ()) | |
| 128 retval(0) = scalar; | |
| 129 | |
| 130 return retval; | |
| 131 } | |
| 132 } | |
| 133 | |
| 2376 | 134 octave_value |
| 5279 | 135 octave_scalar::convert_to_str_internal (bool, bool, char type) const |
| 2376 | 136 { |
| 137 octave_value retval; | |
| 138 | |
| 139 if (xisnan (scalar)) | |
| 140 ::error ("invalid conversion from NaN to character"); | |
| 141 else | |
| 142 { | |
| 4100 | 143 int ival = NINT (scalar); |
| 144 | |
| 145 if (ival < 0 || ival > UCHAR_MAX) | |
| 146 { | |
| 5775 | 147 // FIXME -- is there something better we could do? |
| 4100 | 148 |
| 149 ival = 0; | |
| 150 | |
| 151 ::warning ("range error for conversion to character value"); | |
| 152 } | |
| 153 | |
| 5279 | 154 retval = octave_value (std::string (1, static_cast<char> (ival)), type); |
| 2376 | 155 } |
| 156 | |
| 157 return retval; | |
| 158 } | |
| 159 | |
| 4687 | 160 bool |
| 6974 | 161 octave_scalar::save_ascii (std::ostream& os) |
| 4687 | 162 { |
| 163 double d = double_value (); | |
| 164 | |
| 5958 | 165 octave_write_double (os, d); |
| 166 | |
| 167 os << "\n"; | |
| 4687 | 168 |
| 169 return true; | |
| 170 } | |
| 171 | |
| 172 bool | |
| 173 octave_scalar::load_ascii (std::istream& is) | |
| 174 { | |
| 175 scalar = octave_read_double (is); | |
| 176 if (!is) | |
| 177 { | |
| 178 error ("load: failed to load scalar constant"); | |
| 179 return false; | |
| 180 } | |
| 181 | |
| 182 return true; | |
| 183 } | |
| 184 | |
| 185 bool | |
| 186 octave_scalar::save_binary (std::ostream& os, bool& /* save_as_floats */) | |
| 187 { | |
| 5760 | 188 char tmp = LS_DOUBLE; |
| 189 os.write (reinterpret_cast<char *> (&tmp), 1); | |
| 4687 | 190 double dtmp = double_value (); |
| 5760 | 191 os.write (reinterpret_cast<char *> (&dtmp), 8); |
| 4687 | 192 |
| 193 return true; | |
| 194 } | |
| 195 | |
| 196 bool | |
| 197 octave_scalar::load_binary (std::istream& is, bool swap, | |
| 4815 | 198 oct_mach_info::float_format fmt) |
| 4687 | 199 { |
| 200 char tmp; | |
| 5760 | 201 if (! is.read (reinterpret_cast<char *> (&tmp), 1)) |
| 4687 | 202 return false; |
| 203 | |
| 204 double dtmp; | |
| 5760 | 205 read_doubles (is, &dtmp, static_cast<save_type> (tmp), 1, swap, fmt); |
| 4687 | 206 if (error_state || ! is) |
| 207 return false; | |
| 208 | |
| 209 scalar = dtmp; | |
| 210 return true; | |
| 211 } | |
| 212 | |
| 213 #if defined (HAVE_HDF5) | |
| 4944 | 214 |
| 4687 | 215 bool |
| 216 octave_scalar::save_hdf5 (hid_t loc_id, const char *name, | |
| 217 bool /* save_as_floats */) | |
| 218 { | |
| 219 hsize_t dimens[3]; | |
| 220 hid_t space_hid = -1, data_hid = -1; | |
| 221 bool retval = true; | |
| 222 | |
| 4815 | 223 space_hid = H5Screate_simple (0, dimens, 0); |
| 4687 | 224 if (space_hid < 0) return false; |
| 225 | |
| 226 data_hid = H5Dcreate (loc_id, name, H5T_NATIVE_DOUBLE, space_hid, | |
| 227 H5P_DEFAULT); | |
| 228 if (data_hid < 0) | |
| 229 { | |
| 230 H5Sclose (space_hid); | |
| 231 return false; | |
| 232 } | |
| 233 | |
| 234 double tmp = double_value (); | |
| 235 retval = H5Dwrite (data_hid, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, | |
| 4815 | 236 H5P_DEFAULT, &tmp) >= 0; |
| 4687 | 237 |
| 238 H5Dclose (data_hid); | |
| 239 H5Sclose (space_hid); | |
| 4837 | 240 |
| 4687 | 241 return retval; |
| 242 } | |
| 243 | |
| 244 bool | |
| 245 octave_scalar::load_hdf5 (hid_t loc_id, const char *name, | |
| 246 bool /* have_h5giterate_bug */) | |
| 247 { | |
| 248 hid_t data_hid = H5Dopen (loc_id, name); | |
| 249 hid_t space_id = H5Dget_space (data_hid); | |
| 250 | |
| 251 hsize_t rank = H5Sget_simple_extent_ndims (space_id); | |
| 252 | |
| 253 if (rank != 0) | |
| 254 { | |
| 255 H5Dclose (data_hid); | |
| 256 return false; | |
| 257 } | |
| 258 | |
| 259 double dtmp; | |
| 260 if (H5Dread (data_hid, H5T_NATIVE_DOUBLE, H5S_ALL, H5S_ALL, | |
| 4815 | 261 H5P_DEFAULT, &dtmp) < 0) |
| 4687 | 262 { |
| 263 H5Dclose (data_hid); | |
| 264 return false; | |
| 265 } | |
| 266 | |
| 267 scalar = dtmp; | |
| 4837 | 268 |
| 4687 | 269 H5Dclose (data_hid); |
| 4837 | 270 |
| 4687 | 271 return true; |
| 272 } | |
| 4944 | 273 |
| 4687 | 274 #endif |
| 275 | |
| 5900 | 276 mxArray * |
| 277 octave_scalar::as_mxArray (void) const | |
| 278 { | |
| 279 mxArray *retval = new mxArray (mxDOUBLE_CLASS, 1, 1, mxREAL); | |
| 280 | |
| 281 double *pr = static_cast<double *> (retval->get_data ()); | |
| 282 | |
| 283 pr[0] = scalar; | |
| 284 | |
| 285 return retval; | |
| 286 } | |
| 287 | |
|
7503
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
288 #define SCALAR_MAPPER(MAP, FCN) \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
289 octave_value \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
290 octave_scalar::MAP (void) const \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
291 { \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
292 return octave_value (FCN (scalar)); \ |
|
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 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
295 #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
|
296 octave_value \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
297 octave_scalar::MAP (void) const \ |
|
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 return (scalar < L1 || scalar > L2 \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
300 ? octave_value (CFCN (Complex (scalar))) \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
301 : octave_value (RFCN (scalar))); \ |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
302 } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
303 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
304 static double |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
305 xconj (double x) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
306 { |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
307 return x; |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
308 } |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
309 |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
310 SCALAR_MAPPER (erf, ::erf) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
311 SCALAR_MAPPER (erfc, ::erfc) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
312 SCALAR_MAPPER (gamma, xgamma) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
313 SCALAR_MAPPER (lgamma, xlgamma) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
314 SCALAR_MAPPER (abs, ::fabs) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
315 SCALAR_MAPPER (acos, ::acos) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
316 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
|
317 SCALAR_MAPPER (angle, ::arg) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
318 SCALAR_MAPPER (arg, ::arg) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
319 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
|
320 SCALAR_MAPPER (asinh, ::asinh) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
321 SCALAR_MAPPER (atan, ::atan) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
322 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
|
323 SCALAR_MAPPER (ceil, ::ceil) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
324 SCALAR_MAPPER (conj, xconj) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
325 SCALAR_MAPPER (cos, ::cos) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
326 SCALAR_MAPPER (cosh, ::cosh) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
327 SCALAR_MAPPER (exp, ::exp) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
328 SCALAR_MAPPER (fix, ::fix) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
329 SCALAR_MAPPER (floor, ::floor) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
330 SCALAR_MAPPER (imag, ::imag) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
331 CD_SCALAR_MAPPER (log, ::log, std::log, 0.0, octave_Inf) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
332 CD_SCALAR_MAPPER (log10, ::log10, std::log10, 0.0, octave_Inf) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
333 SCALAR_MAPPER (real, ::real) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
334 SCALAR_MAPPER (round, ::round) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
335 SCALAR_MAPPER (signum, ::signum) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
336 SCALAR_MAPPER (sin, ::sin) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
337 SCALAR_MAPPER (sinh, ::sinh) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
338 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
|
339 SCALAR_MAPPER (tan, ::tan) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
340 SCALAR_MAPPER (tanh, ::tanh) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
341 SCALAR_MAPPER (finite, xfinite) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
342 SCALAR_MAPPER (isinf, xisinf) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
343 SCALAR_MAPPER (isna, octave_is_NA) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
344 SCALAR_MAPPER (isnan, xisnan) |
|
8c32f95c2639
convert mapper functions to new format
David Bateman <dbateman@free.fr>
parents:
7017
diff
changeset
|
345 |
| 2376 | 346 /* |
| 347 ;;; Local Variables: *** | |
| 348 ;;; mode: C++ *** | |
| 349 ;;; End: *** | |
| 350 */ |