AmrSystemCoupler: multi-species SYSTEM coupler on AMR (milestone 2.3). More...

#include <pops/core/model/coupled_system.hpp>
#include <pops/amr/hierarchy/refinement_ratio.hpp>
#include <pops/core/foundation/types.hpp>
#include <pops/coupling/amr/amr_coupler_mp.hpp>
#include <pops/coupling/base/aux_fill.hpp>
#include <pops/coupling/source/coupled_source.hpp>
#include <pops/coupling/base/elliptic_rhs.hpp>
#include <pops/numerics/elliptic/interface/elliptic_problem.hpp>
#include <pops/numerics/elliptic/interface/elliptic_solver.hpp>
#include <pops/numerics/elliptic/mg/geometric_mg.hpp>
#include <pops/numerics/time/amr/reflux/amr_reflux_mf.hpp>
#include <pops/numerics/time/integrators/implicit_stepper.hpp>
#include <pops/numerics/time/schemes/scheduler.hpp>
#include <pops/mesh/index/box2d.hpp>
#include <pops/mesh/layout/box_array.hpp>
#include <pops/mesh/layout/distribution_mapping.hpp>
#include <pops/mesh/boundary/fill_boundary.hpp>
#include <pops/mesh/execution/for_each.hpp>
#include <pops/mesh/geometry/geometry.hpp>
#include <pops/mesh/storage/multifab.hpp>
#include <pops/mesh/boundary/physical_bc.hpp>
#include <pops/parallel/comm.hpp>
#include <cstddef>
#include <functional>
#include <stdexcept>
#include <type_traits>
#include <utility>
#include <vector>

Include dependency graph for amr_system_coupler.hpp:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	pops::AmrHierarchyLayout
	Single source of truth on the GRID shared by all blocks: per level the BoxArray (boxes AND order), the DistributionMapping (rank per box) and dx/dy. More...

class	pops::AmrSystemCoupler< System, RhsAssembler, Elliptic >
	Multi-species system coupler on AMR. More...

struct	pops::AmrImplicitSourceStepper
	Default implicit callback for AmrSystemCoupler::step: backward-Euler (Newton) on the model source, applied to EACH level of the hierarchy, followed by a fine -> coarse cascade (coverage consistency, cf. More...

Namespaces
namespace	pops

namespace	pops::detail

Typedefs
template<CoupledSystemLike System, class RhsAssembler , class Elliptic = GeometricMG>
using	pops::AmrSystemDriver = AmrSystemCoupler< System, RhsAssembler, Elliptic >

Enumerations
enum class	pops::PoissonCadence { pops::OncePerStep , pops::PerSubstep }
	Re-solve frequency of the Poisson on AMR: OncePerStep (phi solved once per macro-step, frozen during the advance; cheapest); PerSubstep (phi re-solved before each species substep, more faithful for a field-driven transport, more expensive). More...

Functions
bool	pops::detail::same_level_layout (const BoxArray &a_ba, const DistributionMapping &a_dm, Real a_dx, Real a_dy, const BoxArray &b_ba, const DistributionMapping &b_dm, Real b_dx, Real b_dy)

void	pops::detail::same_layout_or_throw (const std::vector< std::vector< AmrLevelMP > > &block_levels)

Variables
template<class >
constexpr bool	pops::detail::amr_always_false_v = false

Detailed Description

AmrSystemCoupler: multi-species SYSTEM coupler on AMR (milestone 2.3).

Carries a CoupledSystem on an AMR hierarchy: each block has ITS OWN level hierarchy, all species SHARE the same AMR grid, the same aux field (phi, grad phi [, B_z, ...]) and the same coarse Poisson. Orchestration: sync_down (per block) -> coarse Poisson f = Sum_s q_s n_s -> coarse aux + injection to the fine levels -> each block advances via advance_amr<Disc_block> (with its scheme and its substeps; implicit/IMEX blocks delegated to a callback). STRONG INVARIANT: all blocks live on EXACTLY the same grid per level (the aux is shared); same_layout_or_throw checks this at the ctor. PoissonCadence chooses the re-solve frequency of phi (OncePerStep frozen vs PerSubstep). Single-block = bit-identical path to history (loops over the other blocks are empty).

include/pops/coupling/system/amr_system_coupler.hpp File Reference¶

Classes

Namespaces

Typedefs

Enumerations

Functions

Variables

Detailed Description