Mercurial > hg > octave-lojdl
annotate libinterp/corefcn/ls-mat4.cc @ 17438:3856298f1ff8
eliminate unimplemented vax and cray floating point formats
* mach-info.h
* test/io.tst: Delete vaxd, vaxg, and cray from list of architectures
in fopen test.
# HG changeset patch
# User John W. Eaton <jwe@octave.org>
# Date 1379028062 14400
# Thu Sep 12 19:21:02 2013 -0400
# Node ID ea2763020ac1aa15ad0d43d787b61ff3b8469209
# Parent e89e86e1a5518c11058930fec5e907e60e31251b
eliminate unimplemented vax and cray floating point formats
* mach-info.h (oct_mach_info::float_format): Delete flt_fmt_vaxd,
flt_fmt_vaxg, and flt_fmt_cray enum values.
* ls-mat4.cc (float_format_to_mopt_digit, mopt_digit_to_float_format):
Delete cases for vaxd, vaxg, and cray floating point formats.
* lo-ieee.cc (octave_ieee_init): Likewise.
* mach-info.cc (get_float_format,
oct_mach_info::string_to_float_format,
oct_mach_info::float_format_as_string): Likewise.
* oct-stream.cc: Delete explicit template instantiation for do_write.
(octave_stream::write, do_read): Delete checks for vaxd and vaxg
floating point formats.
* pr-output.cc (pr_any_float): Don't check for cray or unknonwn float
formats.
* data-conv.cc (VAX_D_double_to_IEEE_little_double,
VAX_G_double_to_IEEE_little_double, Cray_to_IEEE_little_double,
IEEE_big_float_to_IEEE_little_float, VAX_D_float_to_IEEE_little_float,
VAX_G_float_to_IEEE_little_float, Cray_to_IEEE_little_float,
IEEE_little_double_to_IEEE_big_double,
VAX_D_double_to_IEEE_big_double, VAX_G_double_to_IEEE_big_double,
Cray_to_IEEE_big_double, IEEE_little_float_to_IEEE_big_float,
VAX_D_float_to_IEEE_big_float, VAX_G_float_to_IEEE_big_float,
Cray_to_IEEE_big_float, IEEE_little_double_to_VAX_D_double,
IEEE_big_double_to_VAX_D_double, VAX_G_double_to_VAX_D_double,
Cray_to_VAX_D_double, IEEE_little_float_to_VAX_D_float,
IEEE_big_float_to_VAX_D_float, VAX_G_float_to_VAX_D_float,
Cray_to_VAX_D_float, IEEE_little_double_to_VAX_G_double,
IEEE_big_double_to_VAX_G_double, VAX_D_double_to_VAX_G_double,
Cray_to_VAX_G_double, IEEE_little_float_to_VAX_G_float,
IEEE_big_float_to_VAX_G_float, VAX_D_float_to_VAX_G_float,
Cray_to_VAX_G_float): Delete.
(do_double_format_conversion): Eliminate unimplemented vaxd, vaxg, and
cray conversions.
* file-io.cc (Ffopen, Ffread): Fix doc strings.
| author | John W. Eaton <jwe@octave.org> |
|---|---|
| date | Thu, 12 Sep 2013 19:21:02 -0400 |
| parents | 68fc671a9339 |
| children |
| rev | line source |
|---|---|
| 4634 | 1 /* |
| 2 | |
|
14138
72c96de7a403
maint: update copyright notices for 2012
John W. Eaton <jwe@octave.org>
parents:
11586
diff
changeset
|
3 Copyright (C) 1996-2012 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, | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
77 int len, bool swap, |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
78 oct_mach_info::float_format flt_fmt) |
| 4634 | 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 | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
112 read_mat_file_header (std::istream& is, bool& swap, int32_t& mopt, |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
113 int32_t& nr, int32_t& nc, |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
114 int32_t& imag, int32_t& len, |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
115 int quiet) |
| 4634 | 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) | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
168 error ("load: can't read binary file"); |
| 4634 | 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 case 3: | |
| 198 case 4: | |
| 199 default: | |
| 200 flt_fmt = oct_mach_info::flt_fmt_unknown; | |
| 201 break; | |
| 202 } | |
| 203 | |
| 204 return flt_fmt; | |
| 205 } | |
| 206 | |
| 207 int | |
| 208 float_format_to_mopt_digit (oct_mach_info::float_format flt_fmt) | |
| 209 { | |
| 210 int retval = -1; | |
| 211 | |
| 212 switch (flt_fmt) | |
| 213 { | |
| 214 case oct_mach_info::flt_fmt_ieee_little_endian: | |
| 215 retval = 0; | |
| 216 break; | |
| 217 | |
| 218 case oct_mach_info::flt_fmt_ieee_big_endian: | |
| 219 retval = 1; | |
| 220 break; | |
| 221 | |
| 222 default: | |
| 223 break; | |
| 224 } | |
| 225 | |
| 226 return retval; | |
| 227 } | |
| 228 | |
| 229 // Extract one value (scalar, matrix, string, etc.) from stream IS and | |
| 230 // place it in TC, returning the name of the variable. | |
| 231 // | |
| 232 // The data is expected to be in Matlab version 4 .mat format, though | |
| 233 // not all the features of that format are supported. | |
| 234 // | |
| 235 // FILENAME is used for error messages. | |
| 236 // | |
| 237 // This format provides no way to tag the data as global. | |
| 238 | |
| 239 std::string | |
| 240 read_mat_binary_data (std::istream& is, const std::string& filename, | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
241 octave_value& tc) |
| 4634 | 242 { |
| 243 std::string retval; | |
| 244 | |
| 245 // These are initialized here instead of closer to where they are | |
| 246 // first used to avoid errors from gcc about goto crossing | |
| 247 // initialization of variable. | |
| 248 | |
| 249 Matrix re; | |
| 250 oct_mach_info::float_format flt_fmt = oct_mach_info::flt_fmt_unknown; | |
| 251 bool swap = false; | |
| 252 int type = 0; | |
| 253 int prec = 0; | |
| 254 int order = 0; | |
| 255 int mach = 0; | |
| 256 int dlen = 0; | |
| 257 | |
| 5828 | 258 int32_t mopt, nr, nc, imag, len; |
| 4634 | 259 |
| 260 int err = read_mat_file_header (is, swap, mopt, nr, nc, imag, len); | |
| 261 if (err) | |
| 262 { | |
| 263 if (err < 0) | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
264 goto data_read_error; |
| 4634 | 265 else |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
266 return retval; |
| 4634 | 267 } |
| 268 | |
| 269 type = mopt % 10; // Full, sparse, etc. | |
| 270 mopt /= 10; // Eliminate first digit. | |
| 271 prec = mopt % 10; // double, float, int, etc. | |
| 272 mopt /= 10; // Eliminate second digit. | |
| 273 order = mopt % 10; // Row or column major ordering. | |
| 274 mopt /= 10; // Eliminate third digit. | |
| 275 mach = mopt % 10; // IEEE, VAX, etc. | |
| 276 | |
| 277 flt_fmt = mopt_digit_to_float_format (mach); | |
| 278 | |
| 279 if (flt_fmt == oct_mach_info::flt_fmt_unknown) | |
| 280 { | |
| 281 error ("load: unrecognized binary format!"); | |
| 282 return retval; | |
| 283 } | |
| 284 | |
| 285 if (imag && type == 1) | |
| 286 { | |
| 287 error ("load: encountered complex matrix with string flag set!"); | |
| 288 return retval; | |
| 289 } | |
| 290 | |
| 291 // LEN includes the terminating character, and the file is also | |
| 292 // supposed to include it, but apparently not all files do. Either | |
| 293 // way, I think this should work. | |
| 294 | |
| 295 { | |
| 296 OCTAVE_LOCAL_BUFFER (char, name, len+1); | |
| 297 name[len] = '\0'; | |
| 5760 | 298 if (! is.read (name, len)) |
| 4634 | 299 goto data_read_error; |
| 300 retval = name; | |
| 301 | |
| 302 dlen = nr * nc; | |
| 303 if (dlen < 0) | |
| 304 goto data_read_error; | |
| 305 | |
| 306 if (order) | |
| 307 { | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
308 octave_idx_type tmp = nr; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
309 nr = nc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
310 nc = tmp; |
| 4634 | 311 } |
| 312 | |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
313 if (type == 2) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
314 { |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
315 if (nc == 4) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
316 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
317 octave_idx_type nr_new, nc_new; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
318 Array<Complex> data (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
319 Array<octave_idx_type> c (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
320 Array<octave_idx_type> r (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
321 OCTAVE_LOCAL_BUFFER (double, dtmp, nr); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
322 OCTAVE_LOCAL_BUFFER (double, ctmp, nr); |
| 4634 | 323 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
324 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
325 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
326 r.xelem (i) = dtmp[i] - 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
327 nr_new = dtmp[nr - 1]; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
328 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
329 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
330 c.xelem (i) = dtmp[i] - 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
331 nc_new = dtmp[nr - 1]; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
332 read_mat_binary_data (is, dtmp, prec, nr - 1, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
333 read_mat_binary_data (is, ctmp, prec, 1, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
334 read_mat_binary_data (is, ctmp, prec, nr - 1, swap, flt_fmt); |
| 4634 | 335 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
336 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
337 data.xelem (i) = Complex (dtmp[i], ctmp[i]); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
338 read_mat_binary_data (is, ctmp, prec, 1, swap, flt_fmt); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
339 |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
340 SparseComplexMatrix smc = SparseComplexMatrix (data, r, c, |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
341 nr_new, nc_new); |
| 4634 | 342 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
343 tc = order ? smc.transpose () : smc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
344 } |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
345 else |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
346 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
347 octave_idx_type nr_new, nc_new; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
348 Array<double> data (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
349 Array<octave_idx_type> c (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
350 Array<octave_idx_type> r (dim_vector (1, nr - 1)); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
351 OCTAVE_LOCAL_BUFFER (double, dtmp, nr); |
| 4634 | 352 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
353 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
354 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
355 r.xelem (i) = dtmp[i] - 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
356 nr_new = dtmp[nr - 1]; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
357 read_mat_binary_data (is, dtmp, prec, nr, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
358 for (octave_idx_type i = 0; i < nr - 1; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
359 c.xelem (i) = dtmp[i] - 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
360 nc_new = dtmp[nr - 1]; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
361 read_mat_binary_data (is, data.fortran_vec (), prec, nr - 1, swap, flt_fmt); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
362 read_mat_binary_data (is, dtmp, prec, 1, swap, flt_fmt); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
363 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
364 SparseMatrix sm = SparseMatrix (data, r, c, nr_new, nc_new); |
| 4634 | 365 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
366 tc = order ? sm.transpose () : sm; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
367 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
368 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
369 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
370 { |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
371 re.resize (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
372 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
373 read_mat_binary_data (is, re.fortran_vec (), prec, dlen, swap, flt_fmt); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
374 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
375 if (! is || error_state) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
376 { |
|
15466
d174210ce1ec
use ' instead of ` in error messages, warnings and most comments
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
377 error ("load: reading matrix data for '%s'", name); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
378 goto data_read_error; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
379 } |
| 4634 | 380 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
381 if (imag) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
382 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
383 Matrix im (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
384 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
385 read_mat_binary_data (is, im.fortran_vec (), prec, dlen, swap, |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
386 flt_fmt); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
387 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
388 if (! is || error_state) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
389 { |
|
15466
d174210ce1ec
use ' instead of ` in error messages, warnings and most comments
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
390 error ("load: reading imaginary matrix data for '%s'", name); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
391 goto data_read_error; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
392 } |
| 4634 | 393 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
394 ComplexMatrix ctmp (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
395 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
396 for (octave_idx_type j = 0; j < nc; j++) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
397 for (octave_idx_type i = 0; i < nr; i++) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
398 ctmp (i, j) = Complex (re (i, j), im (i, j)); |
| 4634 | 399 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
400 tc = order ? ctmp.transpose () : ctmp; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
401 } |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
402 else |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
403 tc = order ? re.transpose () : re; |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
404 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
405 if (type == 1) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
406 tc = tc.convert_to_str (false, true, '\''); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
407 } |
| 4634 | 408 |
| 409 return retval; | |
| 410 } | |
| 411 | |
| 412 data_read_error: | |
|
15466
d174210ce1ec
use ' instead of ` in error messages, warnings and most comments
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
413 error ("load: trouble reading binary file '%s'", filename.c_str ()); |
| 4634 | 414 return retval; |
| 415 } | |
| 416 | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
417 // Save the data from TC along with the corresponding NAME on stream OS |
| 4634 | 418 // in the MatLab version 4 binary format. |
| 419 | |
| 420 bool | |
| 421 save_mat_binary_data (std::ostream& os, const octave_value& tc, | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
422 const std::string& name) |
| 4634 | 423 { |
| 5828 | 424 int32_t mopt = 0; |
| 4634 | 425 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
426 mopt += tc.is_sparse_type () ? 2 : tc.is_string () ? 1 : 0; |
| 4634 | 427 |
| 428 oct_mach_info::float_format flt_fmt = | |
| 429 oct_mach_info::native_float_format ();; | |
| 430 | |
| 431 mopt += 1000 * float_format_to_mopt_digit (flt_fmt); | |
| 432 | |
| 5760 | 433 os.write (reinterpret_cast<char *> (&mopt), 4); |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
434 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
435 octave_idx_type len; |
| 5828 | 436 int32_t nr = tc.rows (); |
| 4634 | 437 |
| 5828 | 438 int32_t nc = tc.columns (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
439 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
440 if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
441 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
442 len = tc.nnz (); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
443 uint32_t nnz = len + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
444 os.write (reinterpret_cast<char *> (&nnz), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
445 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
446 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
|
447 os.write (reinterpret_cast<char *> (&iscmplx), 4); |
| 4634 | 448 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
449 uint32_t tmp = 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
450 os.write (reinterpret_cast<char *> (&tmp), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
451 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
452 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
453 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
454 os.write (reinterpret_cast<char *> (&nr), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
455 os.write (reinterpret_cast<char *> (&nc), 4); |
| 4634 | 456 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
457 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
|
458 os.write (reinterpret_cast<char *> (&imag), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
459 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
460 len = nr * nc; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
461 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
462 |
| 4634 | 463 |
| 464 // LEN includes the terminating character, and the file is also | |
| 465 // supposed to include it. | |
| 466 | |
| 5828 | 467 int32_t name_len = name.length () + 1; |
| 4634 | 468 |
| 5760 | 469 os.write (reinterpret_cast<char *> (&name_len), 4); |
| 4634 | 470 os << name << '\0'; |
| 471 | |
| 472 if (tc.is_string ()) | |
| 473 { | |
|
10066
2cd940306a06
make unwind_protect frames local
Jaroslav Hajek <highegg@gmail.com>
parents:
9377
diff
changeset
|
474 unwind_protect frame; |
| 4634 | 475 |
| 476 charMatrix chm = tc.char_matrix_value (); | |
| 477 | |
| 5275 | 478 octave_idx_type nrow = chm.rows (); |
| 479 octave_idx_type ncol = chm.cols (); | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
480 |
| 4634 | 481 OCTAVE_LOCAL_BUFFER (double, buf, ncol*nrow); |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
482 |
| 5275 | 483 for (octave_idx_type i = 0; i < nrow; i++) |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
484 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
485 std::string tstr = chm.row_as_string (i); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
486 const char *s = tstr.data (); |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
487 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
488 for (octave_idx_type j = 0; j < ncol; j++) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
489 buf[j*nrow+i] = static_cast<double> (*s++ & 0x00FF); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
490 } |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
491 os.write (reinterpret_cast<char *> (buf), nrow*ncol*sizeof (double)); |
| 4634 | 492 } |
| 493 else if (tc.is_range ()) | |
| 494 { | |
| 495 Range r = tc.range_value (); | |
| 496 double base = r.base (); | |
| 497 double inc = r.inc (); | |
| 5275 | 498 octave_idx_type nel = r.nelem (); |
| 499 for (octave_idx_type i = 0; i < nel; i++) | |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
500 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
501 double x = base + i * inc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
502 os.write (reinterpret_cast<char *> (&x), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
503 } |
| 4634 | 504 } |
| 505 else if (tc.is_real_scalar ()) | |
| 506 { | |
| 507 double tmp = tc.double_value (); | |
| 5760 | 508 os.write (reinterpret_cast<char *> (&tmp), 8); |
| 4634 | 509 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
510 else if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
511 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
512 double ds; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
513 OCTAVE_LOCAL_BUFFER (double, dtmp, len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
514 if (tc.is_complex_matrix ()) |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
515 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
516 SparseComplexMatrix m = tc.sparse_complex_matrix_value (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
517 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
518 for (octave_idx_type i = 0; i < len; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
519 dtmp[i] = m.ridx (i) + 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
520 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
521 ds = nr; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
522 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
523 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
524 octave_idx_type ii = 0; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
525 for (octave_idx_type j = 0; j < nc; j++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
526 for (octave_idx_type i = m.cidx (j); i < m.cidx (j+1); i++) |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
527 dtmp[ii++] = j + 1; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
528 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
529 ds = nc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
530 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
531 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
532 for (octave_idx_type i = 0; i < len; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
533 dtmp[i] = std::real (m.data (i)); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
534 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
535 ds = 0.; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
536 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
537 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
538 for (octave_idx_type i = 0; i < len; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
539 dtmp[i] = std::imag (m.data (i)); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
540 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
541 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
542 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
543 else |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
544 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
545 SparseMatrix m = tc.sparse_matrix_value (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
546 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
547 for (octave_idx_type i = 0; i < len; i++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
548 dtmp[i] = m.ridx (i) + 1; |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
549 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
550 ds = nr; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
551 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
552 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
553 octave_idx_type ii = 0; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
554 for (octave_idx_type j = 0; j < nc; j++) |
|
14861
f7afecdd87ef
maint: Use Octave coding conventions for cuddling parentheses in src/ directory
Rik <octave@nomad.inbox5.com>
parents:
14138
diff
changeset
|
555 for (octave_idx_type i = m.cidx (j); i < m.cidx (j+1); i++) |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
556 dtmp[ii++] = j + 1; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
557 os.write (reinterpret_cast<const char *> (dtmp), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
558 ds = nc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
559 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
560 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
561 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
562 ds = 0.; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
563 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
564 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
565 } |
| 4634 | 566 else if (tc.is_real_matrix ()) |
| 567 { | |
| 568 Matrix m = tc.matrix_value (); | |
| 5760 | 569 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 570 } |
| 571 else if (tc.is_complex_scalar ()) | |
| 572 { | |
| 573 Complex tmp = tc.complex_value (); | |
| 5760 | 574 os.write (reinterpret_cast<char *> (&tmp), 16); |
| 4634 | 575 } |
| 576 else if (tc.is_complex_matrix ()) | |
| 577 { | |
| 578 ComplexMatrix m_cmplx = tc.complex_matrix_value (); | |
| 579 Matrix m = ::real (m_cmplx); | |
| 5760 | 580 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 581 m = ::imag (m_cmplx); |
| 5760 | 582 os.write (reinterpret_cast<const char *> (m.data ()), 8 * len); |
| 4634 | 583 } |
| 584 else | |
| 585 gripe_wrong_type_arg ("save", tc, false); | |
| 586 | |
| 587 return os; | |
| 588 } |