Method construct_old_thermo_source

Source/sources/Castro_thermo.cpp:5–47

using namespace amrex;

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

    amrex::ignore_unused(time);
    amrex::ignore_unused(dt);

#ifndef MHD
  if (!(time_integration_method == SpectralDeferredCorrections)) {
      return;
  }
#endif

  const Real strt_time = ParallelDescriptor::second();

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

  thermo_src.setVal(0.0);

  fill_thermo_source(state_in, thermo_src);

  Real mult_factor = 1.0;

  MultiFab::Saxpy(source, mult_factor, thermo_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;

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

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

}