This class extends FEMContext in order to provide extra data required to perform local element residual and Jacobian assembly in the case of a discontinuous Galerkin (DG) discretization. More...

#include <dg_fem_context.h>

Inheritance diagram for libMesh::DGFEMContext:

[legend]

Public Types
enum	AlgebraicType { NONE = 0 , DOFS_ONLY , CURRENT , OLD , OLD_DOFS_ONLY }
	Enum describing what data to use when initializing algebraic structures on each element. More...

typedef const DenseSubVector< Number > &(DiffContext::*	diff_subsolution_getter) (unsigned int) const
	Helper typedef to simplify refactoring.

typedef std::map< constNumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< DenseSubVector< Number > > > >::iterator	localized_vectors_iterator
	Typedef for the localized_vectors iterator.

Public Member Functions
	DGFEMContext (const System &sys)
	Constructor.

virtual	~DGFEMContext ()
	Destructor.

virtual void	side_fe_reinit () override
	Override side_fe_reinit to set a boolean flag so that by default DG terms are assumed to be inactive.

void	neighbor_side_fe_reinit ()
	Initialize neighbor side data needed to assemble DG terms.

const std::vector< dof_id_type > &	get_neighbor_dof_indices () const
	Accessor for neighbor dof indices.

const std::vector< dof_id_type > &	get_neighbor_dof_indices (unsigned int var) const
	Accessor for element dof indices of a particular variable corresponding to the index argument.

const DenseVector< Number > &	get_neighbor_residual () const
	Const accessor for neighbor residual.

DenseVector< Number > &	get_neighbor_residual ()
	Non-const accessor for neighbor residual.

const DenseSubVector< Number > &	get_neighbor_residual (unsigned int var) const
	Const accessor for neighbor residual of a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_neighbor_residual (unsigned int var)
	Non-const accessor for neighbor residual of a particular variable corresponding to the variable index argument.

const DenseMatrix< Number > &	get_elem_elem_jacobian () const
	Const accessor for element-element Jacobian.

DenseMatrix< Number > &	get_elem_elem_jacobian ()
	Non-const accessor for element-element Jacobian.

const DenseSubMatrix< Number > &	get_elem_elem_jacobian (unsigned int var1, unsigned int var2) const
	Const accessor for element-element Jacobian of particular variables corresponding to the variable index arguments.

DenseSubMatrix< Number > &	get_elem_elem_jacobian (unsigned int var1, unsigned int var2)
	Non-const accessor for element-element Jacobian of particular variables corresponding to the variable index arguments.

const DenseMatrix< Number > &	get_elem_neighbor_jacobian () const
	Const accessor for element-neighbor Jacobian.

DenseMatrix< Number > &	get_elem_neighbor_jacobian ()
	Non-const accessor for element -neighborJacobian.

const DenseSubMatrix< Number > &	get_elem_neighbor_jacobian (unsigned int var1, unsigned int var2) const
	Const accessor for element-neighbor Jacobian of particular variables corresponding to the variable index arguments.

DenseSubMatrix< Number > &	get_elem_neighbor_jacobian (unsigned int var1, unsigned int var2)
	Non-const accessor for element-neighbor Jacobian of particular variables corresponding to the variable index arguments.

const DenseMatrix< Number > &	get_neighbor_elem_jacobian () const
	Const accessor for element-neighbor Jacobian.

DenseMatrix< Number > &	get_neighbor_elem_jacobian ()
	Non-const accessor for element Jacobian.

const DenseSubMatrix< Number > &	get_neighbor_elem_jacobian (unsigned int var1, unsigned int var2) const
	Const accessor for neighbor-element Jacobian of particular variables corresponding to the variable index arguments.

DenseSubMatrix< Number > &	get_neighbor_elem_jacobian (unsigned int var1, unsigned int var2)
	Non-const accessor for neighbor-element Jacobian of particular variables corresponding to the variable index arguments.

const DenseMatrix< Number > &	get_neighbor_neighbor_jacobian () const
	Const accessor for element-neighbor Jacobian.

DenseMatrix< Number > &	get_neighbor_neighbor_jacobian ()
	Non-const accessor for element Jacobian.

const DenseSubMatrix< Number > &	get_neighbor_neighbor_jacobian (unsigned int var1, unsigned int var2) const
	Const accessor for neighbor-neighbor Jacobian of particular variables corresponding to the variable index arguments.

DenseSubMatrix< Number > &	get_neighbor_neighbor_jacobian (unsigned int var1, unsigned int var2)
	Non-const accessor for neighbor-neighbor Jacobian of particular variables corresponding to the variable index arguments.

void	set_neighbor (const Elem &neighbor)
	Set the neighbor element which we will use to assemble DG terms.

const Elem &	get_neighbor () const
	Accessor for current neighbor Elem object for assembling DG terms.

bool	dg_terms_are_active () const
	Are the DG terms active, i.e.

template<typename OutputShape >
void	get_neighbor_side_fe (unsigned int var, FEGenericBase< OutputShape > *&fe) const
	Accessor for neighbor edge/face (2D/3D) finite element object for variable var.

void	use_default_quadrature_rules (int extra_quadrature_order=0)
	Use quadrature rules designed to over-integrate a mass matrix, plus `extra_quadrature_order`.

void	use_unweighted_quadrature_rules (int extra_quadrature_order=0)
	Use quadrature rules designed to exactly integrate unweighted undistorted basis functions, plus `extra_quadrature_order`.

void	use_quadrature_rules (QuadratureType qt, int extra_quadrature_order=0)
	Use a specified quadrature rule type `qt` for every element, side, and edge dimension, at the default quadrature order for the hardest FE type plus `extra_quadrature_order`.

bool	has_side_boundary_id (boundary_id_type id) const
	Reports if the boundary id is found on the current side.

void	side_boundary_ids (std::vector< boundary_id_type > &vec_to_fill) const
	As above, but fills in the std::set provided by the user.

