Method construct_new_geom_source

Source/sources/Castro_geom.cpp:67–137

void
Castro::construct_new_geom_source(MultiFab& source, const MultiFab& state_old,
                                  const MultiFab& state_new, Real time, Real dt)
{

  amrex::ignore_unused(source);
  amrex::ignore_unused(state_old);
  amrex::ignore_unused(state_new);
  amrex::ignore_unused(time);
  amrex::ignore_unused(dt);

  if (geom.Coord() == 0) {
      return;
  }

  if (use_geom_source == 0) {
      return;
  }

#if AMREX_SPACEDIM <= 2
  const Real strt_time = ParallelDescriptor::second();

  MultiFab geom_src(grids, dmap, source.nComp(), 0);

  geom_src.setVal(0.0);

  // Subtract off the old-time value first
  Real old_time = time - dt;

  if (geom.Coord() == 1) {
      fill_RZ_geom_source(old_time, dt, state_old, geom_src);
  } else {
      fill_RTheta_geom_source(old_time, dt, state_old, geom_src);
  }

  Real mult_factor = -0.5;

  MultiFab::Saxpy(source, mult_factor, geom_src, 0, 0, source.nComp(), 0);   // NOLINT(readability-suspicious-call-argument)

  // Time center with the new data

  geom_src.setVal(0.0);

  mult_factor = 0.5;

  if (geom.Coord() == 1) {
      fill_RZ_geom_source(time, dt, state_new, geom_src);
  } else {
      fill_RTheta_geom_source(time, dt, state_new, geom_src);
  }

  MultiFab::Saxpy(source, mult_factor, geom_src, 0, 0, source.nComp(), 0);   // NOLINT(readability-suspicious-call-argument)

  if (verbose > 1)
  {
      const int IOProc   = ParallelDescriptor::IOProcessorNumber();
      Real      run_time = ParallelDescriptor::second() - strt_time;