|
604
|
1 /* |
|
|
2 |
|
2847
|
3 Copyright (C) 1996, 1997 John W. Eaton |
|
604
|
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 |
|
|
9 Free Software Foundation; either version 2, or (at your option) any |
|
|
10 later version. |
|
|
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 |
|
|
18 along with Octave; see the file COPYING. If not, write to the Free |
|
1315
|
19 Software Foundation, 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. |
|
604
|
20 |
|
|
21 */ |
|
|
22 |
|
3911
|
23 // Author: John W. Eaton. |
|
|
24 // HDF5 support by Steven G. Johnson <stevenj@alum.mit.edu> |
|
|
25 // Matlab v5 support by James R. Van Zandt <jrv@vanzandt.mv.com> |
|
3687
|
26 |
|
604
|
27 #ifdef HAVE_CONFIG_H |
|
1192
|
28 #include <config.h> |
|
604
|
29 #endif |
|
|
30 |
|
1343
|
31 #include <cfloat> |
|
|
32 #include <cstring> |
|
|
33 #include <cctype> |
|
|
34 |
|
4249
|
35 #include <fstream> |
|
3503
|
36 #include <iomanip> |
|
|
37 #include <iostream> |
|
1728
|
38 #include <string> |
|
|
39 |
|
3687
|
40 #ifdef HAVE_HDF5 |
|
|
41 #include <hdf5.h> |
|
|
42 #endif |
|
|
43 |
|
1961
|
44 #include "byte-swap.h" |
|
|
45 #include "data-conv.h" |
|
2926
|
46 #include "file-ops.h" |
|
|
47 #include "glob-match.h" |
|
2890
|
48 #include "lo-mappers.h" |
|
4051
|
49 #include "lo-sstream.h" |
|
2318
|
50 #include "mach-info.h" |
|
3185
|
51 #include "oct-env.h" |
|
3258
|
52 #include "oct-time.h" |
|
4171
|
53 #include "quit.h" |
|
1755
|
54 #include "str-vec.h" |
|
|
55 |
|
4332
|
56 #include "Cell.h" |
|
1352
|
57 #include "defun.h" |
|
604
|
58 #include "error.h" |
|
777
|
59 #include "gripes.h" |
|
1352
|
60 #include "load-save.h" |
|
1750
|
61 #include "oct-obj.h" |
|
3687
|
62 #include "oct-map.h" |
|
4332
|
63 #include "ov-cell.h" |
|
1352
|
64 #include "pager.h" |
|
1750
|
65 #include "pt-exp.h" |
|
1352
|
66 #include "symtab.h" |
|
|
67 #include "sysdep.h" |
|
|
68 #include "unwind-prot.h" |
|
604
|
69 #include "utils.h" |
|
2371
|
70 #include "variables.h" |
|
3185
|
71 #include "version.h" |
|
3688
|
72 #include "dMatrix.h" |
|
|
73 |
|
4659
|
74 #ifdef HAVE_HDF5 |
|
4633
|
75 #include "ls-hdf5.h" |
|
4659
|
76 #endif |
|
4633
|
77 #include "ls-mat-ascii.h" |
|
|
78 #include "ls-mat4.h" |
|
|
79 #include "ls-mat5.h" |
|
|
80 #include "ls-oct-ascii.h" |
|
|
81 #include "ls-oct-binary.h" |
|
3688
|
82 |
|
3598
|
83 // Write octave-core file if Octave crashes or is killed by a signal. |
|
3189
|
84 static bool Vcrash_dumps_octave_core; |
|
|
85 |
|
4791
|
86 // The maximum amount of memory (in kilobytes) that we will attempt to |
|
|
87 // write to the Octave core file. |
|
|
88 static double Voctave_core_file_limit; |
|
|
89 |
|
|
90 // The name of the Octave core file. |
|
|
91 static std::string Voctave_core_file_name; |
|
|
92 |
|
3687
|
93 // The default output format. May be one of "binary", "text", |
|
|
94 // "mat-binary", or "hdf5". |
|
3523
|
95 static std::string Vdefault_save_format; |
|
2194
|
96 |
|
4791
|
97 // The output format for Octave core files. May be one of "binary", |
|
4788
|
98 // "text", "mat-binary", or "hdf5". |
|
4791
|
99 static std::string Voctave_core_file_format; |
|
4788
|
100 |
|
3709
|
101 // The format string for the comment line at the top of text-format |
|
|
102 // save files. Passed to strftime. Should begin with `#' and contain |
|
|
103 // no newline characters. |
|
|
104 static std::string Vsave_header_format_string; |
|
|
105 |
|
630
|
106 // XXX FIXME XXX -- shouldn't this be implemented in terms of other |
|
|
107 // functions that are already available? |
|
604
|
108 |
|
|
109 // Install a variable with name NAME and the value specified TC in the |
|
3019
|
110 // symbol table. If FORCE is TRUE, replace any existing definition |
|
|
111 // for NAME. If GLOBAL is TRUE, make the variable global. |
|
604
|
112 // |
|
|
113 // Assumes TC is defined. |
|
|
114 |
|
|
115 static void |
|
4171
|
116 install_loaded_variable (int force, const std::string& name, |
|
|
117 const octave_value& val, |
|
|
118 int global, const std::string& doc) |
|
604
|
119 { |
|
1358
|
120 // Is there already a symbol by this name? If so, what is it? |
|
604
|
121 |
|
2856
|
122 symbol_record *lsr = curr_sym_tab->lookup (name); |
|
604
|
123 |
|
3019
|
124 bool is_undefined = true; |
|
|
125 bool is_variable = false; |
|
|
126 bool is_function = false; |
|
|
127 bool is_global = false; |
|
604
|
128 |
|
|
129 if (lsr) |
|
|
130 { |
|
|
131 is_undefined = ! lsr->is_defined (); |
|
|
132 is_variable = lsr->is_variable (); |
|
|
133 is_function = lsr->is_function (); |
|
|
134 is_global = lsr->is_linked_to_global (); |
|
|
135 } |
|
|
136 |
|
|
137 symbol_record *sr = 0; |
|
|
138 |
|
|
139 if (global) |
|
|
140 { |
|
|
141 if (is_global || is_undefined) |
|
|
142 { |
|
|
143 if (force || is_undefined) |
|
|
144 { |
|
2856
|
145 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
146 link_to_global_variable (lsr); |
|
|
147 sr = lsr; |
|
|
148 } |
|
|
149 else |
|
|
150 { |
|
4171
|
151 warning ("load: global variable name `%s' exists", |
|
|
152 name.c_str ()); |
|
604
|
153 warning ("use `load -force' to overwrite"); |
|
|
154 } |
|
|
155 } |
|
|
156 else if (is_function) |
|
|
157 { |
|
|
158 if (force) |
|
|
159 { |
|
2856
|
160 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
161 link_to_global_variable (lsr); |
|
|
162 sr = lsr; |
|
|
163 } |
|
|
164 else |
|
|
165 { |
|
4171
|
166 warning ("load: `%s' is currently a function in this scope", |
|
|
167 name.c_str ()); |
|
604
|
168 warning ("`load -force' will load variable and hide function"); |
|
|
169 } |
|
|
170 } |
|
|
171 else if (is_variable) |
|
|
172 { |
|
|
173 if (force) |
|
|
174 { |
|
2856
|
175 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
176 link_to_global_variable (lsr); |
|
|
177 sr = lsr; |
|
|
178 } |
|
|
179 else |
|
|
180 { |
|
4171
|
181 warning ("load: local variable name `%s' exists", |
|
|
182 name.c_str ()); |
|
604
|
183 warning ("use `load -force' to overwrite"); |
|
|
184 } |
|
|
185 } |
|
|
186 else |
|
774
|
187 error ("load: unable to load data for unknown symbol type"); |
|
604
|
188 } |
|
|
189 else |
|
|
190 { |
|
|
191 if (is_global) |
|
|
192 { |
|
|
193 if (force || is_undefined) |
|
|
194 { |
|
2856
|
195 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
196 link_to_global_variable (lsr); |
|
|
197 sr = lsr; |
|
|
198 } |
|
|
199 else |
|
|
200 { |
|
4171
|
201 warning ("load: global variable name `%s' exists", |
|
|
202 name.c_str ()); |
|
604
|
203 warning ("use `load -force' to overwrite"); |
|
|
204 } |
|
|
205 } |
|
|
206 else if (is_function) |
|
|
207 { |
|
|
208 if (force) |
|
|
209 { |
|
2856
|
210 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
211 link_to_global_variable (lsr); |
|
|
212 sr = lsr; |
|
|
213 } |
|
|
214 else |
|
|
215 { |
|
4171
|
216 warning ("load: `%s' is currently a function in this scope", |
|
|
217 name.c_str ()); |
|
604
|
218 warning ("`load -force' will load variable and hide function"); |
|
|
219 } |
|
|
220 } |
|
|
221 else if (is_variable || is_undefined) |
|
|
222 { |
|
|
223 if (force || is_undefined) |
|
|
224 { |
|
2856
|
225 lsr = curr_sym_tab->lookup (name, true); |
|
604
|
226 sr = lsr; |
|
|
227 } |
|
|
228 else |
|
|
229 { |
|
4171
|
230 warning ("load: local variable name `%s' exists", |
|
|
231 name.c_str ()); |
|
604
|
232 warning ("use `load -force' to overwrite"); |
|
|
233 } |
|
|
234 } |
|
|
235 else |
|
774
|
236 error ("load: unable to load data for unknown symbol type"); |
|
604
|
237 } |
|
|
238 |
|
|
239 if (sr) |
|
|
240 { |
|
2371
|
241 sr->define (val); |
|
4171
|
242 sr->document (doc); |
|
604
|
243 return; |
|
|
244 } |
|
|
245 else |
|
4171
|
246 error ("load: unable to load variable `%s'", name.c_str ()); |
|
604
|
247 |
|
|
248 return; |
|
|
249 } |
|
|
250 |
|
3019
|
251 // Return TRUE if NAME matches one of the given globbing PATTERNS. |
|
604
|
252 |
|
3013
|
253 static bool |
|
3769
|
254 matches_patterns (const string_vector& patterns, int pat_idx, |
|
3523
|
255 int num_pat, const std::string& name) |
|
604
|
256 { |
|
1755
|
257 for (int i = pat_idx; i < num_pat; i++) |
|
604
|
258 { |
|
1792
|
259 glob_match pattern (patterns[i]); |
|
3013
|
260 |
|
1792
|
261 if (pattern.match (name)) |
|
3013
|
262 return true; |
|
604
|
263 } |
|
3688
|
264 |
|
3013
|
265 return false; |
|
604
|
266 } |
|
|
267 |
|
4329
|
268 int |
|
3523
|
269 read_binary_file_header (std::istream& is, bool& swap, |
|
4329
|
270 oct_mach_info::float_format& flt_fmt, bool quiet) |
|
604
|
271 { |
|
3552
|
272 const int magic_len = 10; |
|
|
273 char magic[magic_len+1]; |
|
3557
|
274 is.read (X_CAST (char *, magic), magic_len); |
|
604
|
275 magic[magic_len] = '\0'; |
|
3688
|
276 |
|
604
|
277 if (strncmp (magic, "Octave-1-L", magic_len) == 0) |
|
2318
|
278 swap = oct_mach_info::words_big_endian (); |
|
604
|
279 else if (strncmp (magic, "Octave-1-B", magic_len) == 0) |
|
2318
|
280 swap = ! oct_mach_info::words_big_endian (); |
|
604
|
281 else |
|
|
282 { |
|
|
283 if (! quiet) |
|
|
284 error ("load: can't read binary file"); |
|
|
285 return -1; |
|
|
286 } |
|
|
287 |
|
|
288 char tmp = 0; |
|
3557
|
289 is.read (X_CAST (char *, &tmp), 1); |
|
604
|
290 |
|
2318
|
291 flt_fmt = mopt_digit_to_float_format (tmp); |
|
|
292 |
|
4574
|
293 if (flt_fmt == oct_mach_info::flt_fmt_unknown) |
|
604
|
294 { |
|
|
295 if (! quiet) |
|
|
296 error ("load: unrecognized binary format!"); |
|
3688
|
297 |
|
604
|
298 return -1; |
|
|
299 } |
|
|
300 |
|
|
301 return 0; |
|
|
302 } |
|
|
303 |
|
|
304 static load_save_format |
|
3523
|
305 get_file_format (const std::string& fname, const std::string& orig_fname) |
|
604
|
306 { |
|
|
307 load_save_format retval = LS_UNKNOWN; |
|
|
308 |
|
4844
|
309 // If the file doesn't exist do nothing |
|
|
310 std::ifstream file_exist (fname.c_str ()); |
|
|
311 if (file_exist) |
|
|
312 file_exist.close (); |
|
|
313 else |
|
4845
|
314 { |
|
|
315 error ("load: nonexistent file `%s'", fname.c_str ()); |
|
|
316 return LS_UNKNOWN; |
|
|
317 } |
|
4844
|
318 |
|
3687
|
319 #ifdef HAVE_HDF5 |
|
3688
|
320 // check this before we open the file |
|
3687
|
321 if (H5Fis_hdf5 (fname.c_str ()) > 0) |
|
|
322 return LS_HDF5; |
|
|
323 #endif /* HAVE_HDF5 */ |
|
|
324 |
|
3523
|
325 std::ifstream file (fname.c_str ()); |
|
604
|
326 |
|
|
327 if (! file) |
|
|
328 { |
|
1750
|
329 error ("load: couldn't open input file `%s'", orig_fname.c_str ()); |
|
604
|
330 return retval; |
|
|
331 } |
|
|
332 |
|
4574
|
333 oct_mach_info::float_format flt_fmt = oct_mach_info::flt_fmt_unknown; |
|
604
|
334 |
|
3019
|
335 bool swap = false; |
|
|
336 |
|
|
337 if (read_binary_file_header (file, swap, flt_fmt, true) == 0) |
|
604
|
338 retval = LS_BINARY; |
|
|
339 else |
|
|
340 { |
|
3538
|
341 file.seekg (0, std::ios::beg); |
|
604
|
342 |
|
|
343 FOUR_BYTE_INT mopt, nr, nc, imag, len; |
|
1180
|
344 |
|
|
345 int err = read_mat_file_header (file, swap, mopt, nr, nc, imag, len, 1); |
|
|
346 |
|
|
347 if (! err) |
|
604
|
348 retval = LS_MAT_BINARY; |
|
|
349 else |
|
|
350 { |
|
2511
|
351 file.clear (); |
|
3538
|
352 file.seekg (0, std::ios::beg); |
|
604
|
353 |
|
3697
|
354 err = read_mat5_binary_file_header (file, swap, true); |
|
3688
|
355 |
|
|
356 if (! err) |
|
|
357 { |
|
|
358 file.clear (); |
|
|
359 file.seekg (0, std::ios::beg); |
|
|
360 retval = LS_MAT5_BINARY; |
|
|
361 } |
|
|
362 else |
|
|
363 { |
|
|
364 file.clear (); |
|
|
365 file.seekg (0, std::ios::beg); |
|
|
366 |
|
4171
|
367 std::string tmp = extract_keyword (file, "name"); |
|
|
368 |
|
|
369 if (! tmp.empty ()) |
|
|
370 retval = LS_ASCII; |
|
3688
|
371 else |
|
|
372 { |
|
|
373 // Try reading the file as numbers only, determining the |
|
|
374 // number of rows and columns from the data. We don't |
|
|
375 // even bother to check to see if the first item in the |
|
|
376 // file is a number, so that get_complete_line() can |
|
|
377 // skip any comments that might appear at the top of the |
|
|
378 // file. |
|
|
379 |
|
|
380 retval = LS_MAT_ASCII; |
|
|
381 } |
|
2511
|
382 } |
|
604
|
383 } |
|
|
384 } |
|
|
385 |
|
|
386 file.close (); |
|
|
387 |
|
|
388 if (retval == LS_UNKNOWN) |
|
1750
|
389 error ("load: unable to determine file format for `%s'", |
|
|
390 orig_fname.c_str ()); |
|
604
|
391 |
|
|
392 return retval; |
|
|
393 } |
|
|
394 |
|
4329
|
395 octave_value |
|
3523
|
396 do_load (std::istream& stream, const std::string& orig_fname, bool force, |
|
2318
|
397 load_save_format format, oct_mach_info::float_format flt_fmt, |
|
4687
|
398 bool list_only, bool swap, bool verbose, |
|
3687
|
399 const string_vector& argv, int argv_idx, int argc, int nargout) |
|
604
|
400 { |
|
3727
|
401 octave_value retval; |
|
|
402 |
|
|
403 Octave_map retstruct; |
|
604
|
404 |
|
4051
|
405 OSSTREAM output_buf; |
|
|
406 |
|
604
|
407 int count = 0; |
|
4051
|
408 |
|
604
|
409 for (;;) |
|
|
410 { |
|
3019
|
411 bool global = false; |
|
2086
|
412 octave_value tc; |
|
604
|
413 |
|
4171
|
414 std::string name; |
|
|
415 std::string doc; |
|
604
|
416 |
|
|
417 switch (format) |
|
|
418 { |
|
|
419 case LS_ASCII: |
|
3136
|
420 name = read_ascii_data (stream, orig_fname, global, tc, count); |
|
604
|
421 break; |
|
|
422 |
|
|
423 case LS_BINARY: |
|
|
424 name = read_binary_data (stream, swap, flt_fmt, orig_fname, |
|
|
425 global, tc, doc); |
|
|
426 break; |
|
|
427 |
|
2511
|
428 case LS_MAT_ASCII: |
|
|
429 name = read_mat_ascii_data (stream, orig_fname, tc); |
|
|
430 break; |
|
|
431 |
|
604
|
432 case LS_MAT_BINARY: |
|
|
433 name = read_mat_binary_data (stream, orig_fname, tc); |
|
|
434 break; |
|
|
435 |
|
3687
|
436 #ifdef HAVE_HDF5 |
|
|
437 case LS_HDF5: |
|
4687
|
438 name = read_hdf5_data (stream, orig_fname, global, tc, doc); |
|
3687
|
439 break; |
|
|
440 #endif /* HAVE_HDF5 */ |
|
|
441 |
|
3688
|
442 case LS_MAT5_BINARY: |
|
|
443 name = read_mat5_binary_element (stream, orig_fname, swap, |
|
|
444 global, tc); |
|
|
445 break; |
|
|
446 |
|
604
|
447 default: |
|
775
|
448 gripe_unrecognized_data_fmt ("load"); |
|
604
|
449 break; |
|
|
450 } |
|
|
451 |
|
4171
|
452 if (error_state || stream.eof () || name.empty ()) |
|
|
453 break; |
|
|
454 else if (! error_state && ! name.empty ()) |
|
604
|
455 { |
|
|
456 if (tc.is_defined ()) |
|
|
457 { |
|
3136
|
458 if (format == LS_MAT_ASCII && argv_idx < argc) |
|
|
459 warning ("load: loaded ASCII file `%s' -- ignoring extra args", |
|
3687
|
460 orig_fname.c_str ()); |
|
3136
|
461 |
|
|
462 if (format == LS_MAT_ASCII |
|
|
463 || argv_idx == argc |
|
1755
|
464 || matches_patterns (argv, argv_idx, argc, name)) |
|
604
|
465 { |
|
|
466 count++; |
|
621
|
467 if (list_only) |
|
|
468 { |
|
|
469 if (verbose) |
|
|
470 { |
|
|
471 if (count == 1) |
|
|
472 output_buf |
|
|
473 << "type rows cols name\n" |
|
|
474 << "==== ==== ==== ====\n"; |
|
|
475 |
|
3013
|
476 output_buf |
|
3548
|
477 << std::setiosflags (std::ios::left) |
|
|
478 << std::setw (16) << tc.type_name () . c_str () |
|
|
479 << std::setiosflags (std::ios::right) |
|
|
480 << std::setw (7) << tc.rows () |
|
|
481 << std::setw (7) << tc.columns () |
|
3013
|
482 << " "; |
|
621
|
483 } |
|
|
484 output_buf << name << "\n"; |
|
|
485 } |
|
|
486 else |
|
|
487 { |
|
3727
|
488 if (nargout == 1) |
|
|
489 { |
|
|
490 if (format == LS_MAT_ASCII) |
|
|
491 retval = tc; |
|
|
492 else |
|
4675
|
493 retstruct.assign (name, tc); |
|
3727
|
494 } |
|
|
495 else |
|
|
496 install_loaded_variable (force, name, tc, global, doc); |
|
621
|
497 } |
|
604
|
498 } |
|
2511
|
499 |
|
|
500 // Only attempt to read one item from a headless text file. |
|
|
501 |
|
|
502 if (format == LS_MAT_ASCII) |
|
|
503 break; |
|
604
|
504 } |
|
|
505 else |
|
4171
|
506 error ("load: unable to load variable `%s'", name.c_str ()); |
|
604
|
507 } |
|
|
508 else |
|
|
509 { |
|
|
510 if (count == 0) |
|
|
511 error ("load: are you sure `%s' is an Octave data file?", |
|
1755
|
512 orig_fname.c_str ()); |
|
604
|
513 |
|
|
514 break; |
|
|
515 } |
|
|
516 } |
|
|
517 |
|
621
|
518 if (list_only && count) |
|
|
519 { |
|
4051
|
520 output_buf << OSSTREAM_ENDS; |
|
|
521 std::string msg = OSSTREAM_STR (output_buf); |
|
|
522 OSSTREAM_FREEZE (output_buf); |
|
2095
|
523 |
|
621
|
524 if (nargout > 0) |
|
2095
|
525 retval = msg; |
|
621
|
526 else |
|
2095
|
527 octave_stdout << msg; |
|
621
|
528 } |
|
3727
|
529 else if (! retstruct.empty ()) |
|
|
530 retval = retstruct; |
|
621
|
531 |
|
863
|
532 return retval; |
|
|
533 } |
|
|
534 |
|
3687
|
535 // HDF5 load/save documentation is included in the Octave manual |
|
|
536 // regardless, but if HDF5 is not linked in we also include a |
|
|
537 // sentence noting this, so the user understands that the features |
|
|
538 // aren't available. Define a macro for this sentence: |
|
|
539 |
|
|
540 #ifdef HAVE_HDF5 |
|
|
541 #define HAVE_HDF5_HELP_STRING "" |
|
|
542 #else /* ! HAVE_HDF5 */ |
|
|
543 #define HAVE_HDF5_HELP_STRING "\n\ |
|
|
544 HDF5 load and save are not available, as this Octave executable was\n\ |
|
|
545 not linked with the HDF5 library." |
|
|
546 #endif /* ! HAVE HDF5 */ |
|
|
547 |
|
4208
|
548 DEFCMD (load, args, nargout, |
|
3372
|
549 "-*- texinfo -*-\n\ |
|
|
550 @deffn {Command} load options file v1 v2 @dots{}\n\ |
|
|
551 Load the named variables from the file @var{file}. As with @code{save},\n\ |
|
|
552 you may specify a list of variables and @code{load} will only extract\n\ |
|
|
553 those variables with names that match. For example, to restore the\n\ |
|
|
554 variables saved in the file @file{data}, use the command\n\ |
|
|
555 \n\ |
|
|
556 @example\n\ |
|
|
557 load data\n\ |
|
|
558 @end example\n\ |
|
863
|
559 \n\ |
|
3372
|
560 If a variable that is not marked as global is loaded from a file when a\n\ |
|
|
561 global symbol with the same name already exists, it is loaded in the\n\ |
|
|
562 global symbol table. Also, if a variable is marked as global in a file\n\ |
|
|
563 and a local symbol exists, the local symbol is moved to the global\n\ |
|
|
564 symbol table and given the value from the file. Since it seems that\n\ |
|
|
565 both of these cases are likely to be the result of some sort of error,\n\ |
|
|
566 they will generate warnings.\n\ |
|
863
|
567 \n\ |
|
3727
|
568 If invoked with a single output argument, Octave returns data instead\n\ |
|
|
569 of inserting variables in the symbol table. If the data file contains\n\ |
|
|
570 only numbers (TAB- or space-delimited columns), a matrix of values is\n\ |
|
|
571 returned. Otherwise, @code{load} returns a structure with members\n\ |
|
|
572 corresponding to the names of the variables in the file.\n\ |
|
|
573 \n\ |
|
3372
|
574 The @code{load} command can read data stored in Octave's text and\n\ |
|
|
575 binary formats, and @sc{Matlab}'s binary format. It will automatically\n\ |
|
|
576 detect the type of file and do conversion from different floating point\n\ |
|
|
577 formats (currently only IEEE big and little endian, though other formats\n\ |
|
|
578 may added in the future).\n\ |
|
|
579 \n\ |
|
|
580 Valid options for @code{load} are listed in the following table.\n\ |
|
863
|
581 \n\ |
|
3372
|
582 @table @code\n\ |
|
|
583 @item -force\n\ |
|
4884
|
584 The @samp{-force} option is accepted but ignored for backward\n\ |
|
|
585 compatiability. Octave now overwrites variables currently in memeory with\n\ |
|
|
586 the same name as those found in the file.\n\ |
|
3372
|
587 \n\ |
|
|
588 @item -ascii\n\ |
|
|
589 Force Octave to assume the file is in Octave's text format.\n\ |
|
|
590 \n\ |
|
|
591 @item -binary\n\ |
|
|
592 Force Octave to assume the file is in Octave's binary format.\n\ |
|
|
593 \n\ |
|
4884
|
594 @item -mat\n\ |
|
|
595 @itemx -mat-binary\n\ |
|
3372
|
596 Force Octave to assume the file is in @sc{Matlab}'s binary format.\n\ |
|
3687
|
597 \n\ |
|
3688
|
598 @item -mat4-binary\n\ |
|
|
599 Force Octave to assume the file is in the binary format written by\n\ |
|
|
600 @sc{Matlab} version 4.\n\ |
|
|
601 \n\ |
|
3687
|
602 @item -hdf5\n\ |
|
|
603 Force Octave to assume the file is in HDF5 format.\n\ |
|
|
604 (HDF5 is a free, portable binary format developed by the National\n\ |
|
|
605 Center for Supercomputing Applications at the University of Illinois.)\n\ |
|
|
606 Note that Octave can read HDF5 files not created by itself, but may\n\ |
|
4687
|
607 skip some datasets in formats that it cannot support.\n" |
|
3687
|
608 |
|
|
609 HAVE_HDF5_HELP_STRING |
|
|
610 |
|
|
611 "\n\ |
|
|
612 @item -import\n\ |
|
4884
|
613 The @samp{-import} is accepted but ignored for backward compatiability.\n\ |
|
|
614 Octave can now support multi-dimensional HDF data and automatically\n\ |
|
|
615 modifies variable names if they are invalid Octave identifiers.\n\ |
|
3687
|
616 \n\ |
|
3372
|
617 @end table\n\ |
|
|
618 @end deffn") |
|
863
|
619 { |
|
2086
|
620 octave_value_list retval; |
|
863
|
621 |
|
1755
|
622 int argc = args.length () + 1; |
|
|
623 |
|
1968
|
624 string_vector argv = args.make_argv ("load"); |
|
1755
|
625 |
|
|
626 if (error_state) |
|
|
627 return retval; |
|
863
|
628 |
|
1358
|
629 // It isn't necessary to have the default load format stored in a |
|
|
630 // user preference variable since we can determine the type of file |
|
|
631 // as we are reading. |
|
863
|
632 |
|
|
633 load_save_format format = LS_UNKNOWN; |
|
|
634 |
|
4691
|
635 bool force = true; |
|
3019
|
636 bool list_only = false; |
|
|
637 bool verbose = false; |
|
863
|
638 |
|
1755
|
639 int i; |
|
|
640 for (i = 1; i < argc; i++) |
|
863
|
641 { |
|
1755
|
642 if (argv[i] == "-force" || argv[i] == "-f") |
|
863
|
643 { |
|
4884
|
644 // Silently ignore this |
|
|
645 // warning ("load: -force ignored"); |
|
863
|
646 } |
|
1755
|
647 else if (argv[i] == "-list" || argv[i] == "-l") |
|
863
|
648 { |
|
3019
|
649 list_only = true; |
|
863
|
650 } |
|
1755
|
651 else if (argv[i] == "-verbose" || argv[i] == "-v") |
|
863
|
652 { |
|
3019
|
653 verbose = true; |
|
863
|
654 } |
|
1755
|
655 else if (argv[i] == "-ascii" || argv[i] == "-a") |
|
863
|
656 { |
|
|
657 format = LS_ASCII; |
|
|
658 } |
|
1755
|
659 else if (argv[i] == "-binary" || argv[i] == "-b") |
|
863
|
660 { |
|
|
661 format = LS_BINARY; |
|
|
662 } |
|
4884
|
663 else if (argv[i] == "-mat-binary" || argv[i] == "-mat" || argv[i] == "-m") |
|
863
|
664 { |
|
3688
|
665 format = LS_MAT5_BINARY; |
|
|
666 } |
|
3797
|
667 else if (argv[i] == "-mat4-binary" || argv[i] == "-4" || argv[i] == "-v4") |
|
3688
|
668 { |
|
863
|
669 format = LS_MAT_BINARY; |
|
|
670 } |
|
3687
|
671 else if (argv[i] == "-hdf5" || argv[i] == "-h") |
|
|
672 { |
|
|
673 #ifdef HAVE_HDF5 |
|
|
674 format = LS_HDF5; |
|
|
675 #else /* ! HAVE_HDF5 */ |
|
|
676 error ("load: octave executable was not linked with HDF5 library"); |
|
|
677 return retval; |
|
|
678 #endif /* ! HAVE_HDF5 */ |
|
|
679 } |
|
|
680 else if (argv[i] == "-import" || argv[i] == "-i") |
|
|
681 { |
|
4687
|
682 warning ("load: -import ignored"); |
|
3687
|
683 } |
|
863
|
684 else |
|
|
685 break; |
|
|
686 } |
|
|
687 |
|
1755
|
688 if (i == argc) |
|
863
|
689 { |
|
2057
|
690 print_usage ("load"); |
|
863
|
691 return retval; |
|
|
692 } |
|
|
693 |
|
3523
|
694 std::string orig_fname = argv[i]; |
|
863
|
695 |
|
4574
|
696 oct_mach_info::float_format flt_fmt = oct_mach_info::flt_fmt_unknown; |
|
863
|
697 |
|
3019
|
698 bool swap = false; |
|
863
|
699 |
|
1755
|
700 if (argv[i] == "-") |
|
863
|
701 { |
|
1755
|
702 i++; |
|
863
|
703 |
|
3687
|
704 #ifdef HAVE_HDF5 |
|
|
705 if (format == LS_HDF5) |
|
|
706 error ("load: cannot read HDF5 format from stdin"); |
|
|
707 else |
|
|
708 #endif /* HAVE_HDF5 */ |
|
863
|
709 if (format != LS_UNKNOWN) |
|
|
710 { |
|
1358
|
711 // XXX FIXME XXX -- if we have already seen EOF on a |
|
3531
|
712 // previous call, how do we fix up the state of std::cin so |
|
|
713 // that we can get additional input? I'm afraid that we |
|
|
714 // can't fix this using std::cin only. |
|
|
715 |
|
|
716 retval = do_load (std::cin, orig_fname, force, format, flt_fmt, |
|
4687
|
717 list_only, swap, verbose, argv, i, argc, |
|
863
|
718 nargout); |
|
|
719 } |
|
|
720 else |
|
|
721 error ("load: must specify file format if reading from stdin"); |
|
|
722 } |
|
|
723 else |
|
|
724 { |
|
3523
|
725 std::string fname = file_ops::tilde_expand (argv[i]); |
|
863
|
726 |
|
|
727 if (format == LS_UNKNOWN) |
|
|
728 format = get_file_format (fname, orig_fname); |
|
|
729 |
|
3687
|
730 #ifdef HAVE_HDF5 |
|
|
731 if (format == LS_HDF5) |
|
|
732 { |
|
|
733 i++; |
|
|
734 |
|
4844
|
735 // If the file doesn't exist do nothing |
|
|
736 std::ifstream file (fname.c_str (), std::ios::in); |
|
|
737 if (file) |
|
3687
|
738 { |
|
4844
|
739 file.close (); |
|
|
740 |
|
|
741 hdf5_ifstream hdf5_file (fname.c_str ()); |
|
3687
|
742 |
|
4844
|
743 if (hdf5_file.file_id >= 0) |
|
|
744 { |
|
|
745 retval = do_load (hdf5_file, orig_fname, force, format, |
|
|
746 flt_fmt, list_only, swap, verbose, |
|
|
747 argv, i, argc, nargout); |
|
|
748 |
|
|
749 hdf5_file.close (); |
|
|
750 } |
|
|
751 else |
|
|
752 error ("load: couldn't open input file `%s'", |
|
|
753 orig_fname.c_str ()); |
|
3687
|
754 } |
|
4845
|
755 else |
|
|
756 error ("load: nonexistent file `%s'", orig_fname.c_str ()); |
|
3687
|
757 } |
|
|
758 else |
|
|
759 #endif /* HAVE_HDF5 */ |
|
|
760 // don't insert any statements here; the "else" above has to |
|
|
761 // go with the "if" below!!!!! |
|
863
|
762 if (format != LS_UNKNOWN) |
|
|
763 { |
|
1755
|
764 i++; |
|
863
|
765 |
|
3775
|
766 std::ios::openmode mode = std::ios::in; |
|
4791
|
767 |
|
|
768 if (format == LS_BINARY |
|
|
769 #ifdef HAVE_HDF5 |
|
|
770 || format == LS_HDF5 |
|
|
771 #endif |
|
|
772 || format == LS_MAT_BINARY |
|
|
773 || format == LS_MAT5_BINARY) |
|
3552
|
774 mode |= std::ios::binary; |
|
863
|
775 |
|
3523
|
776 std::ifstream file (fname.c_str (), mode); |
|
863
|
777 |
|
|
778 if (file) |
|
|
779 { |
|
|
780 if (format == LS_BINARY) |
|
|
781 { |
|
|
782 if (read_binary_file_header (file, swap, flt_fmt) < 0) |
|
|
783 { |
|
|
784 file.close (); |
|
|
785 return retval; |
|
|
786 } |
|
|
787 } |
|
3688
|
788 else if (format == LS_MAT5_BINARY) |
|
|
789 { |
|
|
790 if (read_mat5_binary_file_header (file, swap, false) < 0) |
|
|
791 { |
|
|
792 file.close (); |
|
|
793 return retval; |
|
|
794 } |
|
|
795 } |
|
863
|
796 |
|
|
797 retval = do_load (file, orig_fname, force, format, |
|
4687
|
798 flt_fmt, list_only, swap, verbose, |
|
1755
|
799 argv, i, argc, nargout); |
|
863
|
800 file.close (); |
|
|
801 } |
|
|
802 else |
|
1755
|
803 error ("load: couldn't open input file `%s'", |
|
|
804 orig_fname.c_str ()); |
|
863
|
805 } |
|
4845
|
806 else |
|
|
807 error ("load: nonexistent file: `%s'", orig_fname.c_str ()); |
|
863
|
808 } |
|
604
|
809 |
|
|
810 return retval; |
|
|
811 } |
|
|
812 |
|
3019
|
813 // Return TRUE if PATTERN has any special globbing chars in it. |
|
|
814 |
|
|
815 static bool |
|
3523
|
816 glob_pattern_p (const std::string& pattern) |
|
604
|
817 { |
|
|
818 int open = 0; |
|
|
819 |
|
1755
|
820 int len = pattern.length (); |
|
|
821 |
|
|
822 for (int i = 0; i < len; i++) |
|
604
|
823 { |
|
1755
|
824 char c = pattern[i]; |
|
|
825 |
|
604
|
826 switch (c) |
|
|
827 { |
|
|
828 case '?': |
|
|
829 case '*': |
|
3019
|
830 return true; |
|
604
|
831 |
|
|
832 case '[': // Only accept an open brace if there is a close |
|
|
833 open++; // brace to match it. Bracket expressions must be |
|
|
834 continue; // complete, according to Posix.2 |
|
|
835 |
|
|
836 case ']': |
|
|
837 if (open) |
|
3019
|
838 return true; |
|
604
|
839 continue; |
|
4402
|
840 |
|
604
|
841 case '\\': |
|
1755
|
842 if (i == len - 1) |
|
3019
|
843 return false; |
|
604
|
844 |
|
|
845 default: |
|
|
846 continue; |
|
|
847 } |
|
|
848 } |
|
|
849 |
|
3019
|
850 return false; |
|
604
|
851 } |
|
|
852 |
|
4791
|
853 static void |
|
|
854 do_save (std::ostream& os, const octave_value& tc, |
|
|
855 const std::string& name, const std::string& help, |
|
|
856 int global, load_save_format fmt, bool save_as_floats, |
|
|
857 bool& infnan_warned) |
|
604
|
858 { |
|
|
859 switch (fmt) |
|
|
860 { |
|
|
861 case LS_ASCII: |
|
3738
|
862 save_ascii_data (os, tc, name, infnan_warned, false, global, 0); |
|
604
|
863 break; |
|
|
864 |
|
|
865 case LS_BINARY: |
|
630
|
866 save_binary_data (os, tc, name, help, global, save_as_floats); |
|
604
|
867 break; |
|
|
868 |
|
667
|
869 case LS_MAT_BINARY: |
|
|
870 save_mat_binary_data (os, tc, name); |
|
|
871 break; |
|
|
872 |
|
3687
|
873 #ifdef HAVE_HDF5 |
|
|
874 case LS_HDF5: |
|
|
875 save_hdf5_data (os, tc, name, help, global, save_as_floats); |
|
|
876 break; |
|
|
877 #endif /* HAVE_HDF5 */ |
|
|
878 |
|
3688
|
879 case LS_MAT5_BINARY: |
|
|
880 save_mat5_binary_element (os, tc, name, global, save_as_floats); |
|
|
881 break; |
|
|
882 |
|
604
|
883 default: |
|
775
|
884 gripe_unrecognized_data_fmt ("save"); |
|
604
|
885 break; |
|
|
886 } |
|
|
887 } |
|
|
888 |
|
4791
|
889 // Save the info from SR on stream OS in the format specified by FMT. |
|
|
890 |
|
|
891 void |
|
|
892 do_save (std::ostream& os, symbol_record *sr, load_save_format fmt, |
|
|
893 bool save_as_floats, bool& infnan_warned) |
|
|
894 { |
|
|
895 if (! sr->is_variable ()) |
|
|
896 { |
|
|
897 error ("save: can only save variables, not functions"); |
|
|
898 return; |
|
|
899 } |
|
|
900 |
|
|
901 octave_value tc = sr->def (); |
|
|
902 |
|
|
903 if (tc.is_defined ()) |
|
|
904 { |
|
|
905 std::string name = sr->name (); |
|
|
906 std::string help = sr->help (); |
|
|
907 |
|
|
908 int global = sr->is_linked_to_global (); |
|
|
909 |
|
|
910 do_save (os, tc, name, help, global, fmt, save_as_floats, |
|
|
911 infnan_warned); |
|
|
912 } |
|
|
913 } |
|
|
914 |
|
604
|
915 // Save variables with names matching PATTERN on stream OS in the |
|
3019
|
916 // format specified by FMT. If SAVE_BUILTINS is TRUE, also save |
|
604
|
917 // builtin variables with names that match PATTERN. |
|
|
918 |
|
|
919 static int |
|
3523
|
920 save_vars (std::ostream& os, const std::string& pattern, bool save_builtins, |
|
4791
|
921 load_save_format fmt, bool save_as_floats) |
|
604
|
922 { |
|
3355
|
923 Array<symbol_record *> vars = curr_sym_tab->glob |
|
|
924 (pattern, symbol_record::USER_VARIABLE, SYMTAB_ALL_SCOPES); |
|
|
925 |
|
|
926 int saved = vars.length (); |
|
|
927 |
|
3738
|
928 bool infnan_warned = false; |
|
|
929 |
|
3355
|
930 for (int i = 0; i < saved; i++) |
|
620
|
931 { |
|
4791
|
932 do_save (os, vars(i), fmt, save_as_floats, infnan_warned); |
|
620
|
933 |
|
|
934 if (error_state) |
|
|
935 break; |
|
|
936 } |
|
604
|
937 |
|
620
|
938 if (! error_state && save_builtins) |
|
604
|
939 { |
|
4009
|
940 vars = fbi_sym_tab->glob |
|
3355
|
941 (pattern, symbol_record::BUILTIN_VARIABLE, SYMTAB_ALL_SCOPES); |
|
|
942 |
|
|
943 int count = vars.length (); |
|
604
|
944 |
|
|
945 saved += count; |
|
|
946 |
|
3355
|
947 for (int i = 0; i < count; i++) |
|
620
|
948 { |
|
4791
|
949 do_save (os, vars(i), fmt, save_as_floats, infnan_warned); |
|
620
|
950 |
|
|
951 if (error_state) |
|
|
952 break; |
|
|
953 } |
|
604
|
954 } |
|
|
955 |
|
|
956 return saved; |
|
|
957 } |
|
|
958 |
|
|
959 static load_save_format |
|
4788
|
960 get_save_format (const std::string& fmt, |
|
|
961 load_save_format fallback_format = LS_ASCII) |
|
604
|
962 { |
|
4788
|
963 load_save_format retval = fallback_format; |
|
1755
|
964 |
|
|
965 if (fmt == "binary") |
|
604
|
966 retval = LS_BINARY; |
|
1755
|
967 else if (fmt == "mat-binary" || fmt =="mat_binary") |
|
3688
|
968 retval = LS_MAT5_BINARY; |
|
|
969 else if (fmt == "mat4-binary" || fmt =="mat4_binary") |
|
911
|
970 retval = LS_MAT_BINARY; |
|
3687
|
971 #ifdef HAVE_HDF5 |
|
|
972 else if (fmt == "hdf5") |
|
|
973 retval = LS_HDF5; |
|
|
974 #endif /* HAVE_HDF5 */ |
|
604
|
975 |
|
|
976 return retval; |
|
|
977 } |
|
|
978 |
|
4329
|
979 void |
|
3523
|
980 write_header (std::ostream& os, load_save_format format) |
|
863
|
981 { |
|
3185
|
982 switch (format) |
|
863
|
983 { |
|
3185
|
984 case LS_BINARY: |
|
|
985 { |
|
|
986 os << (oct_mach_info::words_big_endian () |
|
|
987 ? "Octave-1-B" : "Octave-1-L"); |
|
|
988 |
|
|
989 oct_mach_info::float_format flt_fmt = |
|
|
990 oct_mach_info::native_float_format (); |
|
|
991 |
|
|
992 char tmp = (char) float_format_to_mopt_digit (flt_fmt); |
|
|
993 |
|
3557
|
994 os.write (X_CAST (char *, &tmp), 1); |
|
3185
|
995 } |
|
3688
|
996 break; |
|
|
997 |
|
|
998 case LS_MAT5_BINARY: |
|
|
999 { |
|
3775
|
1000 char const * versionmagic; |
|
3688
|
1001 TWO_BYTE_INT number = *(TWO_BYTE_INT *)"\x00\x01"; |
|
|
1002 struct tm bdt; |
|
|
1003 time_t now; |
|
|
1004 char headertext[128]; |
|
|
1005 |
|
|
1006 time (&now); |
|
|
1007 bdt = *gmtime (&now); |
|
|
1008 memset (headertext, ' ', 124); |
|
|
1009 // ISO 8601 format date |
|
|
1010 strftime (headertext, 124, "MATLAB 5.0 MAT-file, written by Octave " |
|
|
1011 OCTAVE_VERSION ", %Y-%m-%d %T UTC", &bdt); |
|
|
1012 |
|
|
1013 // The first pair of bytes give the version of the MAT file |
|
|
1014 // format. The second pair of bytes form a magic number which |
|
|
1015 // signals a MAT file. MAT file data are always written in |
|
|
1016 // native byte order. The order of the bytes in the second |
|
|
1017 // pair indicates whether the file was written by a big- or |
|
|
1018 // little-endian machine. However, the version number is |
|
|
1019 // written in the *opposite* byte order from everything else! |
|
|
1020 if (number == 1) |
|
|
1021 versionmagic = "\x01\x00\x4d\x49"; // this machine is big endian |
|
|
1022 else |
|
|
1023 versionmagic = "\x00\x01\x49\x4d"; // this machine is little endian |
|
|
1024 |
|
|
1025 memcpy (headertext+124, versionmagic, 4); |
|
|
1026 os.write (headertext, 128); |
|
|
1027 } |
|
|
1028 |
|
|
1029 break; |
|
3185
|
1030 |
|
3687
|
1031 #ifdef HAVE_HDF5 |
|
|
1032 case LS_HDF5: |
|
|
1033 #endif /* HAVE_HDF5 */ |
|
3185
|
1034 case LS_ASCII: |
|
|
1035 { |
|
3709
|
1036 octave_localtime now; |
|
|
1037 |
|
|
1038 std::string comment_string = now.strftime (Vsave_header_format_string); |
|
|
1039 |
|
|
1040 if (! comment_string.empty ()) |
|
|
1041 { |
|
3687
|
1042 #ifdef HAVE_HDF5 |
|
3709
|
1043 if (format == LS_HDF5) |
|
|
1044 { |
|
|
1045 hdf5_ofstream& hs = (hdf5_ofstream&) os; |
|
|
1046 H5Gset_comment (hs.file_id, "/", comment_string.c_str ()); |
|
|
1047 } |
|
|
1048 else |
|
3687
|
1049 #endif /* HAVE_HDF5 */ |
|
3709
|
1050 os << comment_string << "\n"; |
|
3687
|
1051 } |
|
3185
|
1052 } |
|
|
1053 break; |
|
|
1054 |
|
|
1055 default: |
|
|
1056 break; |
|
863
|
1057 } |
|
|
1058 } |
|
|
1059 |
|
|
1060 static void |
|
3769
|
1061 save_vars (const string_vector& argv, int argv_idx, int argc, |
|
3523
|
1062 std::ostream& os, bool save_builtins, load_save_format fmt, |
|
3185
|
1063 bool save_as_floats, bool write_header_info) |
|
863
|
1064 { |
|
3185
|
1065 if (write_header_info) |
|
|
1066 write_header (os, fmt); |
|
863
|
1067 |
|
1755
|
1068 if (argv_idx == argc) |
|
863
|
1069 { |
|
|
1070 save_vars (os, "*", save_builtins, fmt, save_as_floats); |
|
|
1071 } |
|
|
1072 else |
|
|
1073 { |
|
1755
|
1074 for (int i = argv_idx; i < argc; i++) |
|
863
|
1075 { |
|
1755
|
1076 if (! save_vars (os, argv[i], save_builtins, fmt, save_as_floats)) |
|
863
|
1077 { |
|
1755
|
1078 warning ("save: no such variable `%s'", argv[i].c_str ()); |
|
863
|
1079 } |
|
|
1080 } |
|
|
1081 } |
|
|
1082 } |
|
|
1083 |
|
1380
|
1084 void |
|
4791
|
1085 dump_octave_core (std::ostream& os, const char *fname, load_save_format fmt) |
|
|
1086 { |
|
|
1087 write_header (os, fmt); |
|
|
1088 |
|
|
1089 Array<symbol_record *> vars = curr_sym_tab->glob |
|
|
1090 ("*", symbol_record::USER_VARIABLE, SYMTAB_ALL_SCOPES); |
|
|
1091 |
|
|
1092 int num_to_save = vars.length (); |
|
|
1093 |
|
|
1094 bool infnan_warned = false; |
|
|
1095 |
|
|
1096 double save_mem_size = 0; |
|
|
1097 |
|
|
1098 for (int i = 0; i < num_to_save; i++) |
|
|
1099 { |
|
|
1100 symbol_record *sr = vars(i); |
|
|
1101 |
|
|
1102 if (sr->is_variable ()) |
|
|
1103 { |
|
|
1104 octave_value tc = sr->def (); |
|
|
1105 |
|
|
1106 if (tc.is_defined ()) |
|
|
1107 { |
|
|
1108 double tc_size = tc.byte_size () / 1024; |
|
|
1109 |
|
|
1110 // XXX FIXME XXX -- maybe we should try to throw out hte |
|
|
1111 // largest first... |
|
|
1112 |
|
|
1113 if (Voctave_core_file_limit < 0 |
|
|
1114 || save_mem_size + tc_size < Voctave_core_file_limit) |
|
|
1115 { |
|
|
1116 save_mem_size += tc_size; |
|
|
1117 |
|
|
1118 std::string name = sr->name (); |
|
|
1119 std::string help = sr->help (); |
|
|
1120 |
|
|
1121 int global = sr->is_linked_to_global (); |
|
|
1122 |
|
|
1123 do_save (os, tc, name, help, global, fmt, false, |
|
|
1124 infnan_warned); |
|
|
1125 |
|
|
1126 if (error_state) |
|
|
1127 break; |
|
|
1128 } |
|
|
1129 } |
|
|
1130 } |
|
|
1131 } |
|
|
1132 |
|
|
1133 message (0, "save to `%s' complete", fname); |
|
|
1134 } |
|
|
1135 |
|
|
1136 void |
|
|
1137 dump_octave_core (void) |
|
1380
|
1138 { |
|
3189
|
1139 if (Vcrash_dumps_octave_core) |
|
1380
|
1140 { |
|
3189
|
1141 // XXX FIXME XXX -- should choose better file name? |
|
|
1142 |
|
4791
|
1143 const char *fname = Voctave_core_file_name.c_str (); |
|
3189
|
1144 |
|
|
1145 message (0, "attempting to save variables to `%s'...", fname); |
|
|
1146 |
|
4788
|
1147 load_save_format format |
|
4791
|
1148 = get_save_format (Voctave_core_file_format, LS_BINARY); |
|
3189
|
1149 |
|
3775
|
1150 std::ios::openmode mode = std::ios::out|std::ios::trunc; |
|
4791
|
1151 |
|
|
1152 if (format == LS_BINARY |
|
|
1153 #ifdef HAVE_HDF5 |
|
|
1154 || format == LS_HDF5 |
|
|
1155 #endif |
|
|
1156 || format == LS_MAT_BINARY |
|
|
1157 || format == LS_MAT5_BINARY) |
|
3552
|
1158 mode |= std::ios::binary; |
|
3189
|
1159 |
|
3687
|
1160 #ifdef HAVE_HDF5 |
|
|
1161 if (format == LS_HDF5) |
|
3189
|
1162 { |
|
3687
|
1163 hdf5_ofstream file (fname); |
|
|
1164 |
|
|
1165 if (file.file_id >= 0) |
|
|
1166 { |
|
4791
|
1167 dump_octave_core (file, fname, format); |
|
3687
|
1168 |
|
|
1169 file.close (); |
|
|
1170 } |
|
|
1171 else |
|
|
1172 warning ("unable to open `%s' for writing...", fname); |
|
3189
|
1173 } |
|
|
1174 else |
|
3687
|
1175 #endif /* HAVE_HDF5 */ |
|
|
1176 // don't insert any commands here! The open brace below must |
|
|
1177 // go with the else above! |
|
|
1178 { |
|
|
1179 std::ofstream file (fname, mode); |
|
|
1180 |
|
|
1181 if (file) |
|
|
1182 { |
|
4791
|
1183 dump_octave_core (file, fname, format); |
|
|
1184 |
|
3687
|
1185 file.close (); |
|
|
1186 } |
|
|
1187 else |
|
|
1188 warning ("unable to open `%s' for writing...", fname); |
|
|
1189 } |
|
1380
|
1190 } |
|
|
1191 } |
|
|
1192 |
|
4208
|
1193 DEFCMD (save, args, , |
|
3372
|
1194 "-*- texinfo -*-\n\ |
|
|
1195 @deffn {Command} save options file v1 v2 @dots{}\n\ |
|
|
1196 Save the named variables @var{v1}, @var{v2}, @dots{} in the file\n\ |
|
|
1197 @var{file}. The special filename @samp{-} can be used to write the\n\ |
|
|
1198 output to your terminal. If no variable names are listed, Octave saves\n\ |
|
|
1199 all the variables in the current scope. Valid options for the\n\ |
|
|
1200 @code{save} command are listed in the following table. Options that\n\ |
|
|
1201 modify the output format override the format specified by the built-in\n\ |
|
|
1202 variable @code{default_save_format}.\n\ |
|
|
1203 \n\ |
|
|
1204 @table @code\n\ |
|
|
1205 @item -ascii\n\ |
|
|
1206 Save the data in Octave's text data format.\n\ |
|
|
1207 \n\ |
|
|
1208 @item -binary\n\ |
|
|
1209 Save the data in Octave's binary data format.\n\ |
|
|
1210 \n\ |
|
|
1211 @item -float-binary\n\ |
|
|
1212 Save the data in Octave's binary data format but only using single\n\ |
|
|
1213 precision. You should use this format only if you know that all the\n\ |
|
|
1214 values to be saved can be represented in single precision.\n\ |
|
|
1215 \n\ |
|
4884
|
1216 @item -mat\n\ |
|
|
1217 @itemx -mat-binary\n\ |
|
3372
|
1218 Save the data in @sc{Matlab}'s binary data format.\n\ |
|
|
1219 \n\ |
|
3688
|
1220 @item -mat4-binary\n\ |
|
|
1221 Save the data in the binary format written by @sc{Matlab} version 4.\n\ |
|
|
1222 \n\ |
|
3687
|
1223 @item -hdf5\n\ |
|
|
1224 Save the data in HDF5 format.\n\ |
|
|
1225 (HDF5 is a free, portable binary format developed by the National\n\ |
|
|
1226 Center for Supercomputing Applications at the University of Illinois.)\n" |
|
|
1227 |
|
|
1228 HAVE_HDF5_HELP_STRING |
|
|
1229 |
|
|
1230 "\n\ |
|
|
1231 @item -float-hdf5\n\ |
|
|
1232 Save the data in HDF5 format but only using single precision.\n\ |
|
|
1233 You should use this format only if you know that all the\n\ |
|
|
1234 values to be saved can be represented in single precision.\n\ |
|
|
1235 \n\ |
|
3372
|
1236 @item -save-builtins\n\ |
|
|
1237 Force Octave to save the values of built-in variables too. By default,\n\ |
|
|
1238 Octave does not save built-in variables.\n\ |
|
|
1239 @end table\n\ |
|
604
|
1240 \n\ |
|
3372
|
1241 The list of variables to save may include wildcard patterns containing\n\ |
|
|
1242 the following special characters:\n\ |
|
|
1243 @table @code\n\ |
|
|
1244 @item ?\n\ |
|
|
1245 Match any single character.\n\ |
|
|
1246 \n\ |
|
|
1247 @item *\n\ |
|
|
1248 Match zero or more characters.\n\ |
|
|
1249 \n\ |
|
|
1250 @item [ @var{list} ]\n\ |
|
|
1251 Match the list of characters specified by @var{list}. If the first\n\ |
|
|
1252 character is @code{!} or @code{^}, match all characters except those\n\ |
|
|
1253 specified by @var{list}. For example, the pattern @samp{[a-zA-Z]} will\n\ |
|
|
1254 match all lower and upper case alphabetic characters. \n\ |
|
|
1255 @end table\n\ |
|
|
1256 \n\ |
|
|
1257 Except when using the @sc{Matlab} binary data file format, saving global\n\ |
|
|
1258 variables also saves the global status of the variable, so that if it is\n\ |
|
|
1259 restored at a later time using @samp{load}, it will be restored as a\n\ |
|
|
1260 global variable.\n\ |
|
|
1261 \n\ |
|
|
1262 The command\n\ |
|
|
1263 \n\ |
|
|
1264 @example\n\ |
|
|
1265 save -binary data a b*\n\ |
|
|
1266 @end example\n\ |
|
|
1267 \n\ |
|
|
1268 @noindent\n\ |
|
|
1269 saves the variable @samp{a} and all variables beginning with @samp{b} to\n\ |
|
|
1270 the file @file{data} in Octave's binary format.\n\ |
|
|
1271 @end deffn") |
|
604
|
1272 { |
|
2086
|
1273 octave_value_list retval; |
|
604
|
1274 |
|
1755
|
1275 int argc = args.length () + 1; |
|
|
1276 |
|
1968
|
1277 string_vector argv = args.make_argv ("save"); |
|
1755
|
1278 |
|
|
1279 if (error_state) |
|
|
1280 return retval; |
|
604
|
1281 |
|
1358
|
1282 // Here is where we would get the default save format if it were |
|
|
1283 // stored in a user preference variable. |
|
604
|
1284 |
|
3019
|
1285 bool save_builtins = false; |
|
|
1286 |
|
|
1287 bool save_as_floats = false; |
|
630
|
1288 |
|
4788
|
1289 load_save_format format = get_save_format (Vdefault_save_format); |
|
604
|
1290 |
|
3185
|
1291 bool append = false; |
|
|
1292 |
|
1755
|
1293 int i; |
|
|
1294 for (i = 1; i < argc; i++) |
|
604
|
1295 { |
|
3185
|
1296 if (argv[i] == "-append") |
|
|
1297 { |
|
|
1298 append = true; |
|
|
1299 } |
|
3465
|
1300 else if (argv[i] == "-ascii" || argv[i] == "-a") |
|
604
|
1301 { |
|
|
1302 format = LS_ASCII; |
|
|
1303 } |
|
1755
|
1304 else if (argv[i] == "-binary" || argv[i] == "-b") |
|
604
|
1305 { |
|
|
1306 format = LS_BINARY; |
|
|
1307 } |
|
3687
|
1308 else if (argv[i] == "-hdf5" || argv[i] == "-h") |
|
|
1309 { |
|
|
1310 #ifdef HAVE_HDF5 |
|
|
1311 format = LS_HDF5; |
|
|
1312 #else /* ! HAVE_HDF5 */ |
|
|
1313 error ("save: octave executable was not linked with HDF5 library"); |
|
|
1314 return retval; |
|
|
1315 #endif /* ! HAVE_HDF5 */ |
|
|
1316 } |
|
4884
|
1317 else if (argv[i] == "-mat-binary" || argv[i] == "-mat" || argv[i] == "-m") |
|
667
|
1318 { |
|
3688
|
1319 format = LS_MAT5_BINARY; |
|
|
1320 } |
|
3797
|
1321 else if (argv[i] == "-mat4-binary" || argv[i] == "-4" || argv[i] == "-v4") |
|
3688
|
1322 { |
|
667
|
1323 format = LS_MAT_BINARY; |
|
|
1324 } |
|
1755
|
1325 else if (argv[i] == "-float-binary" || argv[i] == "-f") |
|
630
|
1326 { |
|
|
1327 format = LS_BINARY; |
|
3019
|
1328 save_as_floats = true; |
|
630
|
1329 } |
|
3687
|
1330 else if (argv[i] == "-float-hdf5") |
|
|
1331 { |
|
|
1332 #ifdef HAVE_HDF5 |
|
|
1333 format = LS_HDF5; |
|
|
1334 save_as_floats = true; |
|
|
1335 #else /* ! HAVE_HDF5 */ |
|
|
1336 error ("save: octave executable was not linked with HDF5 library"); |
|
|
1337 return retval; |
|
|
1338 #endif /* ! HAVE_HDF5 */ |
|
|
1339 } |
|
1755
|
1340 else if (argv[i] == "-save-builtins") |
|
604
|
1341 { |
|
3019
|
1342 save_builtins = true; |
|
604
|
1343 } |
|
|
1344 else |
|
|
1345 break; |
|
|
1346 } |
|
|
1347 |
|
2057
|
1348 if (i == argc) |
|
604
|
1349 { |
|
|
1350 print_usage ("save"); |
|
|
1351 return retval; |
|
|
1352 } |
|
|
1353 |
|
630
|
1354 if (save_as_floats && format == LS_ASCII) |
|
|
1355 { |
|
|
1356 error ("save: cannot specify both -ascii and -float-binary"); |
|
|
1357 return retval; |
|
|
1358 } |
|
|
1359 |
|
1755
|
1360 if (argv[i] == "-") |
|
604
|
1361 { |
|
1755
|
1362 i++; |
|
863
|
1363 |
|
3687
|
1364 #ifdef HAVE_HDF5 |
|
|
1365 if (format == LS_HDF5) |
|
4687
|
1366 error ("save: cannot write HDF5 format to stdout"); |
|
3687
|
1367 else |
|
|
1368 #endif /* HAVE_HDF5 */ |
|
|
1369 // don't insert any commands here! the brace below must go |
|
|
1370 // with the "else" above! |
|
|
1371 { |
|
|
1372 // XXX FIXME XXX -- should things intended for the screen end up |
|
|
1373 // in a octave_value (string)? |
|
|
1374 |
|
|
1375 save_vars (argv, i, argc, octave_stdout, save_builtins, format, |
|
|
1376 save_as_floats, true); |
|
|
1377 } |
|
604
|
1378 } |
|
1755
|
1379 |
|
|
1380 // Guard against things like `save a*', which are probably mistakes... |
|
|
1381 |
|
|
1382 else if (i == argc - 1 && glob_pattern_p (argv[i])) |
|
|
1383 { |
|
|
1384 print_usage ("save"); |
|
604
|
1385 return retval; |
|
|
1386 } |
|
|
1387 else |
|
|
1388 { |
|
3523
|
1389 std::string fname = file_ops::tilde_expand (argv[i]); |
|
1755
|
1390 |
|
|
1391 i++; |
|
604
|
1392 |
|
3775
|
1393 std::ios::openmode mode = std::ios::out; |
|
4791
|
1394 |
|
|
1395 if (format == LS_BINARY |
|
|
1396 #ifdef HAVE_HDF5 |
|
|
1397 || format == LS_HDF5 |
|
|
1398 #endif |
|
|
1399 || format == LS_MAT_BINARY |
|
|
1400 || format == LS_MAT5_BINARY) |
|
3552
|
1401 mode |= std::ios::binary; |
|
3538
|
1402 |
|
|
1403 mode |= append ? std::ios::ate : std::ios::trunc; |
|
3185
|
1404 |
|
3687
|
1405 #ifdef HAVE_HDF5 |
|
|
1406 if (format == LS_HDF5) |
|
863
|
1407 { |
|
3687
|
1408 hdf5_ofstream hdf5_file (fname.c_str ()); |
|
|
1409 |
|
4687
|
1410 if (hdf5_file.file_id >= 0) |
|
|
1411 { |
|
|
1412 save_vars (argv, i, argc, hdf5_file, save_builtins, format, |
|
|
1413 save_as_floats, true); |
|
3687
|
1414 |
|
4687
|
1415 hdf5_file.close (); |
|
3687
|
1416 } |
|
|
1417 else |
|
|
1418 { |
|
|
1419 error ("save: couldn't open output file `%s'", fname.c_str ()); |
|
|
1420 return retval; |
|
|
1421 } |
|
863
|
1422 } |
|
|
1423 else |
|
3687
|
1424 #endif /* HAVE_HDF5 */ |
|
|
1425 // don't insert any statements here! The brace below must go |
|
|
1426 // with the "else" above! |
|
604
|
1427 { |
|
3687
|
1428 std::ofstream file (fname.c_str (), mode); |
|
|
1429 |
|
|
1430 if (file) |
|
|
1431 { |
|
|
1432 bool write_header_info |
|
3775
|
1433 = ((file.rdbuf ())->pubseekoff (0, std::ios::cur) |
|
3769
|
1434 == static_cast<std::streampos> (0)); |
|
3687
|
1435 |
|
|
1436 save_vars (argv, i, argc, file, save_builtins, format, |
|
|
1437 save_as_floats, write_header_info); |
|
|
1438 } |
|
|
1439 else |
|
|
1440 { |
|
|
1441 error ("save: couldn't open output file `%s'", fname.c_str ()); |
|
|
1442 return retval; |
|
|
1443 } |
|
604
|
1444 } |
|
|
1445 } |
|
|
1446 |
|
|
1447 return retval; |
|
|
1448 } |
|
|
1449 |
|
2194
|
1450 static int |
|
3189
|
1451 crash_dumps_octave_core (void) |
|
|
1452 { |
|
|
1453 Vcrash_dumps_octave_core = check_preference ("crash_dumps_octave_core"); |
|
4449
|
1454 |
|
3189
|
1455 return 0; |
|
|
1456 } |
|
|
1457 |
|
|
1458 static int |
|
2194
|
1459 default_save_format (void) |
|
|
1460 { |
|
|
1461 int status = 0; |
|
|
1462 |
|
3523
|
1463 std::string s = builtin_string_variable ("default_save_format"); |
|
2194
|
1464 |
|
|
1465 if (s.empty ()) |
|
|
1466 { |
|
|
1467 gripe_invalid_value_specified ("default_save_format"); |
|
|
1468 status = -1; |
|
|
1469 } |
|
|
1470 else |
|
|
1471 Vdefault_save_format = s; |
|
|
1472 |
|
|
1473 return status; |
|
|
1474 } |
|
|
1475 |
|
4788
|
1476 static int |
|
4791
|
1477 octave_core_file_limit (void) |
|
|
1478 { |
|
|
1479 double val; |
|
|
1480 |
|
|
1481 if (builtin_real_scalar_variable ("octave_core_file_limit", val)) |
|
|
1482 { |
|
|
1483 Voctave_core_file_limit = val; |
|
|
1484 return 0; |
|
|
1485 } |
|
|
1486 else |
|
|
1487 gripe_invalid_value_specified ("octave_core_file_limit"); |
|
|
1488 |
|
|
1489 return -1; |
|
|
1490 } |
|
|
1491 |
|
|
1492 static int |
|
|
1493 octave_core_file_name (void) |
|
4788
|
1494 { |
|
|
1495 int status = 0; |
|
|
1496 |
|
4791
|
1497 std::string s = builtin_string_variable ("octave_core_file_name"); |
|
4788
|
1498 |
|
|
1499 if (s.empty ()) |
|
|
1500 { |
|
4791
|
1501 gripe_invalid_value_specified ("octave_core_file_name"); |
|
4788
|
1502 status = -1; |
|
|
1503 } |
|
|
1504 else |
|
4791
|
1505 Voctave_core_file_name = s; |
|
|
1506 |
|
|
1507 return status; |
|
|
1508 } |
|
|
1509 |
|
|
1510 static int |
|
|
1511 octave_core_file_format (void) |
|
|
1512 { |
|
|
1513 int status = 0; |
|
|
1514 |
|
|
1515 std::string s = builtin_string_variable ("octave_core_file_format"); |
|
|
1516 |
|
|
1517 if (s.empty ()) |
|
|
1518 { |
|
|
1519 gripe_invalid_value_specified ("octave_core_file_format"); |
|
|
1520 status = -1; |
|
|
1521 } |
|
|
1522 else |
|
|
1523 Voctave_core_file_format = s; |
|
4788
|
1524 |
|
|
1525 return status; |
|
|
1526 } |
|
|
1527 |
|
3769
|
1528 static std::string |
|
3709
|
1529 default_save_header_format (void) |
|
|
1530 { |
|
|
1531 return |
|
4633
|
1532 std::string ("# Created by Octave " OCTAVE_VERSION |
|
|
1533 ", %a %b %d %H:%M:%S %Y %Z <") |
|
3709
|
1534 + octave_env::get_user_name () |
|
|
1535 + std::string ("@") |
|
|
1536 + octave_env::get_host_name () |
|
|
1537 + std::string (">"); |
|
|
1538 } |
|
|
1539 |
|
|
1540 static int |
|
|
1541 save_header_format_string (void) |
|
|
1542 { |
|
|
1543 int status = 0; |
|
|
1544 |
|
|
1545 octave_value v = builtin_any_variable ("save_header_format_string"); |
|
|
1546 |
|
|
1547 if (v.is_string ()) |
|
|
1548 Vsave_header_format_string = v.string_value (); |
|
|
1549 else |
|
|
1550 { |
|
|
1551 gripe_invalid_value_specified ("save_header_format_string"); |
|
|
1552 status = -1; |
|
|
1553 } |
|
|
1554 |
|
|
1555 return status; |
|
|
1556 } |
|
|
1557 |
|
2194
|
1558 void |
|
|
1559 symbols_of_load_save (void) |
|
|
1560 { |
|
4233
|
1561 DEFVAR (crash_dumps_octave_core, true, crash_dumps_octave_core, |
|
3372
|
1562 "-*- texinfo -*-\n\ |
|
|
1563 @defvr {Built-in Variable} crash_dumps_octave_core\n\ |
|
|
1564 If this variable is set to a nonzero value, Octave tries to save all\n\ |
|
|
1565 current variables the the file \"octave-core\" if it crashes or receives a\n\ |
|
|
1566 hangup, terminate or similar signal. The default value is 1.\n\ |
|
4791
|
1567 @seealso{octave_core_file_limit, octave_core_file_name, and octave_core_file_format}\n\ |
|
3372
|
1568 @end defvr"); |
|
3189
|
1569 |
|
3258
|
1570 DEFVAR (default_save_format, "ascii", default_save_format, |
|
3372
|
1571 "-*- texinfo -*-\n\ |
|
|
1572 @defvr {Built-in Variable} default_save_format\n\ |
|
|
1573 This variable specifies the default format for the @code{save} command.\n\ |
|
|
1574 It should have one of the following values: @code{\"ascii\"},\n\ |
|
|
1575 @code{\"binary\"}, @code{float-binary}, or @code{\"mat-binary\"}. The\n\ |
|
|
1576 initial default save format is Octave's text format.\n\ |
|
4791
|
1577 @end defvr"); |
|
|
1578 |
|
|
1579 DEFVAR (octave_core_file_limit, -1.0, octave_core_file_limit, |
|
|
1580 "-*- texinfo -*-\n\ |
|
|
1581 @defvr {Built-in Variable} octave_core_file_limit\n\ |
|
|
1582 The maximum amount of memory (in kilobytes) of the top-level workspace\n\ |
|
|
1583 that Octave will attempt to write when saving data to the\n\ |
|
|
1584 @var{octave_core_file_name}. If @var{octave_core_file_format} is a\n\ |
|
|
1585 binary format, then @var{octave_core_file_limit} will be approximately\n\ |
|
|
1586 the maximum size of the file. If a text file format is used, then the\n\ |
|
|
1587 file could be much larger than the limit.\n\ |
|
|
1588 The default value is -1 (unlimited)\n\ |
|
|
1589 @seealso{crash_dumps_octave_core, octave_core_file_name, and octave_core_file_format}\n\ |
|
4788
|
1590 @end defvr"); |
|
|
1591 |
|
4791
|
1592 DEFVAR (octave_core_file_name, "octave-core", octave_core_file_name, |
|
4788
|
1593 "-*- texinfo -*-\n\ |
|
4791
|
1594 @defvr {Built-in Variable} octave_core_file_name\n\ |
|
|
1595 The name of the file used for saving data from the top-level workspace\n\ |
|
|
1596 when Octave aborts. The default value is @code{\"octave-core\"}\n\ |
|
|
1597 @seealso{crash_dumps_octave_core, octave_core_file_name, and octave_core_file_format}\n\ |
|
|
1598 @end defvr"); |
|
|
1599 |
|
|
1600 DEFVAR (octave_core_file_format, "binary", octave_core_file_format, |
|
|
1601 "-*- texinfo -*-\n\ |
|
|
1602 @defvr {Built-in Variable} octave_core_file_format\n\ |
|
4788
|
1603 If Octave aborts, it attempts to save the contents of the top-level\n\ |
|
|
1604 workspace in a file using this format. The value of\n\ |
|
4791
|
1605 @code{octave_core_file_format} should have one of the following values:\n\ |
|
4788
|
1606 @code{\"ascii\"}, @code{\"binary\"}, @code{float-binary}, or\n\ |
|
|
1607 @code{\"mat-binary\"}. The default value is Octave's binary format.\n\ |
|
4791
|
1608 @seealso{crash_dumps_octave_core, octave_core_file_name, and octave_core_file_limit}\n\ |
|
3372
|
1609 @end defvr"); |
|
2194
|
1610 |
|
3709
|
1611 DEFVAR (save_header_format_string, default_save_header_format (), |
|
|
1612 save_header_format_string, |
|
|
1613 "-*- texinfo -*-\n\ |
|
|
1614 @defvr {Built-in Variable} save_header_format_string\n\ |
|
|
1615 This variable specifies the the format string for the comment line\n\ |
|
|
1616 that is written at the beginning of text-format data files saved by\n\ |
|
|
1617 Octave. The format string is passed to @code{strftime} and should\n\ |
|
|
1618 begin with the character @samp{#} and contain no newline characters.\n\ |
|
|
1619 If the value of @code{save_header_format_string} is the empty string,\n\ |
|
|
1620 the header comment is omitted from text-format data files. The\n\ |
|
|
1621 default value is\n\ |
|
|
1622 \n\ |
|
|
1623 @example\n\ |
|
4060
|
1624 \"# Created by Octave VERSION, %a %b %d %H:%M:%S %Y %Z <USER@@HOST>\"\n\ |
|
3709
|
1625 @end example\n\ |
|
|
1626 @seealso{strftime}\n\ |
|
|
1627 @end defvr"); |
|
2194
|
1628 } |
|
|
1629 |
|
604
|
1630 /* |
|
|
1631 ;;; Local Variables: *** |
|
|
1632 ;;; mode: C++ *** |
|
|
1633 ;;; End: *** |
|
|
1634 */ |
