Mercurial > hg > octave-nkf
annotate libinterp/corefcn/ls-mat4.cc @ 20750:3339c9bdfe6a
Activate FSAL property in dorpri timestepper
* scripts/ode/private/runge_kutta_45_dorpri.m: don't compute
first stage if values from previous iteration are passed.
* scripts/ode/private/integrate_adaptive.m: do not update
cmputed stages if timestep is rejected.
| author | Carlo de Falco <carlo.defalco@polimi.it> |
|---|---|
| date | Sat, 03 Oct 2015 07:32:50 +0200 |
| parents | 00cf2847355d |
| children | f90c8372b7ba |
| rev | line source |
|---|---|
| 4634 | 1 /* |
| 2 | |
|
19898
4197fc428c7d
maint: Update copyright notices for 2015.
John W. Eaton <jwe@octave.org>
parents:
18846
diff
changeset
|
3 Copyright (C) 1996-2015 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 | |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
174 data_read_error: |
| 4634 | 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]; |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
361 read_mat_binary_data (is, data.fortran_vec (), prec, nr - 1, |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
362 swap, flt_fmt); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
363 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
|
364 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
365 SparseMatrix sm = SparseMatrix (data, r, c, nr_new, nc_new); |
| 4634 | 366 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
367 tc = order ? sm.transpose () : sm; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
368 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
369 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
370 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
371 { |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
372 re.resize (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
373 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
374 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
|
375 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
376 if (! is || error_state) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
377 { |
|
15466
d174210ce1ec
use ' instead of ` in error messages, warnings and most comments
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
378 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
|
379 goto data_read_error; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
380 } |
| 4634 | 381 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
382 if (imag) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
383 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
384 Matrix im (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
385 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
386 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
|
387 flt_fmt); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
388 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
389 if (! is || error_state) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
390 { |
|
15466
d174210ce1ec
use ' instead of ` in error messages, warnings and most comments
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
391 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
|
392 goto data_read_error; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
393 } |
| 4634 | 394 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
395 ComplexMatrix ctmp (nr, nc); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
396 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
397 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
|
398 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
|
399 ctmp (i, j) = Complex (re (i, j), im (i, j)); |
| 4634 | 400 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
401 tc = order ? ctmp.transpose () : ctmp; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
402 } |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
403 else |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
404 tc = order ? re.transpose () : re; |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
405 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
406 if (type == 1) |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
407 tc = tc.convert_to_str (false, true, '\''); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
408 } |
| 4634 | 409 |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
410 return retval; |
|
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
411 } |
| 4634 | 412 |
|
17787
175b392e91fe
Use GNU style coding conventions for code in libinterp/
Rik <rik@octave.org>
parents:
17744
diff
changeset
|
413 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
|
414 error ("load: trouble reading binary file '%s'", filename.c_str ()); |
| 4634 | 415 return retval; |
| 416 } | |
| 417 | |
|
11586
12df7854fa7c
strip trailing whitespace from source files
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
418 // Save the data from TC along with the corresponding NAME on stream OS |
| 4634 | 419 // in the MatLab version 4 binary format. |
| 420 | |
| 421 bool | |
| 422 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
|
423 const std::string& name) |
| 4634 | 424 { |
| 5828 | 425 int32_t mopt = 0; |
| 4634 | 426 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
427 mopt += tc.is_sparse_type () ? 2 : tc.is_string () ? 1 : 0; |
| 4634 | 428 |
| 429 oct_mach_info::float_format flt_fmt = | |
| 430 oct_mach_info::native_float_format ();; | |
| 431 | |
| 432 mopt += 1000 * float_format_to_mopt_digit (flt_fmt); | |
| 433 | |
| 5760 | 434 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
|
435 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
436 octave_idx_type len; |
| 5828 | 437 int32_t nr = tc.rows (); |
| 4634 | 438 |
| 5828 | 439 int32_t nc = tc.columns (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
440 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
441 if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
442 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
443 len = tc.nnz (); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
444 uint32_t nnz = len + 1; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
445 os.write (reinterpret_cast<char *> (&nnz), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
446 |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
447 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
|
448 os.write (reinterpret_cast<char *> (&iscmplx), 4); |
| 4634 | 449 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
450 uint32_t tmp = 0; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
451 os.write (reinterpret_cast<char *> (&tmp), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
452 } |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
453 else |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
454 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
455 os.write (reinterpret_cast<char *> (&nr), 4); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
456 os.write (reinterpret_cast<char *> (&nc), 4); |
| 4634 | 457 |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
458 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
|
459 os.write (reinterpret_cast<char *> (&imag), 4); |
|
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 len = nr * nc; |
|
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 |
| 4634 | 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 } |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
491 std::streamsize n_bytes = static_cast<std::streamsize> (nrow) * |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
492 static_cast<std::streamsize> (ncol) * |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
493 sizeof (double); |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
494 os.write (reinterpret_cast<char *> (buf), n_bytes); |
| 4634 | 495 } |
| 496 else if (tc.is_range ()) | |
| 497 { | |
| 498 Range r = tc.range_value (); | |
| 499 double base = r.base (); | |
| 500 double inc = r.inc (); | |
|
20438
00cf2847355d
Deprecate Array::nelem() and Range::nelem() in favour of ::numel().
Carnë Draug <carandraug@octave.org>
parents:
19898
diff
changeset
|
501 octave_idx_type nel = r.numel (); |
| 5275 | 502 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
|
503 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
504 double x = base + i * inc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
505 os.write (reinterpret_cast<char *> (&x), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
506 } |
| 4634 | 507 } |
| 508 else if (tc.is_real_scalar ()) | |
| 509 { | |
| 510 double tmp = tc.double_value (); | |
| 5760 | 511 os.write (reinterpret_cast<char *> (&tmp), 8); |
| 4634 | 512 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
513 else if (tc.is_sparse_type ()) |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
514 { |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
515 double ds; |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
516 OCTAVE_LOCAL_BUFFER (double, dtmp, len); |
|
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
517 if (tc.is_complex_matrix ()) |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
518 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
519 SparseComplexMatrix m = tc.sparse_complex_matrix_value (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
520 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
521 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
|
522 dtmp[i] = m.ridx (i) + 1; |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
523 std::streamsize n_bytes = 8 * static_cast<std::streamsize> (len); |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
524 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
525 ds = nr; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
526 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
|
527 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
528 octave_idx_type ii = 0; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
529 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
|
530 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
|
531 dtmp[ii++] = j + 1; |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
532 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
533 ds = nc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
534 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
|
535 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
536 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
|
537 dtmp[i] = std::real (m.data (i)); |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
538 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
539 ds = 0.; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
540 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
|
541 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
542 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
|
543 dtmp[i] = std::imag (m.data (i)); |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
544 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
545 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
546 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
547 else |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
548 { |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
549 SparseMatrix m = tc.sparse_matrix_value (); |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
550 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
551 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
|
552 dtmp[i] = m.ridx (i) + 1; |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
553 std::streamsize n_bytes = 8 * static_cast<std::streamsize> (len); |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
554 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
555 ds = nr; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
556 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
|
557 |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
558 octave_idx_type ii = 0; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
559 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
|
560 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
|
561 dtmp[ii++] = j + 1; |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
562 os.write (reinterpret_cast<const char *> (dtmp), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
563 ds = nc; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
564 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
|
565 |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
566 os.write (reinterpret_cast<const char *> (m.data ()), n_bytes); |
|
10315
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
567 ds = 0.; |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
568 os.write (reinterpret_cast<const char *> (&ds), 8); |
|
57a59eae83cc
untabify src C++ source files
John W. Eaton <jwe@octave.org>
parents:
10160
diff
changeset
|
569 } |
|
7918
78eef61f75d5
Add matlab V4 sparse matrixload/save
David Bateman <dbateman@free.fr>
parents:
7336
diff
changeset
|
570 } |
| 4634 | 571 else if (tc.is_real_matrix ()) |
| 572 { | |
| 573 Matrix m = tc.matrix_value (); | |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
574 std::streamsize n_bytes = 8 * static_cast<std::streamsize> (len); |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
575 os.write (reinterpret_cast<const char *> (m.data ()), n_bytes); |
| 4634 | 576 } |
| 577 else if (tc.is_complex_scalar ()) | |
| 578 { | |
| 579 Complex tmp = tc.complex_value (); | |
| 5760 | 580 os.write (reinterpret_cast<char *> (&tmp), 16); |
| 4634 | 581 } |
| 582 else if (tc.is_complex_matrix ()) | |
| 583 { | |
| 584 ComplexMatrix m_cmplx = tc.complex_matrix_value (); | |
| 585 Matrix m = ::real (m_cmplx); | |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
586 std::streamsize n_bytes = 8 * static_cast<std::streamsize> (len); |
|
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
587 os.write (reinterpret_cast<const char *> (m.data ()), n_bytes); |
| 4634 | 588 m = ::imag (m_cmplx); |
|
18846
9a6646cc7c3e
Fix load/save to handle up to 2^32-1 elements.
Rik <rik@octave.org>
parents:
18452
diff
changeset
|
589 os.write (reinterpret_cast<const char *> (m.data ()), n_bytes); |
| 4634 | 590 } |
| 591 else | |
| 592 gripe_wrong_type_arg ("save", tc, false); | |
| 593 | |
|
18452
bd9d34f28b0f
Use std::ostream::fail instead of unsafe implicit bool conversion (bug #41335)
Mike Miller <mtmiller@ieee.org>
parents:
17787
diff
changeset
|
594 return ! os.fail (); |
| 4634 | 595 } |