Mercurial > octave-nkf
annotate src/ls-mat4.cc @ 8874:bd1b1fe9c6e9 ss-3-1-53
bump version info for snapshot
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Wed, 25 Feb 2009 18:35:47 -0500 |
| parents | 25bc2d31e1bf |
| children | eb63fbe60fab |
| rev | line source |
|---|---|
| 4634 | 1 /* |
| 2 | |
| 7017 | 3 Copyright (C) 1996, 1997, 2003, 2004, 2005, 2006, 2007 John W. Eaton |
| 4634 | 4 |
| 5 This file is part of Octave. | |
| 6 | |
| 7 Octave is free software; you can redistribute it and/or modify it | |
| 8 under the terms of the GNU General Public License as published by the | |
| 7016 | 9 Free Software Foundation; either version 3 of the License, or (at your |
| 10 option) any later version. | |
| 4634 | 11 |
| 12 Octave is distributed in the hope that it will be useful, but WITHOUT | |
| 13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
| 14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License | |
| 15 for more details. | |
| 16 | |
| 17 You should have received a copy of the GNU General Public License | |
| 7016 | 18 along with Octave; see the file COPYING. If not, see |
| 19 <http://www.gnu.org/licenses/>. | |
| 4634 | 20 |
| 21 */ | |
| 22 | |
| 23 #ifdef HAVE_CONFIG_H | |
| 24 #include <config.h> | |
| 25 #endif | |
| 26 | |
| 27 #include <cfloat> | |
| 28 #include <cstring> | |
| 29 #include <cctype> | |
| 30 | |
| 31 #include <fstream> | |
| 32 #include <iomanip> | |
| 33 #include <iostream> | |
| 34 #include <string> | |
| 4726 | 35 #include <vector> |
| 4634 | 36 |
| 37 #include "byte-swap.h" | |
| 38 #include "data-conv.h" | |
| 39 #include "file-ops.h" | |
| 40 #include "glob-match.h" | |
| 41 #include "lo-mappers.h" | |
| 42 #include "mach-info.h" | |
| 43 #include "oct-env.h" | |
| 44 #include "oct-time.h" | |
| 45 #include "quit.h" | |
| 46 #include "str-vec.h" | |
|
8377
25bc2d31e1bf
improve OCTAVE_LOCAL_BUFFER
Jaroslav Hajek <highegg@gmail.com>
parents:
7918
diff
changeset
|
47 #include "oct-locbuf.h" |
| 4634 | 48 |
| 49 #include "Cell.h" | |
| 50 #include "defun.h" | |
| 51 #include "error.h" | |
| 52 #include "gripes.h" | |
| 53 #include "load-save.h" | |
| 54 #include "oct-obj.h" | |
| 55 #include "oct-map.h" | |
| 56 #include "ov-cell.h" | |
| 57 #include "pager.h" | |
| 58 #include "pt-exp.h" | |
| 59 #include "sysdep.h" | |
| 60 #include "unwind-prot.h" | |
| 61 #include "utils.h" | |
| 62 #include "variables.h" | |
| 63 #include "version.h" | |
| 64 #include "dMatrix.h" | |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
65 #include "dSparse.h" |
| 4634 | 66 |
| 67 #include "ls-mat4.h" | |
| 68 | |
| 69 // Read LEN elements of data from IS in the format specified by | |
| 70 // PRECISION, placing the result in DATA. If SWAP is TRUE, swap | |
| 71 // the bytes of each element before copying to DATA. FLT_FMT | |
| 72 // specifies the format of the data if we are reading floating point | |
| 73 // numbers. | |
| 74 | |
| 75 static void | |
| 76 read_mat_binary_data (std::istream& is, double *data, int precision, | |
| 77 int len, bool swap, | |
| 78 oct_mach_info::float_format flt_fmt) | |
| 79 { | |
| 80 switch (precision) | |
| 81 { | |
| 82 case 0: | |
| 83 read_doubles (is, data, LS_DOUBLE, len, swap, flt_fmt); | |
| 84 break; | |
| 85 | |
| 86 case 1: | |
| 87 read_doubles (is, data, LS_FLOAT, len, swap, flt_fmt); | |
| 88 break; | |
| 89 | |
| 90 case 2: | |
| 91 read_doubles (is, data, LS_INT, len, swap, flt_fmt); | |
| 92 break; | |
| 93 | |
| 94 case 3: | |
| 95 read_doubles (is, data, LS_SHORT, len, swap, flt_fmt); | |
| 96 break; | |
| 97 | |
| 98 case 4: | |
| 99 read_doubles (is, data, LS_U_SHORT, len, swap, flt_fmt); | |
| 100 break; | |
| 101 | |
| 102 case 5: | |
| 103 read_doubles (is, data, LS_U_CHAR, len, swap, flt_fmt); | |
| 104 break; | |
| 105 | |
| 106 default: | |
| 107 break; | |
| 108 } | |
| 109 } | |
| 110 | |
| 111 int | |
| 5828 | 112 read_mat_file_header (std::istream& is, bool& swap, int32_t& mopt, |
| 113 int32_t& nr, int32_t& nc, | |
| 114 int32_t& imag, int32_t& len, | |
| 4634 | 115 int quiet) |
| 116 { | |
| 117 swap = false; | |
| 118 | |
| 119 // We expect to fail here, at the beginning of a record, so not | |
| 120 // being able to read another mopt value should not result in an | |
| 121 // error. | |
| 122 | |
| 5760 | 123 is.read (reinterpret_cast<char *> (&mopt), 4); |
| 4634 | 124 if (! is) |
| 125 return 1; | |
| 126 | |
| 5760 | 127 if (! is.read (reinterpret_cast<char *> (&nr), 4)) |
| 4634 | 128 goto data_read_error; |
| 129 | |
| 5760 | 130 if (! is.read (reinterpret_cast<char *> (&nc), 4)) |
| 4634 | 131 goto data_read_error; |
| 132 | |
| 5760 | 133 if (! is.read (reinterpret_cast<char *> (&imag), 4)) |
| 4634 | 134 goto data_read_error; |
| 135 | |
| 5760 | 136 if (! is.read (reinterpret_cast<char *> (&len), 4)) |
| 4634 | 137 goto data_read_error; |
| 138 | |
| 139 // If mopt is nonzero and the byte order is swapped, mopt will be | |
| 140 // bigger than we expect, so we swap bytes. | |
| 141 // | |
| 142 // If mopt is zero, it means the file was written on a little endian | |
| 143 // machine, and we only need to swap if we are running on a big endian | |
| 144 // machine. | |
| 145 // | |
| 146 // Gag me. | |
| 147 | |
| 148 if (oct_mach_info::words_big_endian () && mopt == 0) | |
| 149 swap = true; | |
| 150 | |
| 151 // mopt is signed, therefore byte swap may result in negative value. | |
| 152 | |
| 153 if (mopt > 9999 || mopt < 0) | |
| 154 swap = true; | |
| 155 | |
| 156 if (swap) | |
| 157 { | |
| 4944 | 158 swap_bytes<4> (&mopt); |
| 159 swap_bytes<4> (&nr); | |
| 160 swap_bytes<4> (&nc); | |
| 161 swap_bytes<4> (&imag); | |
| 162 swap_bytes<4> (&len); | |
| 4634 | 163 } |
| 164 | |
| 165 if (mopt > 9999 || mopt < 0 || imag > 1 || imag < 0) | |
| 166 { | |
| 167 if (! quiet) | |
| 168 error ("load: can't read binary file"); | |
| 169 return -1; | |
| 170 } | |
| 171 | |
| 172 return 0; | |
| 173 | |
| 174 data_read_error: | |
| 175 return -1; | |
| 176 } | |
| 177 | |
| 178 // We don't just use a cast here, because we need to be able to detect | |
| 179 // possible errors. | |
| 180 | |
| 181 oct_mach_info::float_format | |
| 182 mopt_digit_to_float_format (int mach) | |
| 183 { | |
| 184 oct_mach_info::float_format flt_fmt = oct_mach_info::flt_fmt_unknown; | |
| 185 | |
| 186 switch (mach) | |
| 187 { | |
| 188 case 0: | |
| 189 flt_fmt = oct_mach_info::flt_fmt_ieee_little_endian; | |
| 190 break; | |
| 191 | |
| 192 case 1: | |
| 193 flt_fmt = oct_mach_info::flt_fmt_ieee_big_endian; | |
| 194 break; | |
| 195 | |
| 196 case 2: | |
| 197 flt_fmt = oct_mach_info::flt_fmt_vax_d; | |
| 198 break; | |
| 199 | |
| 200 case 3: | |
| 201 flt_fmt = oct_mach_info::flt_fmt_vax_g; | |
| 202 break; | |
| 203 | |
| 204 case 4: | |
| 205 flt_fmt = oct_mach_info::flt_fmt_cray; | |
| 206 break; | |
| 207 | |
| 208 default: | |
| 209 flt_fmt = oct_mach_info::flt_fmt_unknown; | |
| 210 break; | |
| 211 } | |
| 212 | |
| 213 return flt_fmt; | |
| 214 } | |
| 215 | |
| 216 int | |
| 217 float_format_to_mopt_digit (oct_mach_info::float_format flt_fmt) | |
| 218 { | |
| 219 int retval = -1; | |
| 220 | |
| 221 switch (flt_fmt) | |
| 222 { | |
| 223 case oct_mach_info::flt_fmt_ieee_little_endian: | |
| 224 retval = 0; | |
| 225 break; | |
| 226 | |
| 227 case oct_mach_info::flt_fmt_ieee_big_endian: | |
| 228 retval = 1; | |
| 229 break; | |
| 230 | |
| 231 case oct_mach_info::flt_fmt_vax_d: | |
| 232 retval = 2; | |
| 233 break; | |
| 234 | |
| 235 case oct_mach_info::flt_fmt_vax_g: | |
| 236 retval = 3; | |
| 237 break; | |
| 238 | |
| 239 case oct_mach_info::flt_fmt_cray: | |
| 240 retval = 4; | |
| 241 break; | |
| 242 | |
| 243 default: | |
| 244 break; | |
| 245 } | |
| 246 | |
| 247 return retval; | |
| 248 } | |
| 249 | |
| 250 // Extract one value (scalar, matrix, string, etc.) from stream IS and | |
| 251 // place it in TC, returning the name of the variable. | |
| 252 // | |
| 253 // The data is expected to be in Matlab version 4 .mat format, though | |
| 254 // not all the features of that format are supported. | |
| 255 // | |
| 256 // FILENAME is used for error messages. | |
| 257 // | |
| 258 // This format provides no way to tag the data as global. | |
| 259 | |
| 260 std::string | |
| 261 read_mat_binary_data (std::istream& is, const std::string& filename, | |
| 262 octave_value& tc) | |
| 263 { | |
| 264 std::string retval; | |
| 265 | |
| 266 // These are initialized here instead of closer to where they are | |
| 267 // first used to avoid errors from gcc about goto crossing | |
| 268 // initialization of variable. | |
| 269 | |
| 270 Matrix re; | |
| 271 oct_mach_info::float_format flt_fmt = oct_mach_info::flt_fmt_unknown; | |
| 272 bool swap = false; | |
| 273 int type = 0; | |
| 274 int prec = 0; | |
| 275 int order = 0; | |
| 276 int mach = 0; | |
| 277 int dlen = 0; | |
| 278 | |
| 5828 | 279 int32_t mopt, nr, nc, imag, len; |
| 4634 | 280 |
| 281 int err = read_mat_file_header (is, swap, mopt, nr, nc, imag, len); | |
| 282 if (err) | |
| 283 { | |
| 284 if (err < 0) | |
| 285 goto data_read_error; | |
| 286 else | |
| 287 return retval; | |
| 288 } | |
| 289 | |
| 290 type = mopt % 10; // Full, sparse, etc. | |
| 291 mopt /= 10; // Eliminate first digit. | |
| 292 prec = mopt % 10; // double, float, int, etc. | |
| 293 mopt /= 10; // Eliminate second digit. | |
| 294 order = mopt % 10; // Row or column major ordering. | |
| 295 mopt /= 10; // Eliminate third digit. | |
| 296 mach = mopt % 10; // IEEE, VAX, etc. | |
| 297 | |
| 298 flt_fmt = mopt_digit_to_float_format (mach); | |
| 299 | |
| 300 if (flt_fmt == oct_mach_info::flt_fmt_unknown) | |
| 301 { | |
| 302 error ("load: unrecognized binary format!"); | |
| 303 return retval; | |
| 304 } | |
| 305 | |
| 306 if (imag && type == 1) | |
| 307 { | |
| 308 error ("load: encountered complex matrix with string flag set!"); | |
| 309 return retval; | |
| 310 } | |
| 311 | |
| 312 // LEN includes the terminating character, and the file is also | |
| 313 // supposed to include it, but apparently not all files do. Either | |
| 314 // way, I think this should work. | |
| 315 | |
| 316 { | |
| 317 OCTAVE_LOCAL_BUFFER (char, name, len+1); | |
| 318 name[len] = '\0'; | |
| 5760 | 319 if (! is.read (name, len)) |
| 4634 | 320 goto data_read_error; |
| 321 retval = name; | |
| 322 | |
| 323 dlen = nr * nc; | |
| 324 if (dlen < 0) | |
| 325 goto data_read_error; | |
| 326 | |
| 327 if (order) | |
| 328 { | |
| 5275 | 329 octave_idx_type tmp = nr; |
| 4634 | 330 nr = nc; |
| 331 nc = tmp; | |
| 332 } | |
| 333 | |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
334 if (type == 2) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
335 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
336 if (nc == 4) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
337 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
338 octave_idx_type nr_new, nc_new; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
339 Array<Complex> data (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
340 Array<octave_idx_type> c (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
341 Array<octave_idx_type> r (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
342 OCTAVE_LOCAL_BUFFER (double, dtmp, nr); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
343 OCTAVE_LOCAL_BUFFER (double, ctmp, nr); |
| 4634 | 344 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
345 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
346 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
347 r.xelem(i) = dtmp[i] - 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
348 nr_new = dtmp[nr - 1]; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
349 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
350 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
351 c.xelem(i) = dtmp[i] - 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
352 nc_new = dtmp[nr - 1]; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
353 read_mat_binary_data (is, dtmp, prec, nr - 1, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
354 read_mat_binary_data (is, ctmp, prec, 1, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
355 read_mat_binary_data (is, ctmp, prec, nr - 1, swap, flt_fmt); |
| 4634 | 356 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
357 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
358 data.xelem(i) = Complex (dtmp[i], ctmp[i]); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
359 read_mat_binary_data (is, ctmp, prec, 1, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
360 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
361 SparseComplexMatrix smc = SparseComplexMatrix (data, r, c, |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
362 nr_new, nc_new); |
| 4634 | 363 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
364 tc = order ? smc.transpose () : smc; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
365 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
366 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
367 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
368 octave_idx_type nr_new, nc_new; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
369 Array<double> data (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
370 Array<octave_idx_type> c (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
371 Array<octave_idx_type> r (dim_vector (1, nr - 1)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
372 OCTAVE_LOCAL_BUFFER (double, dtmp, nr); |
| 4634 | 373 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
374 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
375 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
376 r.xelem(i) = dtmp[i] - 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
377 nr_new = dtmp[nr - 1]; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
378 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
379 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
380 c.xelem(i) = dtmp[i] - 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
381 nc_new = dtmp[nr - 1]; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
382 read_mat_binary_data (is, data.fortran_vec (), prec, nr - 1, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
383 read_mat_binary_data (is, dtmp, prec, 1, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
384 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
385 SparseMatrix sm = SparseMatrix (data, r, c, nr_new, nc_new); |
| 4634 | 386 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
387 tc = order ? sm.transpose () : sm; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
388 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
389 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
390 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
391 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
392 re.resize (nr, nc); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
393 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
394 read_mat_binary_data (is, re.fortran_vec (), prec, dlen, swap, flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
395 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
396 if (! is || error_state) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
397 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
398 error ("load: reading matrix data for `%s'", name); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
399 goto data_read_error; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
400 } |
| 4634 | 401 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
402 if (imag) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
403 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
404 Matrix im (nr, nc); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
405 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
406 read_mat_binary_data (is, im.fortran_vec (), prec, dlen, swap, |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
407 flt_fmt); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
408 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
409 if (! is || error_state) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
410 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
411 error ("load: reading imaginary matrix data for `%s'", name); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
412 goto data_read_error; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
413 } |
| 4634 | 414 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
415 ComplexMatrix ctmp (nr, nc); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
416 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
417 for (octave_idx_type j = 0; j < nc; j++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
418 for (octave_idx_type i = 0; i < nr; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
419 ctmp (i, j) = Complex (re (i, j), im (i, j)); |
| 4634 | 420 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
421 tc = order ? ctmp.transpose () : ctmp; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
422 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
423 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
424 tc = order ? re.transpose () : re; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
425 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
426 if (type == 1) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
427 tc = tc.convert_to_str (false, true, '\''); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
428 } |
| 4634 | 429 |
| 430 return retval; | |
| 431 } | |
| 432 | |
| 433 data_read_error: | |
| 434 error ("load: trouble reading binary file `%s'", filename.c_str ()); | |
| 435 return retval; | |
| 436 } | |
| 437 | |
| 438 // Save the data from TC along with the corresponding NAME on stream OS | |
| 439 // in the MatLab version 4 binary format. | |
| 440 | |
| 441 bool | |
| 442 save_mat_binary_data (std::ostream& os, const octave_value& tc, | |
| 443 const std::string& name) | |
| 444 { | |
| 5828 | 445 int32_t mopt = 0; |
| 4634 | 446 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
447 mopt += tc.is_sparse_type () ? 2 : tc.is_string () ? 1 : 0; |
| 4634 | 448 |
| 449 oct_mach_info::float_format flt_fmt = | |
| 450 oct_mach_info::native_float_format ();; | |
| 451 | |
| 452 mopt += 1000 * float_format_to_mopt_digit (flt_fmt); | |
| 453 | |
| 5760 | 454 os.write (reinterpret_cast<char *> (&mopt), 4); |
| 4634 | 455 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
456 octave_idx_type len; |
| 5828 | 457 int32_t nr = tc.rows (); |
| 4634 | 458 |
| 5828 | 459 int32_t nc = tc.columns (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
460 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
461 if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
462 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
463 len = tc.nnz (); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
464 uint32_t nnz = len + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
465 os.write (reinterpret_cast<char *> (&nnz), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
466 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
467 uint32_t iscmplx = tc.is_complex_type () ? 4 : 3; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
468 os.write (reinterpret_cast<char *> (&iscmplx), 4); |
| 4634 | 469 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
470 uint32_t tmp = 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
471 os.write (reinterpret_cast<char *> (&tmp), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
472 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
473 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
474 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
475 os.write (reinterpret_cast<char *> (&nr), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
476 os.write (reinterpret_cast<char *> (&nc), 4); |
| 4634 | 477 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
478 int32_t imag = tc.is_complex_type () ? 1 : 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
479 os.write (reinterpret_cast<char *> (&imag), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
480 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
481 len = nr * nc; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
482 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
483 |
| 4634 | 484 |
| 485 // LEN includes the terminating character, and the file is also | |
| 486 // supposed to include it. | |
| 487 | |
| 5828 | 488 int32_t name_len = name.length () + 1; |
| 4634 | 489 |
| 5760 | 490 os.write (reinterpret_cast<char *> (&name_len), 4); |
| 4634 | 491 os << name << '\0'; |
| 492 | |
| 493 if (tc.is_string ()) | |
| 494 { | |
| 495 unwind_protect::begin_frame ("save_mat_binary_data"); | |
| 496 | |
| 497 charMatrix chm = tc.char_matrix_value (); | |
| 498 | |
| 5275 | 499 octave_idx_type nrow = chm.rows (); |
| 500 octave_idx_type ncol = chm.cols (); | |
| 4634 | 501 |
| 502 OCTAVE_LOCAL_BUFFER (double, buf, ncol*nrow); | |
| 503 | |
| 5275 | 504 for (octave_idx_type i = 0; i < nrow; i++) |
| 4634 | 505 { |
| 506 std::string tstr = chm.row_as_string (i); | |
| 507 const char *s = tstr.data (); | |
| 508 | |
| 5275 | 509 for (octave_idx_type j = 0; j < ncol; j++) |
| 4634 | 510 buf[j*nrow+i] = static_cast<double> (*s++ & 0x00FF); |
| 511 } | |
| 5760 | 512 os.write (reinterpret_cast<char *> (buf), nrow*ncol*sizeof(double)); |
| 4634 | 513 |
| 514 unwind_protect::run_frame ("save_mat_binary_data"); | |
| 515 } | |
| 516 else if (tc.is_range ()) | |
| 517 { | |
| 518 Range r = tc.range_value (); | |
| 519 double base = r.base (); | |
| 520 double inc = r.inc (); | |
| 5275 | 521 octave_idx_type nel = r.nelem (); |
| 522 for (octave_idx_type i = 0; i < nel; i++) | |
| 4634 | 523 { |
| 524 double x = base + i * inc; | |
| 5760 | 525 os.write (reinterpret_cast<char *> (&x), 8); |
| 4634 | 526 } |
| 527 } | |
| 528 else if (tc.is_real_scalar ()) | |
| 529 { | |
| 530 double tmp = tc.double_value (); | |
| 5760 | 531 os.write (reinterpret_cast<char *> (&tmp), 8); |
| 4634 | 532 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
533 else if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
534 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
535 double ds; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
536 OCTAVE_LOCAL_BUFFER (double, dtmp, len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
537 if (tc.is_complex_matrix ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
538 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
539 SparseComplexMatrix m = tc.sparse_complex_matrix_value (); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
540 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
541 for (octave_idx_type i = 0; i < len; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
542 dtmp [i] = m.ridx(i) + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
543 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
544 ds = nr; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
545 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
546 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
547 octave_idx_type ii = 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
548 for (octave_idx_type j = 0; j < nc; j++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
549 for (octave_idx_type i = m.cidx(j); i < m.cidx(j+1); i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
550 dtmp[ii++] = j + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
551 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
552 ds = nc; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
553 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
554 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
555 for (octave_idx_type i = 0; i < len; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
556 dtmp [i] = std::real (m.data(i)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
557 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
558 ds = 0.; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
559 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
560 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
561 for (octave_idx_type i = 0; i < len; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
562 dtmp [i] = std::imag (m.data(i)); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
563 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
564 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
565 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
566 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
567 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
568 SparseMatrix m = tc.sparse_matrix_value (); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
569 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
570 for (octave_idx_type i = 0; i < len; i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
571 dtmp [i] = m.ridx(i) + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
572 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
573 ds = nr; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
574 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
575 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
576 octave_idx_type ii = 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
577 for (octave_idx_type j = 0; j < nc; j++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
578 for (octave_idx_type i = m.cidx(j); i < m.cidx(j+1); i++) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
579 dtmp[ii++] = j + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
580 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
581 ds = nc; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
582 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
583 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
584 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
585 ds = 0.; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
586 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
587 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
588 } |
| 4634 | 589 else if (tc.is_real_matrix ()) |
| 590 { | |
| 591 Matrix m = tc.matrix_value (); | |
| 5760 | 592 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 593 } |
| 594 else if (tc.is_complex_scalar ()) | |
| 595 { | |
| 596 Complex tmp = tc.complex_value (); | |
| 5760 | 597 os.write (reinterpret_cast<char *> (&tmp), 16); |
| 4634 | 598 } |
| 599 else if (tc.is_complex_matrix ()) | |
| 600 { | |
| 601 ComplexMatrix m_cmplx = tc.complex_matrix_value (); | |
| 602 Matrix m = ::real (m_cmplx); | |
| 5760 | 603 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 604 m = ::imag (m_cmplx); |
| 5760 | 605 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 606 } |
| 607 else | |
| 608 gripe_wrong_type_arg ("save", tc, false); | |
| 609 | |
| 610 return os; | |
| 611 } | |
| 612 | |
| 613 /* | |
| 614 ;;; Local Variables: *** | |
| 615 ;;; mode: C++ *** | |
| 616 ;;; End: *** | |
| 617 */ |