Method construct_old_geom_source

Source/sources/Castro_geom.cpp:12–63

///

void
Castro::construct_old_geom_source(MultiFab& source, const MultiFab& state_in, Real time, Real dt)
{

  amrex::ignore_unused(source);
  amrex::ignore_unused(state_in);
  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);

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

  Real mult_factor = 1.0;

  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();
      amrex::Real run_time = ParallelDescriptor::second() - strt_time;

#ifdef BL_LAZY
      Lazy::QueueReduction( [=] () mutable {
#endif
      ParallelDescriptor::ReduceRealMax(run_time,IOProc);

      amrex::Print() << "Castro::construct_old_geom_source() time = " << run_time << "\n" << "\n";
#ifdef BL_LAZY
      });
#endif
    }

#endif
}