Number	interior_value (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	interior_value (unsigned int var, unsigned int qp, OutputType &u) const

Number	side_value (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	side_value (unsigned int var, unsigned int qp, OutputType &u) const

Number	point_value (unsigned int var, const Point &p) const

template<typename OutputType >
void	point_value (unsigned int var, const Point &p, OutputType &u, const Real tolerance=TOLERANCE) const

Gradient	interior_gradient (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	interior_gradient (unsigned int var, unsigned int qp, OutputType &du) const

Gradient	side_gradient (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	side_gradient (unsigned int var, unsigned int qp, OutputType &du) const

Gradient	point_gradient (unsigned int var, const Point &p) const

template<typename OutputType >
void	point_gradient (unsigned int var, const Point &p, OutputType &grad_u, const Real tolerance=TOLERANCE) const

Tensor	interior_hessian (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	interior_hessian (unsigned int var, unsigned int qp, OutputType &d2u) const

Tensor	side_hessian (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	side_hessian (unsigned int var, unsigned int qp, OutputType &d2u) const

Tensor	point_hessian (unsigned int var, const Point &p) const

template<typename OutputType >
void	point_hessian (unsigned int var, const Point &p, OutputType &hess_u, const Real tolerance=TOLERANCE) const

Number	fixed_interior_value (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_interior_value (unsigned int var, unsigned int qp, OutputType &u) const

Number	fixed_side_value (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_side_value (unsigned int var, unsigned int qp, OutputType &u) const

Number	fixed_point_value (unsigned int var, const Point &p) const

template<typename OutputType >
void	fixed_point_value (unsigned int var, const Point &p, OutputType &u, const Real tolerance=TOLERANCE) const

Gradient	fixed_interior_gradient (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_interior_gradient (unsigned int var, unsigned int qp, OutputType &grad_u) const

Gradient	fixed_side_gradient (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_side_gradient (unsigned int var, unsigned int qp, OutputType &grad_u) const

Gradient	fixed_point_gradient (unsigned int var, const Point &p) const

template<typename OutputType >
void	fixed_point_gradient (unsigned int var, const Point &p, OutputType &grad_u, const Real tolerance=TOLERANCE) const

Tensor	fixed_interior_hessian (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_interior_hessian (unsigned int var, unsigned int qp, OutputType &hess_u) const

Tensor	fixed_side_hessian (unsigned int var, unsigned int qp) const

template<typename OutputType >
void	fixed_side_hessian (unsigned int var, unsigned int qp, OutputType &hess_u) const

Tensor	fixed_point_hessian (unsigned int var, const Point &p) const

template<typename OutputType >
void	fixed_point_hessian (unsigned int var, const Point &p, OutputType &hess_u, const Real tolerance=TOLERANCE) const

template<typename OutputShape >
void	get_element_fe (unsigned int var, FEGenericBase< OutputShape > *&fe) const
	Accessor for interior finite element object for variable var for the largest dimension in the mesh.

FEBase *	get_element_fe (unsigned int var) const
	Accessor for interior finite element object for scalar-valued variable var for the largest dimension in the mesh.

template<typename OutputShape >
void	get_element_fe (unsigned int var, FEGenericBase< OutputShape > *&fe, unsigned short dim) const
	Accessor for interior finite element object for variable var for dimension dim.

void	get_element_fe (unsigned int var, FEAbstract *&fe, unsigned short dim) const
	Accessor for interior finite element object for variable var for dimension dim.

FEBase *	get_element_fe (unsigned int var, unsigned short dim) const
	Accessor for interior finite element object for scalar-valued variable var for dimension dim.

template<typename OutputShape >
void	get_side_fe (unsigned int var, FEGenericBase< OutputShape > *&fe) const
	Accessor for edge/face (2D/3D) finite element object for variable var for the largest dimension in the mesh.

FEBase *	get_side_fe (unsigned int var) const
	Accessor for side finite element object for scalar-valued variable var for the largest dimension in the mesh.

template<typename OutputShape >
void	get_side_fe (unsigned int var, FEGenericBase< OutputShape > *&fe, unsigned short dim) const
	Accessor for edge/face (2D/3D) finite element object for variable var for dimension dim.

void	get_side_fe (unsigned int var, FEAbstract *&fe, unsigned short dim) const
	Accessor for edge/face (2D/3D) finite element object for variable var for dimension dim.

FEBase *	get_side_fe (unsigned int var, unsigned short dim) const
	Accessor for side finite element object for scalar-valued variable var for dimension dim.

template<typename OutputShape >
void	get_edge_fe (unsigned int var, FEGenericBase< OutputShape > *&fe) const
	Accessor for edge (3D only!) finite element object for variable var.

void	get_edge_fe (unsigned int var, FEAbstract *&fe) const

FEBase *	get_edge_fe (unsigned int var) const
	Accessor for edge (3D only!) finite element object for scalar-valued variable var.

template<typename OutputType >
void	interior_values (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &interior_values_vector) const
	Fills a vector of values of the _system_vector at the all the quadrature points in the current element interior.

template<typename OutputType >
void	side_values (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &side_values_vector) const
	Fills a vector of values of the _system_vector at the all the quadrature points on the current element side.

template<typename OutputType >
void	interior_gradients (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &interior_gradients_vector) const
	Fills a vector with the gradient of the solution variable `var` at all the quadrature points in the current element interior.

template<typename OutputType >
void	side_gradients (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &side_gradients_vector) const
	Fills a vector with the gradient of the solution variable `var` at all the quadrature points on the current element side.

template<typename OutputType >
void	interior_hessians (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &d2u_vals) const
	Fills a vector of hessians of the _system_vector at the all the quadrature points in the current element interior.

template<typename OutputType >
void	side_hessians (unsigned int var, const NumericVector< Number > &_system_vector, std::vector< OutputType > &d2u_vals) const
	Fills a vector of hessians of the _system_vector at the all the quadrature points on the current element side.

template<typename OutputType >
void	interior_rate (unsigned int var, unsigned int qp, OutputType &u) const

template<typename OutputType >
void	interior_rate_gradient (unsigned int var, unsigned int qp, OutputType &u) const

template<typename OutputType >
void	side_rate (unsigned int var, unsigned int qp, OutputType &u) const

template<typename OutputType >
void	point_rate (unsigned int var, const Point &p, OutputType &u) const

template<typename OutputType >
void	interior_accel (unsigned int var, unsigned int qp, OutputType &u) const

template<typename OutputType >
void	side_accel (unsigned int var, unsigned int qp, OutputType &u) const

template<typename OutputType >
void	point_accel (unsigned int var, const Point &p, OutputType &u) const

template<typename OutputType >
void	interior_curl (unsigned int var, unsigned int qp, OutputType &curl_u) const

template<typename OutputType >
void	point_curl (unsigned int var, const Point &p, OutputType &curl_u, const Real tolerance=TOLERANCE) const

template<typename OutputType >
void	interior_div (unsigned int var, unsigned int qp, OutputType &div_u) const

virtual void	elem_reinit (Real theta) override
	Resets the current time in the context.

virtual void	elem_side_reinit (Real theta) override
	Resets the current time in the context.

virtual void	elem_edge_reinit (Real theta) override
	Resets the current time in the context.

virtual void	nonlocal_reinit (Real theta) override
	Gives derived classes the opportunity to reinitialize data needed for nonlocal calculations at a new point within a timestep.

virtual void	pre_fe_reinit (const System &, const Elem *e)
	Reinitializes local data vectors/matrices on the current geometric element.

virtual void	elem_fe_reinit (const std::vector< Point > *const pts=nullptr)
	Reinitializes interior FE objects on the current geometric element.

virtual void	edge_fe_reinit ()
	Reinitializes edge FE objects on the current geometric element.

const QBase &	get_element_qrule () const
	Accessor for element interior quadrature rule for the dimension of the current _elem.

const QBase &	get_element_qrule (unsigned short dim) const
	Accessor for element interior quadrature rule.

const QBase &	get_side_qrule () const
	Accessor for element side quadrature rule for the dimension of the current _elem.

const QBase &	get_side_qrule (unsigned short dim) const
	Accessor for element side quadrature rule.

const QBase &	get_edge_qrule () const
	Accessor for element edge quadrature rule.

virtual void	set_mesh_system (System *sys)
	Tells the FEMContext that system `sys` contains the isoparametric Lagrangian variables which correspond to the coordinates of mesh nodes, in problems where the mesh itself is expected to move in time.

const System *	get_mesh_system () const
	Accessor for moving mesh System.

System *	get_mesh_system ()
	Accessor for moving mesh System.

unsigned int	get_mesh_x_var () const
	Accessor for x-variable of moving mesh System.

void	set_mesh_x_var (unsigned int x_var)
	Accessor for x-variable of moving mesh System.

unsigned int	get_mesh_y_var () const
	Accessor for y-variable of moving mesh System.

void	set_mesh_y_var (unsigned int y_var)
	Accessor for y-variable of moving mesh System.

unsigned int	get_mesh_z_var () const
	Accessor for z-variable of moving mesh System.

void	set_mesh_z_var (unsigned int z_var)
	Accessor for z-variable of moving mesh System.

bool	has_elem () const
	Test for current Elem object.

const Elem &	get_elem () const
	Accessor for current Elem object.

Elem &	get_elem ()
	Accessor for current Elem object.

unsigned char	get_side () const
	Accessor for current side of Elem object.

unsigned char	get_edge () const
	Accessor for current edge of Elem object.

unsigned char	get_dim () const
	Accessor for cached mesh dimension.

unsigned char	get_elem_dim () const

const std::set< unsigned char > &	elem_dimensions () const

void	elem_position_set (Real theta)
	Uses the coordinate data specified by mesh_*_position configuration to set the geometry of `elem` to the value it would take after a fraction `theta` of a timestep.

void	elem_position_get ()
	Uses the geometry of `elem` to set the coordinate data specified by mesh_*_position configuration.

void	set_algebraic_type (const AlgebraicType atype)
	Setting which determines whether to initialize algebraic structures (elem_*) on each element and set their values from current_local_solution.

AlgebraicType	algebraic_type () const

void	set_custom_solution (const NumericVector< Number > *custom_sol)
	Set a NumericVector to be used in place of current_local_solution for calculating elem_solution.

void	set_jacobian_tolerance (Real tol)
	Calls set_jacobian_tolerance() on all the FE objects controlled by this class.

FEType	find_hardest_fe_type ()
	Helper function for creating quadrature rules.

void	attach_quadrature_rules ()
	Helper function for attaching quadrature rules.

const std::vector< unsigned int > *	active_vars () const
	Return a pointer to the vector of active variables being computed for, or a null pointer if all variables in the system are active.

template<typename OutputShape >
FEGenericBase< OutputShape > *	build_new_fe (const FEGenericBase< OutputShape > *fe, const Point &p, const Real tolerance=TOLERANCE, const int get_derivative_level=-1) const
	Helper function to reduce some code duplication in the _point_ methods.

template<>
FEGenericBase< Real > *	cached_fe (const unsigned int elem_dim, const FEType fe_type, const int get_derivative_level) const

template<>
FEGenericBase< RealGradient > *	cached_fe (const unsigned int elem_dim, const FEType fe_type, const int get_derivative_level) const

unsigned int	n_vars () const
	Number of variables in solution.

const System &	get_system () const
	Accessor for associated system.

const DenseVector< Number > &	get_elem_solution () const
	Accessor for element solution.

DenseVector< Number > &	get_elem_solution ()
	Non-const accessor for element solution.

const DenseSubVector< Number > &	get_elem_solution (unsigned int var) const
	Accessor for element solution of a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_elem_solution (unsigned int var)
	Accessor for element solution of a particular variable corresponding to the variable index argument.

const DenseVector< Number > &	get_elem_solution_rate () const
	Accessor for element solution rate of change w.r.t.

DenseVector< Number > &	get_elem_solution_rate ()
	Non-const accessor for element solution rate of change w.r.t.

const DenseSubVector< Number > &	get_elem_solution_rate (unsigned int var) const
	Accessor for element solution rate for a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_elem_solution_rate (unsigned int var)
	Accessor for element solution rate for a particular variable corresponding to the variable index argument.

const DenseVector< Number > &	get_elem_solution_accel () const
	Accessor for element solution accel of change w.r.t.

DenseVector< Number > &	get_elem_solution_accel ()
	Non-const accessor for element solution accel of change w.r.t.

const DenseSubVector< Number > &	get_elem_solution_accel (unsigned int var) const
	Accessor for element solution accel for a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_elem_solution_accel (unsigned int var)
	Accessor for element solution accel for a particular variable corresponding to the variable index argument.

const DenseVector< Number > &	get_elem_fixed_solution () const
	Accessor for element fixed solution.

DenseVector< Number > &	get_elem_fixed_solution ()
	Non-const accessor for element fixed solution.

const DenseSubVector< Number > &	get_elem_fixed_solution (unsigned int var) const
	Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_elem_fixed_solution (unsigned int var)
	Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

const DenseVector< Number > &	get_elem_residual () const
	Const accessor for element residual.

DenseVector< Number > &	get_elem_residual ()
	Non-const accessor for element residual.

const DenseSubVector< Number > &	get_elem_residual (unsigned int var) const
	Const accessor for element residual of a particular variable corresponding to the variable index argument.

DenseSubVector< Number > &	get_elem_residual (unsigned int var)
	Non-const accessor for element residual of a particular variable corresponding to the variable index argument.

const DenseMatrix< Number > &	get_elem_jacobian () const
	Const accessor for element Jacobian.

DenseMatrix< Number > &	get_elem_jacobian ()
	Non-const accessor for element Jacobian.

const DenseSubMatrix< Number > &	get_elem_jacobian (unsigned int var1, unsigned int var2) const
	Const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

DenseSubMatrix< Number > &	get_elem_jacobian (unsigned int var1, unsigned int var2)
	Non-const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

const std::vector< Number > &	get_qois () const
	Const accessor for QoI vector.

std::vector< Number > &	get_qois ()
	Non-const accessor for QoI vector.

const std::vector< DenseVector< Number > > &	get_qoi_derivatives () const
	Const accessor for QoI derivatives.

std::vector< DenseVector< Number > > &	get_qoi_derivatives ()
	Non-const accessor for QoI derivatives.

const DenseSubVector< Number > &	get_qoi_derivatives (std::size_t qoi, unsigned int var) const
	Const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

DenseSubVector< Number > &	get_qoi_derivatives (std::size_t qoi, unsigned int var)
	Non-const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

const std::vector< dof_id_type > &	get_dof_indices () const
	Accessor for element dof indices.

std::vector< dof_id_type > &	get_dof_indices ()
	Non-const accessor for element dof indices.

const std::vector< dof_id_type > &	get_dof_indices (unsigned int var) const
	Accessor for element dof indices of a particular variable corresponding to the index argument.

std::vector< dof_id_type > &	get_dof_indices (unsigned int var)
	Accessor for element dof indices of a particular variable corresponding to the index argument.

unsigned int	n_dof_indices () const
	Total number of dof indices on the element.

unsigned int	n_dof_indices (unsigned int var) const
	Total number of dof indices of the particular variable corresponding to the index argument.

Real	get_system_time () const
	Accessor for the time variable stored in the system class.

Real	get_time () const
	Accessor for the time for which the current nonlinear_solution is defined.

void	set_time (Real time_in)
	Set the time for which the current nonlinear_solution is defined.

Real	get_elem_solution_derivative () const
	The derivative of the current elem_solution w.r.t.

Real	get_elem_solution_rate_derivative () const
	The derivative of the current elem_solution_rate w.r.t.

Real	get_elem_solution_accel_derivative () const
	The derivative of the current elem_solution_accel w.r.t.

Real	get_fixed_solution_derivative () const
	The derivative of the current fixed_elem_solution w.r.t.

bool	is_adjoint () const
	Accessor for querying whether we need to do a primal or adjoint solve.

bool &	is_adjoint ()
	Accessor for setting whether we need to do a primal or adjoint solve.

void	set_deltat_pointer (Real *dt)
	Points the _deltat member of this class at a timestep value stored in the creating System, for example DiffSystem::deltat.

Real	get_deltat_value ()

void	add_localized_vector (NumericVector< Number > &localized_vector, const System &sys)
	Adds a vector to the map of localized vectors.

DenseVector< Number > &	get_localized_vector (const NumericVector< Number > &localized_vector)
	Return a reference to DenseVector localization of localized_vector contained in the _localized_vectors map.

const DenseVector< Number > &	get_localized_vector (const NumericVector< Number > &localized_vector) const
	const accessible version of get_localized_vector function

DenseSubVector< Number > &	get_localized_subvector (const NumericVector< Number > &localized_vector, unsigned int var)
	Return a reference to DenseSubVector localization of localized_vector at variable var contained in the _localized_vectors map.

const DenseSubVector< Number > &	get_localized_subvector (const NumericVector< Number > &localized_vector, unsigned int var) const
	const accessible version of get_localized_subvector function

Public Attributes
unsigned char	side
	Current side for side_* to examine.

unsigned char	edge
	Current edge for edge_* to examine.

Real	time
	For time-dependent problems, this is the time t for which the current nonlinear_solution is defined.

const Real	system_time
	This is the time stored in the System class at the time this context was created, i.e.

Real	elem_solution_derivative
	The derivative of elem_solution with respect to the current nonlinear solution.

Real	elem_solution_rate_derivative
	The derivative of elem_solution_rate with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Real	elem_solution_accel_derivative
	The derivative of elem_solution_accel with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Real	fixed_solution_derivative
	The derivative of elem_fixed_solution with respect to the nonlinear solution, for use by systems constructing jacobians with elem_fixed_solution based methods.

Protected Member Functions
template<typename OutputShape >
FEGenericBase< OutputShape > *	cached_fe (const unsigned int elem_dim, const FEType fe_type, const int get_derivative_level) const

void	set_elem (const Elem *e)
	Helper function to promote accessor usage.

template<typename OutputType , typename FENeeded< OutputType >::value_getter fe_getter, diff_subsolution_getter subsolution_getter>
void	some_value (unsigned int var, unsigned int qp, OutputType &u) const
	Helper function to reduce some code duplication in the *interior_value methods.

template<typename OutputType , typename FENeeded< OutputType >::grad_getter fe_getter, diff_subsolution_getter subsolution_getter>
void	some_gradient (unsigned int var, unsigned int qp, OutputType &u) const
	Helper function to reduce some code duplication in the *interior_gradient methods.

template<typename OutputType , typename FENeeded< OutputType >::hess_getter fe_getter, diff_subsolution_getter subsolution_getter>
void	some_hessian (unsigned int var, unsigned int qp, OutputType &u) const
	Helper function to reduce some code duplication in the *interior_hessian methods.

Protected Attributes
std::unique_ptr< const std::vector< unsigned int > >	_active_vars
	Variables on which to enable calculations, or nullptr if all variables in the System are to be enabled.

System *	_mesh_sys
	System from which to acquire moving mesh information.

unsigned int	_mesh_x_var
	Variables from which to acquire moving mesh information.

unsigned int	_mesh_y_var

unsigned int	_mesh_z_var

AlgebraicType	_atype
	Keep track of what type of algebra reinitialization is to be done.

const NumericVector< Number > *	_custom_solution
	Data with which to do algebra reinitialization.

std::unique_ptr< FEGenericBase< Real > >	_real_fe

std::unique_ptr< FEGenericBase< RealGradient > >	_real_grad_fe

int	_real_fe_derivative_level

int	_real_grad_fe_derivative_level

bool	_real_fe_is_inf

bool	_real_grad_fe_is_inf

std::vector< std::map< FEType, std::unique_ptr< FEAbstract > > >	_element_fe
	Finite element objects for each variable's interior, sides and edges.

std::vector< std::map< FEType, std::unique_ptr< FEAbstract > > >	_side_fe

std::map< FEType, std::unique_ptr< FEAbstract > >	_edge_fe

std::vector< std::vector< FEAbstract * > >	_element_fe_var
	Pointers to the same finite element objects, but indexed by variable number.

std::vector< std::vector< FEAbstract * > >	_side_fe_var

std::vector< FEAbstract * >	_edge_fe_var

const BoundaryInfo &	_boundary_info
	Saved reference to BoundaryInfo on the mesh for this System.

const Elem *	_elem
	Current element for element_* to examine.

unsigned char	_dim
	Cached dimension of largest dimension element in this mesh.

unsigned char	_elem_dim
	Cached dimension of this->_elem.

std::set< unsigned char >	_elem_dims
	Cached dimensions of elements in the mesh, plus dimension 0 if SCALAR variables are in use.

std::vector< std::unique_ptr< QBase > >	_element_qrule
	Quadrature rule for element interior.

std::vector< std::unique_ptr< QBase > >	_side_qrule
	Quadrature rules for element sides The FEM context will try to find a quadrature rule that correctly integrates all variables.

std::unique_ptr< QBase >	_edge_qrule
	Quadrature rules for element edges.

int	_extra_quadrature_order
	The extra quadrature order for this context.

std::map< const NumericVector< Number > *, std::pair< DenseVector< Number >, std::vector< DenseSubVector< Number > > > >	_localized_vectors
	Contains pointers to vectors the user has asked to be localized, keyed with pairs of element localized versions of that vector and per variable views.

const bool	_have_local_matrices
	Whether we have local matrices allocated/initialized.

DenseVector< Number >	_elem_solution
	Element by element components of nonlinear_solution as adjusted by a time_solver.

std::vector< DenseSubVector< Number > >	_elem_subsolutions

DenseVector< Number >	_elem_solution_rate
	Element by element components of du/dt as adjusted by a time_solver.

std::vector< DenseSubVector< Number > >	_elem_subsolution_rates

DenseVector< Number >	_elem_solution_accel
	Element by element components of du/dt as adjusted by a time_solver.

std::vector< DenseSubVector< Number > >	_elem_subsolution_accels

DenseVector< Number >	_elem_fixed_solution
	Element by element components of nonlinear_solution at a fixed point in a timestep, for optional use by e.g.

std::vector< DenseSubVector< Number > >	_elem_fixed_subsolutions

DenseVector< Number >	_elem_residual
	Element residual vector.

DenseMatrix< Number >	_elem_jacobian
	Element jacobian: derivatives of elem_residual with respect to elem_solution.

std::vector< Number >	_elem_qoi
	Element quantity of interest contributions.

std::vector< DenseVector< Number > >	_elem_qoi_derivative
	Element quantity of interest derivative contributions.

std::vector< std::vector< DenseSubVector< Number > > >	_elem_qoi_subderivatives

std::vector< DenseSubVector< Number > >	_elem_subresiduals
	Element residual subvectors and (if `_have_local_matrices`) Jacobian submatrices.

std::vector< std::vector< DenseSubMatrix< Number > > >	_elem_subjacobians

std::vector< dof_id_type >	_dof_indices
	Global Degree of freedom index lists.

std::vector< std::vector< dof_id_type > >	_dof_indices_var

Private Member Functions
void	init_internal_data (const System &sys)
	Helper function used in constructors to set up internal data.

void	_do_elem_position_set (Real theta)
	Uses the coordinate data specified by mesh_*_position configuration to set the geometry of `elem` to the value it would take after a fraction `theta` of a timestep.

void	_update_time_from_system (Real theta)
	Update the time in the context object for the given value of theta, based on the values of "time" and "deltat" stored in the system which created this context.

Private Attributes
const Elem *	_neighbor
	Current neighbor element for assembling DG terms.

DenseVector< Number >	_neighbor_residual
	Residual vector of the neighbor component.

DenseMatrix< Number >	_elem_elem_jacobian
	The DG Jacobian terms.

DenseMatrix< Number >	_elem_neighbor_jacobian

DenseMatrix< Number >	_neighbor_elem_jacobian

DenseMatrix< Number >	_neighbor_neighbor_jacobian

std::vector< std::unique_ptr< DenseSubVector< Number > > >	_neighbor_subresiduals
	Element residual subvectors and Jacobian submatrices.

std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > >	_elem_elem_subjacobians

std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > >	_elem_neighbor_subjacobians

std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > >	_neighbor_elem_subjacobians

std::vector< std::vector< std::unique_ptr< DenseSubMatrix< Number > > > >	_neighbor_neighbor_subjacobians

std::vector< dof_id_type >	_neighbor_dof_indices
	Global Degree of freedom index lists for the neighbor element.

std::vector< std::vector< dof_id_type > >	_neighbor_dof_indices_var

std::map< FEType, std::unique_ptr< FEAbstract > >	_neighbor_side_fe
	Finite element objects for each variable's sides on the neighbor element.

std::vector< FEAbstract * >	_neighbor_side_fe_var
	Pointers to the same finite element objects on the neighbor element, but indexed by variable number.

bool	_dg_terms_active
	Boolean flag to indicate whether or not the DG terms have been assembled and should be used in the global matrix assembly.

Real *	_deltat
	Defaults to nullptr, can optionally be used to point to a timestep value in the System-derived class responsible for creating this DiffContext.

const System &	_system
	A reference to the system this context is constructed with.

bool	_is_adjoint
	Is this context to be used for a primal or adjoint solve?

Detailed Description

This class extends FEMContext in order to provide extra data required to perform local element residual and Jacobian assembly in the case of a discontinuous Galerkin (DG) discretization.

Author: David Knezevic

Date: 2015

Extends FEMContext to work for DG problems.

Definition at line 39 of file dg_fem_context.h.

Member Typedef Documentation

◆ diff_subsolution_getter

typedef const DenseSubVector< Number > &(DiffContext::* libMesh::FEMContext::diff_subsolution_getter) (unsigned int) const

inherited

Helper typedef to simplify refactoring.

Definition at line 1121 of file fem_context.h.

◆ localized_vectors_iterator

typedef std::map<constNumericVector<Number>*,std::pair<DenseVector<Number>,std::vector<DenseSubVector<Number>>>>::iterator libMesh::DiffContext::localized_vectors_iterator

inherited

Typedef for the localized_vectors iterator.

Definition at line 540 of file diff_context.h.

Member Enumeration Documentation

◆ AlgebraicType

enum libMesh::FEMContext::AlgebraicType

inherited

Enum describing what data to use when initializing algebraic structures on each element.

Enumerator
NONE
DOFS_ONLY
CURRENT
OLD
OLD_DOFS_ONLY

Definition at line 989 of file fem_context.h.

                     { NONE = 0,  // Do not reinitialize dof_indices
                       DOFS_ONLY, // Reinitialize dof_indices, not
                                  // algebraic structures
                       CURRENT,   // Use dof_indices, current solution
                       OLD,       // Use old_dof_indices, custom solution
                       OLD_DOFS_ONLY}; // Reinitialize old_dof_indices, not

Constructor & Destructor Documentation

◆ DGFEMContext()

libMesh::DGFEMContext::DGFEMContext ( const System & sys )

explicit

Constructor.

Allocates some but fills no data structures.

Definition at line 29 of file dg_fem_context.C.

  : FEMContext(sys),
    _neighbor(nullptr),
    _neighbor_dof_indices_var(sys.n_vars()),
    _dg_terms_active(false)
{
  unsigned int nv = sys.n_vars();
  libmesh_assert (nv);
 
  _neighbor_subresiduals.reserve(nv);
  _elem_elem_subjacobians.resize(nv);
  _elem_neighbor_subjacobians.resize(nv);
  _neighbor_elem_subjacobians.resize(nv);
  _neighbor_neighbor_subjacobians.resize(nv);
 
  for (unsigned int i=0; i != nv; ++i)
    {
      _neighbor_subresiduals.emplace_back(std::make_unique<DenseSubVector<Number>>(_neighbor_residual));
      _elem_elem_subjacobians[i].reserve(nv);
      _elem_neighbor_subjacobians[i].reserve(nv);
      _neighbor_elem_subjacobians[i].reserve(nv);
      _neighbor_neighbor_subjacobians[i].reserve(nv);
 
      for (unsigned int j=0; j != nv; ++j)
        {
          _elem_elem_subjacobians[i].emplace_back(std::make_unique<DenseSubMatrix<Number>>(_elem_elem_jacobian));
          _elem_neighbor_subjacobians[i].emplace_back(std::make_unique<DenseSubMatrix<Number>>(_elem_neighbor_jacobian));
          _neighbor_elem_subjacobians[i].emplace_back(std::make_unique<DenseSubMatrix<Number>>(_neighbor_elem_jacobian));
          _neighbor_neighbor_subjacobians[i].emplace_back(std::make_unique<DenseSubMatrix<Number>>(_neighbor_neighbor_jacobian));
        }
    }
 
  _neighbor_side_fe_var.resize(nv);
  for (unsigned int i=0; i != nv; ++i)
    {
      FEType fe_type = sys.variable_type(i);
 
      if (_neighbor_side_fe[fe_type] == nullptr)
        _neighbor_side_fe[fe_type] = FEAbstract::build(this->_dim, fe_type);
 
      _neighbor_side_fe_var[i] = _neighbor_side_fe[fe_type].get();
    }
}

References libMesh::FEMContext::_dim, _elem_elem_jacobian, _elem_elem_subjacobians, _elem_neighbor_jacobian, _elem_neighbor_subjacobians, _neighbor_elem_jacobian, _neighbor_elem_subjacobians, _neighbor_neighbor_jacobian, _neighbor_neighbor_subjacobians, _neighbor_residual, _neighbor_side_fe, _neighbor_side_fe_var, _neighbor_subresiduals, libMesh::FEAbstract::build(), libMesh::libmesh_assert(), libMesh::System::n_vars(), and libMesh::System::variable_type().

◆ ~DGFEMContext()

libMesh::DGFEMContext::~DGFEMContext ( )

virtualdefault

Destructor.

Member Function Documentation

◆ _do_elem_position_set()

void libMesh::FEMContext::_do_elem_position_set ( Real theta )

privateinherited

Uses the coordinate data specified by mesh_*_position configuration to set the geometry of elem to the value it would take after a fraction theta of a timestep.

This does the work of elem_position_set, but isn't safe to call without _mesh_sys/etc. defined first.

Definition at line 1632 of file fem_context.C.

{
  // This is too expensive to call unless we've been asked to move the mesh
  libmesh_assert (_mesh_sys);
 
  // This will probably break with threading when two contexts are
  // operating on elements which share a node
  libmesh_assert_equal_to (libMesh::n_threads(), 1);
 
  // If the coordinate data is in our own system, it's already
  // been set up for us, and we can ignore our input parameter theta
  //  if (_mesh_sys == this->number())
  //    {
  unsigned int n_nodes = this->get_elem().n_nodes();
 
#ifndef NDEBUG
  const unsigned char dim = this->get_elem_dim();
 
  // For simplicity we demand that mesh coordinates be stored
  // in a format that allows a direct copy
  libmesh_assert(this->get_mesh_x_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_x_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_elem_solution(this->get_mesh_x_var()).size() == n_nodes));
  libmesh_assert(this->get_mesh_y_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_y_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_elem_solution(this->get_mesh_y_var()).size() == n_nodes));
  libmesh_assert(this->get_mesh_z_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_z_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_elem_solution(this->get_mesh_z_var()).size() == n_nodes));
#endif
 
  // Set the new point coordinates
  if (this->get_mesh_x_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      const_cast<Elem &>(this->get_elem()).point(i)(0) =
        libmesh_real(this->get_elem_solution(this->get_mesh_x_var())(i));
 
  if (this->get_mesh_y_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      const_cast<Elem &>(this->get_elem()).point(i)(1) =
        libmesh_real(this->get_elem_solution(this->get_mesh_y_var())(i));
 
  if (this->get_mesh_z_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      const_cast<Elem &>(this->get_elem()).point(i)(2) =
        libmesh_real(this->get_elem_solution(this->get_mesh_z_var())(i));
  //    }
  // FIXME - If the coordinate data is not in our own system, someone
  // had better get around to implementing that... - RHS
  //  else
  //    {
  //      libmesh_not_implemented();
  //    }
}

References libMesh::FEMContext::_mesh_sys, dim, libMesh::FEMContext::get_elem(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::get_element_fe(), libMesh::FEMContext::get_mesh_x_var(), libMesh::FEMContext::get_mesh_y_var(), libMesh::FEMContext::get_mesh_z_var(), libMesh::invalid_uint, libMesh::libmesh_assert(), libMesh::libmesh_real(), libMesh::FEMap::map_fe_type(), libMesh::Elem::n_nodes(), n_nodes, libMesh::n_threads(), and libMesh::Elem::point().

Referenced by libMesh::FEMContext::elem_position_set().

◆ _update_time_from_system()

void libMesh::FEMContext::_update_time_from_system ( Real theta )

privateinherited

Update the time in the context object for the given value of theta, based on the values of "time" and "deltat" stored in the system which created this context.

Definition at line 1948 of file fem_context.C.

{
  // Update the "time" variable based on the value of theta.  For this
  // to work, we need to know the value of deltat, a pointer to which is now
  // stored by our parent DiffContext class.  Note: get_deltat_value() will
  // assert in debug mode if the requested pointer is nullptr.
  const Real deltat = this->get_deltat_value();
 
  this->set_time(theta*(this->get_system_time() + deltat) + (1.-theta)*this->get_system_time());
}

References libMesh::DiffContext::get_deltat_value(), libMesh::DiffContext::get_system_time(), libMesh::Real, and libMesh::DiffContext::set_time().

Referenced by libMesh::FEMContext::elem_edge_reinit(), libMesh::FEMContext::elem_reinit(), libMesh::FEMContext::elem_side_reinit(), and libMesh::FEMContext::nonlocal_reinit().

◆ active_vars()

const std::vector< unsigned int > * libMesh::FEMContext::active_vars ( ) const

inlineinherited

Return a pointer to the vector of active variables being computed for, or a null pointer if all variables in the system are active.

Definition at line 1052 of file fem_context.h.

1052{ return _active_vars.get(); }

libMesh::FEMContext::_active_vars

std::unique_ptr< const std::vector< unsigned int > > _active_vars

Variables on which to enable calculations, or nullptr if all variables in the System are to be enable...

Definition fem_context.h:1072

References libMesh::FEMContext::_active_vars.

Referenced by libMesh::FEMContext::FEMContext(), and libMesh::WrappedFunctor< Output >::init_context().

◆ add_localized_vector()

void libMesh::DiffContext::add_localized_vector	(	NumericVector< Number > &	localized_vector,
		const System &	sys
	)

inherited

Adds a vector to the map of localized vectors.

We can later evaluate interior_values, interior_gradients and side_values for these fields these vectors represent.

Definition at line 125 of file diff_context.C.

{
  // Make an empty pair keyed with a reference to this _localized_vector
  _localized_vectors[&localized_vector] = std::make_pair(DenseVector<Number>(), std::vector<DenseSubVector<Number>>());
 
  unsigned int nv = sys.n_vars();
 
  _localized_vectors[&localized_vector].second.reserve(nv);
 
  // Fill the DenseSubVector with nv copies of DenseVector
  for (unsigned int i=0; i != nv; ++i)
    _localized_vectors[&localized_vector].second.emplace_back(_localized_vectors[&localized_vector].first);
}

References libMesh::DiffContext::_localized_vectors, and libMesh::System::n_vars().

Referenced by SigmaPhysics::init_context().

◆ algebraic_type()

AlgebraicType libMesh::FEMContext::algebraic_type ( ) const

inlineinherited

Definition at line 1010 of file fem_context.h.

1010{ return _atype; }

libMesh::FEMContext::_atype

AlgebraicType _atype

Keep track of what type of algebra reinitialization is to be done.

Definition fem_context.h:1087

References libMesh::FEMContext::_atype.

Referenced by libMesh::FEMContext::build_new_fe(), and libMesh::FEMContext::pre_fe_reinit().

◆ attach_quadrature_rules()

void libMesh::FEMContext::attach_quadrature_rules ( )

inherited

Helper function for attaching quadrature rules.

Definition at line 126 of file fem_context.C.

{
  const System & sys = this->get_system();
 
  auto attach_rules = [this, &sys](unsigned int v)
    {
      for (const auto & dim : _elem_dims)
        {
          FEType fe_type = sys.variable_type(v);
 
          _element_fe[dim][fe_type]->attach_quadrature_rule(_element_qrule[dim].get());
          if (dim)
            _side_fe[dim][fe_type]->attach_quadrature_rule(_side_qrule[dim].get());
          if (dim == 3)
            _edge_fe[fe_type]->attach_quadrature_rule(_edge_qrule.get());
        };
    };
 
  if (_active_vars)
    for (auto v : *_active_vars)
      attach_rules(v);
  else
    for (auto v : make_range(sys.n_vars()))
      attach_rules(v);
}

References libMesh::FEMContext::_active_vars, libMesh::FEMContext::_edge_fe, libMesh::FEMContext::_edge_qrule, libMesh::FEMContext::_elem_dims, libMesh::FEMContext::_element_fe, libMesh::FEMContext::_element_qrule, libMesh::FEMContext::_side_fe, libMesh::FEMContext::_side_qrule, dim, libMesh::DiffContext::get_system(), libMesh::make_range(), libMesh::System::n_vars(), and libMesh::System::variable_type().

Referenced by libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ build_new_fe()

template<typename OutputShape >

template LIBMESH_EXPORT FEGenericBase< RealGradient > * libMesh::FEMContext::build_new_fe	(	const FEGenericBase< OutputShape > *	fe,
		const Point &	p,
		const Real	tolerance = `TOLERANCE`,
		const int	get_derivative_level = `-1`
	)		const

inherited

Helper function to reduce some code duplication in the *_point_* methods.

get_derivative_level should be -1 to get_ everything, 0 to get_phi, 1 to get_dphi, 2 to get_d2phi, or 3 to get_curl_phi

Definition at line 2048 of file fem_context.C.

{
  FEType fe_type = fe->get_fe_type();
 
  // If we don't have an Elem to evaluate on, then the only functions
  // we can sensibly evaluate are the scalar dofs which are the same
  // everywhere.
  libmesh_assert(this->has_elem() || fe_type.family == SCALAR);
 
#ifdef LIBMESH_ENABLE_AMR
  const int add_p_level = fe->add_p_level_in_reinit();
  if ((algebraic_type() == OLD) &&
      this->has_elem())
    {
      if (this->get_elem().p_refinement_flag() == Elem::JUST_REFINED)
        fe_type.order -= add_p_level;
      else if (this->get_elem().p_refinement_flag() == Elem::JUST_COARSENED)
        fe_type.order += add_p_level;
    }
#endif // LIBMESH_ENABLE_AMR
 
  const unsigned int elem_dim = this->has_elem() ? this->get_elem().dim() : 0;
 
  FEGenericBase<OutputShape>* fe_new =
    cached_fe<OutputShape>(elem_dim, fe_type, get_derivative_level);
#ifdef LIBMESH_ENABLE_AMR
  fe_new->add_p_level_in_reinit(add_p_level);
#endif // LIBMESH_ENABLE_AMR
 
  // Map the physical co-ordinates to the master co-ordinates using the inverse_map from fe_interface.h
  // Build a vector of point co-ordinates to send to reinit
  Point master_point = this->has_elem() ?
    FEMap::inverse_map (elem_dim, &this->get_elem(), p, tolerance) :
    Point(0);
 
  std::vector<Point> coor(1, master_point);
 
  switch (get_derivative_level)
    {
    case -1:
      fe_new->get_phi();
      fe_new->get_dphi();
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
      fe_new->get_d2phi();
#endif
      fe_new->get_curl_phi();
      break;
    case 0:
      fe_new->get_phi();
      break;
    case 1:
      fe_new->get_dphi();
      break;
    case 2:
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
      fe_new->get_d2phi();
#else
      // here a different configuration is required.
      libmesh_not_implemented();
#endif
      break;
    case 3:
      fe_new->get_curl_phi();
      break;
    default:
      libmesh_error();
    }
 
  // Reinitialize the element and compute the shape function values at coor
  if (this->has_elem())
    fe_new->reinit (&this->get_elem(), &coor);
  else
    // If !this->has_elem(), then we assume we are dealing with a SCALAR variable
    fe_new->reinit (nullptr, &coor);
 
  return fe_new;
}

References libMesh::FEAbstract::add_p_level_in_reinit(), libMesh::FEMContext::algebraic_type(), libMesh::Elem::dim(), libMesh::FEType::family, libMesh::FEGenericBase< OutputType >::get_curl_phi(), libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem(), libMesh::FEAbstract::get_fe_type(), libMesh::FEGenericBase< OutputType >::get_phi(), libMesh::FEMContext::has_elem(), libMesh::FEMap::inverse_map(), libMesh::Elem::JUST_COARSENED, libMesh::Elem::JUST_REFINED, libMesh::libmesh_assert(), libMesh::FEMContext::OLD, libMesh::FEType::order, libMesh::FEAbstract::reinit(), and libMesh::SCALAR.

Referenced by libMesh::FEMContext::fixed_point_gradient(), libMesh::FEMContext::fixed_point_hessian(), libMesh::FEMContext::fixed_point_value(), libMesh::FEMContext::point_curl(), libMesh::FEMContext::point_gradient(), libMesh::FEMContext::point_hessian(), and libMesh::FEMContext::point_value().

◆ cached_fe() [1/3]

template<typename OutputShape >

FEGenericBase< OutputShape > * libMesh::FEMContext::cached_fe	(	const unsigned int	elem_dim,
		const FEType	fe_type,
		const int	get_derivative_level
	)		const

protectedinherited

◆ cached_fe() [2/3]

template<>

FEGenericBase< Real > * libMesh::FEMContext::cached_fe	(	const unsigned int	elem_dim,
		const FEType	fe_type,
		const int	get_derivative_level
	)		const

inherited

Definition at line 1963 of file fem_context.C.

{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
  const bool fe_needs_inf =
    this->has_elem() && this->get_elem().infinite();
#endif
 
  if (!_real_fe ||
      elem_dim != _real_fe->get_dim() ||
      fe_type != _real_fe->get_fe_type() ||
      get_derivative_level != _real_fe_derivative_level)
    {
      _real_fe_derivative_level = get_derivative_level;
 
      _real_fe =
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
        fe_needs_inf ?
        FEGenericBase<Real>::build_InfFE(elem_dim, fe_type) :
#endif
        FEGenericBase<Real>::build(elem_dim, fe_type);
    }
 
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
  else if (fe_needs_inf && !_real_fe_is_inf)
    _real_fe =
      FEGenericBase<Real>::build_InfFE(elem_dim, fe_type);
  else if (!fe_needs_inf && _real_fe_is_inf)
    _real_fe =
      FEGenericBase<Real>::build(elem_dim, fe_type);
 
  _real_fe_is_inf =
    (this->has_elem() && this->get_elem().infinite());
#endif
 
  return _real_fe.get();
}

References libMesh::FEMContext::_real_fe, libMesh::FEMContext::_real_fe_derivative_level, libMesh::FEMContext::_real_fe_is_inf, libMesh::FEGenericBase< OutputType >::build(), libMesh::FEGenericBase< OutputType >::build_InfFE(), libMesh::FEMContext::get_elem(), libMesh::FEMContext::has_elem(), and libMesh::Elem::infinite().

◆ cached_fe() [3/3]

template<>

FEGenericBase< RealGradient > * libMesh::FEMContext::cached_fe	(	const unsigned int	elem_dim,
		const FEType	fe_type,
		const int	get_derivative_level
	)		const

inherited

Definition at line 2005 of file fem_context.C.

{
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
  const bool fe_needs_inf =
    this->has_elem() && this->get_elem().infinite();
#endif
 
  if (!_real_grad_fe ||
      elem_dim != _real_grad_fe->get_dim() ||
      fe_type != _real_grad_fe->get_fe_type() ||
      get_derivative_level != _real_grad_fe_derivative_level)
    {
      _real_grad_fe_derivative_level = get_derivative_level;
 
      _real_grad_fe =
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
        fe_needs_inf ?
        FEGenericBase<RealGradient>::build_InfFE(elem_dim, fe_type) :
#endif
        FEGenericBase<RealGradient>::build(elem_dim, fe_type);
    }
 
#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
  else if (fe_needs_inf && !_real_grad_fe_is_inf)
    _real_grad_fe =
      FEGenericBase<RealGradient>::build_InfFE(elem_dim, fe_type);
  else if (!fe_needs_inf && _real_grad_fe_is_inf)
    _real_grad_fe =
      FEGenericBase<RealGradient>::build(elem_dim, fe_type);
 
  _real_grad_fe_is_inf =
    (this->has_elem() && this->get_elem().infinite());
#endif
 
  return _real_grad_fe.get();
}

References libMesh::FEMContext::_real_grad_fe, libMesh::FEMContext::_real_grad_fe_derivative_level, libMesh::FEMContext::_real_grad_fe_is_inf, libMesh::FEGenericBase< OutputType >::build(), libMesh::FEGenericBase< OutputType >::build_InfFE(), libMesh::FEMContext::get_elem(), libMesh::FEMContext::has_elem(), and libMesh::Elem::infinite().

◆ dg_terms_are_active()

bool libMesh::DGFEMContext::dg_terms_are_active ( ) const

inline

Are the DG terms active, i.e.

have they been assembled?

Definition at line 232 of file dg_fem_context.h.

233 { return _dg_terms_active; }

References _dg_terms_active.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector().

◆ edge_fe_reinit()

void libMesh::FEMContext::edge_fe_reinit ( )

virtualinherited

Reinitializes edge FE objects on the current geometric element.

Definition at line 1539 of file fem_context.C.

{
  libmesh_assert_equal_to (this->get_elem_dim(), 3);
 
  // Initialize all the interior FE objects on elem/edge.
  // Logging of FE::reinit is done in the FE functions
  for (auto & pr : _edge_fe)
    pr.second->edge_reinit(&(this->get_elem()), this->get_edge());
}

References libMesh::FEMContext::_edge_fe, libMesh::FEMContext::get_edge(), and libMesh::FEMContext::get_elem_dim().

Referenced by libMesh::FEMContext::elem_edge_reinit().

◆ elem_dimensions()

const std::set< unsigned char > & libMesh::FEMContext::elem_dimensions ( ) const

inlineinherited

Returns: Set of dimensions of elements present in the mesh at context initialization.

Definition at line 969 of file fem_context.h.

970 { return _elem_dims; }

References libMesh::FEMContext::_elem_dims.

Referenced by libMesh::RBEIMConstruction::enrich_eim_approximation_on_interiors(), libMesh::WrappedFunctor< Output >::init_context(), libMesh::FEMSystem::init_context(), HeatSystem::init_context(), ElasticitySystem::init_context(), libMesh::DiscontinuityMeasure::init_context(), libMesh::LaplacianErrorEstimator::init_context(), libMesh::KellyErrorEstimator::init_context(), HilbertSystem::init_context(), libMesh::VariationalSmootherSystem::init_context(), libMesh::RBParametrizedFunction::preevaluate_parametrized_function_on_mesh(), libMesh::RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_sides(), and libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::SubFunctor().

◆ elem_edge_reinit()

void libMesh::FEMContext::elem_edge_reinit ( Real theta )

overridevirtualinherited

Resets the current time in the context.

Additionally, reinitialize Elem and FE objects if there's a moving mesh present in the system such that the mesh is deformed to its position at \( t_{\theta} \).

Reimplemented from libMesh::DiffContext.

Definition at line 1476 of file fem_context.C.

{
  // Update the "time" variable of this context object
  this->_update_time_from_system(theta);
 
  // Handle a moving element if necessary
  if (_mesh_sys)
    {
      // FIXME - not threadsafe yet!
      elem_position_set(theta);
      edge_fe_reinit();
    }
}

References libMesh::FEMContext::_mesh_sys, libMesh::FEMContext::_update_time_from_system(), libMesh::FEMContext::edge_fe_reinit(), and libMesh::FEMContext::elem_position_set().

◆ elem_fe_reinit()

void libMesh::FEMContext::elem_fe_reinit ( const std::vector< Point > *const pts = nullptr )

virtualinherited

Reinitializes interior FE objects on the current geometric element.

Definition at line 1501 of file fem_context.C.

{
  // Initialize all the interior FE objects on elem.
  // Logging of FE::reinit is done in the FE functions
  // We only reinit the FE objects for the current element
  // dimension
  const unsigned char dim = this->get_elem_dim();
 
  libmesh_assert( !_element_fe[dim].empty() );
 
  for (const auto & pr : _element_fe[dim])
    {
      if (this->has_elem())
        pr.second->reinit(&(this->get_elem()), pts);
        // If !this->has_elem(), then still might need to reinit for a
        // SCALAR variable; everything else will depend on an elem
      else if (pr.first.family == SCALAR)
        pr.second->reinit(nullptr);
    }
}

References libMesh::FEMContext::_element_fe, dim, libMesh::FEMContext::get_elem_dim(), libMesh::FEMContext::has_elem(), libMesh::libmesh_assert(), and libMesh::SCALAR.

Referenced by libMesh::ExactSolution::_compute_error(), libMesh::RBConstruction::add_scaled_matrix_and_vector(), assembly_with_dg_fem_context(), libMesh::VariationalSmootherSystem::compute_mesh_quality_info(), libMesh::FEMContext::elem_reinit(), HilbertSystem::element_time_derivative(), libMesh::RBEIMConstruction::enrich_eim_approximation_on_sides(), libMesh::RBEIMConstruction::initialize_qp_data(), libMesh::FEMSystem::mesh_position_set(), libMesh::FEMContext::nonlocal_reinit(), Integrate::operator()(), libMesh::VariationalSmootherSystem::prepare_for_smoothing(), OverlappingAlgebraicGhostingTest::run_ghosting_test(), OverlappingCouplingGhostingTest::run_sparsity_pattern_test(), SlitMeshRefinedSystemTest::testRestart(), and SlitMeshRefinedSystemTest::testSystem().

◆ elem_position_get()

void libMesh::FEMContext::elem_position_get ( )

inherited

Uses the geometry of elem to set the coordinate data specified by mesh_*_position configuration.

Definition at line 1551 of file fem_context.C.

{
  // This is too expensive to call unless we've been asked to move the mesh
  libmesh_assert (_mesh_sys);
 
  // This will probably break with threading when two contexts are
  // operating on elements which share a node
  libmesh_assert_equal_to (libMesh::n_threads(), 1);
 
  // If the coordinate data is in our own system, it's already
  // been set up for us
  //  if (_mesh_sys == this->number())
  //    {
  unsigned int n_nodes = this->get_elem().n_nodes();
 
#ifndef NDEBUG
  const unsigned char dim = this->get_elem_dim();
 
  // For simplicity we demand that mesh coordinates be stored
  // in a format that allows a direct copy
  libmesh_assert(this->get_mesh_x_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_x_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_element_fe(this->get_mesh_x_var(), dim)->get_fe_type().order.get_order()
                  == this->get_elem().default_order()));
  libmesh_assert(this->get_mesh_y_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_y_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_element_fe(this->get_mesh_y_var(), dim)->get_fe_type().order.get_order()
                  == this->get_elem().default_order()));
  libmesh_assert(this->get_mesh_z_var() == libMesh::invalid_uint ||
                 (this->get_element_fe(this->get_mesh_z_var(), dim)->get_fe_type().family
                  == FEMap::map_fe_type(this->get_elem()) &&
                  this->get_element_fe(this->get_mesh_z_var(), dim)->get_fe_type().order.get_order()
                  == this->get_elem().default_order()));
#endif
 
  // Get degree of freedom coefficients from point coordinates
  if (this->get_mesh_x_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      (this->get_elem_solution(this->get_mesh_x_var()))(i) = this->get_elem().point(i)(0);
 
  if (this->get_mesh_y_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      (this->get_elem_solution(this->get_mesh_y_var()))(i) = this->get_elem().point(i)(1);
 
  if (this->get_mesh_z_var() != libMesh::invalid_uint)
    for (unsigned int i=0; i != n_nodes; ++i)
      (this->get_elem_solution(this->get_mesh_z_var()))(i) = this->get_elem().point(i)(2);
  //    }
  // FIXME - If the coordinate data is not in our own system, someone
  // had better get around to implementing that... - RHS
  //  else
  //    {
  //      libmesh_not_implemented();
  //    }
}

References libMesh::FEMContext::_mesh_sys, dim, libMesh::FEMContext::get_elem(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::get_element_fe(), libMesh::FEMContext::get_mesh_x_var(), libMesh::FEMContext::get_mesh_y_var(), libMesh::FEMContext::get_mesh_z_var(), libMesh::invalid_uint, libMesh::libmesh_assert(), libMesh::FEMap::map_fe_type(), libMesh::Elem::n_nodes(), n_nodes, libMesh::n_threads(), and libMesh::Elem::point().

Referenced by libMesh::FEMSystem::mesh_position_get().

◆ elem_position_set()

void libMesh::FEMContext::elem_position_set ( Real theta )

inlineinherited

Uses the coordinate data specified by mesh_*_position configuration to set the geometry of elem to the value it would take after a fraction theta of a timestep.

Definition at line 1288 of file fem_context.h.

{
  if (_mesh_sys)
    this->_do_elem_position_set(theta);
}

References libMesh::FEMContext::_do_elem_position_set(), and libMesh::FEMContext::_mesh_sys.

Referenced by libMesh::FEMContext::elem_edge_reinit(), libMesh::FEMContext::elem_reinit(), libMesh::FEMContext::elem_side_reinit(), and libMesh::FEMSystem::mesh_position_set().

◆ elem_reinit()

void libMesh::FEMContext::elem_reinit ( Real theta )

overridevirtualinherited

Resets the current time in the context.

Additionally, reinitialize Elem and FE objects if there's a moving mesh present in the system such that the mesh is deformed to its position at \( t_{\theta} \).

Reimplemented from libMesh::DiffContext.

Definition at line 1437 of file fem_context.C.

{
  // Update the "time" variable of this context object
  this->_update_time_from_system(theta);
 
  // Handle a moving element if necessary.
  if (_mesh_sys)
    {
      // We assume that the ``default'' state
      // of the mesh is its final, theta=1.0
      // position, so we don't bother with
      // mesh motion in that case.
      if (theta != 1.0)
        {
          // FIXME - ALE is not threadsafe yet!
          libmesh_assert_equal_to (libMesh::n_threads(), 1);
 
          elem_position_set(theta);
        }
      elem_fe_reinit();
    }
}

References libMesh::FEMContext::_mesh_sys, libMesh::FEMContext::_update_time_from_system(), libMesh::FEMContext::elem_fe_reinit(), libMesh::FEMContext::elem_position_set(), and libMesh::n_threads().

◆ elem_side_reinit()

void libMesh::FEMContext::elem_side_reinit ( Real theta )

overridevirtualinherited

Resets the current time in the context.

Additionally, reinitialize Elem and FE objects if there's a moving mesh present in the system such that the mesh is deformed to its position at \( t_{\theta} \).

Reimplemented from libMesh::DiffContext.

Definition at line 1461 of file fem_context.C.

{
  // Update the "time" variable of this context object
  this->_update_time_from_system(theta);
 
  // Handle a moving element if necessary
  if (_mesh_sys)
    {
      // FIXME - not threadsafe yet!
      elem_position_set(theta);
      side_fe_reinit();
    }
}

References libMesh::FEMContext::_mesh_sys, libMesh::FEMContext::_update_time_from_system(), libMesh::FEMContext::elem_position_set(), and libMesh::FEMContext::side_fe_reinit().

◆ find_hardest_fe_type()

FEType libMesh::FEMContext::find_hardest_fe_type ( )

inherited

Helper function for creating quadrature rules.

Definition at line 85 of file fem_context.C.

{
  const System & sys = this->get_system();
  FEType hardest_fe_type =
    sys.variable_type(_active_vars ?
                      (*_active_vars)[0] : 0);
 
  auto check_var = [&hardest_fe_type, &sys](unsigned int v)
    {
      FEType fe_type = sys.variable_type(v);
 
      // Make sure we find a non-SCALAR FE family, even in the case
      // where the first variable(s) weren't
      if (hardest_fe_type.family == SCALAR)
        {
          hardest_fe_type.family = fe_type.family;
          hardest_fe_type.order = fe_type.order;
        }
 
      // FIXME - we don't yet handle mixed finite elements from
      // different families which require different quadrature rules
      // libmesh_assert_equal_to (fe_type.family, hardest_fe_type.family);
 
      // We need to detect SCALAR's so we can prepare FE objects for
      // them, and so we don't mistake high order scalars as a reason
      // to crank up the quadrature order on other types.
      if (fe_type.family != SCALAR && fe_type.order > hardest_fe_type.order)
        hardest_fe_type = fe_type;
    };
 
  if (_active_vars)
    for (auto v : *_active_vars)
      check_var(v);
  else
    for (auto v : make_range(sys.n_vars()))
      check_var(v);
 
  return hardest_fe_type;
}

References libMesh::FEMContext::_active_vars, libMesh::FEType::family, libMesh::DiffContext::get_system(), libMesh::make_range(), libMesh::System::n_vars(), libMesh::FEType::order, libMesh::SCALAR, and libMesh::System::variable_type().

Referenced by libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ fixed_interior_gradient() [1/2]

Gradient libMesh::FEMContext::fixed_interior_gradient	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 1106 of file fem_context.C.

{
  Gradient du;
 
  this->fixed_interior_gradient( var, qp, du );
 
  return du;
}

References libMesh::FEMContext::fixed_interior_gradient().

Referenced by libMesh::FEMContext::fixed_interior_gradient().

◆ fixed_interior_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_interior_gradient	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	grad_u
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 1117 of file fem_context.C.

{
  this->some_gradient
    <OutputType,
     &FEMContext::get_element_fe
     <typename TensorTools::MakeReal
      <typename TensorTools::DecrementRank
       <OutputType>::type>::type>,
     &DiffContext::get_elem_fixed_solution>
    (var, qp, du);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_gradient().

◆ fixed_interior_hessian() [1/2]

Tensor libMesh::FEMContext::fixed_interior_hessian	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 1133 of file fem_context.C.

{
  Tensor d2u;
 
  this->fixed_interior_hessian( var, qp, d2u );
 
  return d2u;
}

References libMesh::FEMContext::fixed_interior_hessian().

Referenced by libMesh::FEMContext::fixed_interior_hessian().

◆ fixed_interior_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_interior_hessian	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	hess_u
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 1144 of file fem_context.C.

{
  this->some_hessian<OutputType,
                     &FEMContext::get_element_fe
                     <typename TensorTools::MakeReal
                      <typename TensorTools::DecrementRank
                       <typename TensorTools::DecrementRank
                        <OutputType>::type>::type>::type>,
                     &DiffContext::get_elem_fixed_solution>(var, qp, d2u);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_hessian().

◆ fixed_interior_value() [1/2]

Number libMesh::FEMContext::fixed_interior_value	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The value of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 1083 of file fem_context.C.

{
  Number u = 0.;
 
  this->fixed_interior_value( var, qp, u );
 
  return u;
}

References libMesh::FEMContext::fixed_interior_value().

Referenced by libMesh::FEMContext::fixed_interior_value().

◆ fixed_interior_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_interior_value	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The value of the fixed_solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 1095 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_element_fe
                   <typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_fixed_solution>(var, qp, u);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_value().

◆ fixed_point_gradient() [1/2]

Gradient libMesh::FEMContext::fixed_point_gradient	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the physical point p on the current element.

Note: This API is currently present for backward compatibility.

Definition at line 1279 of file fem_context.C.

{
  Gradient grad_u;
 
  this->fixed_point_gradient( var, p, grad_u );
 
  return grad_u;
}

References libMesh::FEMContext::fixed_point_gradient().

Referenced by libMesh::FEMContext::fixed_point_gradient().

◆ fixed_point_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_point_gradient	(	unsigned int	var,
		const Point &	p,
		OutputType &	grad_u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 1291 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal
    <typename TensorTools::DecrementRank<OutputType>::type>::type
    OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (_elem_fixed_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_fixed_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 1 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputGradient>> &  dphi = fe_new->get_dphi();
 
  grad_u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    grad_u.add_scaled(dphi[l][0], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_fixed_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_fixed_solution().

◆ fixed_point_hessian() [1/2]

Tensor libMesh::FEMContext::fixed_point_hessian	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the physical point p on the current element.

Note: This API is currently present for backward compatibility.

Definition at line 1330 of file fem_context.C.

{
  Tensor hess_u;
 
  this->fixed_point_hessian( var, p, hess_u );
 
  return hess_u;
}

References libMesh::FEMContext::fixed_point_hessian().

Referenced by libMesh::FEMContext::fixed_point_hessian().

◆ fixed_point_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_point_hessian	(	unsigned int	var,
		const Point &	p,
		OutputType &	hess_u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 1342 of file fem_context.C.

{
  typedef typename TensorTools::DecrementRank<OutputType>::type Rank1Decrement;
  typedef typename TensorTools::DecrementRank<Rank1Decrement>::type Rank2Decrement;
  typedef typename TensorTools::MakeReal<Rank2Decrement>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (_elem_fixed_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_fixed_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 2 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputTensor>> &  d2phi = fe_new->get_d2phi();
 
  hess_u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    hess_u.add_scaled(d2phi[l][0], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_fixed_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_fixed_solution().

◆ fixed_point_value() [1/2]

Number libMesh::FEMContext::fixed_point_value	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The value of the fixed_solution variable var at the physical point p on the current element.

Note: This API is currently present for backward compatibility.

Definition at line 1233 of file fem_context.C.

{
  Number u = 0.;
 
  this->fixed_point_value( var, p, u );
 
  return u;
}

References libMesh::FEMContext::fixed_point_value().

Referenced by libMesh::FEMContext::fixed_point_value().

◆ fixed_point_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_point_value	(	unsigned int	var,
		const Point &	p,
		OutputType &	u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The value of the fixed_solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 1243 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (_elem_fixed_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_fixed_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 0 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<OutputShape>> &  phi = fe_new->get_phi();
 
  u = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    u += phi[l][0] * coef(l);
 
  return;
}

References libMesh::DiffContext::_elem_fixed_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_elem_fixed_solution(), and libMesh::FEGenericBase< OutputType >::get_phi().

◆ fixed_side_gradient() [1/2]

Gradient libMesh::FEMContext::fixed_side_gradient	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This API is currently present for backward compatibility.

Definition at line 1183 of file fem_context.C.

{
  Gradient du;
 
  this->fixed_side_gradient( var, qp, du );
 
  return du;
}

References libMesh::FEMContext::fixed_side_gradient().

Referenced by libMesh::FEMContext::fixed_side_gradient().

◆ fixed_side_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_side_gradient	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	grad_u
	)		const

inherited

Returns: The gradient of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 1194 of file fem_context.C.

{
  this->some_gradient<OutputType,
                      &FEMContext::get_side_fe
                      <typename TensorTools::MakeReal
                       <typename TensorTools::DecrementRank
                        <OutputType>::type>::type>,
                      &DiffContext::get_elem_fixed_solution>(var, qp, du);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_gradient().

◆ fixed_side_hessian() [1/2]

Tensor libMesh::FEMContext::fixed_side_hessian	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This API is currently present for backward compatibility.

Definition at line 1208 of file fem_context.C.

{
  Tensor d2u;
 
  this->fixed_side_hessian( var, qp, d2u );
 
  return d2u;
}

References libMesh::FEMContext::fixed_side_hessian().

Referenced by libMesh::FEMContext::fixed_side_hessian().

◆ fixed_side_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_side_hessian	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	hess_u
	)		const

inherited

Returns: The hessian of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 1218 of file fem_context.C.

{
  this->some_hessian<OutputType,
                     &FEMContext::get_side_fe
                     <typename TensorTools::MakeReal
                      <typename TensorTools::DecrementRank
                       <typename TensorTools::DecrementRank
                        <OutputType>::type>::type>::type>,
                     &DiffContext::get_elem_fixed_solution>(var, qp, d2u);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_hessian().

◆ fixed_side_value() [1/2]

Number libMesh::FEMContext::fixed_side_value	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The value of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This API is currently present for backward compatibility.

Definition at line 1159 of file fem_context.C.

{
  Number u = 0.;
 
  this->fixed_side_value( var, qp, u );
 
  return u;
}

References libMesh::FEMContext::fixed_side_value().

Referenced by libMesh::FEMContext::fixed_side_value().

◆ fixed_side_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::fixed_side_value	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The value of the fixed_solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 1170 of file fem_context.C.

{
  this->some_value
    <OutputType,
     &FEMContext::get_side_fe
     <typename TensorTools::MakeReal<OutputType>::type>,
     &DiffContext::get_elem_fixed_solution>
    (var, qp, u);
}

References libMesh::DiffContext::get_elem_fixed_solution(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_value().

◆ get_deltat_value()

Real libMesh::DiffContext::get_deltat_value ( )

inherited

Returns: The value currently pointed to by this class's _deltat member

Definition at line 117 of file diff_context.C.

{
  libmesh_assert(_deltat);
 
  return *_deltat;
}

References libMesh::DiffContext::_deltat, and libMesh::libmesh_assert().

Referenced by libMesh::FEMContext::_update_time_from_system().

◆ get_dim()

unsigned char libMesh::FEMContext::get_dim ( ) const

inlineinherited

Accessor for cached mesh dimension.

This is the largest dimension of the elements in the mesh. For the dimension of this->_elem, use get_elem_dim();

Definition at line 954 of file fem_context.h.

955 { return this->_dim; }

References libMesh::FEMContext::_dim.

Referenced by libMesh::VariationalSmootherSystem::element_time_derivative(), neighbor_side_fe_reinit(), and SolidSystem::side_time_derivative().

◆ get_dof_indices() [1/4]

std::vector< dof_id_type > & libMesh::DiffContext::get_dof_indices ( )

inlineinherited

Non-const accessor for element dof indices.

Definition at line 369 of file diff_context.h.

370 { return _dof_indices; }

libMesh::DiffContext::_dof_indices

std::vector< dof_id_type > _dof_indices

Global Degree of freedom index lists.

Definition diff_context.h:639

References libMesh::DiffContext::_dof_indices.

◆ get_dof_indices() [2/4]

const std::vector< dof_id_type > & libMesh::DiffContext::get_dof_indices ( ) const

inlineinherited

Accessor for element dof indices.

Definition at line 363 of file diff_context.h.

364 { return _dof_indices; }

References libMesh::DiffContext::_dof_indices.

◆ get_dof_indices() [3/4]

std::vector< dof_id_type > & libMesh::DiffContext::get_dof_indices ( unsigned int var )

inlineinherited

Accessor for element dof indices of a particular variable corresponding to the index argument.

Definition at line 386 of file diff_context.h.

  {
    libmesh_assert_greater(_dof_indices_var.size(), var);
    return _dof_indices_var[var];
  }

References libMesh::DiffContext::_dof_indices_var.

◆ get_dof_indices() [4/4]

const std::vector< dof_id_type > & libMesh::DiffContext::get_dof_indices ( unsigned int var ) const

inlineinherited

Accessor for element dof indices of a particular variable corresponding to the index argument.

Definition at line 376 of file diff_context.h.

  {
    libmesh_assert_greater(_dof_indices_var.size(), var);
    return _dof_indices_var[var];
  }

References libMesh::DiffContext::_dof_indices_var.

◆ get_edge()

unsigned char libMesh::FEMContext::get_edge ( ) const

inlineinherited

Accessor for current edge of Elem object.

Definition at line 946 of file fem_context.h.

947 { return edge; }

libMesh::FEMContext::edge

unsigned char edge

Current edge for edge_* to examine.

Definition fem_context.h:1036

References libMesh::FEMContext::edge.

Referenced by libMesh::FEMContext::edge_fe_reinit().

◆ get_edge_fe() [1/3]

FEBase * libMesh::FEMContext::get_edge_fe ( unsigned int var ) const

inlineinherited

Accessor for edge (3D only!) finite element object for scalar-valued variable var.

Definition at line 1364 of file fem_context.h.

{
  libmesh_assert_less ( var, _edge_fe_var.size() );
  return cast_ptr<FEBase *>( _edge_fe_var[var] );
}

References libMesh::FEMContext::_edge_fe_var.

◆ get_edge_fe() [2/3]

void libMesh::FEMContext::get_edge_fe	(	unsigned int	var,
		FEAbstract *&	fe
	)		const

inlineinherited

Definition at line 1357 of file fem_context.h.

{
  libmesh_assert_less ( var, _edge_fe_var.size() );
  fe = _edge_fe_var[var];
}

References libMesh::FEMContext::_edge_fe_var.

◆ get_edge_fe() [3/3]

template<typename OutputShape >

void libMesh::FEMContext::get_edge_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe
	)		const

inlineinherited

Accessor for edge (3D only!) finite element object for variable var.

Definition at line 1350 of file fem_context.h.

{
  libmesh_assert_less ( var, _edge_fe_var.size() );
  fe = cast_ptr<FEGenericBase<OutputShape> *>( _edge_fe_var[var] );
}

References libMesh::FEMContext::_edge_fe_var.

Referenced by libMesh::GenericProjector< FFunctor, GFunctor, FValue, ProjectionAction >::SubFunctor::SubFunctor().

◆ get_edge_qrule()

const QBase & libMesh::FEMContext::get_edge_qrule ( ) const

inlineinherited

Accessor for element edge quadrature rule.

Definition at line 849 of file fem_context.h.

850 { return *(this->_edge_qrule); }

References libMesh::FEMContext::_edge_qrule.

◆ get_elem() [1/2]

Elem & libMesh::FEMContext::get_elem ( )

inlineinherited

Accessor for current Elem object.

Definition at line 933 of file fem_context.h.

934 { libmesh_assert(this->_elem);

935 return *(const_cast<Elem *>(this->_elem)); }

libMesh::FEMContext::_elem

const Elem * _elem

Current element for element_* to examine.

Definition fem_context.h:1210

References libMesh::FEMContext::_elem, and libMesh::libmesh_assert().

◆ get_elem() [2/2]

const Elem & libMesh::FEMContext::get_elem ( ) const

inlineinherited

Accessor for current Elem object.

Definition at line 926 of file fem_context.h.

927 { libmesh_assert(this->_elem);

928 return *(this->_elem); }

References libMesh::FEMContext::_elem, and libMesh::libmesh_assert().

Referenced by libMesh::FEMContext::_do_elem_position_set(), libMesh::RBConstruction::add_scaled_matrix_and_vector(), AssemblyF0::boundary_assembly(), AssemblyF1::boundary_assembly(), AssemblyF2::boundary_assembly(), AssemblyA0::boundary_assembly(), AssemblyA1::boundary_assembly(), AssemblyA2::boundary_assembly(), libMesh::FEMContext::build_new_fe(), libMesh::FEMContext::cached_fe(), libMesh::OldSolutionBase< Output, point_output >::check_old_context(), libMesh::OldSolutionBase< Output, point_output >::check_old_context(), libMesh::FEMContext::elem_position_get(), HeatSystem::element_time_derivative(), SolidSystem::element_time_derivative(), HilbertSystem::element_time_derivative(), libMesh::VariationalSmootherSystem::element_time_derivative(), libMesh::OldSolutionValue< Output, point_output >::eval_at_node(), libMesh::OldSolutionCoefs< Output, point_output >::eval_at_node(), libMesh::OldSolutionValue< Output, point_output >::eval_at_node(), libMesh::OldSolutionCoefs< Output, point_output >::eval_at_node(), libMesh::FEMContext::has_side_boundary_id(), libMesh::RBEIMConstruction::initialize_qp_data(), OutputAssembly::interior_assembly(), A0::interior_assembly(), A1::interior_assembly(), EIM_F::interior_assembly(), AssemblyEIM::interior_assembly(), libMesh::FEMSystem::mesh_position_set(), libMesh::FEMSystem::numerical_jacobian(), SlitFunc::operator()(), libMesh::FEMContext::pre_fe_reinit(), libMesh::FEMContext::side_boundary_ids(), NavierSystem::side_constraint(), and SolidSystem::side_time_derivative().

◆ get_elem_dim()

unsigned char libMesh::FEMContext::get_elem_dim ( ) const

inlineinherited

Returns: The dimension of this->_elem. For mixed dimension meshes, this may be different from get_dim(). If no element init has happened yet, fall back on get_dim().

Definition at line 962 of file fem_context.h.

963 { return this->_elem ? this->_elem_dim : this->_dim; }

libMesh::FEMContext::_elem_dim

unsigned char _elem_dim

Cached dimension of this->_elem.

Definition fem_context.h:1220

References libMesh::FEMContext::_dim, libMesh::FEMContext::_elem, and libMesh::FEMContext::_elem_dim.

◆ get_elem_elem_jacobian() [1/4]

DenseMatrix< Number > & libMesh::DGFEMContext::get_elem_elem_jacobian ( )

inline

Non-const accessor for element-element Jacobian.

Definition at line 116 of file dg_fem_context.h.

117 { return _elem_elem_jacobian; }

References _elem_elem_jacobian.

◆ get_elem_elem_jacobian() [2/4]

const DenseMatrix< Number > & libMesh::DGFEMContext::get_elem_elem_jacobian ( ) const

inline

Const accessor for element-element Jacobian.

Definition at line 110 of file dg_fem_context.h.

111 { return _elem_elem_jacobian; }

References _elem_elem_jacobian.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ get_elem_elem_jacobian() [3/4]

DenseSubMatrix< Number > & libMesh::DGFEMContext::get_elem_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)

inline

Non-const accessor for element-element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 130 of file dg_fem_context.h.

131 { return *(_elem_elem_subjacobians[var1][var2]); }

References _elem_elem_subjacobians.

◆ get_elem_elem_jacobian() [4/4]

const DenseSubMatrix< Number > & libMesh::DGFEMContext::get_elem_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)		const

inline

Const accessor for element-element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 123 of file dg_fem_context.h.

124 { return *(_elem_elem_subjacobians[var1][var2]); }

References _elem_elem_subjacobians.

◆ get_elem_fixed_solution() [1/4]

DenseVector< Number > & libMesh::DiffContext::get_elem_fixed_solution ( )

inlineinherited

Non-const accessor for element fixed solution.

Definition at line 216 of file diff_context.h.

217 { return _elem_fixed_solution; }

libMesh::DiffContext::_elem_fixed_solution

DenseVector< Number > _elem_fixed_solution

Element by element components of nonlinear_solution at a fixed point in a timestep,...

Definition diff_context.h:604

References libMesh::DiffContext::_elem_fixed_solution.

◆ get_elem_fixed_solution() [2/4]

const DenseVector< Number > & libMesh::DiffContext::get_elem_fixed_solution ( ) const

inlineinherited

Accessor for element fixed solution.

Definition at line 210 of file diff_context.h.

211 { return _elem_fixed_solution; }

References libMesh::DiffContext::_elem_fixed_solution.

◆ get_elem_fixed_solution() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_elem_fixed_solution ( unsigned int var )

inlineinherited

Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

Definition at line 233 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_fixed_subsolutions.size(), var);
    return _elem_fixed_subsolutions[var];
  }

References libMesh::DiffContext::_elem_fixed_subsolutions.

◆ get_elem_fixed_solution() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_elem_fixed_solution ( unsigned int var ) const

inlineinherited

Accessor for element fixed solution of a particular variable corresponding to the variable index argument.

Definition at line 223 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_fixed_subsolutions.size(), var);
    return _elem_fixed_subsolutions[var];
  }

References libMesh::DiffContext::_elem_fixed_subsolutions.

◆ get_elem_jacobian() [1/4]

DenseMatrix< Number > & libMesh::DiffContext::get_elem_jacobian ( )

inlineinherited

Non-const accessor for element Jacobian.

Definition at line 283 of file diff_context.h.

  {
    libmesh_assert(_have_local_matrices);
    return _elem_jacobian;
  }

References libMesh::DiffContext::_elem_jacobian, libMesh::DiffContext::_have_local_matrices, and libMesh::libmesh_assert().

◆ get_elem_jacobian() [2/4]

const DenseMatrix< Number > & libMesh::DiffContext::get_elem_jacobian ( ) const

inlineinherited

Const accessor for element Jacobian.

Definition at line 274 of file diff_context.h.

  {
    libmesh_assert(_have_local_matrices);
    return _elem_jacobian;
  }

References libMesh::DiffContext::_elem_jacobian, libMesh::DiffContext::_have_local_matrices, and libMesh::libmesh_assert().

◆ get_elem_jacobian() [3/4]

DenseSubMatrix< Number > & libMesh::DiffContext::get_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)

inlineinherited

Non-const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

Only available if _have_local_matrices

Definition at line 306 of file diff_context.h.

  {
    libmesh_assert(_have_local_matrices);
    libmesh_assert_greater(_elem_subjacobians.size(), var1);
    libmesh_assert_greater(_elem_subjacobians[var1].size(), var2);
    return _elem_subjacobians[var1][var2];
  }

References libMesh::DiffContext::_elem_subjacobians, libMesh::DiffContext::_have_local_matrices, and libMesh::libmesh_assert().

◆ get_elem_jacobian() [4/4]

const DenseSubMatrix< Number > & libMesh::DiffContext::get_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)		const

inlineinherited

Const accessor for element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 293 of file diff_context.h.

  {
    libmesh_assert(_have_local_matrices);
    libmesh_assert_greater(_elem_subjacobians.size(), var1);
    libmesh_assert_greater(_elem_subjacobians[var1].size(), var2);
    return _elem_subjacobians[var1][var2];
  }

References libMesh::DiffContext::_elem_subjacobians, libMesh::DiffContext::_have_local_matrices, and libMesh::libmesh_assert().

◆ get_elem_neighbor_jacobian() [1/4]

DenseMatrix< Number > & libMesh::DGFEMContext::get_elem_neighbor_jacobian ( )

inline

Non-const accessor for element -neighborJacobian.

Definition at line 142 of file dg_fem_context.h.

143 { return _elem_neighbor_jacobian; }

References _elem_neighbor_jacobian.

◆ get_elem_neighbor_jacobian() [2/4]

const DenseMatrix< Number > & libMesh::DGFEMContext::get_elem_neighbor_jacobian ( ) const

inline

Const accessor for element-neighbor Jacobian.

Definition at line 136 of file dg_fem_context.h.

137 { return _elem_neighbor_jacobian; }

References _elem_neighbor_jacobian.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ get_elem_neighbor_jacobian() [3/4]

DenseSubMatrix< Number > & libMesh::DGFEMContext::get_elem_neighbor_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)

inline

Non-const accessor for element-neighbor Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 156 of file dg_fem_context.h.

157 { return *(_elem_neighbor_subjacobians[var1][var2]); }

References _elem_neighbor_subjacobians.

◆ get_elem_neighbor_jacobian() [4/4]

const DenseSubMatrix< Number > & libMesh::DGFEMContext::get_elem_neighbor_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)		const

inline

Const accessor for element-neighbor Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 149 of file dg_fem_context.h.

150 { return *(_elem_neighbor_subjacobians[var1][var2]); }

References _elem_neighbor_subjacobians.

◆ get_elem_residual() [1/4]

DenseVector< Number > & libMesh::DiffContext::get_elem_residual ( )

inlineinherited

Non-const accessor for element residual.

Definition at line 248 of file diff_context.h.

249 { return _elem_residual; }

libMesh::DiffContext::_elem_residual

DenseVector< Number > _elem_residual

Element residual vector.

Definition diff_context.h:610

References libMesh::DiffContext::_elem_residual.

◆ get_elem_residual() [2/4]

const DenseVector< Number > & libMesh::DiffContext::get_elem_residual ( ) const

inlineinherited

Const accessor for element residual.

Definition at line 242 of file diff_context.h.

243 { return _elem_residual; }

References libMesh::DiffContext::_elem_residual.

◆ get_elem_residual() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_elem_residual ( unsigned int var )

inlineinherited

Non-const accessor for element residual of a particular variable corresponding to the variable index argument.

Definition at line 265 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subresiduals.size(), var);
    return _elem_subresiduals[var];
  }

References libMesh::DiffContext::_elem_subresiduals.

◆ get_elem_residual() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_elem_residual ( unsigned int var ) const

inlineinherited

Const accessor for element residual of a particular variable corresponding to the variable index argument.

Definition at line 255 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subresiduals.size(), var);
    return _elem_subresiduals[var];
  }

References libMesh::DiffContext::_elem_subresiduals.

◆ get_elem_solution() [1/4]

DenseVector< Number > & libMesh::DiffContext::get_elem_solution ( )

inlineinherited

Non-const accessor for element solution.

Definition at line 118 of file diff_context.h.

119 { return _elem_solution; }

libMesh::DiffContext::_elem_solution

DenseVector< Number > _elem_solution

Element by element components of nonlinear_solution as adjusted by a time_solver.

Definition diff_context.h:582

References libMesh::DiffContext::_elem_solution.

◆ get_elem_solution() [2/4]

const DenseVector< Number > & libMesh::DiffContext::get_elem_solution ( ) const

inlineinherited

Accessor for element solution.

Definition at line 112 of file diff_context.h.

113 { return _elem_solution; }

References libMesh::DiffContext::_elem_solution.

◆ get_elem_solution() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution ( unsigned int var )

inlineinherited

Accessor for element solution of a particular variable corresponding to the variable index argument.

Definition at line 135 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolutions.size(), var);
    return _elem_subsolutions[var];
  }

References libMesh::DiffContext::_elem_subsolutions.

◆ get_elem_solution() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution ( unsigned int var ) const

inlineinherited

Accessor for element solution of a particular variable corresponding to the variable index argument.

Definition at line 125 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolutions.size(), var);
    return _elem_subsolutions[var];
  }

References libMesh::DiffContext::_elem_subsolutions.

◆ get_elem_solution_accel() [1/4]

DenseVector< Number > & libMesh::DiffContext::get_elem_solution_accel ( )

inlineinherited

Non-const accessor for element solution accel of change w.r.t.

time.

Definition at line 184 of file diff_context.h.

185 { return _elem_solution_accel; }

libMesh::DiffContext::_elem_solution_accel

DenseVector< Number > _elem_solution_accel

Element by element components of du/dt as adjusted by a time_solver.

Definition diff_context.h:596

References libMesh::DiffContext::_elem_solution_accel.

◆ get_elem_solution_accel() [2/4]

const DenseVector< Number > & libMesh::DiffContext::get_elem_solution_accel ( ) const

inlineinherited

Accessor for element solution accel of change w.r.t.

time.

Definition at line 177 of file diff_context.h.

178 { return _elem_solution_accel; }

References libMesh::DiffContext::_elem_solution_accel.

Referenced by libMesh::NewmarkSolver::_general_residual(), libMesh::FEMContext::interior_accel(), libMesh::FEMContext::pre_fe_reinit(), libMesh::FirstOrderUnsteadySolver::prepare_accel(), and libMesh::FEMContext::side_accel().

◆ get_elem_solution_accel() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution_accel ( unsigned int var )

inlineinherited

Accessor for element solution accel for a particular variable corresponding to the variable index argument.

Definition at line 201 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolution_accels.size(), var);
    return _elem_subsolution_accels[var];
  }

References libMesh::DiffContext::_elem_subsolution_accels.

◆ get_elem_solution_accel() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution_accel ( unsigned int var ) const

inlineinherited

Accessor for element solution accel for a particular variable corresponding to the variable index argument.

Definition at line 191 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolution_accels.size(), var);
    return _elem_subsolution_accels[var];
  }

References libMesh::DiffContext::_elem_subsolution_accels.

◆ get_elem_solution_accel_derivative()

Real libMesh::DiffContext::get_elem_solution_accel_derivative ( ) const

inlineinherited

The derivative of the current elem_solution_accel w.r.t.

the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_accel_derivative() is 0.

Definition at line 450 of file diff_context.h.

451 { return elem_solution_accel_derivative; }

libMesh::DiffContext::elem_solution_accel_derivative

Real elem_solution_accel_derivative

The derivative of elem_solution_accel with respect to the current nonlinear solution,...

Definition diff_context.h:510

References libMesh::DiffContext::elem_solution_accel_derivative.

Referenced by ElasticitySystem::mass_residual(), SecondOrderScalarSystemSecondOrderTimeSolverBase::mass_residual(), and SecondOrderScalarSystemFirstOrderTimeSolverBase::mass_residual().

◆ get_elem_solution_derivative()

Real libMesh::DiffContext::get_elem_solution_derivative ( ) const

inlineinherited

The derivative of the current elem_solution w.r.t.

the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_derivative() is 0.

Definition at line 432 of file diff_context.h.

433 { return elem_solution_derivative; }

libMesh::DiffContext::elem_solution_derivative

Real elem_solution_derivative

The derivative of elem_solution with respect to the current nonlinear solution.

Definition diff_context.h:496

References libMesh::DiffContext::elem_solution_derivative.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), NavierSystem::element_constraint(), HeatSystem::element_time_derivative(), ElasticitySystem::element_time_derivative(), NonlinearDiffusionSystem::element_time_derivative(), HilbertSystem::element_time_derivative(), NavierSystem::mass_residual(), and NavierSystem::side_constraint().

