ryujin 2.1.1 revision ee5cbcbf2346c1299c942d0e1f13b46449973c18
Loading...
Searching...
No Matches
List of all members
ryujin::Skeleton::LimiterView< dim, Number > Class Template Reference
Minimal equation interface

#include <source/skeleton/limiter.h>

Public Types

Typedefs and constexpr constants
using View = HyperbolicSystemView< dim, Number >
 
using ScalarNumber = typename View::ScalarNumber
 
using state_type = typename View::state_type
 
using flux_contribution_type = typename View::flux_contribution_type
 
using PrecomputedVectorView = typename View::PrecomputedVectorView
 

Public Member Functions

Stencil-based computation of bounds

Intended usage:

LimiterView<dim, Number> limiter_view;
for (unsigned int i = n_internal; i < n_owned; ++i) {
// ...
limiter_view.reset(pv, i, U_i, flux_i);
for (unsigned int col_idx = 1; col_idx < row_length; ++col_idx) {
// ...
limiter_view.accumulate(pv, js, U_j, flux_j, scaled_c_ij,
affine_shift);
}
limiter_view.bounds(hd_i);
}
ryujin::Skeleton::LimiterView::bounds
Bounds bounds(const Number hd_i) const
Definition limiter.h:243
ryujin::Skeleton::LimiterView::reset
void reset(const PrecomputedVectorView &, const unsigned int, const state_type &, const flux_contribution_type &)
Definition limiter.h:218
ryujin::Skeleton::LimiterView::accumulate
void accumulate(const PrecomputedVectorView &, const unsigned int *, const state_type &, const flux_contribution_type &, const dealii::Tensor< 1, dim, Number > &, const state_type &)
Definition limiter.h:230
void reset (const PrecomputedVectorView &, const unsigned int, const state_type &, const flux_contribution_type &)
 
void accumulate (const PrecomputedVectorView &, const unsigned int *, const state_type &, const flux_contribution_type &, const dealii::Tensor< 1, dim, Number > &, const state_type &)
 
Bounds bounds (const Number hd_i) const
 
Convex limiter
std::tuple< Number, bool > limit (const Bounds &, const state_type &, const state_type &, const Number=Number(0.), const Number t_max=Number(1.)) const
 

Computation and manipulation of bounds

using Bounds = std::array< Number, n_bounds >
 
static constexpr unsigned int n_bounds = 0
 
 LimiterView (const View &view, const Limiter< ScalarNumber > &limiter)
 
unsigned int iterations () const
 
Bounds projection_bounds_from_state (const PrecomputedVectorView &, const unsigned int, const state_type &) const
 
Bounds combine_bounds (const Bounds &, const Bounds &) const
 
Bounds fully_relax_bounds (const Bounds &, const Number &) const
 

Detailed Description

template<int dim, typename Number>
class ryujin::Skeleton::LimiterView< dim, Number >

A view of the Limiter that makes the interface available for a given dimension dim and choice of number type Number (which can be a scalar float, or double, as well as a VectorizedArray holding packed scalars).

Definition at line 110 of file limiter.h.

Member Typedef Documentation

◆ View

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::View = HyperbolicSystemView<dim, Number>

Definition at line 118 of file limiter.h.

◆ ScalarNumber

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::ScalarNumber = typename View::ScalarNumber

Definition at line 120 of file limiter.h.

◆ state_type

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::state_type = typename View::state_type

Definition at line 122 of file limiter.h.

◆ flux_contribution_type

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::flux_contribution_type = typename View::flux_contribution_type

Definition at line 124 of file limiter.h.

◆ PrecomputedVectorView

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::PrecomputedVectorView = typename View::PrecomputedVectorView

Definition at line 126 of file limiter.h.

◆ Bounds

template<int dim, typename Number >
using ryujin::Skeleton::LimiterView< dim, Number >::Bounds = std::array<Number, n_bounds>

Array type used to store accumulated bounds.

Definition at line 141 of file limiter.h.

