Method add_srcdata

src/sources.cpp:314–338

}

void fields::add_srcdata(struct sourcedata cur_data, src_time *src, size_t n,
                         std::complex<double> *amp_arr, bool needs_boundary_fix) {
  if (n == 0) {
    n = cur_data.idx_arr.size();
    assert(amp_arr == NULL);
    amp_arr = cur_data.amp_arr.data();
  }
  assert(n == cur_data.idx_arr.size());
  sources = src->add_to(sources, &src);
  std::vector<ptrdiff_t> index_arr(cur_data.idx_arr);
  std::vector<std::complex<double> > amplitudes(amp_arr, amp_arr + n);
  component c = cur_data.near_fd_comp;
  field_type ft = is_H_or_B(c) ? B_stuff : D_stuff;
  if (0 > cur_data.fc_idx or cur_data.fc_idx >= num_chunks)
    meep::abort("fields chunk index out of range");
  fields_chunk *fc = chunks[cur_data.fc_idx];
  if (!fc->is_mine()) meep::abort("wrong fields chunk");

  fc->add_source(ft,
                 src_vol(c, src, std::move(index_arr), std::move(amplitudes), needs_boundary_fix));
  // We can't do require_component(c) since that only works if all processes are adding
  // srcdata for the same components in the same order, which may not be true.
  // ... instead, the caller should call fields::require_source_components()
  //     after all add_srcdata calls are complete.
}

static double *amp_func_data_re = NULL;
static double *amp_func_data_im = NULL;