◆ get_elem_solution_rate() [1/4]

DenseVector< Number > & libMesh::DiffContext::get_elem_solution_rate ( )

inlineinherited

Non-const accessor for element solution rate of change w.r.t.

time.

Definition at line 151 of file diff_context.h.

152 { return _elem_solution_rate; }

libMesh::DiffContext::_elem_solution_rate

DenseVector< Number > _elem_solution_rate

Element by element components of du/dt as adjusted by a time_solver.

Definition diff_context.h:589

References libMesh::DiffContext::_elem_solution_rate.

◆ get_elem_solution_rate() [2/4]

const DenseVector< Number > & libMesh::DiffContext::get_elem_solution_rate ( ) const

inlineinherited

Accessor for element solution rate of change w.r.t.

time.

Definition at line 144 of file diff_context.h.

145 { return _elem_solution_rate; }

References libMesh::DiffContext::_elem_solution_rate.

Referenced by libMesh::NewmarkSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::EulerSolver::_general_residual(), libMesh::FEMContext::interior_rate(), libMesh::FEMContext::interior_rate_gradient(), libMesh::FEMContext::pre_fe_reinit(), libMesh::FirstOrderUnsteadySolver::prepare_accel(), and libMesh::FEMContext::side_rate().

◆ get_elem_solution_rate() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution_rate ( unsigned int var )

