Method get_source_slice

src/array_slice.cpp:748–779

}

complex<realnum> *fields::get_source_slice(const volume &where, component source_slice_component,
                                           complex<realnum> *slice) {
  size_t dims[3];
  direction dirs[3];
  vec min_max_loc[2];
  int rank = get_array_slice_dimensions(where, dims, dirs, false, false, min_max_loc);
  size_t slice_size = dims[0] * (rank >= 2 ? dims[1] : 1) * (rank == 3 ? dims[2] : 1);

  source_slice_data data;
  data.source_component = source_slice_component;
  data.slice_imin = gv.round_vec(min_max_loc[0]);
  data.slice_imax = gv.round_vec(min_max_loc[1]);
  data.slice = new complex<realnum>[slice_size];
  if (!data.slice) meep::abort("%s:%i: out of memory (%zu)", __FILE__, __LINE__, slice_size);

  loop_in_chunks(get_source_slice_chunkloop, (void *)&data, where, Centered, true, false);

  complex<realnum> *slice_collapsed = collapse_array(data.slice, &rank, dims, dirs, where);
  rank = get_array_slice_dimensions(where, dims, dirs, true, false);
  slice_size = dims[0] * (rank >= 2 ? dims[1] : 1) * (rank == 3 ? dims[2] : 1);

  if (slice) {
    memcpy(slice, slice_collapsed, 2 * slice_size * sizeof(realnum));
    delete[] (complex<realnum> *)slice_collapsed;
  }
  else
    slice = slice_collapsed;

  array_to_all(slice, slice_size);

  return slice;
}

/***************************************************************/
/***************************************************************/