Method dump

src/fields_dump.cpp:108–145

}

void fields::dump(const char *filename, bool single_parallel_file) {
  if (verbosity > 0) {
    printf("creating fields output file \"%s\" (%d)...\n", filename, single_parallel_file);
  }

  h5file file(filename, h5file::WRITE, single_parallel_file, !single_parallel_file);

  // Write out the current time 't'
  size_t dims[1] = {1};
  size_t start[1] = {0};
  size_t _t[1] = {(size_t)t};
  file.create_data("t", 1, dims);
  if (am_master() || !single_parallel_file) file.write_chunk(1, start, dims, _t);

  dump_fields_chunk_field(&file, single_parallel_file, "f",
                          [](fields_chunk *chunk, int c, int d) { return &(chunk->f[c][d]); });
  dump_fields_chunk_field(&file, single_parallel_file, "f_u",
                          [](fields_chunk *chunk, int c, int d) { return &(chunk->f_u[c][d]); });
  dump_fields_chunk_field(&file, single_parallel_file, "f_w",
                          [](fields_chunk *chunk, int c, int d) { return &(chunk->f_w[c][d]); });
  dump_fields_chunk_field(&file, single_parallel_file, "f_cond",
                          [](fields_chunk *chunk, int c, int d) { return &(chunk->f_cond[c][d]); });
  dump_fields_chunk_field(
      &file, single_parallel_file, "f_bfast",
      [](fields_chunk *chunk, int c, int d) { return &(chunk->f_bfast[c][d]); });
  dump_fields_chunk_field(
      &file, single_parallel_file, "f_w_prev",
      [](fields_chunk *chunk, int c, int d) { return &(chunk->f_w_prev[c][d]); });

  // Dump DFT chunks.
  for (int i = 0; i < num_chunks; i++) {
    if (single_parallel_file || chunks[i]->is_mine()) {
      char dataname[1024];
      snprintf(dataname, 1024, "chunk%02d", i);
      save_dft_hdf5(chunks[i]->dft_chunks, dataname, &file, 0, single_parallel_file);
    }
  }
}

void fields::load_fields_chunk_field(h5file *h5f, bool single_parallel_file,
                                     const std::string &field_name,