inlineinherited

Accessor for element solution rate for a particular variable corresponding to the variable index argument.

Definition at line 168 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolution_rates.size(), var);
    return _elem_subsolution_rates[var];
  }

References libMesh::DiffContext::_elem_subsolution_rates.

◆ get_elem_solution_rate() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_elem_solution_rate ( unsigned int var ) const

inlineinherited

Accessor for element solution rate for a particular variable corresponding to the variable index argument.

Definition at line 158 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_subsolution_rates.size(), var);
    return _elem_subsolution_rates[var];
  }

References libMesh::DiffContext::_elem_subsolution_rates.

◆ get_elem_solution_rate_derivative()

Real libMesh::DiffContext::get_elem_solution_rate_derivative ( ) const

inlineinherited

The derivative of the current elem_solution_rate w.r.t.

the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_elem_solution_rate_derivative() is 0.

Definition at line 441 of file diff_context.h.

442 { return elem_solution_rate_derivative; }

libMesh::DiffContext::elem_solution_rate_derivative

Real elem_solution_rate_derivative

The derivative of elem_solution_rate with respect to the current nonlinear solution,...

Definition diff_context.h:503

References libMesh::DiffContext::elem_solution_rate_derivative.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), SecondOrderScalarSystemSecondOrderTimeSolverBase::damping_residual(), SecondOrderScalarSystemFirstOrderTimeSolverBase::damping_residual(), FirstOrderScalarSystemBase::mass_residual(), and libMesh::FirstOrderUnsteadySolver::prepare_accel().

◆ get_element_fe() [1/5]

FEBase * libMesh::FEMContext::get_element_fe ( unsigned int var ) const

inlineinherited

Accessor for interior finite element object for scalar-valued variable var for the largest dimension in the mesh.

We default to the largest mesh dim if this method is called before the Elem * is set in the FEMContext, e.g. in FEMSystem::init_context (or a subclass).

Definition at line 304 of file fem_context.h.

305 { return this->get_element_fe(var,this->get_elem_dim()); }

References libMesh::FEMContext::get_elem_dim(), and libMesh::FEMContext::get_element_fe().

◆ get_element_fe() [2/5]

void libMesh::FEMContext::get_element_fe	(	unsigned int	var,
		FEAbstract *&	fe,
		unsigned short	dim
	)		const

inlineinherited

Accessor for interior finite element object for variable var for dimension dim.

Definition at line 1305 of file fem_context.h.

{
  libmesh_assert( !_element_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_element_fe_var[dim].size() ) );
  fe = _element_fe_var[dim][var];
}

References libMesh::FEMContext::_element_fe_var, dim, and libMesh::libmesh_assert().

◆ get_element_fe() [3/5]

template<typename OutputShape >

void libMesh::FEMContext::get_element_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe
	)		const

inlineinherited

Accessor for interior finite element object for variable var for the largest dimension in the mesh.

We default to the largest mesh dim if this method is called before the Elem * is set in the FEMContext, e.g. in FEMSystem::init_context (or a subclass).

Definition at line 295 of file fem_context.h.

296 { this->get_element_fe<OutputShape>(var,fe,this->get_elem_dim()); }

References libMesh::FEMContext::get_elem_dim().

◆ get_element_fe() [4/5]

template<typename OutputShape >

void libMesh::FEMContext::get_element_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe,
		unsigned short	dim
	)		const

inlineinherited

Accessor for interior finite element object for variable var for dimension dim.

Definition at line 1296 of file fem_context.h.

{
  libmesh_assert( !_element_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_element_fe_var[dim].size() ) );
  fe = cast_ptr<FEGenericBase<OutputShape> *>( (_element_fe_var[dim][var] ) );
}

References libMesh::FEMContext::_element_fe_var, dim, and libMesh::libmesh_assert().

◆ get_element_fe() [5/5]

FEBase * libMesh::FEMContext::get_element_fe	(	unsigned int	var,
		unsigned short	dim
	)		const

inlineinherited

Accessor for interior finite element object for scalar-valued variable var for dimension dim.

Definition at line 1314 of file fem_context.h.

{
  libmesh_assert( !_element_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_element_fe_var[dim].size() ) );
  return cast_ptr<FEBase *>( (_element_fe_var[dim][var] ) );
}

References libMesh::FEMContext::_element_fe_var, dim, and libMesh::libmesh_assert().

◆ get_element_qrule() [1/2]

const QBase & libMesh::FEMContext::get_element_qrule ( ) const

inlineinherited

Accessor for element interior quadrature rule for the dimension of the current _elem.

Definition at line 820 of file fem_context.h.

821 { return this->get_element_qrule(this->get_elem_dim()); }

libMesh::FEMContext::get_element_qrule

const QBase & get_element_qrule() const

Accessor for element interior quadrature rule for the dimension of the current _elem.

Definition fem_context.h:820

References libMesh::FEMContext::get_elem_dim(), and libMesh::FEMContext::get_element_qrule().

◆ get_element_qrule() [2/2]

const QBase & libMesh::FEMContext::get_element_qrule ( unsigned short dim ) const

inlineinherited

Accessor for element interior quadrature rule.

Definition at line 833 of file fem_context.h.

834 { libmesh_assert(_element_qrule[dim]);

835 return *(this->_element_qrule[dim]); }

References libMesh::FEMContext::_element_qrule, dim, and libMesh::libmesh_assert().

◆ get_fixed_solution_derivative()

Real libMesh::DiffContext::get_fixed_solution_derivative ( ) const

inlineinherited

The derivative of the current fixed_elem_solution w.r.t.

the unknown solution. Corresponding Jacobian contributions should be multiplied by this amount, or may be skipped if get_fixed_elem_solution_derivative() is 0.

Definition at line 459 of file diff_context.h.

460 { return fixed_solution_derivative; }

libMesh::DiffContext::fixed_solution_derivative

Real fixed_solution_derivative

The derivative of elem_fixed_solution with respect to the nonlinear solution, for use by systems cons...

Definition diff_context.h:517

References libMesh::DiffContext::fixed_solution_derivative.

◆ get_localized_subvector() [1/2]

DenseSubVector< Number > & libMesh::DiffContext::get_localized_subvector	(	const NumericVector< Number > &	localized_vector,
		unsigned int	var
	)

inherited

Return a reference to DenseSubVector localization of localized_vector at variable var contained in the _localized_vectors map.

Definition at line 154 of file diff_context.C.

{
  return _localized_vectors[&localized_vector].second[var];
}

References libMesh::DiffContext::_localized_vectors.

Referenced by libMesh::FEMContext::interior_gradients(), libMesh::FEMContext::interior_hessians(), libMesh::FEMContext::interior_values(), libMesh::FEMContext::side_gradients(), libMesh::FEMContext::side_hessians(), and libMesh::FEMContext::side_values().

◆ get_localized_subvector() [2/2]

const DenseSubVector< Number > & libMesh::DiffContext::get_localized_subvector	(	const NumericVector< Number > &	localized_vector,
		unsigned int	var
	)		const

inherited

const accessible version of get_localized_subvector function

Definition at line 160 of file diff_context.C.

{
  auto localized_vectors_it = _localized_vectors.find(&localized_vector);
  libmesh_assert(localized_vectors_it != _localized_vectors.end());
  return localized_vectors_it->second.second[var];
}

References libMesh::DiffContext::_localized_vectors, and libMesh::libmesh_assert().

◆ get_localized_vector() [1/2]

DenseVector< Number > & libMesh::DiffContext::get_localized_vector ( const NumericVector< Number > & localized_vector )

inherited

Return a reference to DenseVector localization of localized_vector contained in the _localized_vectors map.

Definition at line 140 of file diff_context.C.

{
  return _localized_vectors[&localized_vector].first;
}

References libMesh::DiffContext::_localized_vectors.

◆ get_localized_vector() [2/2]

const DenseVector< Number > & libMesh::DiffContext::get_localized_vector ( const NumericVector< Number > & localized_vector ) const

inherited

const accessible version of get_localized_vector function

Definition at line 146 of file diff_context.C.

{
  auto localized_vectors_it = _localized_vectors.find(&localized_vector);
  libmesh_assert(localized_vectors_it != _localized_vectors.end());
  return localized_vectors_it->second.first;
}

References libMesh::DiffContext::_localized_vectors, and libMesh::libmesh_assert().

◆ get_mesh_system() [1/2]

System * libMesh::FEMContext::get_mesh_system ( )

inlineinherited

Accessor for moving mesh System.

Definition at line 872 of file fem_context.h.

873 { return this->_mesh_sys; }

References libMesh::FEMContext::_mesh_sys.

◆ get_mesh_system() [2/2]

const System * libMesh::FEMContext::get_mesh_system ( ) const

inlineinherited

Accessor for moving mesh System.

Definition at line 866 of file fem_context.h.

867 { return this->_mesh_sys; }

References libMesh::FEMContext::_mesh_sys.

◆ get_mesh_x_var()

unsigned int libMesh::FEMContext::get_mesh_x_var ( ) const

inlineinherited

Accessor for x-variable of moving mesh System.

Definition at line 878 of file fem_context.h.

879 { return _mesh_x_var; }

libMesh::FEMContext::_mesh_x_var

unsigned int _mesh_x_var

Variables from which to acquire moving mesh information.

Definition fem_context.h:1082

References libMesh::FEMContext::_mesh_x_var.

Referenced by libMesh::FEMContext::_do_elem_position_set(), and libMesh::FEMContext::elem_position_get().

◆ get_mesh_y_var()

unsigned int libMesh::FEMContext::get_mesh_y_var ( ) const

inlineinherited

Accessor for y-variable of moving mesh System.

Definition at line 892 of file fem_context.h.

893 { return _mesh_y_var; }

libMesh::FEMContext::_mesh_y_var

unsigned int _mesh_y_var

Definition fem_context.h:1082

References libMesh::FEMContext::_mesh_y_var.

Referenced by libMesh::FEMContext::_do_elem_position_set(), and libMesh::FEMContext::elem_position_get().

◆ get_mesh_z_var()

unsigned int libMesh::FEMContext::get_mesh_z_var ( ) const

inlineinherited

Accessor for z-variable of moving mesh System.

Definition at line 906 of file fem_context.h.

907 { return _mesh_z_var; }

libMesh::FEMContext::_mesh_z_var

unsigned int _mesh_z_var

Definition fem_context.h:1082

References libMesh::FEMContext::_mesh_z_var.

Referenced by libMesh::FEMContext::_do_elem_position_set(), and libMesh::FEMContext::elem_position_get().

◆ get_neighbor()

const Elem & libMesh::DGFEMContext::get_neighbor ( ) const

inline

Accessor for current neighbor Elem object for assembling DG terms.

Definition at line 226 of file dg_fem_context.h.

227 { return *_neighbor; }

References _neighbor.

Referenced by neighbor_side_fe_reinit().

◆ get_neighbor_dof_indices() [1/2]

const std::vector< dof_id_type > & libMesh::DGFEMContext::get_neighbor_dof_indices ( ) const

inline

