grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System > Class Template Reference

Positivity-preserving adaptive order reconstruction. See fd::reconstruction::positivity_preserving_adaptive_order() for details. The rest mass density, electron fraction, and the pressure are kept positive. Use unlimited 5th order (degree 4 polynomial) reconstruction on the metric variables. More...

#include <PositivityPreservingAdaptiveOrder.hpp>

Classes
struct	Alpha5
struct	Alpha7
struct	Alpha9
struct	AtmosphereTreatment
struct	LowOrderReconstructor
struct	ReconstructRhoTimesTemperature

Public Types
using	options
using	reconstruction_argument_tags
Public Types inherited from grmhd::GhValenciaDivClean::fd::Reconstructor< System >
using	system = System
using	creatable_classes

Public Member Functions
	PositivityPreservingAdaptiveOrderPrim (PositivityPreservingAdaptiveOrderPrim &&)=default
PositivityPreservingAdaptiveOrderPrim &	operator= (PositivityPreservingAdaptiveOrderPrim &&)=default
	PositivityPreservingAdaptiveOrderPrim (const PositivityPreservingAdaptiveOrderPrim &)=default
PositivityPreservingAdaptiveOrderPrim &	operator= (const PositivityPreservingAdaptiveOrderPrim &)=default
	PositivityPreservingAdaptiveOrderPrim (double alpha_5, std::optional< double > alpha_7, std::optional< double > alpha_9, FallbackReconstructorType low_order_reconstructor, ::VariableFixing::FixReconstructedStateToAtmosphere fix_reconstructed_state_to_atmosphere, bool reconstruct_rho_times_temperature, const Options::Context &context={})
	PositivityPreservingAdaptiveOrderPrim (CkMigrateMessage *msg)
	WRAPPED_PUPable_decl_base_template (Reconstructor< System >, PositivityPreservingAdaptiveOrderPrim< System >)
auto	get_clone () const -> std::unique_ptr< Reconstructor< System > > override
bool	supports_adaptive_order () const override
void	pup (PUP::er &p) override
size_t	ghost_zone_size () const override
template<size_t ThermodynamicDim, typename TagsList>
void	reconstruct (gsl::not_null< std::array< Variables< TagsList >, dim > * > vars_on_lower_face, gsl::not_null< std::array< Variables< TagsList >, dim > * > vars_on_upper_face, gsl::not_null< std::optional< std::array< gsl::span< std::uint8_t >, dim > > * > reconstruction_order, const Variables< hydro::grmhd_tags< DataVector > > &volume_prims, const Variables< typename System::variables_tag::type::tags_list > &volume_spacetime_and_cons_vars, const EquationsOfState::EquationOfState< true, ThermodynamicDim > &eos, const Element< dim > &element, const DirectionalIdMap< dim, evolution::dg::subcell::GhostData > &ghost_data, const Mesh< dim > &subcell_mesh, const VariableFixing::FixToAtmosphere< dim > &fix_to_atmosphere) const
template<size_t ThermodynamicDim, typename TagsList>
void	reconstruct_fd_neighbor (gsl::not_null< Variables< TagsList > * > vars_on_face, const Variables< hydro::grmhd_tags< DataVector > > &subcell_volume_prims, const Variables< grmhd::GhValenciaDivClean::Tags::spacetime_reconstruction_tags > &subcell_volume_spacetime_metric, const EquationsOfState::EquationOfState< true, ThermodynamicDim > &eos, const Element< dim > &element, const DirectionalIdMap< dim, evolution::dg::subcell::GhostData > &ghost_data, const Mesh< dim > &subcell_mesh, const VariableFixing::FixToAtmosphere< dim > &fix_to_atmosphere, Direction< dim > direction_to_reconstruct) const
	Called by an element doing DG when the neighbor is doing subcell.
bool	reconstruct_rho_times_temperature () const override
Public Member Functions inherited from grmhd::GhValenciaDivClean::fd::Reconstructor< System >
	Reconstructor (const Reconstructor &)=default
Reconstructor &	operator= (const Reconstructor &)=default
	Reconstructor (Reconstructor &&)=default
Reconstructor &	operator= (Reconstructor &&)=default
void	pup (PUP::er &p) override

Static Public Attributes
static constexpr size_t	dim = 3
static constexpr Options::String	help
static constexpr bool	use_adaptive_order = true

Friends
template<typename LocalSystem>
bool	operator== (const PositivityPreservingAdaptiveOrderPrim< LocalSystem > &lhs, const PositivityPreservingAdaptiveOrderPrim< LocalSystem > &rhs)
template<typename LocalSystem>
bool	operator!= (const PositivityPreservingAdaptiveOrderPrim< LocalSystem > &lhs, const PositivityPreservingAdaptiveOrderPrim< LocalSystem > &rhs)

Detailed Description

template<typename System>
class grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >

Positivity-preserving adaptive order reconstruction. See fd::reconstruction::positivity_preserving_adaptive_order() for details. The rest mass density, electron fraction, and the pressure are kept positive. Use unlimited 5th order (degree 4 polynomial) reconstruction on the metric variables.

Only the spacetime metric is reconstructed when we and the neighboring element in the direction are doing FD. If we are doing DG and a neighboring element is doing FD, then the spacetime metric, \(\Phi_{iab}\), and \(\Pi_{ab}\) are all reconstructed since the Riemann solver on the DG element also needs to solve for the metric variables.

Member Typedef Documentation

options

template<typename System>

using grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::options

Initial value:

 
      tmpl::list<Alpha5, Alpha7, Alpha9, LowOrderReconstructor,
                 AtmosphereTreatment, ReconstructRhoTimesTemperature>

reconstruction_argument_tags

template<typename System>

using grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::reconstruction_argument_tags

Initial value:

 
      tmpl::list<::Tags::Variables<hydro::grmhd_tags<DataVector>>,
                 typename System::variables_tag,
                 hydro::Tags::GrmhdEquationOfState, domain::Tags::Element<dim>,
                 evolution::dg::subcell::Tags::GhostDataForReconstruction<dim>,
                 evolution::dg::subcell::Tags::Mesh<dim>,
                 ::Tags::VariableFixer<VariableFixing::FixToAtmosphere<dim>>>

Member Function Documentation

get_clone()

template<typename System>

auto grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::get_clone ( ) const -> std::unique_ptr< Reconstructor< System > >

overridevirtual

Implements grmhd::GhValenciaDivClean::fd::Reconstructor< System >.

ghost_zone_size()

template<typename System>

size_t grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::ghost_zone_size ( ) const

inlineoverridevirtual

Implements grmhd::GhValenciaDivClean::fd::Reconstructor< System >.

reconstruct_rho_times_temperature()

template<typename System>

bool grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::reconstruct_rho_times_temperature ( ) const

overridevirtual

Implements grmhd::GhValenciaDivClean::fd::Reconstructor< System >.

supports_adaptive_order()

template<typename System>

bool grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::supports_adaptive_order ( ) const

inlineoverridevirtual

Reimplemented from grmhd::GhValenciaDivClean::fd::Reconstructor< System >.

Member Data Documentation

help

template<typename System>

Options::String grmhd::GhValenciaDivClean::fd::PositivityPreservingAdaptiveOrderPrim< System >::help

staticconstexpr

Initial value:

{
      "Positivity-preserving adaptive-order reconstruction."}

---

The documentation for this class was generated from the following file:

• src/Evolution/Systems/GrMhd/GhValenciaDivClean/FiniteDifference/PositivityPreservingAdaptiveOrder.hpp