Method load

src/fields_dump.cpp:232–278

}

void fields::load(const char *filename, bool single_parallel_file) {
  if (verbosity > 0)
    printf("reading fields from file \"%s\" (%d)...\n", filename, single_parallel_file);

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

  // Read in the current time 't'
  int rank;
  size_t dims[1] = {1};
  size_t start[1] = {0};
  size_t _t[1];
  file.read_size("t", &rank, dims, 1);
  if (rank != 1 || dims[0] != 1) meep::abort("time size mismatch in fields::load");
  if (am_master() || !single_parallel_file) file.read_chunk(1, start, dims, _t);

  if (single_parallel_file) {
    file.prevent_deadlock();
    broadcast(0, _t, dims[0]);
  }

  t = static_cast<int>(_t[0]);
  calc_sources(time());

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

  // Load 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);
      load_dft_hdf5(chunks[i]->dft_chunks, dataname, &file, 0, single_parallel_file);
    }
  }
}

} // namespace meep