Accessor for neighbor dof indices.

Definition at line 71 of file dg_fem_context.h.

72 { return _neighbor_dof_indices; }

libMesh::DGFEMContext::_neighbor_dof_indices

std::vector< dof_id_type > _neighbor_dof_indices

Global Degree of freedom index lists for the neighbor element.

Definition dg_fem_context.h:274

References _neighbor_dof_indices.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ get_neighbor_dof_indices() [2/2]

const std::vector< dof_id_type > & libMesh::DGFEMContext::get_neighbor_dof_indices ( unsigned int var ) const

inline

Accessor for element dof indices of a particular variable corresponding to the index argument.

Definition at line 78 of file dg_fem_context.h.

79 { return _neighbor_dof_indices_var[var]; }

References _neighbor_dof_indices_var.

◆ get_neighbor_elem_jacobian() [1/4]

DenseMatrix< Number > & libMesh::DGFEMContext::get_neighbor_elem_jacobian ( )

inline

Non-const accessor for element Jacobian.

Definition at line 168 of file dg_fem_context.h.

169 { return _neighbor_elem_jacobian; }

References _neighbor_elem_jacobian.

◆ get_neighbor_elem_jacobian() [2/4]

const DenseMatrix< Number > & libMesh::DGFEMContext::get_neighbor_elem_jacobian ( ) const

inline

Const accessor for element-neighbor Jacobian.

Definition at line 162 of file dg_fem_context.h.

163 { return _neighbor_elem_jacobian; }

References _neighbor_elem_jacobian.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ get_neighbor_elem_jacobian() [3/4]

DenseSubMatrix< Number > & libMesh::DGFEMContext::get_neighbor_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)

inline

Non-const accessor for neighbor-element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 182 of file dg_fem_context.h.

183 { return *(_neighbor_elem_subjacobians[var1][var2]); }

References _neighbor_elem_subjacobians.

◆ get_neighbor_elem_jacobian() [4/4]

const DenseSubMatrix< Number > & libMesh::DGFEMContext::get_neighbor_elem_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)		const

inline

Const accessor for neighbor-element Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 175 of file dg_fem_context.h.

176 { return *(_neighbor_elem_subjacobians[var1][var2]); }

References _neighbor_elem_subjacobians.

◆ get_neighbor_neighbor_jacobian() [1/4]

DenseMatrix< Number > & libMesh::DGFEMContext::get_neighbor_neighbor_jacobian ( )

inline

Non-const accessor for element Jacobian.

Definition at line 194 of file dg_fem_context.h.

195 { return _neighbor_neighbor_jacobian; }

References _neighbor_neighbor_jacobian.

◆ get_neighbor_neighbor_jacobian() [2/4]

const DenseMatrix< Number > & libMesh::DGFEMContext::get_neighbor_neighbor_jacobian ( ) const

inline

Const accessor for element-neighbor Jacobian.

Definition at line 188 of file dg_fem_context.h.

189 { return _neighbor_neighbor_jacobian; }

References _neighbor_neighbor_jacobian.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ get_neighbor_neighbor_jacobian() [3/4]

DenseSubMatrix< Number > & libMesh::DGFEMContext::get_neighbor_neighbor_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)

inline

Non-const accessor for neighbor-neighbor Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 208 of file dg_fem_context.h.

209 { return *(_neighbor_neighbor_subjacobians[var1][var2]); }

References _neighbor_neighbor_subjacobians.

◆ get_neighbor_neighbor_jacobian() [4/4]

const DenseSubMatrix< Number > & libMesh::DGFEMContext::get_neighbor_neighbor_jacobian	(	unsigned int	var1,
		unsigned int	var2
	)		const

inline

Const accessor for neighbor-neighbor Jacobian of particular variables corresponding to the variable index arguments.

Definition at line 201 of file dg_fem_context.h.

202 { return *(_neighbor_neighbor_subjacobians[var1][var2]); }

References _neighbor_neighbor_subjacobians.

◆ get_neighbor_residual() [1/4]

DenseVector< Number > & libMesh::DGFEMContext::get_neighbor_residual ( )

inline

Non-const accessor for neighbor residual.

Definition at line 90 of file dg_fem_context.h.

91 { return _neighbor_residual; }

References _neighbor_residual.

◆ get_neighbor_residual() [2/4]

const DenseVector< Number > & libMesh::DGFEMContext::get_neighbor_residual ( ) const

inline

Const accessor for neighbor residual.

Definition at line 84 of file dg_fem_context.h.

85 { return _neighbor_residual; }

References _neighbor_residual.

◆ get_neighbor_residual() [3/4]

DenseSubVector< Number > & libMesh::DGFEMContext::get_neighbor_residual ( unsigned int var )

inline

Non-const accessor for neighbor residual of a particular variable corresponding to the variable index argument.

Definition at line 104 of file dg_fem_context.h.

105 { return *(_neighbor_subresiduals[var]); }

References _neighbor_subresiduals.

◆ get_neighbor_residual() [4/4]

const DenseSubVector< Number > & libMesh::DGFEMContext::get_neighbor_residual ( unsigned int var ) const

inline

Const accessor for neighbor residual of a particular variable corresponding to the variable index argument.

Definition at line 97 of file dg_fem_context.h.

98 { return *(_neighbor_subresiduals[var]); }

References _neighbor_subresiduals.

◆ get_neighbor_side_fe()

template<typename OutputShape >

void libMesh::DGFEMContext::get_neighbor_side_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe
	)		const

inline

Accessor for neighbor edge/face (2D/3D) finite element object for variable var.

Definition at line 301 of file dg_fem_context.h.

{
  libmesh_assert_less ( var, _neighbor_side_fe_var.size() );
  fe = cast_ptr<FEGenericBase<OutputShape> *>( _neighbor_side_fe_var[var] );
}

References _neighbor_side_fe_var.

Referenced by assembly_with_dg_fem_context().

◆ get_qoi_derivatives() [1/4]

std::vector< DenseVector< Number > > & libMesh::DiffContext::get_qoi_derivatives ( )

inlineinherited

Non-const accessor for QoI derivatives.

Definition at line 335 of file diff_context.h.

336 { return _elem_qoi_derivative; }

libMesh::DiffContext::_elem_qoi_derivative

std::vector< DenseVector< Number > > _elem_qoi_derivative

Element quantity of interest derivative contributions.

Definition diff_context.h:626

References libMesh::DiffContext::_elem_qoi_derivative.

◆ get_qoi_derivatives() [2/4]

const std::vector< DenseVector< Number > > & libMesh::DiffContext::get_qoi_derivatives ( ) const

inlineinherited

Const accessor for QoI derivatives.

Definition at line 329 of file diff_context.h.

330 { return _elem_qoi_derivative; }

References libMesh::DiffContext::_elem_qoi_derivative.

Referenced by HeatSystem::element_qoi_derivative(), LaplaceSystem::element_qoi_derivative(), LaplaceQoI::element_qoi_derivative(), libMesh::FEMContext::pre_fe_reinit(), LaplaceSystem::side_qoi_derivative(), and CoupledSystemQoI::side_qoi_derivative().

◆ get_qoi_derivatives() [3/4]

DenseSubVector< Number > & libMesh::DiffContext::get_qoi_derivatives	(	std::size_t	qoi,
		unsigned int	var
	)

inlineinherited

Non-const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

Definition at line 353 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_qoi_subderivatives.size(), qoi);
    libmesh_assert_greater(_elem_qoi_subderivatives[qoi].size(), var);
    return _elem_qoi_subderivatives[qoi][var];
  }

References libMesh::DiffContext::_elem_qoi_subderivatives.

◆ get_qoi_derivatives() [4/4]

const DenseSubVector< Number > & libMesh::DiffContext::get_qoi_derivatives	(	std::size_t	qoi,
		unsigned int	var
	)		const

inlineinherited

Const accessor for QoI derivative of a particular qoi and variable corresponding to the index arguments.

Definition at line 342 of file diff_context.h.

  {
    libmesh_assert_greater(_elem_qoi_subderivatives.size(), qoi);
    libmesh_assert_greater(_elem_qoi_subderivatives[qoi].size(), var);
    return _elem_qoi_subderivatives[qoi][var];
  }

References libMesh::DiffContext::_elem_qoi_subderivatives.

◆ get_qois() [1/2]

std::vector< Number > & libMesh::DiffContext::get_qois ( )

inlineinherited

Non-const accessor for QoI vector.

Definition at line 323 of file diff_context.h.

324 { return _elem_qoi; }

libMesh::DiffContext::_elem_qoi

std::vector< Number > _elem_qoi

Element quantity of interest contributions.

Definition diff_context.h:621

References libMesh::DiffContext::_elem_qoi.

◆ get_qois() [2/2]

const std::vector< Number > & libMesh::DiffContext::get_qois ( ) const

inlineinherited

Const accessor for QoI vector.

Definition at line 317 of file diff_context.h.

318 { return _elem_qoi; }

References libMesh::DiffContext::_elem_qoi.

Referenced by LaplaceQoI::element_qoi(), HeatSystem::element_qoi(), and CoupledSystemQoI::side_qoi().

◆ get_side()

unsigned char libMesh::FEMContext::get_side ( ) const

inlineinherited

Accessor for current side of Elem object.

Definition at line 940 of file fem_context.h.

941 { return side; }

libMesh::FEMContext::side

unsigned char side

Current side for side_* to examine.

Definition fem_context.h:1031

References libMesh::FEMContext::side.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), assembly_with_dg_fem_context(), libMesh::FEMContext::side_fe_reinit(), and SolidSystem::side_time_derivative().

◆ get_side_fe() [1/5]

FEBase * libMesh::FEMContext::get_side_fe ( unsigned int var ) const

inlineinherited

Accessor for side finite element object for scalar-valued variable var for the largest dimension in the mesh.

We default to the largest mesh dim if this method is called before the Elem * is set in the FEMContext, e.g. in FEMSystem::init_context (or a subclass).

Definition at line 344 of file fem_context.h.

345 { return this->get_side_fe(var,this->get_elem_dim()); }

References libMesh::FEMContext::get_elem_dim(), and libMesh::FEMContext::get_side_fe().

◆ get_side_fe() [2/5]

void libMesh::FEMContext::get_side_fe	(	unsigned int	var,
		FEAbstract *&	fe,
		unsigned short	dim
	)		const

inlineinherited

Accessor for edge/face (2D/3D) finite element object for variable var for dimension dim.

Definition at line 1332 of file fem_context.h.

{
  libmesh_assert( !_side_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_side_fe_var[dim].size() ) );
  fe = _side_fe_var[dim][var];
}

References libMesh::FEMContext::_side_fe_var, dim, and libMesh::libmesh_assert().

◆ get_side_fe() [3/5]

template<typename OutputShape >

void libMesh::FEMContext::get_side_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe
	)		const

inlineinherited

Accessor for edge/face (2D/3D) finite element object for variable var for the largest dimension in the mesh.

We default to the largest mesh dim if this method is called before the Elem * is set in the FEMContext, e.g. in FEMSystem::init_context (or a subclass).

Definition at line 335 of file fem_context.h.

336 { this->get_side_fe<OutputShape>(var,fe,this->get_elem_dim()); }

References libMesh::FEMContext::get_elem_dim().

◆ get_side_fe() [4/5]

template<typename OutputShape >

void libMesh::FEMContext::get_side_fe	(	unsigned int	var,
		FEGenericBase< OutputShape > *&	fe,
		unsigned short	dim
	)		const

inlineinherited

Accessor for edge/face (2D/3D) finite element object for variable var for dimension dim.

Definition at line 1323 of file fem_context.h.

{
  libmesh_assert( !_side_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_side_fe_var[dim].size() ) );
  fe = cast_ptr<FEGenericBase<OutputShape> *>( (_side_fe_var[dim][var] ) );
}

References libMesh::FEMContext::_side_fe_var, dim, and libMesh::libmesh_assert().

◆ get_side_fe() [5/5]

FEBase * libMesh::FEMContext::get_side_fe	(	unsigned int	var,
		unsigned short	dim
	)		const

inlineinherited

Accessor for side finite element object for scalar-valued variable var for dimension dim.

Definition at line 1341 of file fem_context.h.

{
  libmesh_assert( !_side_fe_var[dim].empty() );
  libmesh_assert_less ( var, (_side_fe_var[dim].size() ) );
  return cast_ptr<FEBase *>( (_side_fe_var[dim][var] ) );
}

References libMesh::FEMContext::_side_fe_var, dim, and libMesh::libmesh_assert().

◆ get_side_qrule() [1/2]

const QBase & libMesh::FEMContext::get_side_qrule ( ) const

inlineinherited

Accessor for element side quadrature rule for the dimension of the current _elem.

Definition at line 827 of file fem_context.h.

828 { return this->get_side_qrule(this->get_elem_dim()); }

libMesh::FEMContext::get_side_qrule

const QBase & get_side_qrule() const

Accessor for element side quadrature rule for the dimension of the current _elem.

Definition fem_context.h:827

References libMesh::FEMContext::get_elem_dim(), and libMesh::FEMContext::get_side_qrule().

Referenced by assembly_with_dg_fem_context(), AssemblyF0::boundary_assembly(), AssemblyF1::boundary_assembly(), AssemblyF2::boundary_assembly(), AssemblyA0::boundary_assembly(), AssemblyA1::boundary_assembly(), AssemblyA2::boundary_assembly(), A2::boundary_assembly(), A3::boundary_assembly(), F0::boundary_assembly(), Output0::boundary_assembly(), libMesh::FEMContext::get_side_qrule(), LaplaceSystem::side_constraint(), LaplaceSystem::side_postprocess(), CoupledSystemQoI::side_qoi(), LaplaceSystem::side_qoi_derivative(), CoupledSystemQoI::side_qoi_derivative(), SolidSystem::side_time_derivative(), ElasticitySystem::side_time_derivative(), and CurlCurlSystem::side_time_derivative().

◆ get_side_qrule() [2/2]

const QBase & libMesh::FEMContext::get_side_qrule ( unsigned short dim ) const

inlineinherited

Accessor for element side quadrature rule.

Definition at line 840 of file fem_context.h.

  {
    libmesh_assert(_side_qrule[dim]);
    return *(this->_side_qrule[dim]);
  }

References libMesh::FEMContext::_side_qrule, dim, and libMesh::libmesh_assert().

◆ get_system()

const System & libMesh::DiffContext::get_system ( ) const

inlineinherited

Accessor for associated system.

Definition at line 106 of file diff_context.h.

107 { return _system; }

libMesh::DiffContext::_system

const System & _system

A reference to the system this context is constructed with.

Definition diff_context.h:660

References libMesh::DiffContext::_system.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), AssemblyA0::boundary_assembly(), AssemblyA1::boundary_assembly(), AssemblyA2::boundary_assembly(), libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), HeatSystem::element_qoi(), HeatSystem::element_qoi_derivative(), HeatSystem::element_time_derivative(), libMesh::FEMContext::find_hardest_fe_type(), CoupledSystemQoI::init_context(), SigmaPhysics::init_context(), libMesh::RBEIMConstruction::init_context(), neighbor_side_fe_reinit(), and libMesh::DifferentiablePhysics::nonlocal_mass_residual().

◆ get_system_time()

Real libMesh::DiffContext::get_system_time ( ) const

inlineinherited

Accessor for the time variable stored in the system class.

Definition at line 411 of file diff_context.h.

412 { return system_time; }

libMesh::DiffContext::system_time

const Real system_time

This is the time stored in the System class at the time this context was created, i....

Definition diff_context.h:490

References libMesh::DiffContext::system_time.

Referenced by libMesh::FEMContext::_update_time_from_system().

◆ get_time()

Real libMesh::DiffContext::get_time ( ) const

inlineinherited

Accessor for the time for which the current nonlinear_solution is defined.

Definition at line 417 of file diff_context.h.

418 { return time; }

libMesh::DiffContext::time

Real time

For time-dependent problems, this is the time t for which the current nonlinear_solution is defined.

Definition diff_context.h:481

References libMesh::DiffContext::time.

◆ has_elem()

bool libMesh::FEMContext::has_elem ( ) const

inlineinherited

Test for current Elem object.

Definition at line 920 of file fem_context.h.

921 { return (this->_elem != nullptr); }

References libMesh::FEMContext::_elem.

Referenced by libMesh::FEMContext::build_new_fe(), libMesh::FEMContext::cached_fe(), libMesh::FEMContext::elem_fe_reinit(), and libMesh::FEMContext::pre_fe_reinit().

◆ has_side_boundary_id()

bool libMesh::FEMContext::has_side_boundary_id ( boundary_id_type id ) const

inherited

Reports if the boundary id is found on the current side.

Definition at line 322 of file fem_context.C.

{
  return _boundary_info.has_boundary_id(&(this->get_elem()), side, id);
}

References libMesh::FEMContext::_boundary_info, libMesh::FEMContext::get_elem(), libMesh::BoundaryInfo::has_boundary_id(), and libMesh::FEMContext::side.

Referenced by A2::boundary_assembly(), A3::boundary_assembly(), F0::boundary_assembly(), Output0::boundary_assembly(), CoupledSystemQoI::side_qoi(), CoupledSystemQoI::side_qoi_derivative(), and ElasticitySystem::side_time_derivative().

◆ init_internal_data()

void libMesh::FEMContext::init_internal_data ( const System & sys )

privateinherited

Helper function used in constructors to set up internal data.

Definition at line 223 of file fem_context.C.

{
  // Reserve space for the FEAbstract and QBase objects for each
  // element dimension possibility (0,1,2,3)
 
  // Note: we would simply resize() all four of these vectors, but
  // some compilers (ICC 19, MSVC) generate a diagnostic about copying
  // a std::unique_ptr in this case, so the following two lines are a
  // workaround.
  _element_fe = std::vector<std::map<FEType, std::unique_ptr<FEAbstract>>>(4);
  _side_fe = std::vector<std::map<FEType, std::unique_ptr<FEAbstract>>>(4);
  _element_fe_var.resize(4);
  _side_fe_var.resize(4);
 
  // We need to know which of our variables has the hardest
  // shape functions to numerically integrate.
 
  unsigned int nv = sys.n_vars();
  libmesh_assert (nv);
 
  bool have_scalar = false;
 
  if (_active_vars)
    {
      for (auto v : *_active_vars)
        if (sys.variable_type(v).family == SCALAR)
          {
            have_scalar = true;
            break;
          }
    }
  else
    {
      for (auto v : make_range(sys.n_vars()))
        if (sys.variable_type(v).family == SCALAR)
          {
            have_scalar = true;
            break;
          }
    }
 
  if (have_scalar)
    // SCALAR FEs have dimension 0 by assumption
    _elem_dims.insert(0);
 
  auto build_var_fe = [this, &sys](unsigned int dim,
                                   unsigned int i)
    {
      FEType fe_type = sys.variable_type(i);
      const bool add_p_level = fe_type.p_refinement;
 
      auto & element_fe = _element_fe[dim][fe_type];
      auto & side_fe = _side_fe[dim][fe_type];
      if (!element_fe)
        {
          element_fe = FEAbstract::build(dim, fe_type);
          element_fe->add_p_level_in_reinit(add_p_level);
          side_fe = FEAbstract::build(dim, fe_type);
          side_fe->add_p_level_in_reinit(add_p_level);
 
          if (dim == 3)
          {
            auto & edge_fe = _edge_fe[fe_type];
            edge_fe = FEAbstract::build(dim, fe_type);
            edge_fe->add_p_level_in_reinit(add_p_level);
          }
        }
 
      _element_fe_var[dim][i] = element_fe.get();
      _side_fe_var[dim][i] = side_fe.get();
      if ((dim) == 3)
        _edge_fe_var[i] = _edge_fe[fe_type].get();
    };
 
  for (const auto & dim : _elem_dims)
    {
      // Create finite element objects
      _element_fe_var[dim].resize(nv);
      _side_fe_var[dim].resize(nv);
      if (dim == 3)
        _edge_fe_var.resize(nv);
 
      if (_active_vars)
        for (auto v : *_active_vars)
          build_var_fe(dim, v);
      else
        for (auto v : make_range(nv))
          build_var_fe(dim, v);
    }
 
  this->use_default_quadrature_rules(_extra_quadrature_order);
}

Referenced by libMesh::FEMContext::FEMContext().

◆ interior_accel()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_accel< Gradient >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The second time derivative (acceleration) of the solution variable var at the quadrature point qp on the current element interior.

Definition at line 1414 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_element_fe
                   <typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution_accel>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution_accel(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_value().

Referenced by ElasticitySystem::mass_residual(), SecondOrderScalarSystemSecondOrderTimeSolverBase::mass_residual(), and SecondOrderScalarSystemFirstOrderTimeSolverBase::mass_residual().

◆ interior_curl()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_curl< Gradient >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	curl_u
	)		const

inherited

Returns: The curl of the solution variable var at the physical point p on the current element.

Definition at line 622 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputShape>> & curl_phi = fe->get_curl_phi();
 
  // Accumulate solution curl
  curl_u = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    curl_u.add_scaled(curl_phi[l][qp], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEGenericBase< OutputType >::get_curl_phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

Referenced by CurlCurlSystem::element_time_derivative().

◆ interior_div()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_div< Number >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	div_u
	)		const

inherited

Returns: The divergence of the solution variable var at the physical point p on the current element.

Definition at line 653 of file fem_context.C.

{
  typedef typename
    TensorTools::IncrementRank
    <typename TensorTools::MakeReal<OutputType>::type>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputDivergence>> & div_phi = fe->get_div_phi();
 
  // Accumulate solution curl
  div_u = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    div_u += div_phi[l][qp] * coef(l);
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEGenericBase< OutputType >::get_div_phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

◆ interior_gradient() [1/2]

Gradient libMesh::FEMContext::interior_gradient	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The gradient of the solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 491 of file fem_context.C.

{
  Gradient du;
 
  this->interior_gradient( var, qp, du );
 
  return du;
}

References libMesh::FEMContext::interior_gradient().

◆ interior_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::interior_gradient	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	du
	)		const

inherited

Returns: The gradient of the solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 504 of file fem_context.C.

{
  this->some_gradient<OutputType,
                      &FEMContext::get_element_fe<typename TensorTools::MakeReal
                                                  <typename TensorTools::DecrementRank
                                                   <OutputType>::type>::type>,
                      &DiffContext::get_elem_solution>(var, qp, du);
}

References libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_gradient().

◆ interior_gradients()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_gradients< Tensor >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	interior_gradients_vector
	)		const

inherited

Fills a vector with the gradient of the solution variable var at all the quadrature points in the current element interior.

Note: This is the preferred API.

Definition at line 518 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal
    <typename TensorTools::DecrementRank<OutputType>::type>::type
    OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputGradient>> & dphi = fe->get_dphi();
 
  // Loop over all the q_points in this finite element
  for (auto qp : index_range(du_vals))
    {
      OutputType & du = du_vals[qp];
 
      // Compute the gradient at this q_point
      du = 0;
 
      for (unsigned int l=0; l != n_dofs; l++)
        du.add_scaled(dphi[l][qp], coef(l));
    }
 
  return;
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), and libMesh::index_range().

◆ interior_hessian() [1/2]

Tensor libMesh::FEMContext::interior_hessian	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The hessian of the solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 556 of file fem_context.C.

{
  Tensor d2u;
 
  this->interior_hessian( var, qp, d2u );
 
  return d2u;
}

References libMesh::FEMContext::interior_hessian().

Referenced by libMesh::FEMContext::interior_hessian().

◆ interior_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::interior_hessian	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	d2u
	)		const

inherited

Returns: The hessian of the solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 566 of file fem_context.C.

{
  this->some_hessian<OutputType,
                     &FEMContext::get_element_fe
                     <typename TensorTools::MakeReal
                      <typename TensorTools::DecrementRank
                       <typename TensorTools::DecrementRank
                        <OutputType>::type>::type>::type>,
                     &DiffContext::get_elem_solution>(var, qp, d2u);
}

References libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_hessian().

◆ interior_hessians()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_hessians< Tensor >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	d2u_vals
	)		const

inherited

Fills a vector of hessians of the _system_vector at the all the quadrature points in the current element interior.

This is the preferred API.

Definition at line 580 of file fem_context.C.

{
  typedef typename TensorTools::DecrementRank<OutputType>::type Rank1Decrement;
  typedef typename TensorTools::DecrementRank<Rank1Decrement>::type Rank2Decrement;
  typedef typename TensorTools::MakeReal<Rank2Decrement>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputTensor>> & d2phi = fe->get_d2phi();
 
  // Loop over all the q_points in this finite element
  for (auto qp : index_range(d2u_vals))
    {
      OutputType & d2u = d2u_vals[qp];
 
      // Compute the gradient at this q_point
      d2u = 0;
 
      for (unsigned int l=0; l != n_dofs; l++)
        d2u.add_scaled(d2phi[l][qp], coef(l));
    }
 
  return;
}

References libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), and libMesh::index_range().

◆ interior_rate()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_rate< Gradient >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The time derivative (rate) of the solution variable var at the quadrature point qp on the current element interior.

Definition at line 1383 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_element_fe
                   <typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution_rate>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution_rate(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_value().

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), SecondOrderScalarSystemSecondOrderTimeSolverBase::damping_residual(), SecondOrderScalarSystemFirstOrderTimeSolverBase::damping_residual(), libMesh::FEMPhysics::mass_residual(), NavierSystem::mass_residual(), and FirstOrderScalarSystemBase::mass_residual().

◆ interior_rate_gradient()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_rate_gradient< Tensor >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The time derivative (rate) of the solution gradient of variable var at the quadrature point qp on the current element interior.

Definition at line 1393 of file fem_context.C.