Constructor & Destructor Documentation

◆ LimiterView()

template<int dim, typename Number >
ryujin::Skeleton::LimiterView< dim, Number >::LimiterView ( const View &  view,
const Limiter< ScalarNumber > &  limiter 
)
inline

Constructor taking a HyperbolicSystemView and a Limiter object as arguments

Definition at line 147 of file limiter.h.

Member Function Documentation

◆ iterations()

template<int dim, typename Number >
unsigned int ryujin::Skeleton::LimiterView< dim, Number >::iterations ( ) const
inline

Return the number of limiter iterations.

Definition at line 156 of file limiter.h.

◆ projection_bounds_from_state()

template<int dim, typename Number >
Bounds ryujin::Skeleton::LimiterView< dim, Number >::projection_bounds_from_state ( const PrecomputedVectorView &  ,
const unsigned int  ,
const state_type &   
) const
inline

Given a state U_i and an index i return "strict" bounds, i.e., a minimal convex set containing the state.

Definition at line 165 of file limiter.h.

◆ combine_bounds()

template<int dim, typename Number >
Bounds ryujin::Skeleton::LimiterView< dim, Number >::combine_bounds ( const Bounds &  ,
const Bounds &   
) const
inline

Given two bounds bounds_left, bounds_right, this function computes a larger, combined set of bounds that this is a (convex) superset of the two.

Definition at line 177 of file limiter.h.

◆ fully_relax_bounds()

template<int dim, typename Number >
Bounds ryujin::Skeleton::LimiterView< dim, Number >::fully_relax_bounds ( const Bounds &  ,
const Number &   
) const
inline

This function applies a relaxation to a given a (strict) bound bounds using a non dimensionalized measure hd (that should scale as $h^d$, where $h$ is the local mesh size).

Definition at line 188 of file limiter.h.

Referenced by ryujin::Skeleton::LimiterView< dim, Number >::bounds().

◆ reset()

template<int dim, typename Number >
void ryujin::Skeleton::LimiterView< dim, Number >::reset ( const PrecomputedVectorView &  ,
const unsigned int  ,
const state_type &  ,
const flux_contribution_type &   
)
inline

Reset temporary storage

Definition at line 218 of file limiter.h.

◆ accumulate()

template<int dim, typename Number >
void ryujin::Skeleton::LimiterView< dim, Number >::accumulate ( const PrecomputedVectorView &  ,
const unsigned int *  ,
const state_type &  ,
const flux_contribution_type &  ,
const dealii::Tensor< 1, dim, Number > &  ,
const state_type &   
)
inline

When looping over the sparsity row, add the contribution associated with the neighboring state U_j.

Definition at line 230 of file limiter.h.

◆ bounds()

template<int dim, typename Number >
Bounds ryujin::Skeleton::LimiterView< dim, Number >::bounds ( const Number  hd_i) const
inline

Return the computed bounds (with relaxation applied).

Definition at line 243 of file limiter.h.

References ryujin::Skeleton::LimiterView< dim, Number >::fully_relax_bounds().

◆ limit()

template<int dim, typename Number >
std::tuple< Number, bool > ryujin::Skeleton::LimiterView< dim, Number >::limit ( const Bounds &  ,
const state_type &  ,
const state_type &  ,
const Number  = Number(0.),
const Number  t_max = Number(1.) 
) const
inline

Given a state \(\mathbf U\) and an update \(\mathbf P\) this function computes and returns the maximal coefficient \(t\), obeying \(t_{\text{min}} < t < t_{\text{max}}\), such that the selected local minimum principles are obeyed.

Definition at line 262 of file limiter.h.

Member Data Documentation

◆ n_bounds

template<int dim, typename Number >
constexpr unsigned int ryujin::Skeleton::LimiterView< dim, Number >::n_bounds = 0
staticconstexpr

The number of stored entries in the bounds array.

Definition at line 136 of file limiter.h.

The documentation for this class was generated from the following file:
Generated by doxygen 1.9.8