{
  this->some_gradient<OutputType,
                      &FEMContext::get_element_fe<typename TensorTools::MakeReal
                                                  <typename TensorTools::DecrementRank
                                                   <OutputType>::type>::type>,
                      &DiffContext::get_elem_solution_rate>(var, qp, dudot);
}

References libMesh::DiffContext::get_elem_solution_rate(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::some_gradient().

◆ interior_value() [1/2]

Number libMesh::FEMContext::interior_value	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The value of the solution variable var at the quadrature point qp on the current element interior.

Note: This API is currently present for backward compatibility.

Definition at line 436 of file fem_context.C.

{
  Number u;
 
  this->interior_value( var, qp, u );
 
  return u;
}

References libMesh::FEMContext::interior_value().

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), LaplaceSystem::element_postprocess(), LaplaceQoI::element_qoi(), HeatSystem::element_qoi(), CoupledSystem::element_time_derivative(), NavierSystem::element_time_derivative(), CurlCurlSystem::element_time_derivative(), HilbertSystem::element_time_derivative(), libMesh::FEMContext::interior_value(), SlitMeshRefinedSystemTest::testRestart(), and SlitMeshRefinedSystemTest::testSystem().

◆ interior_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::interior_value	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The value of the solution variable var at the quadrature point qp on the current element interior.

Note: This is the preferred API.

Definition at line 446 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_element_fe<typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution(), and libMesh::FEMContext::some_value().

◆ interior_values()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::interior_values< Gradient >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	interior_values_vector
	)		const

inherited

Fills a vector of values of the _system_vector at the all the quadrature points in the current element interior.

Definition at line 456 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get the finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<OutputShape>> & phi = fe->get_phi();
 
  // Loop over all the q_points on this element
  for (auto qp : index_range(u_vals))
    {
      OutputType & u = u_vals[qp];
 
      // Compute the value at this q_point
      u = 0.;
 
      for (unsigned int l=0; l != n_dofs; l++)
        u += phi[l][qp] * coef(l);
    }
 
  return;
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), libMesh::FEGenericBase< OutputType >::get_phi(), and libMesh::index_range().

Referenced by HeatSystem::element_qoi_derivative().

◆ is_adjoint() [1/2]

bool & libMesh::DiffContext::is_adjoint ( )

inlineinherited

Accessor for setting whether we need to do a primal or adjoint solve.

Definition at line 473 of file diff_context.h.

474 { return _is_adjoint; }

libMesh::DiffContext::_is_adjoint

bool _is_adjoint

Is this context to be used for a primal or adjoint solve?

Definition diff_context.h:665

References libMesh::DiffContext::_is_adjoint.

◆ is_adjoint() [2/2]

bool libMesh::DiffContext::is_adjoint ( ) const

inlineinherited

Accessor for querying whether we need to do a primal or adjoint solve.

Definition at line 466 of file diff_context.h.

467 { return _is_adjoint; }

References libMesh::DiffContext::_is_adjoint.

Referenced by SigmaPhysics::init_context().

◆ n_dof_indices() [1/2]

unsigned int libMesh::DiffContext::n_dof_indices ( ) const

inlineinherited

Total number of dof indices on the element.

Definition at line 395 of file diff_context.h.

396 { return cast_int<unsigned int>(_dof_indices.size()); }

References libMesh::DiffContext::_dof_indices.

◆ n_dof_indices() [2/2]

unsigned int libMesh::DiffContext::n_dof_indices ( unsigned int var ) const

inlineinherited

Total number of dof indices of the particular variable corresponding to the index argument.

Definition at line 402 of file diff_context.h.

  {
    libmesh_assert_greater(_dof_indices_var.size(), var);
    return cast_int<unsigned int>(_dof_indices_var[var].size());
  }

References libMesh::DiffContext::_dof_indices_var.

◆ n_vars()

unsigned int libMesh::DiffContext::n_vars ( ) const

inlineinherited

Number of variables in solution.

Definition at line 100 of file diff_context.h.

101 { return cast_int<unsigned int>(_dof_indices_var.size()); }

References libMesh::DiffContext::_dof_indices_var.

Referenced by libMesh::FirstOrderUnsteadySolver::compute_second_order_eqns(), libMesh::FEMPhysics::eulerian_residual(), libMesh::WrappedFunctor< Output >::init_context(), libMesh::DiscontinuityMeasure::init_context(), libMesh::LaplacianErrorEstimator::init_context(), libMesh::KellyErrorEstimator::init_context(), libMesh::FEMPhysics::mass_residual(), neighbor_side_fe_reinit(), libMesh::DifferentiablePhysics::nonlocal_mass_residual(), and libMesh::FEMSystem::numerical_jacobian().

◆ neighbor_side_fe_reinit()

void libMesh::DGFEMContext::neighbor_side_fe_reinit ( )

Initialize neighbor side data needed to assemble DG terms.

The neighbor element is determined by the current value of get_neighbor().

Definition at line 84 of file dg_fem_context.C.

{
  // Call this *after* side_fe_reinit
 
  // Initialize all the neighbor side FE objects based on inverse mapping
  // the quadrature points on the current side
  std::vector<Point> qface_side_points;
  std::vector<Point> qface_neighbor_points;
  for (auto & [neighbor_side_fe_type, fe] : _neighbor_side_fe)
    {
      FEAbstract * side_fe = _side_fe[this->get_dim()][neighbor_side_fe_type].get();
      qface_side_points = side_fe->get_xyz();
 
      FEMap::inverse_map (this->get_dim(), &get_neighbor(),
                          qface_side_points, qface_neighbor_points);
 
      fe->reinit(&get_neighbor(), &qface_neighbor_points);
    }
 
  // Set boolean flag to indicate that the DG terms are active on this element
  _dg_terms_active = true;
 
  // Also, initialize data required for DG assembly on the current side,
  // analogously to FEMContext::pre_fe_reinit
 
  // Initialize the per-element data for elem.
  get_system().get_dof_map().dof_indices (&get_neighbor(), _neighbor_dof_indices);
 
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices().size());
  const unsigned int n_neighbor_dofs = cast_int<unsigned int>
    (_neighbor_dof_indices.size());
 
  // These resize calls also zero out the residual and jacobian
  _neighbor_residual.resize(n_neighbor_dofs);
  _elem_elem_jacobian.resize(n_dofs, n_dofs);
  _elem_neighbor_jacobian.resize(n_dofs, n_neighbor_dofs);
  _neighbor_elem_jacobian.resize(n_neighbor_dofs, n_dofs);
  _neighbor_neighbor_jacobian.resize(n_neighbor_dofs, n_neighbor_dofs);
 
  // Initialize the per-variable data for elem.
  {
    unsigned int sub_dofs = 0;
    for (auto i : make_range(get_system().n_vars()))
      {
        get_system().get_dof_map().dof_indices (&get_neighbor(), _neighbor_dof_indices_var[i], i);
 
        const unsigned int n_dofs_var = cast_int<unsigned int>
          (_neighbor_dof_indices_var[i].size());
 
        _neighbor_subresiduals[i]->reposition
          (sub_dofs, n_dofs_var);
 
        for (unsigned int j=0; j != i; ++j)
          {
            const unsigned int n_dofs_var_j =
              cast_int<unsigned int>
              (this->get_dof_indices(j).size());
 
            _elem_elem_subjacobians[i][j]->reposition
              (sub_dofs, _neighbor_subresiduals[j]->i_off(),
               n_dofs_var, n_dofs_var_j);
            _elem_elem_subjacobians[j][i]->reposition
              (_neighbor_subresiduals[j]->i_off(), sub_dofs,
               n_dofs_var_j, n_dofs_var);
 
            _elem_neighbor_subjacobians[i][j]->reposition
              (sub_dofs, _neighbor_subresiduals[j]->i_off(),
               n_dofs_var, n_dofs_var_j);
            _elem_neighbor_subjacobians[j][i]->reposition
              (_neighbor_subresiduals[j]->i_off(), sub_dofs,
               n_dofs_var_j, n_dofs_var);
 
            _neighbor_elem_subjacobians[i][j]->reposition
              (sub_dofs, _neighbor_subresiduals[j]->i_off(),
               n_dofs_var, n_dofs_var_j);
            _neighbor_elem_subjacobians[j][i]->reposition
              (_neighbor_subresiduals[j]->i_off(), sub_dofs,
               n_dofs_var_j, n_dofs_var);
 
            _neighbor_neighbor_subjacobians[i][j]->reposition
              (sub_dofs, _neighbor_subresiduals[j]->i_off(),
               n_dofs_var, n_dofs_var_j);
            _neighbor_neighbor_subjacobians[j][i]->reposition
              (_neighbor_subresiduals[j]->i_off(), sub_dofs,
               n_dofs_var_j, n_dofs_var);
          }
        _elem_elem_subjacobians[i][i]->reposition
          (sub_dofs, sub_dofs,
           n_dofs_var,
           n_dofs_var);
        _elem_neighbor_subjacobians[i][i]->reposition
          (sub_dofs, sub_dofs,
           n_dofs_var,
           n_dofs_var);
        _neighbor_elem_subjacobians[i][i]->reposition
          (sub_dofs, sub_dofs,
           n_dofs_var,
           n_dofs_var);
        _neighbor_neighbor_subjacobians[i][i]->reposition
          (sub_dofs, sub_dofs,
           n_dofs_var,
           n_dofs_var);
        sub_dofs += n_dofs_var;
      }
    libmesh_assert_equal_to (sub_dofs, n_dofs);
  }
 
}

References _dg_terms_active, _elem_elem_jacobian, _elem_elem_subjacobians, _elem_neighbor_jacobian, _elem_neighbor_subjacobians, _neighbor_dof_indices, _neighbor_dof_indices_var, _neighbor_elem_jacobian, _neighbor_elem_subjacobians, _neighbor_neighbor_jacobian, _neighbor_neighbor_subjacobians, _neighbor_residual, _neighbor_side_fe, _neighbor_subresiduals, libMesh::FEMContext::_side_fe, libMesh::DofMap::dof_indices(), libMesh::FEMContext::get_dim(), libMesh::DiffContext::get_dof_indices(), libMesh::System::get_dof_map(), get_neighbor(), libMesh::DiffContext::get_system(), libMesh::FEAbstract::get_xyz(), libMesh::FEMap::inverse_map(), libMesh::make_range(), libMesh::DiffContext::n_vars(), libMesh::DenseVector< T >::resize(), and libMesh::DenseMatrix< T >::resize().

Referenced by assembly_with_dg_fem_context().

◆ nonlocal_reinit()

void libMesh::FEMContext::nonlocal_reinit ( Real theta )

overridevirtualinherited

Gives derived classes the opportunity to reinitialize data needed for nonlocal calculations at a new point within a timestep.

Reimplemented from libMesh::DiffContext.

Definition at line 1491 of file fem_context.C.

{
  // Update the "time" variable of this context object
  this->_update_time_from_system(theta);
 
  // We can reuse the Elem FE safely here.
  elem_fe_reinit();
}

References libMesh::FEMContext::_update_time_from_system(), and libMesh::FEMContext::elem_fe_reinit().

◆ point_accel()

template<typename OutputType >

void libMesh::FEMContext::point_accel	(	unsigned int	var,
		const Point &	p,
		OutputType &	u
	)		const

inherited

Returns: The second time derivative (acceleration) of the solution variable var at the physical point p on the current element.

◆ point_curl()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::point_curl< Gradient >	(	unsigned int	var,
		const Point &	p,
		OutputType &	curl_u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The curl of the solution variable var at the physical point p on the current element.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 1047 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 3 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputShape>> &  curl_phi = fe_new->get_curl_phi();
 
  curl_u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    curl_u.add_scaled(curl_phi[l][0], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::FEGenericBase< OutputType >::get_curl_phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

◆ point_gradient() [1/2]

Gradient libMesh::FEMContext::point_gradient	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The gradient of the solution variable var at the physical point p on the current element.

Note: This API is currently present for backward compatibility.

Definition at line 944 of file fem_context.C.

{
  Gradient grad_u;
 
  this->point_gradient( var, p, grad_u );
 
  return grad_u;
}

References libMesh::FEMContext::point_gradient().

Referenced by libMesh::ParsedFEMFunction< Output >::eval_args(), CoupledFEMFunctionsx::operator()(), CoupledFEMFunctionsy::operator()(), and libMesh::FEMContext::point_gradient().

◆ point_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::point_gradient	(	unsigned int	var,
		const Point &	p,
		OutputType &	grad_u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The gradient of the solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 956 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal
    <typename TensorTools::DecrementRank<OutputType>::type>::type
    OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 1 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputGradient>> &  dphi = fe_new->get_dphi();
 
  grad_u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    grad_u.add_scaled(dphi[l][0], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

◆ point_hessian() [1/2]

Tensor libMesh::FEMContext::point_hessian	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The hessian of the solution variable var at the physical point p on the current element.

Note: This API currently present for backward compatibility.

Definition at line 996 of file fem_context.C.

{
  Tensor hess_u;
 
  this->point_hessian( var, p, hess_u );
 
  return hess_u;
}

References libMesh::FEMContext::point_hessian().

Referenced by libMesh::ParsedFEMFunction< Output >::eval_args(), and libMesh::FEMContext::point_hessian().

◆ point_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::point_hessian	(	unsigned int	var,
		const Point &	p,
		OutputType &	hess_u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The hessian of the solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 1007 of file fem_context.C.

{
  typedef typename TensorTools::DecrementRank<OutputType>::type Rank1Decrement;
  typedef typename TensorTools::DecrementRank<Rank1Decrement>::type Rank2Decrement;
  typedef typename TensorTools::MakeReal<Rank2Decrement>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 2 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputTensor>> &  d2phi = fe_new->get_d2phi();
 
  hess_u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    hess_u.add_scaled(d2phi[l][0], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

◆ point_rate()

template<typename OutputType >

void libMesh::FEMContext::point_rate	(	unsigned int	var,
		const Point &	p,
		OutputType &	u
	)		const

inherited

Returns: The time derivative (rate) of the solution variable var at the physical point p on the current element.

◆ point_value() [1/2]

Number libMesh::FEMContext::point_value	(	unsigned int	var,
		const Point &	p
	)		const

inherited

Returns: The value of the solution variable var at the physical point p on the current element.

Note: This API is currently present for backward compatibility.

Definition at line 898 of file fem_context.C.

{
  Number u = 0.;
 
  this->point_value( var, p, u );
 
  return u;
}

References libMesh::FEMContext::point_value().

Referenced by libMesh::ParsedFEMFunction< Output >::eval_args(), CoupledFEMFunctionsx::operator()(), CoupledFEMFunctionsy::operator()(), libMesh::FEMContext::point_value(), and NavierSystem::side_constraint().

◆ point_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::point_value	(	unsigned int	var,
		const Point &	p,
		OutputType &	u,
		const Real	tolerance = `TOLERANCE`
	)		const

inherited

Returns: The value of the solution variable var at the physical point p on the current element.

Note: This is the preferred API.

Allows evaluation of points within a relative tolerance outside the element.

Definition at line 908 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * fe = nullptr;
  this->get_element_fe<OutputShape>( var, fe, this->get_elem_dim() );
 
  // Build a FE for calculating u(p)
  FEGenericBase<OutputShape> * fe_new =
    this->build_new_fe( fe, p, tolerance, 0 );
 
  // Get the values of the shape function derivatives
  const std::vector<std::vector<OutputShape>> &  phi = fe_new->get_phi();
 
  u = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    u += phi[l][0] * coef(l);
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::FEMContext::build_new_fe(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_elem_solution(), and libMesh::FEGenericBase< OutputType >::get_phi().

◆ pre_fe_reinit()

void libMesh::FEMContext::pre_fe_reinit	(	const System &	sys,
		const Elem *	e
	)

virtualinherited

Reinitializes local data vectors/matrices on the current geometric element.

Definition at line 1708 of file fem_context.C.

{
  this->set_elem(e);
 
  if (algebraic_type() == CURRENT ||
      algebraic_type() == DOFS_ONLY)
    {
      // Initialize the per-element data for elem.
      if (this->has_elem())
        sys.get_dof_map().dof_indices (&(this->get_elem()), this->get_dof_indices());
      else
        // If !this->has_elem(), then we assume we are dealing with a SCALAR variable
        sys.get_dof_map().dof_indices
          (static_cast<Elem*>(nullptr), this->get_dof_indices());
    }
#ifdef LIBMESH_ENABLE_AMR
  else if (algebraic_type() == OLD ||
           algebraic_type() == OLD_DOFS_ONLY)
    {
      // Initialize the per-element data for elem.
      if (this->has_elem())
        sys.get_dof_map().old_dof_indices (&(this->get_elem()), this->get_dof_indices());
      else
        // If !this->has_elem(), then we assume we are dealing with a SCALAR variable
        sys.get_dof_map().old_dof_indices
          (static_cast<Elem*>(nullptr), this->get_dof_indices());
    }
#endif // LIBMESH_ENABLE_AMR
 
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices().size());
  const unsigned int n_qoi = sys.n_qois();
 
  if (this->algebraic_type() != NONE &&
      this->algebraic_type() != DOFS_ONLY &&
      this->algebraic_type() != OLD_DOFS_ONLY)
    {
      // This also resizes elem_solution
      if (_custom_solution == nullptr)
        sys.current_local_solution->get(this->get_dof_indices(), this->get_elem_solution().get_values());
      else
        _custom_solution->get(this->get_dof_indices(), this->get_elem_solution().get_values());
 
      if (sys.use_fixed_solution)
        this->get_elem_fixed_solution().resize(n_dofs);
 
      // Only make space for these if we're using DiffSystem
      // This is assuming *only* DiffSystem is using elem_solution_rate/accel
      const DifferentiableSystem * diff_system = dynamic_cast<const DifferentiableSystem *>(&sys);
      if (diff_system)
        {
          // Now, we only need these if the solver is unsteady
          if (!diff_system->get_time_solver().is_steady())
            {
              this->get_elem_solution_rate().resize(n_dofs);
 
              // We only need accel space if the TimeSolver is second order
              const UnsteadySolver & time_solver = cast_ref<const UnsteadySolver &>(diff_system->get_time_solver());
 
              if (time_solver.time_order() >= 2 || !diff_system->get_second_order_vars().empty())
                this->get_elem_solution_accel().resize(n_dofs);
            }
        }
 
      if (algebraic_type() != OLD)
        {
          // These resize calls also zero out the residual and jacobian
          this->get_elem_residual().resize(n_dofs);
          if (this->_have_local_matrices)
            this->get_elem_jacobian().resize(n_dofs, n_dofs);
 
          this->get_qoi_derivatives().resize(n_qoi);
          this->_elem_qoi_subderivatives.resize(n_qoi);
          for (std::size_t q=0; q != n_qoi; ++q)
            (this->get_qoi_derivatives())[q].resize(n_dofs);
        }
    }
 
  // Initialize the per-variable data for elem.
  {
    unsigned int sub_dofs = 0;
    for (auto i : make_range(sys.n_vars()))
      {
        if (algebraic_type() == CURRENT ||
            algebraic_type() == DOFS_ONLY)
          {
            if (this->has_elem())
              sys.get_dof_map().dof_indices (&(this->get_elem()), this->get_dof_indices(i), i);
            else
              // If !this->has_elem(), then we assume we are dealing with a SCALAR variable
              sys.get_dof_map().dof_indices
                (static_cast<Elem*>(nullptr), this->get_dof_indices(i), i);
          }
#ifdef LIBMESH_ENABLE_AMR
        else if (algebraic_type() == OLD ||
                 algebraic_type() == OLD_DOFS_ONLY)
          {
            if (this->has_elem())
              sys.get_dof_map().old_dof_indices (&(this->get_elem()), this->get_dof_indices(i), i);
            else
              // If !this->has_elem(), then we assume we are dealing with a SCALAR variable
              sys.get_dof_map().old_dof_indices
                (static_cast<Elem*>(nullptr), this->get_dof_indices(i), i);
          }
#endif // LIBMESH_ENABLE_AMR
 
        if (this->algebraic_type() != NONE &&
            this->algebraic_type() != DOFS_ONLY &&
            this->algebraic_type() != OLD_DOFS_ONLY)
          {
            const unsigned int n_dofs_var = cast_int<unsigned int>
              (this->get_dof_indices(i).size());
 
 
            if (!_active_vars ||
                std::binary_search(_active_vars->begin(),
                                   _active_vars->end(), i))
              {
                this->get_elem_solution(i).reposition
                  (sub_dofs, n_dofs_var);
 
                // Only make space for these if we're using DiffSystem
                // This is assuming *only* DiffSystem is using elem_solution_rate/accel
                const DifferentiableSystem * diff_system = dynamic_cast<const DifferentiableSystem *>(&sys);
                if (diff_system)
                  {
                    // Now, we only need these if the solver is unsteady
                    if (!diff_system->get_time_solver().is_steady())
                      {
                        this->get_elem_solution_rate(i).reposition
                          (sub_dofs, n_dofs_var);
 
                        // We only need accel space if the TimeSolver is second order
                        const UnsteadySolver & time_solver = cast_ref<const UnsteadySolver &>(diff_system->get_time_solver());
 
                        if (time_solver.time_order() >= 2 || !diff_system->get_second_order_vars().empty())
                          this->get_elem_solution_accel(i).reposition
                            (sub_dofs, n_dofs_var);
                      }
                  }
 
                if (sys.use_fixed_solution)
                  this->get_elem_fixed_solution(i).reposition
                    (sub_dofs, n_dofs_var);
 
                if (algebraic_type() != OLD)
                  {
                    this->get_elem_residual(i).reposition
                      (sub_dofs, n_dofs_var);
 
                    for (std::size_t q=0; q != n_qoi; ++q)
                      this->get_qoi_derivatives(q,i).reposition
                        (sub_dofs, n_dofs_var);
 
                    if (this->_have_local_matrices)
                      {
                        for (unsigned int j=0; j != i; ++j)
                          {
                            const unsigned int n_dofs_var_j =
                              cast_int<unsigned int>
                              (this->get_dof_indices(j).size());
 
                            this->get_elem_jacobian(i,j).reposition
                              (sub_dofs, this->get_elem_residual(j).i_off(),
                               n_dofs_var, n_dofs_var_j);
                            this->get_elem_jacobian(j,i).reposition
                              (this->get_elem_residual(j).i_off(), sub_dofs,
                               n_dofs_var_j, n_dofs_var);
                          }
                        this->get_elem_jacobian(i,i).reposition
                          (sub_dofs, sub_dofs,
                           n_dofs_var,
                           n_dofs_var);
                      }
                  }
              }
 
            sub_dofs += n_dofs_var;
          }
      }
 
    if (this->algebraic_type() != NONE &&
        this->algebraic_type() != DOFS_ONLY &&
        this->algebraic_type() != OLD &&
        this->algebraic_type() != OLD_DOFS_ONLY)
      libmesh_assert_equal_to (sub_dofs, n_dofs);
  }
 
  // Now do the localization for the user requested vectors
  if (this->algebraic_type() != NONE &&
      this->algebraic_type() != DOFS_ONLY &&
      this->algebraic_type() != OLD_DOFS_ONLY)
    {
      DiffContext::localized_vectors_iterator localized_vec_it = this->_localized_vectors.begin();
      const DiffContext::localized_vectors_iterator localized_vec_end = this->_localized_vectors.end();
 
      for (; localized_vec_it != localized_vec_end; ++localized_vec_it)
        {
          const NumericVector<Number> & current_localized_vector = *localized_vec_it->first;
          DenseVector<Number> & target_vector = localized_vec_it->second.first;
 
          current_localized_vector.get(this->get_dof_indices(), target_vector.get_values());
 
          // Initialize the per-variable data for elem.
          unsigned int sub_dofs = 0;
          auto init_localized_var_data = [this, localized_vec_it, &sub_dofs](unsigned int i)
            {
              const unsigned int n_dofs_var = cast_int<unsigned int>
                (this->get_dof_indices(i).size());
 
              // This is redundant with earlier initialization, isn't it? - RHS
              // sys.get_dof_map().dof_indices (&(this->get_elem()), this->get_dof_indices(i), i);
 
              localized_vec_it->second.second[i].reposition
                (sub_dofs, n_dofs_var);
 
              sub_dofs += n_dofs_var;
            };
 
          if (_active_vars)
            for (auto v : *_active_vars)
              init_localized_var_data(v);
          else
            for (auto v : make_range(sys.n_vars()))
              init_localized_var_data(v);
 
          libmesh_assert_equal_to (sub_dofs, n_dofs);
        }
    }
}

Referenced by libMesh::ExactSolution::_compute_error(), libMesh::RBConstruction::add_scaled_matrix_and_vector(), libMesh::FEMSystem::assembly(), assembly_with_dg_fem_context(), libMesh::VariationalSmootherSystem::compute_mesh_quality_info(), HilbertSystem::element_time_derivative(), libMesh::RBEIMConstruction::enrich_eim_approximation_on_interiors(), libMesh::RBEIMConstruction::enrich_eim_approximation_on_sides(), libMesh::RBEIMConstruction::initialize_qp_data(), libMesh::FEMSystem::mesh_position_get(), libMesh::FEMSystem::mesh_position_set(), Integrate::operator()(), libMesh::WeightedPatchRecoveryErrorEstimator::EstimateError::operator()(), libMesh::RBParametrizedFunction::preevaluate_parametrized_function_on_mesh(), libMesh::RBParametrizedFunction::preevaluate_parametrized_function_on_mesh_sides(), libMesh::VariationalSmootherSystem::prepare_for_smoothing(), libMesh::System::project_vector(), OverlappingAlgebraicGhostingTest::run_ghosting_test(), OverlappingCouplingGhostingTest::run_sparsity_pattern_test(), SlitMeshRefinedSystemTest::testRestart(), and SlitMeshRefinedSystemTest::testSystem().

◆ set_algebraic_type()

void libMesh::FEMContext::set_algebraic_type ( const AlgebraicType atype )

inlineinherited

Setting which determines whether to initialize algebraic structures (elem_*) on each element and set their values from current_local_solution.

Algebraic initialization may be disabled for efficiency in cases where FEMContext is only used as a convenient container of FE objects.

Definition at line 1004 of file fem_context.h.

1005 { _atype = atype; }

References libMesh::FEMContext::_atype.

Referenced by libMesh::OldSolutionBase< Output, point_output >::init_context().

◆ set_custom_solution()

void libMesh::FEMContext::set_custom_solution ( const NumericVector< Number > * custom_sol )

inlineinherited

Set a NumericVector to be used in place of current_local_solution for calculating elem_solution.

Set to nullptr to restore the current_local_solution behavior. Advanced DifferentiableSystem specific capabilities will only be enabled in the current_local_solution case.

Definition at line 1019 of file fem_context.h.

1020 { _custom_solution = custom_sol; }

References libMesh::FEMContext::_custom_solution.

◆ set_deltat_pointer()

void libMesh::DiffContext::set_deltat_pointer ( Real * dt )

inherited

Points the _deltat member of this class at a timestep value stored in the creating System, for example DiffSystem::deltat.

Definition at line 109 of file diff_context.C.

{
  // We may actually want to be able to set this pointer to nullptr, so
  // don't report an error for that.
  _deltat = dt;
}

References libMesh::DiffContext::_deltat.

Referenced by libMesh::FEMSystem::init_context().

◆ set_elem()

void libMesh::FEMContext::set_elem ( const Elem * e )

protectedinherited

Helper function to promote accessor usage.

Definition at line 1939 of file fem_context.C.

{
  this->_elem = e;
 
  // If e is nullptr, we assume it's SCALAR and set _elem_dim to 0.
  this->_elem_dim =
    cast_int<unsigned char>(this->_elem ? this->_elem->dim() : 0);
}

References libMesh::FEMContext::_elem, libMesh::FEMContext::_elem_dim, and libMesh::Elem::dim().

Referenced by libMesh::FEMContext::pre_fe_reinit().

◆ set_jacobian_tolerance()

void libMesh::FEMContext::set_jacobian_tolerance ( Real tol )

inherited

Calls set_jacobian_tolerance() on all the FE objects controlled by this class.

(Actually, it calls this on the underlying)

Definition at line 1611 of file fem_context.C.

{
  for (auto & m : _element_fe)
    for (auto & pr : m)
      pr.second->get_fe_map().set_jacobian_tolerance(tol);
 
  for (auto & m : _side_fe)
    for (auto & pr : m)
      pr.second->get_fe_map().set_jacobian_tolerance(tol);
 
  for (auto & pr : _edge_fe)
    pr.second->get_fe_map().set_jacobian_tolerance(tol);
}

References libMesh::FEMContext::_edge_fe, libMesh::FEMContext::_element_fe, and libMesh::FEMContext::_side_fe.

◆ set_mesh_system()

virtual void libMesh::FEMContext::set_mesh_system ( System * sys )

inlinevirtualinherited

Tells the FEMContext that system sys contains the isoparametric Lagrangian variables which correspond to the coordinates of mesh nodes, in problems where the mesh itself is expected to move in time.

This should be set automatically if the FEMPhysics requires it.

Definition at line 860 of file fem_context.h.

861 { this->_mesh_sys = sys; }

References libMesh::FEMContext::_mesh_sys.

◆ set_mesh_x_var()

void libMesh::FEMContext::set_mesh_x_var ( unsigned int x_var )

inlineinherited

Accessor for x-variable of moving mesh System.

This should be set automatically if the FEMPhysics requires it.

Definition at line 886 of file fem_context.h.

887 { _mesh_x_var = x_var; }

References libMesh::FEMContext::_mesh_x_var.

◆ set_mesh_y_var()

void libMesh::FEMContext::set_mesh_y_var ( unsigned int y_var )

inlineinherited

Accessor for y-variable of moving mesh System.

This should be set automatically if the FEMPhysics requires it.

Definition at line 900 of file fem_context.h.

901 { _mesh_y_var = y_var; }

References libMesh::FEMContext::_mesh_y_var.

◆ set_mesh_z_var()

void libMesh::FEMContext::set_mesh_z_var ( unsigned int z_var )

inlineinherited

Accessor for z-variable of moving mesh System.

This should be set automatically if the FEMPhysics requires it.

Definition at line 914 of file fem_context.h.

915 { _mesh_z_var = z_var; }

References libMesh::FEMContext::_mesh_z_var.

◆ set_neighbor()

void libMesh::DGFEMContext::set_neighbor ( const Elem & neighbor )

inline

Set the neighbor element which we will use to assemble DG terms.

Note: We do not assume that this element is get_elem().neighbor(side) because we also need to be able to handle the special case of DG terms on "cracks" in a mesh to model certain types of interface conditions. In this case, we need to be able to specify the neighbor element manually. Also, this should give us more flexibility to handle non-conforming meshes.

Definition at line 220 of file dg_fem_context.h.

221 { _neighbor = &neighbor; }

References _neighbor.

Referenced by assembly_with_dg_fem_context().

◆ set_time()

void libMesh::DiffContext::set_time ( Real time_in )

inlineinherited

Set the time for which the current nonlinear_solution is defined.

Definition at line 423 of file diff_context.h.

424 { time = time_in; }

References libMesh::DiffContext::time.

Referenced by libMesh::FEMContext::_update_time_from_system().

◆ side_accel()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::side_accel< Gradient >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The second time derivative (acceleration) of the solution variable var at the quadrature point qp on the current element side.

Definition at line 1426 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_side_fe
                   <typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution_accel>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution_accel(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_value().

◆ side_boundary_ids()

void libMesh::FEMContext::side_boundary_ids ( std::vector< boundary_id_type > & vec_to_fill ) const

inherited

As above, but fills in the std::set provided by the user.

Definition at line 329 of file fem_context.C.

{
  _boundary_info.boundary_ids(&(this->get_elem()), side, vec_to_fill);
}

References libMesh::FEMContext::_boundary_info, libMesh::BoundaryInfo::boundary_ids(), libMesh::FEMContext::get_elem(), and libMesh::FEMContext::side.

◆ side_fe_reinit()

void libMesh::DGFEMContext::side_fe_reinit ( )

overridevirtual

Override side_fe_reinit to set a boolean flag so that by default DG terms are assumed to be inactive.

DG terms are only active if neighbor_side_fe_reinit is called.

Reimplemented from libMesh::FEMContext.

Definition at line 75 of file dg_fem_context.C.

{
  FEMContext::side_fe_reinit();
 
  // By default we assume that the DG terms are inactive
  // They are only active if neighbor_side_fe_reinit is called
  _dg_terms_active = false;
}

References _dg_terms_active, and libMesh::FEMContext::side_fe_reinit().

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), and assembly_with_dg_fem_context().

◆ side_gradient() [1/2]

Gradient libMesh::FEMContext::side_gradient	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The gradient of the solution variable var at the quadrature point qp on the current element side.

Note: This API is currently present for backward compatibility.

Definition at line 743 of file fem_context.C.

{
  Gradient du;
 
  this->side_gradient( var, qp, du );
 
  return du;
}

References libMesh::FEMContext::side_gradient().

Referenced by libMesh::FEMContext::side_gradient(), and LaplaceSystem::side_postprocess().

◆ side_gradient() [2/2]

template<typename OutputType >

void libMesh::FEMContext::side_gradient	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	du
	)		const

inherited

Returns: The gradient of the solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 754 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal
    <typename TensorTools::DecrementRank<OutputType>::type>::type
    OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  libmesh_assert_greater (this->_elem_subsolutions.size(), var);
  const DenseSubVector<Number> & coef = this->get_elem_solution(var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * the_side_fe = nullptr;
  this->get_side_fe<OutputShape>( var, the_side_fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputGradient>> & dphi = the_side_fe->get_dphi();
 
  // Accumulate solution derivatives
  du = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    du.add_scaled(dphi[l][qp], coef(l));
 
  return;
}

References libMesh::DiffContext::_elem_subsolutions, libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem_dim(), and libMesh::DiffContext::get_elem_solution().

◆ side_gradients()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::side_gradients< Tensor >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	side_gradients_vector
	)		const

inherited

Fills a vector with the gradient of the solution variable var at all the quadrature points on the current element side.

Note: This is the preferred API.

Definition at line 788 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal
    <typename TensorTools::DecrementRank<OutputType>::type>::type
    OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * the_side_fe = nullptr;
  this->get_side_fe<OutputShape>( var, the_side_fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputGradient>> & dphi = the_side_fe->get_dphi();
 
  // Loop over all the q_points in this finite element
  for (auto qp : index_range(du_vals))
    {
      OutputType & du = du_vals[qp];
 
      du = 0;
 
      // Compute the gradient at this q_point
      for (unsigned int l=0; l != n_dofs; l++)
        du.add_scaled(dphi[l][qp], coef(l));
    }
 
  return;
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), and libMesh::index_range().

◆ side_hessian() [1/2]

Tensor libMesh::FEMContext::side_hessian	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The hessian of the solution variable var at the quadrature point qp on the current element side.

Note: This API is currently present for backward compatibility.

Definition at line 826 of file fem_context.C.

{
  Tensor d2u;
 
  this->side_hessian( var, qp, d2u );
 
  return d2u;
}

References libMesh::FEMContext::side_hessian().

Referenced by libMesh::FEMContext::side_hessian().

◆ side_hessian() [2/2]

template<typename OutputType >

void libMesh::FEMContext::side_hessian	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	d2u
	)		const

inherited

Returns: The hessian of the solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 839 of file fem_context.C.

{
  this->some_hessian<OutputType,
                     &FEMContext::get_side_fe
                     <typename TensorTools::MakeReal
                      <typename TensorTools::DecrementRank
                       <typename TensorTools::DecrementRank
                        <OutputType>::type>::type>::type>,
                     &DiffContext::get_elem_solution>(var, qp, d2u);
}

References libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_hessian().

◆ side_hessians()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::side_hessians< Tensor >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	d2u_vals
	)		const

inherited

Fills a vector of hessians of the _system_vector at the all the quadrature points on the current element side.

This is the preferred API.

Definition at line 855 of file fem_context.C.

{
  typedef typename TensorTools::DecrementRank<OutputType>::type Rank1Decrement;
  typedef typename TensorTools::DecrementRank<Rank1Decrement>::type Rank2Decrement;
  typedef typename TensorTools::MakeReal<Rank2Decrement>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get finite element object
  FEGenericBase<OutputShape> * the_side_fe = nullptr;
  this->get_side_fe<OutputShape>( var, the_side_fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<typename FEGenericBase<OutputShape>::OutputTensor>> & d2phi = the_side_fe->get_d2phi();
 
  // Loop over all the q_points in this finite element
  for (auto qp : index_range(d2u_vals))
    {
      OutputType & d2u = d2u_vals[qp];
 
      // Compute the gradient at this q_point
      d2u = 0;
 
      for (unsigned int l=0; l != n_dofs; l++)
        d2u.add_scaled(d2phi[l][qp], coef(l));
    }
 
  return;
}

References libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), and libMesh::index_range().

◆ side_rate()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::side_rate< Gradient >	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The time derivative (rate) of the solution variable var at the quadrature point qp on the current element side.

Definition at line 1404 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_side_fe
                   <typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution_rate>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution_rate(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::some_value().

◆ side_value() [1/2]

Number libMesh::FEMContext::side_value	(	unsigned int	var,
		unsigned int	qp
	)		const

inherited

Returns: The value of the solution variable var at the quadrature point qp on the current element side.

Note: This API currently is present for backward compatibility.

Definition at line 685 of file fem_context.C.

{
  Number u = 0.;
 
  this->side_value( var, qp, u );
 
  return u;
}

References libMesh::FEMContext::side_value().

Referenced by LaplaceSystem::side_constraint(), CoupledSystemQoI::side_qoi(), CoupledSystemQoI::side_qoi_derivative(), CurlCurlSystem::side_time_derivative(), and libMesh::FEMContext::side_value().

◆ side_value() [2/2]

template<typename OutputType >

void libMesh::FEMContext::side_value	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

inherited

Returns: The value of the solution variable var at the quadrature point qp on the current element side.

Note: This is the preferred API.

Definition at line 697 of file fem_context.C.

{
  this->some_value<OutputType,
                   &FEMContext::get_side_fe<typename TensorTools::MakeReal<OutputType>::type>,
                   &DiffContext::get_elem_solution>(var, qp, u);
}

References libMesh::DiffContext::get_elem_solution(), and libMesh::FEMContext::some_value().

◆ side_values()

template<typename OutputType >

template LIBMESH_EXPORT void libMesh::FEMContext::side_values< Gradient >	(	unsigned int	var,
		const NumericVector< Number > &	_system_vector,
		std::vector< OutputType > &	side_values_vector
	)		const

inherited

Fills a vector of values of the _system_vector at the all the quadrature points on the current element side.

Definition at line 708 of file fem_context.C.

{
  typedef typename TensorTools::MakeReal<OutputType>::type OutputShape;
 
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = get_localized_subvector(_system_vector, var);
 
  // Get the finite element object
  FEGenericBase<OutputShape> * the_side_fe = nullptr;
  this->get_side_fe<OutputShape>( var, the_side_fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector<OutputShape>> & phi = the_side_fe->get_phi();
 
  // Loop over all the q_points on this element
  for (auto qp : index_range(u_vals))
    {
      OutputType & u = u_vals[qp];
 
      // Compute the value at this q_point
      u = 0.;
 
      for (unsigned int l=0; l != n_dofs; l++)
        u += phi[l][qp] * coef(l);
    }
 
  return;
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::DiffContext::get_localized_subvector(), libMesh::FEGenericBase< OutputType >::get_phi(), and libMesh::index_range().

◆ some_gradient()

template<typename OutputType , typename FEMContext::FENeeded< OutputType >::grad_getter fe_getter, FEMContext::diff_subsolution_getter subsolution_getter>

void libMesh::FEMContext::some_gradient	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

protectedinherited

Helper function to reduce some code duplication in the *interior_gradient methods.

Definition at line 373 of file fem_context.C.

{
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = (this->*subsolution_getter)(var);
 
  // Get finite element object
  typename FENeeded<OutputType>::grad_base * fe = nullptr;
  (this->*fe_getter)( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector
                    <typename FENeeded<OutputType>::grad_base::OutputGradient>>
    & dphi = fe->get_dphi();
 
  // Accumulate solution derivatives
  du = 0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    du.add_scaled(dphi[l][qp], coef(l));
 
  return;
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEGenericBase< OutputType >::get_dphi(), and libMesh::FEMContext::get_elem_dim().

Referenced by libMesh::FEMContext::fixed_interior_gradient(), libMesh::FEMContext::fixed_side_gradient(), libMesh::FEMContext::interior_gradient(), and libMesh::FEMContext::interior_rate_gradient().

◆ some_hessian()

template<typename OutputType , typename FEMContext::FENeeded< OutputType >::hess_getter fe_getter, FEMContext::diff_subsolution_getter subsolution_getter>

void libMesh::FEMContext::some_hessian	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

protectedinherited

Helper function to reduce some code duplication in the *interior_hessian methods.

Definition at line 406 of file fem_context.C.

{
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = (this->*subsolution_getter)(var);
 
  // Get finite element object
  typename FENeeded<OutputType>::hess_base * fe = nullptr;
  (this->*fe_getter)( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector
                    <typename FENeeded<OutputType>::hess_base::OutputTensor>>
    & d2phi = fe->get_d2phi();
 
  // Accumulate solution second derivatives
  d2u = 0.0;
 
  for (unsigned int l=0; l != n_dofs; l++)
    d2u.add_scaled(d2phi[l][qp], coef(l));
 
  return;
}

References libMesh::FEGenericBase< OutputType >::get_d2phi(), libMesh::DiffContext::get_dof_indices(), and libMesh::FEMContext::get_elem_dim().

Referenced by libMesh::FEMContext::fixed_interior_hessian(), libMesh::FEMContext::fixed_side_hessian(), libMesh::FEMContext::interior_hessian(), and libMesh::FEMContext::side_hessian().

◆ some_value()

template<typename OutputType , typename FEMContext::FENeeded< OutputType >::value_getter fe_getter, FEMContext::diff_subsolution_getter subsolution_getter>

void libMesh::FEMContext::some_value	(	unsigned int	var,
		unsigned int	qp,
		OutputType &	u
	)		const

protectedinherited

Helper function to reduce some code duplication in the *interior_value methods.

Definition at line 339 of file fem_context.C.

{
  // Get local-to-global dof index lookup
  const unsigned int n_dofs = cast_int<unsigned int>
    (this->get_dof_indices(var).size());
 
  // Get current local coefficients
  const DenseSubVector<Number> & coef = (this->*subsolution_getter)(var);
  libmesh_assert_equal_to(coef.size(), n_dofs);
 
  // Get finite element object
  typename FENeeded<OutputType>::value_base * fe = nullptr;
  (this->*fe_getter)( var, fe, this->get_elem_dim() );
 
  // Get shape function values at quadrature point
  const std::vector<std::vector
                    <typename FENeeded<OutputType>::value_shape>> & phi = fe->get_phi();
  libmesh_assert_equal_to(phi.size(), n_dofs);
 
  // Accumulate solution value
  u = 0.;
 
  for (unsigned int l=0; l != n_dofs; l++)
    {
      libmesh_assert_less(qp, phi[l].size());
      u += phi[l][qp] * coef(l);
    }
}

References libMesh::DiffContext::get_dof_indices(), libMesh::FEMContext::get_elem_dim(), libMesh::FEGenericBase< OutputType >::get_phi(), and libMesh::DenseSubVector< T >::size().

Referenced by libMesh::FEMContext::fixed_interior_value(), libMesh::FEMContext::fixed_side_value(), libMesh::FEMContext::interior_accel(), libMesh::FEMContext::interior_rate(), libMesh::FEMContext::interior_value(), libMesh::FEMContext::side_accel(), libMesh::FEMContext::side_rate(), and libMesh::FEMContext::side_value().

◆ use_default_quadrature_rules()

void libMesh::FEMContext::use_default_quadrature_rules ( int extra_quadrature_order = 0 )

inherited

Use quadrature rules designed to over-integrate a mass matrix, plus extra_quadrature_order.

Definition at line 154 of file fem_context.C.

{
  _extra_quadrature_order = extra_quadrature_order;
 
  FEType hardest_fe_type = this->find_hardest_fe_type();
 
  for (const auto & dim : _elem_dims)
    {
      // Create an adequate quadrature rule
      _element_qrule[dim] =
        hardest_fe_type.default_quadrature_rule(dim, _extra_quadrature_order);
      if (dim)
        _side_qrule[dim] =
          hardest_fe_type.default_quadrature_rule(dim-1, _extra_quadrature_order);
      if (dim == 3)
        _edge_qrule = hardest_fe_type.default_quadrature_rule(1, _extra_quadrature_order);
    }
 
  this->attach_quadrature_rules();
}

References libMesh::FEMContext::_edge_qrule, libMesh::FEMContext::_elem_dims, libMesh::FEMContext::_element_qrule, libMesh::FEMContext::_extra_quadrature_order, libMesh::FEMContext::_side_qrule, libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEType::default_quadrature_rule(), dim, and libMesh::FEMContext::find_hardest_fe_type().

Referenced by libMesh::FEMContext::init_internal_data().

◆ use_quadrature_rules()

void libMesh::FEMContext::use_quadrature_rules	(	QuadratureType	qt,
		int	extra_quadrature_order = `0`
	)

inherited

Use a specified quadrature rule type qt for every element, side, and edge dimension, at the default quadrature order for the hardest FE type plus extra_quadrature_order.

This is useful when the default Gauss rules (whose points lie in the element interior) are inadequate. For example, rules whose points include the element vertices (QTRAP, QSIMPSON, QNODAL, QGAUSS_LOBATTO) can detect integrand behavior that is localized at element corners — a feature interior Gauss points can miss entirely. Note that not every quadrature type is defined for every element type (e.g. QGAUSS_LOBATTO is only available for tensor-product elements); it is the caller's responsibility to pick a type compatible with the mesh.

Definition at line 197 of file fem_context.C.

{
  _extra_quadrature_order = extra_quadrature_order;
 
  const FEType hardest_fe_type = this->find_hardest_fe_type();
 
  // Match the accuracy of the default rule (which over-integrates a mass
  // matrix), then add any requested extra order on top.
  const Order order = static_cast<Order>(hardest_fe_type.default_quadrature_order()
                                         + extra_quadrature_order);
 
  for (const auto & dim : _elem_dims)
    {
      // Create the requested quadrature rule for the element and its
      // lower-dimensional sides/edges
      _element_qrule[dim] = QBase::build(qt, dim, order);
      if (dim)
        _side_qrule[dim] = QBase::build(qt, cast_int<unsigned int>(dim - 1), order);
      if (dim == 3)
        _edge_qrule = QBase::build(qt, 1, order);
    }
 
  this->attach_quadrature_rules();
}

References libMesh::FEMContext::_edge_qrule, libMesh::FEMContext::_elem_dims, libMesh::FEMContext::_element_qrule, libMesh::FEMContext::_extra_quadrature_order, libMesh::FEMContext::_side_qrule, libMesh::FEMContext::attach_quadrature_rules(), libMesh::QBase::build(), libMesh::FEType::default_quadrature_order(), dim, and libMesh::FEMContext::find_hardest_fe_type().

Referenced by libMesh::VariationalSmootherSystem::init_context().

◆ use_unweighted_quadrature_rules()

void libMesh::FEMContext::use_unweighted_quadrature_rules ( int extra_quadrature_order = 0 )

inherited

Use quadrature rules designed to exactly integrate unweighted undistorted basis functions, plus extra_quadrature_order.

Definition at line 176 of file fem_context.C.

{
  _extra_quadrature_order = extra_quadrature_order;
 
  FEType hardest_fe_type = this->find_hardest_fe_type();
 
  for (const auto & dim : _elem_dims)
    {
      // Create an adequate quadrature rule
      _element_qrule[dim] =
        hardest_fe_type.unweighted_quadrature_rule(dim, _extra_quadrature_order);
      _side_qrule[dim] =
        hardest_fe_type.unweighted_quadrature_rule(dim-1, _extra_quadrature_order);
      if (dim == 3)
        _edge_qrule = hardest_fe_type.unweighted_quadrature_rule(1, _extra_quadrature_order);
    }
 
  this->attach_quadrature_rules();
}

References libMesh::FEMContext::_edge_qrule, libMesh::FEMContext::_elem_dims, libMesh::FEMContext::_element_qrule, libMesh::FEMContext::_extra_quadrature_order, libMesh::FEMContext::_side_qrule, libMesh::FEMContext::attach_quadrature_rules(), dim, libMesh::FEMContext::find_hardest_fe_type(), and libMesh::FEType::unweighted_quadrature_rule().

Member Data Documentation

◆ _active_vars

std::unique_ptr<const std::vector<unsigned int> > libMesh::FEMContext::_active_vars

protectedinherited

Variables on which to enable calculations, or nullptr if all variables in the System are to be enabled.

Definition at line 1072 of file fem_context.h.

Referenced by libMesh::FEMContext::active_vars(), libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::FEMContext(), libMesh::FEMContext::find_hardest_fe_type(), libMesh::FEMContext::init_internal_data(), and libMesh::FEMContext::pre_fe_reinit().

◆ _atype

AlgebraicType libMesh::FEMContext::_atype

protectedinherited

Keep track of what type of algebra reinitialization is to be done.

Definition at line 1087 of file fem_context.h.

Referenced by libMesh::FEMContext::algebraic_type(), and libMesh::FEMContext::set_algebraic_type().

◆ _boundary_info

const BoundaryInfo& libMesh::FEMContext::_boundary_info

protectedinherited

Saved reference to BoundaryInfo on the mesh for this System.

Used to answer boundary id requests.

Definition at line 1205 of file fem_context.h.

Referenced by libMesh::FEMContext::has_side_boundary_id(), and libMesh::FEMContext::side_boundary_ids().

◆ _custom_solution

const NumericVector<Number>* libMesh::FEMContext::_custom_solution

protectedinherited

Data with which to do algebra reinitialization.

Definition at line 1092 of file fem_context.h.

Referenced by libMesh::FEMContext::pre_fe_reinit(), and libMesh::FEMContext::set_custom_solution().

◆ _deltat

Real* libMesh::DiffContext::_deltat

privateinherited

Defaults to nullptr, can optionally be used to point to a timestep value in the System-derived class responsible for creating this DiffContext.

In DiffSystem's build_context() function, is assigned to point to the deltat member of that class.

Accessible via public get_deltat()/set_deltat() methods of this class.

Always test for nullptr before using!

Definition at line 655 of file diff_context.h.

Referenced by libMesh::DiffContext::get_deltat_value(), and libMesh::DiffContext::set_deltat_pointer().

◆ _dg_terms_active

bool libMesh::DGFEMContext::_dg_terms_active

private

Boolean flag to indicate whether or not the DG terms have been assembled and should be used in the global matrix assembly.

Definition at line 296 of file dg_fem_context.h.

Referenced by dg_terms_are_active(), neighbor_side_fe_reinit(), and side_fe_reinit().

◆ _dim

unsigned char libMesh::FEMContext::_dim

protectedinherited

Cached dimension of largest dimension element in this mesh.

Definition at line 1215 of file fem_context.h.

Referenced by DGFEMContext(), libMesh::FEMContext::get_dim(), and libMesh::FEMContext::get_elem_dim().

◆ _dof_indices

std::vector<dof_id_type> libMesh::DiffContext::_dof_indices

protectedinherited

Global Degree of freedom index lists.

Definition at line 639 of file diff_context.h.

Referenced by libMesh::DiffContext::get_dof_indices(), libMesh::DiffContext::get_dof_indices(), and libMesh::DiffContext::n_dof_indices().

◆ _dof_indices_var

std::vector<std::vector<dof_id_type> > libMesh::DiffContext::_dof_indices_var

protectedinherited

Definition at line 640 of file diff_context.h.

Referenced by libMesh::DiffContext::get_dof_indices(), libMesh::DiffContext::get_dof_indices(), libMesh::DiffContext::n_dof_indices(), and libMesh::DiffContext::n_vars().

◆ _edge_fe

std::map<FEType, std::unique_ptr<FEAbstract> > libMesh::FEMContext::_edge_fe

protectedinherited

Definition at line 1188 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::edge_fe_reinit(), libMesh::FEMContext::init_internal_data(), and libMesh::FEMContext::set_jacobian_tolerance().

◆ _edge_fe_var

std::vector<FEAbstract *> libMesh::FEMContext::_edge_fe_var

protectedinherited

Definition at line 1199 of file fem_context.h.

Referenced by libMesh::FEMContext::get_edge_fe(), libMesh::FEMContext::get_edge_fe(), libMesh::FEMContext::get_edge_fe(), and libMesh::FEMContext::init_internal_data().

◆ _edge_qrule

std::unique_ptr<QBase> libMesh::FEMContext::_edge_qrule

protectedinherited

Quadrature rules for element edges.

If the FEM context is told to prepare for edge integration on 3D elements, it will try to find a quadrature rule that correctly integrates all variables. Because edge rules only apply to 3D elements, we don't need to worry about multiple dimensions

Definition at line 1251 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::get_edge_qrule(), libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ _elem

const Elem* libMesh::FEMContext::_elem

protectedinherited

Current element for element_* to examine.

Definition at line 1210 of file fem_context.h.

Referenced by libMesh::FEMContext::get_elem(), libMesh::FEMContext::get_elem(), libMesh::FEMContext::get_elem_dim(), libMesh::FEMContext::has_elem(), and libMesh::FEMContext::set_elem().

◆ _elem_dim

unsigned char libMesh::FEMContext::_elem_dim

protectedinherited

Cached dimension of this->_elem.

Definition at line 1220 of file fem_context.h.

Referenced by libMesh::FEMContext::get_elem_dim(), and libMesh::FEMContext::set_elem().

◆ _elem_dims

std::set<unsigned char> libMesh::FEMContext::_elem_dims

protectedinherited

Cached dimensions of elements in the mesh, plus dimension 0 if SCALAR variables are in use.

Definition at line 1226 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::elem_dimensions(), libMesh::FEMContext::init_internal_data(), libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ _elem_elem_jacobian

DenseMatrix<Number> libMesh::DGFEMContext::_elem_elem_jacobian

private

The DG Jacobian terms.

Trial and test functions come from either element or neighbor.

Definition at line 257 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_elem_elem_jacobian(), get_elem_elem_jacobian(), and neighbor_side_fe_reinit().

◆ _elem_elem_subjacobians

std::vector<std::vector<std::unique_ptr<DenseSubMatrix<Number> > > > libMesh::DGFEMContext::_elem_elem_subjacobians

private

Definition at line 266 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_elem_elem_jacobian(), get_elem_elem_jacobian(), and neighbor_side_fe_reinit().

◆ _elem_fixed_solution

DenseVector<Number> libMesh::DiffContext::_elem_fixed_solution

protectedinherited

Element by element components of nonlinear_solution at a fixed point in a timestep, for optional use by e.g.

stabilized methods

Definition at line 604 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_fixed_solution(), and libMesh::DiffContext::get_elem_fixed_solution().

◆ _elem_fixed_subsolutions

std::vector<DenseSubVector<Number> > libMesh::DiffContext::_elem_fixed_subsolutions

protectedinherited

Definition at line 605 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::FEMContext::fixed_point_gradient(), libMesh::FEMContext::fixed_point_hessian(), libMesh::FEMContext::fixed_point_value(), libMesh::DiffContext::get_elem_fixed_solution(), and libMesh::DiffContext::get_elem_fixed_solution().

◆ _elem_jacobian

DenseMatrix<Number> libMesh::DiffContext::_elem_jacobian

protectedinherited

Element jacobian: derivatives of elem_residual with respect to elem_solution.

Only initialized if _have_local_matrices

Definition at line 616 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_jacobian(), and libMesh::DiffContext::get_elem_jacobian().

◆ _elem_neighbor_jacobian

DenseMatrix<Number> libMesh::DGFEMContext::_elem_neighbor_jacobian

private

Definition at line 258 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_elem_neighbor_jacobian(), get_elem_neighbor_jacobian(), and neighbor_side_fe_reinit().

◆ _elem_neighbor_subjacobians

std::vector<std::vector<std::unique_ptr<DenseSubMatrix<Number> > > > libMesh::DGFEMContext::_elem_neighbor_subjacobians

private

Definition at line 267 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_elem_neighbor_jacobian(), get_elem_neighbor_jacobian(), and neighbor_side_fe_reinit().

◆ _elem_qoi

std::vector<Number> libMesh::DiffContext::_elem_qoi

protectedinherited

Element quantity of interest contributions.

Definition at line 621 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_qois(), and libMesh::DiffContext::get_qois().

◆ _elem_qoi_derivative

std::vector<DenseVector<Number> > libMesh::DiffContext::_elem_qoi_derivative

protectedinherited

Element quantity of interest derivative contributions.

Definition at line 626 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_qoi_derivatives(), and libMesh::DiffContext::get_qoi_derivatives().

◆ _elem_qoi_subderivatives

std::vector<std::vector<DenseSubVector<Number> > > libMesh::DiffContext::_elem_qoi_subderivatives

protectedinherited

Definition at line 627 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_qoi_derivatives(), libMesh::DiffContext::get_qoi_derivatives(), and libMesh::FEMContext::pre_fe_reinit().

◆ _elem_residual

DenseVector<Number> libMesh::DiffContext::_elem_residual

protectedinherited

Element residual vector.

Definition at line 610 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_residual(), and libMesh::DiffContext::get_elem_residual().

◆ _elem_solution

DenseVector<Number> libMesh::DiffContext::_elem_solution

protectedinherited

Element by element components of nonlinear_solution as adjusted by a time_solver.

Definition at line 582 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution(), and libMesh::DiffContext::get_elem_solution().

◆ _elem_solution_accel

DenseVector<Number> libMesh::DiffContext::_elem_solution_accel

protectedinherited

Element by element components of du/dt as adjusted by a time_solver.

Definition at line 596 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution_accel(), and libMesh::DiffContext::get_elem_solution_accel().

◆ _elem_solution_rate

DenseVector<Number> libMesh::DiffContext::_elem_solution_rate

protectedinherited

Element by element components of du/dt as adjusted by a time_solver.

Definition at line 589 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution_rate(), and libMesh::DiffContext::get_elem_solution_rate().

◆ _elem_subjacobians

std::vector<std::vector<DenseSubMatrix<Number> > > libMesh::DiffContext::_elem_subjacobians

protectedinherited

Definition at line 634 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_jacobian(), and libMesh::DiffContext::get_elem_jacobian().

◆ _elem_subresiduals

std::vector<DenseSubVector<Number> > libMesh::DiffContext::_elem_subresiduals

protectedinherited

Element residual subvectors and (if _have_local_matrices) Jacobian submatrices.

Definition at line 633 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_residual(), and libMesh::DiffContext::get_elem_residual().

◆ _elem_subsolution_accels

std::vector<DenseSubVector<Number> > libMesh::DiffContext::_elem_subsolution_accels

protectedinherited

Definition at line 597 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution_accel(), and libMesh::DiffContext::get_elem_solution_accel().

◆ _elem_subsolution_rates

std::vector<DenseSubVector<Number> > libMesh::DiffContext::_elem_subsolution_rates

protectedinherited

Definition at line 590 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution_rate(), and libMesh::DiffContext::get_elem_solution_rate().

◆ _elem_subsolutions

std::vector<DenseSubVector<Number> > libMesh::DiffContext::_elem_subsolutions

protectedinherited

Definition at line 583 of file diff_context.h.

Referenced by libMesh::DiffContext::DiffContext(), libMesh::DiffContext::get_elem_solution(), libMesh::DiffContext::get_elem_solution(), libMesh::FEMContext::interior_curl(), libMesh::FEMContext::interior_div(), libMesh::FEMContext::point_curl(), libMesh::FEMContext::point_gradient(), libMesh::FEMContext::point_hessian(), libMesh::FEMContext::point_value(), and libMesh::FEMContext::side_gradient().

◆ _element_fe

std::vector<std::map<FEType, std::unique_ptr<FEAbstract> > > libMesh::FEMContext::_element_fe

protectedinherited

Finite element objects for each variable's interior, sides and edges.

We store FE objects for each element dimension present in the mesh, except for edge_fe which only applies to 3D elements.

Definition at line 1186 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::elem_fe_reinit(), libMesh::FEMContext::init_internal_data(), and libMesh::FEMContext::set_jacobian_tolerance().

◆ _element_fe_var

std::vector<std::vector<FEAbstract *> > libMesh::FEMContext::_element_fe_var

protectedinherited

Pointers to the same finite element objects, but indexed by variable number.

We store FE objects for each element dimension present in the mesh, except for edge_fe_var which only applies for 3D elements.

Definition at line 1197 of file fem_context.h.

Referenced by libMesh::FEMContext::get_element_fe(), libMesh::FEMContext::get_element_fe(), libMesh::FEMContext::get_element_fe(), and libMesh::FEMContext::init_internal_data().

◆ _element_qrule

std::vector<std::unique_ptr<QBase> > libMesh::FEMContext::_element_qrule

protectedinherited

Quadrature rule for element interior.

The FEM context will try to find a quadrature rule that correctly integrates all variables. We prepare quadrature rules for each element dimension in the mesh.

Definition at line 1234 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::get_element_qrule(), libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ _extra_quadrature_order

int libMesh::FEMContext::_extra_quadrature_order

protectedinherited

The extra quadrature order for this context.

Definition at line 1256 of file fem_context.h.

Referenced by libMesh::FEMContext::init_internal_data(), libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ _have_local_matrices

const bool libMesh::DiffContext::_have_local_matrices

protectedinherited

Whether we have local matrices allocated/initialized.

Definition at line 576 of file diff_context.h.

Referenced by libMesh::DiffContext::get_elem_jacobian(), libMesh::DiffContext::get_elem_jacobian(), libMesh::DiffContext::get_elem_jacobian(), libMesh::DiffContext::get_elem_jacobian(), and libMesh::FEMContext::pre_fe_reinit().

◆ _is_adjoint

bool libMesh::DiffContext::_is_adjoint

privateinherited

Is this context to be used for a primal or adjoint solve?

Definition at line 665 of file diff_context.h.

Referenced by libMesh::DiffContext::is_adjoint(), and libMesh::DiffContext::is_adjoint().

◆ _localized_vectors

std::map<const NumericVector<Number> *, std::pair<DenseVector<Number>, std::vector<DenseSubVector<Number> > > > libMesh::DiffContext::_localized_vectors

protectedinherited

Contains pointers to vectors the user has asked to be localized, keyed with pairs of element localized versions of that vector and per variable views.

Definition at line 571 of file diff_context.h.

Referenced by libMesh::DiffContext::add_localized_vector(), libMesh::DiffContext::get_localized_subvector(), libMesh::DiffContext::get_localized_subvector(), libMesh::DiffContext::get_localized_vector(), libMesh::DiffContext::get_localized_vector(), and libMesh::FEMContext::pre_fe_reinit().

◆ _mesh_sys

System* libMesh::FEMContext::_mesh_sys

protectedinherited

System from which to acquire moving mesh information.

Definition at line 1077 of file fem_context.h.

Referenced by libMesh::FEMContext::_do_elem_position_set(), libMesh::FEMContext::elem_edge_reinit(), libMesh::FEMContext::elem_position_get(), libMesh::FEMContext::elem_position_set(), libMesh::FEMContext::elem_reinit(), libMesh::FEMContext::elem_side_reinit(), libMesh::FEMContext::get_mesh_system(), libMesh::FEMContext::get_mesh_system(), and libMesh::FEMContext::set_mesh_system().

◆ _mesh_x_var

unsigned int libMesh::FEMContext::_mesh_x_var

protectedinherited

Variables from which to acquire moving mesh information.

Definition at line 1082 of file fem_context.h.

Referenced by libMesh::FEMContext::get_mesh_x_var(), and libMesh::FEMContext::set_mesh_x_var().

◆ _mesh_y_var

unsigned int libMesh::FEMContext::_mesh_y_var

protectedinherited

Definition at line 1082 of file fem_context.h.

Referenced by libMesh::FEMContext::get_mesh_y_var(), and libMesh::FEMContext::set_mesh_y_var().

◆ _mesh_z_var

unsigned int libMesh::FEMContext::_mesh_z_var

protectedinherited

Definition at line 1082 of file fem_context.h.

Referenced by libMesh::FEMContext::get_mesh_z_var(), and libMesh::FEMContext::set_mesh_z_var().

◆ _neighbor

const Elem* libMesh::DGFEMContext::_neighbor

private

Current neighbor element for assembling DG terms.

Definition at line 246 of file dg_fem_context.h.

Referenced by get_neighbor(), and set_neighbor().

◆ _neighbor_dof_indices

std::vector<dof_id_type> libMesh::DGFEMContext::_neighbor_dof_indices

private

Global Degree of freedom index lists for the neighbor element.

Definition at line 274 of file dg_fem_context.h.

Referenced by get_neighbor_dof_indices(), and neighbor_side_fe_reinit().

◆ _neighbor_dof_indices_var

std::vector<std::vector<dof_id_type> > libMesh::DGFEMContext::_neighbor_dof_indices_var

private

Definition at line 275 of file dg_fem_context.h.

Referenced by get_neighbor_dof_indices(), and neighbor_side_fe_reinit().

◆ _neighbor_elem_jacobian

DenseMatrix<Number> libMesh::DGFEMContext::_neighbor_elem_jacobian

private

Definition at line 259 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_elem_jacobian(), get_neighbor_elem_jacobian(), and neighbor_side_fe_reinit().

◆ _neighbor_elem_subjacobians

std::vector<std::vector<std::unique_ptr<DenseSubMatrix<Number> > > > libMesh::DGFEMContext::_neighbor_elem_subjacobians

private

Definition at line 268 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_elem_jacobian(), get_neighbor_elem_jacobian(), and neighbor_side_fe_reinit().

◆ _neighbor_neighbor_jacobian

DenseMatrix<Number> libMesh::DGFEMContext::_neighbor_neighbor_jacobian

private

Definition at line 260 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_neighbor_jacobian(), get_neighbor_neighbor_jacobian(), and neighbor_side_fe_reinit().

◆ _neighbor_neighbor_subjacobians

std::vector<std::vector<std::unique_ptr<DenseSubMatrix<Number> > > > libMesh::DGFEMContext::_neighbor_neighbor_subjacobians

private

Definition at line 269 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_neighbor_jacobian(), get_neighbor_neighbor_jacobian(), and neighbor_side_fe_reinit().

◆ _neighbor_residual

DenseVector<Number> libMesh::DGFEMContext::_neighbor_residual

private

Residual vector of the neighbor component.

Definition at line 251 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_residual(), get_neighbor_residual(), and neighbor_side_fe_reinit().

◆ _neighbor_side_fe

std::map<FEType, std::unique_ptr<FEAbstract> > libMesh::DGFEMContext::_neighbor_side_fe

private

Finite element objects for each variable's sides on the neighbor element.

We do not need FE objects for neighbor element interior since we just need to handle DG interface terms here.

Definition at line 284 of file dg_fem_context.h.

Referenced by DGFEMContext(), and neighbor_side_fe_reinit().

◆ _neighbor_side_fe_var

std::vector<FEAbstract *> libMesh::DGFEMContext::_neighbor_side_fe_var

private

Pointers to the same finite element objects on the neighbor element, but indexed by variable number.

Definition at line 290 of file dg_fem_context.h.

Referenced by DGFEMContext(), and get_neighbor_side_fe().

◆ _neighbor_subresiduals

std::vector<std::unique_ptr<DenseSubVector<Number> > > libMesh::DGFEMContext::_neighbor_subresiduals

private

Element residual subvectors and Jacobian submatrices.

Definition at line 265 of file dg_fem_context.h.

Referenced by DGFEMContext(), get_neighbor_residual(), get_neighbor_residual(), and neighbor_side_fe_reinit().

◆ _real_fe

std::unique_ptr<FEGenericBase<Real> > libMesh::FEMContext::_real_fe

mutableprotectedinherited

Definition at line 1094 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _real_fe_derivative_level

int libMesh::FEMContext::_real_fe_derivative_level

mutableprotectedinherited

Definition at line 1096 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _real_fe_is_inf

bool libMesh::FEMContext::_real_fe_is_inf

mutableprotectedinherited

Definition at line 1100 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _real_grad_fe

std::unique_ptr<FEGenericBase<RealGradient> > libMesh::FEMContext::_real_grad_fe

mutableprotectedinherited

Definition at line 1095 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _real_grad_fe_derivative_level

int libMesh::FEMContext::_real_grad_fe_derivative_level

mutableprotectedinherited

Definition at line 1097 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _real_grad_fe_is_inf

bool libMesh::FEMContext::_real_grad_fe_is_inf

mutableprotectedinherited

Definition at line 1101 of file fem_context.h.

Referenced by libMesh::FEMContext::cached_fe().

◆ _side_fe

std::vector<std::map<FEType, std::unique_ptr<FEAbstract> > > libMesh::FEMContext::_side_fe

protectedinherited

Definition at line 1187 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::init_internal_data(), neighbor_side_fe_reinit(), libMesh::FEMContext::set_jacobian_tolerance(), and libMesh::FEMContext::side_fe_reinit().

◆ _side_fe_var

std::vector<std::vector<FEAbstract *> > libMesh::FEMContext::_side_fe_var

protectedinherited

Definition at line 1198 of file fem_context.h.

Referenced by libMesh::FEMContext::get_side_fe(), libMesh::FEMContext::get_side_fe(), libMesh::FEMContext::get_side_fe(), and libMesh::FEMContext::init_internal_data().

◆ _side_qrule

std::vector<std::unique_ptr<QBase> > libMesh::FEMContext::_side_qrule

protectedinherited

Quadrature rules for element sides The FEM context will try to find a quadrature rule that correctly integrates all variables.

We prepare quadrature rules for each element dimension in the mesh.

Definition at line 1242 of file fem_context.h.

Referenced by libMesh::FEMContext::attach_quadrature_rules(), libMesh::FEMContext::get_side_qrule(), libMesh::FEMContext::use_default_quadrature_rules(), libMesh::FEMContext::use_quadrature_rules(), and libMesh::FEMContext::use_unweighted_quadrature_rules().

◆ _system

const System& libMesh::DiffContext::_system

privateinherited

A reference to the system this context is constructed with.

Definition at line 660 of file diff_context.h.

Referenced by libMesh::DiffContext::get_system().

◆ edge

unsigned char libMesh::FEMContext::edge

inherited

Current edge for edge_* to examine.

Definition at line 1036 of file fem_context.h.

Referenced by libMesh::FEMContext::get_edge().

◆ elem_solution_accel_derivative

Real libMesh::DiffContext::elem_solution_accel_derivative

inherited

The derivative of elem_solution_accel with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Definition at line 510 of file diff_context.h.

Referenced by libMesh::NewmarkSolver::_general_residual(), libMesh::DiffContext::get_elem_solution_accel_derivative(), and libMesh::FirstOrderUnsteadySolver::prepare_accel().

◆ elem_solution_derivative

Real libMesh::DiffContext::elem_solution_derivative

inherited

The derivative of elem_solution with respect to the current nonlinear solution.

Definition at line 496 of file diff_context.h.

Referenced by libMesh::NewmarkSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::EulerSolver::_general_residual(), CoupledSystem::element_constraint(), libMesh::EigenTimeSolver::element_residual(), CoupledSystem::element_time_derivative(), NavierSystem::element_time_derivative(), SolidSystem::element_time_derivative(), libMesh::DiffContext::get_elem_solution_derivative(), and libMesh::EigenTimeSolver::side_residual().

◆ elem_solution_rate_derivative

Real libMesh::DiffContext::elem_solution_rate_derivative

inherited

The derivative of elem_solution_rate with respect to the current nonlinear solution, for use by systems with non default mass_residual terms.

Definition at line 503 of file diff_context.h.

Referenced by libMesh::NewmarkSolver::_general_residual(), libMesh::Euler2Solver::_general_residual(), libMesh::EulerSolver::_general_residual(), libMesh::EigenTimeSolver::element_residual(), libMesh::DiffContext::get_elem_solution_rate_derivative(), libMesh::FEMPhysics::mass_residual(), libMesh::DifferentiablePhysics::nonlocal_mass_residual(), and libMesh::EigenTimeSolver::side_residual().

◆ fixed_solution_derivative

Real libMesh::DiffContext::fixed_solution_derivative

inherited

The derivative of elem_fixed_solution with respect to the nonlinear solution, for use by systems constructing jacobians with elem_fixed_solution based methods.

Definition at line 517 of file diff_context.h.

Referenced by libMesh::Euler2Solver::_general_residual(), libMesh::EulerSolver::_general_residual(), libMesh::SteadySolver::_general_residual(), and libMesh::DiffContext::get_fixed_solution_derivative().

◆ side

unsigned char libMesh::FEMContext::side

inherited

Current side for side_* to examine.

Definition at line 1031 of file fem_context.h.

Referenced by libMesh::RBConstruction::add_scaled_matrix_and_vector(), assembly_with_dg_fem_context(), AssemblyF0::boundary_assembly(), AssemblyF1::boundary_assembly(), AssemblyF2::boundary_assembly(), AssemblyA0::boundary_assembly(), AssemblyA1::boundary_assembly(), AssemblyA2::boundary_assembly(), libMesh::FEMContext::get_side(), libMesh::FEMContext::has_side_boundary_id(), libMesh::RBEIMConstruction::initialize_qp_data(), and libMesh::FEMContext::side_boundary_ids().

◆ system_time

const Real libMesh::DiffContext::system_time

inherited

This is the time stored in the System class at the time this context was created, i.e.

the time at the beginning of the current timestep. This value gets set in the constructor and unlike DiffContext::time, is not tweaked mid-timestep by transient solvers: it remains equal to the value it was assigned at construction.

Definition at line 490 of file diff_context.h.

Referenced by libMesh::DiffContext::get_system_time().

◆ time

Real libMesh::DiffContext::time

inherited

For time-dependent problems, this is the time t for which the current nonlinear_solution is defined.

FIXME - this needs to be tweaked mid-timestep by all transient solvers!

Definition at line 481 of file diff_context.h.

Referenced by libMesh::DiffContext::get_time(), and libMesh::DiffContext::set_time().

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

include/systems/dg_fem_context.h
src/systems/dg_fem_context.C

Public Types

Public Member Functions

Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ diff_subsolution_getter

◆ localized_vectors_iterator

Member Enumeration Documentation

◆ AlgebraicType

Constructor & Destructor Documentation

◆ DGFEMContext()

◆ ~DGFEMContext()

Member Function Documentation

◆ _do_elem_position_set()

◆ _update_time_from_system()

◆ active_vars()

◆ add_localized_vector()

◆ algebraic_type()

◆ attach_quadrature_rules()

◆ build_new_fe()

◆ cached_fe() [1/3]

◆ cached_fe() [2/3]

◆ cached_fe() [3/3]

◆ dg_terms_are_active()

◆ edge_fe_reinit()

◆ elem_dimensions()

◆ elem_edge_reinit()

◆ elem_fe_reinit()

◆ elem_position_get()

◆ elem_position_set()

◆ elem_reinit()

◆ elem_side_reinit()

◆ find_hardest_fe_type()

◆ fixed_interior_gradient() [1/2]

◆ fixed_interior_gradient() [2/2]

◆ fixed_interior_hessian() [1/2]

◆ fixed_interior_hessian() [2/2]

◆ fixed_interior_value() [1/2]

◆ fixed_interior_value() [2/2]

◆ fixed_point_gradient() [1/2]

◆ fixed_point_gradient() [2/2]

◆ fixed_point_hessian() [1/2]

◆ fixed_point_hessian() [2/2]

◆ fixed_point_value() [1/2]

◆ fixed_point_value() [2/2]

◆ fixed_side_gradient() [1/2]

◆ fixed_side_gradient() [2/2]

◆ fixed_side_hessian() [1/2]

◆ fixed_side_hessian() [2/2]

◆ fixed_side_value() [1/2]

◆ fixed_side_value() [2/2]

◆ get_deltat_value()

◆ get_dim()

◆ get_dof_indices() [1/4]

◆ get_dof_indices() [2/4]

◆ get_dof_indices() [3/4]

◆ get_dof_indices() [4/4]

◆ get_edge()

◆ get_edge_fe() [1/3]

◆ get_edge_fe() [2/3]

◆ get_edge_fe() [3/3]

◆ get_edge_qrule()

◆ get_elem() [1/2]

◆ get_elem() [2/2]

◆ get_elem_dim()

◆ get_elem_elem_jacobian() [1/4]

◆ get_elem_elem_jacobian() [2/4]

◆ get_elem_elem_jacobian() [3/4]

◆ get_elem_elem_jacobian() [4/4]

◆ get_elem_fixed_solution() [1/4]

◆ get_elem_fixed_solution() [2/4]

◆ get_elem_fixed_solution() [3/4]

◆ get_elem_fixed_solution() [4/4]

◆ get_elem_jacobian() [1/4]

◆ get_elem_jacobian() [2/4]

◆ get_elem_jacobian() [3/4]