libMesh::FE< Dim, T > Class Template Reference

#include <fe.h>

Inheritance diagram for libMesh::FE< Dim, T >:

[legend]

List of all members.

Public Types
typedef FEGenericBase < typename FEOutputType< T > ::type >::OutputShape	OutputShape
typedef TensorTools::IncrementRank < OutputShape >::type	OutputGradient
typedef TensorTools::IncrementRank < OutputGradient >::type	OutputTensor
typedef TensorTools::DecrementRank < OutputShape >::type	OutputDivergence
typedef TensorTools::MakeNumber < OutputShape >::type	OutputNumber
typedef TensorTools::IncrementRank < OutputNumber >::type	OutputNumberGradient
typedef TensorTools::IncrementRank < OutputNumberGradient >::type	OutputNumberTensor
typedef TensorTools::DecrementRank < OutputNumber >::type	OutputNumberDivergence
Public Member Functions
	FE (const FEType &fet)
virtual unsigned int	n_shape_functions () const
virtual FEContinuity	get_continuity () const
virtual bool	is_hierarchic () const
virtual void	reinit (const Elem *elem, const std::vector< Point > *const pts=NULL, const std::vector< Real > *const weights=NULL)
virtual void	reinit (const Elem *elem, const unsigned int side, const Real tolerance=TOLERANCE, const std::vector< Point > *const pts=NULL, const std::vector< Real > *const weights=NULL)
virtual void	edge_reinit (const Elem *elem, const unsigned int edge, const Real tolerance=TOLERANCE, const std::vector< Point > *const pts=NULL, const std::vector< Real > *const weights=NULL)
virtual void	side_map (const Elem *elem, const Elem *side, const unsigned int s, const std::vector< Point > &reference_side_points, std::vector< Point > &reference_points)
virtual void	attach_quadrature_rule (QBase *q)
virtual unsigned int	n_quadrature_points () const
virtual bool	shapes_need_reinit () const
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order, const unsigned int, const unsigned int, const Point &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int i, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int i, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order order, const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
RealGradient	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const ElemType type, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const ElemType type, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
void	nodal_soln (const Elem *, const Order, const std::vector< Number > &, std::vector< Number > &)
template<>
void	nodal_soln (const Elem *, const Order, const std::vector< Number > &, std::vector< Number > &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType, const Order)
template<>
unsigned int	n_dofs (const ElemType, const Order)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
template<>
RealGradient	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
RealGradient	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &)
template<>
RealGradient	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &)
template<>
RealGradient	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
RealGradient	shape (const Elem *elem, const Order order, const unsigned int i, const Point &)
template<>
RealGradient	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &)
template<>
RealGradient	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
RealGradient	shape_second_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType, const Order o)
template<>
unsigned int	n_dofs (const ElemType, const Order o)
template<>
unsigned int	n_dofs (const ElemType, const Order o)
template<>
unsigned int	n_dofs (const ElemType, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
unsigned int	n_dofs_per_elem (const ElemType, const Order)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const Point &p)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &p)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order order, const unsigned int i, const Point &p)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order order, const unsigned int i, const unsigned int j, const Point &p)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int, const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
void	nodal_soln (const Elem *elem, const Order order, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_at_node (const ElemType, const Order, const unsigned int)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
FEContinuity	get_continuity () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
bool	is_hierarchic () const
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
void	compute_constraints (DofConstraints &, DofMap &, const unsigned int, const Elem *)
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
bool	shapes_need_reinit () const
template<>
Real	shape (const ElemType, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape (const Elem *, const Order, const unsigned int libmesh_dbg_var(i), const Point &)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const Point &point_in)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &point_in)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int libmesh_dbg_var(j), const Point &point_in)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const Point &point_in)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &point_in)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &point_in)
template<>
Real	shape (const ElemType, const Order, const unsigned int, const Point &)
template<>
Real	shape (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const Point &point_in)
template<>
Real	shape_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &point_in)
template<>
Real	shape_second_deriv (const ElemType, const Order, const unsigned int, const unsigned int, const Point &)
template<>
Real	shape_second_deriv (const Elem *elem, const Order libmesh_dbg_var(order), const unsigned int i, const unsigned int j, const Point &point_in)
const std::vector< std::vector < OutputShape > > &	get_phi () const
const std::vector< std::vector < OutputGradient > > &	get_dphi () const
const std::vector< std::vector < OutputShape > > &	get_curl_phi () const
const std::vector< std::vector < OutputDivergence > > &	get_div_phi () const
const std::vector< std::vector < OutputShape > > &	get_dphidx () const
const std::vector< std::vector < OutputShape > > &	get_dphidy () const
const std::vector< std::vector < OutputShape > > &	get_dphidz () const
const std::vector< std::vector < OutputShape > > &	get_dphidxi () const
const std::vector< std::vector < OutputShape > > &	get_dphideta () const
const std::vector< std::vector < OutputShape > > &	get_dphidzeta () const
const std::vector< std::vector < OutputTensor > > &	get_d2phi () const
const std::vector< std::vector < OutputShape > > &	get_d2phidx2 () const
const std::vector< std::vector < OutputShape > > &	get_d2phidxdy () const
const std::vector< std::vector < OutputShape > > &	get_d2phidxdz () const
const std::vector< std::vector < OutputShape > > &	get_d2phidy2 () const
const std::vector< std::vector < OutputShape > > &	get_d2phidydz () const
const std::vector< std::vector < OutputShape > > &	get_d2phidz2 () const
const std::vector< std::vector < OutputShape > > &	get_d2phidxi2 () const
const std::vector< std::vector < OutputShape > > &	get_d2phidxideta () const
const std::vector< std::vector < OutputShape > > &	get_d2phidxidzeta () const
const std::vector< std::vector < OutputShape > > &	get_d2phideta2 () const
const std::vector< std::vector < OutputShape > > &	get_d2phidetadzeta () const
const std::vector< std::vector < OutputShape > > &	get_d2phidzeta2 () const
const std::vector < OutputGradient > &	get_dphase () const
const std::vector< Real > &	get_Sobolev_weight () const
const std::vector< RealGradient > &	get_Sobolev_dweight () const
void	print_phi (std::ostream &os) const
void	print_dphi (std::ostream &os) const
void	print_d2phi (std::ostream &os) const
const std::vector< Point > &	get_xyz () const
const std::vector< Real > &	get_JxW () const
const std::vector< RealGradient > &	get_dxyzdxi () const
const std::vector< RealGradient > &	get_dxyzdeta () const
const std::vector< RealGradient > &	get_dxyzdzeta () const
const std::vector< RealGradient > &	get_d2xyzdxi2 () const
const std::vector< RealGradient > &	get_d2xyzdeta2 () const
const std::vector< RealGradient > &	get_d2xyzdzeta2 () const
const std::vector< RealGradient > &	get_d2xyzdxideta () const
const std::vector< RealGradient > &	get_d2xyzdxidzeta () const
const std::vector< RealGradient > &	get_d2xyzdetadzeta () const
const std::vector< Real > &	get_dxidx () const
const std::vector< Real > &	get_dxidy () const
const std::vector< Real > &	get_dxidz () const
const std::vector< Real > &	get_detadx () const
const std::vector< Real > &	get_detady () const
const std::vector< Real > &	get_detadz () const
const std::vector< Real > &	get_dzetadx () const
const std::vector< Real > &	get_dzetady () const
const std::vector< Real > &	get_dzetadz () const
const std::vector< std::vector < Point > > &	get_tangents () const
const std::vector< Point > &	get_normals () const
const std::vector< Real > &	get_curvatures () const
ElemType	get_type () const
unsigned int	get_p_level () const
FEType	get_fe_type () const
Order	get_order () const
FEFamily	get_family () const
const FEMap &	get_fe_map () const
void	print_JxW (std::ostream &os) const
void	print_xyz (std::ostream &os) const
void	print_info (std::ostream &os) const
Static Public Member Functions
static OutputShape	shape (const ElemType t, const Order o, const unsigned int i, const Point &p)
static OutputShape	shape (const Elem *elem, const Order o, const unsigned int i, const Point &p)
static OutputShape	shape_deriv (const ElemType t, const Order o, const unsigned int i, const unsigned int j, const Point &p)
static OutputShape	shape_deriv (const Elem *elem, const Order o, const unsigned int i, const unsigned int j, const Point &p)
static OutputShape	shape_second_deriv (const ElemType t, const Order o, const unsigned int i, const unsigned int j, const Point &p)
static OutputShape	shape_second_deriv (const Elem *elem, const Order o, const unsigned int i, const unsigned int j, const Point &p)
static void	nodal_soln (const Elem *elem, const Order o, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln)
static unsigned int	n_shape_functions (const ElemType t, const Order o)
static unsigned int	n_dofs (const ElemType t, const Order o)
static unsigned int	n_dofs_at_node (const ElemType t, const Order o, const unsigned int n)
static unsigned int	n_dofs_per_elem (const ElemType t, const Order o)
static void	dofs_on_side (const Elem *const elem, const Order o, unsigned int s, std::vector< unsigned int > &di)
static void	dofs_on_edge (const Elem *const elem, const Order o, unsigned int e, std::vector< unsigned int > &di)
static Point	inverse_map (const Elem *elem, const Point &p, const Real tolerance=TOLERANCE, const bool secure=true)
static void	inverse_map (const Elem *elem, const std::vector< Point > &physical_points, std::vector< Point > &reference_points, const Real tolerance=TOLERANCE, const bool secure=true)
static void	compute_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
static Point	map (const Elem *elem, const Point &reference_point)
static Point	map_xi (const Elem *elem, const Point &reference_point)
static Point	map_eta (const Elem *elem, const Point &reference_point)
static Point	map_zeta (const Elem *elem, const Point &reference_point)
static AutoPtr< FEGenericBase >	build (const unsigned int dim, const FEType &type)
static AutoPtr< FEGenericBase >	build_InfFE (const unsigned int dim, const FEType &type)
static void	compute_proj_constraints (DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
static void	coarsened_dof_values (const NumericVector< Number > &global_vector, const DofMap &dof_map, const Elem *coarse_elem, DenseVector< Number > &coarse_dofs, const unsigned int var, const bool use_old_dof_indices=false)
static void	compute_periodic_constraints (DofConstraints &constraints, DofMap &dof_map, const PeriodicBoundaries &boundaries, const MeshBase &mesh, const PointLocatorBase *point_locator, const unsigned int variable_number, const Elem *elem)
static bool	on_reference_element (const Point &p, const ElemType t, const Real eps=TOLERANCE)
static void	get_refspace_nodes (const ElemType t, std::vector< Point > &nodes)
static void	compute_node_constraints (NodeConstraints &constraints, const Elem *elem)
static void	compute_periodic_node_constraints (NodeConstraints &constraints, const PeriodicBoundaries &boundaries, const MeshBase &mesh, const PointLocatorBase *point_locator, const Elem *elem)
static void	print_info (std::ostream &out=libMesh::out)
static std::string	get_info ()
static unsigned int	n_objects ()
static void	enable_print_counter_info ()
static void	disable_print_counter_info ()
Protected Types
typedef std::map< std::string, std::pair< unsigned int, unsigned int > >	Counts
Protected Member Functions
virtual void	init_shape_functions (const std::vector< Point > &qp, const Elem *e)
virtual void	init_base_shape_functions (const std::vector< Point > &qp, const Elem *e)
virtual void	compute_shape_functions (const Elem *elem, const std::vector< Point > &qp)
void	increment_constructor_count (const std::string &name)
void	increment_destructor_count (const std::string &name)
Protected Attributes
std::vector< Point >	cached_nodes
ElemType	last_side
unsigned int	last_edge
AutoPtr< FETransformationBase < FEOutputType< T >::type > >	_fe_trans
std::vector< std::vector < OutputShape > >	phi
std::vector< std::vector < OutputGradient > >	dphi
std::vector< std::vector < OutputShape > >	curl_phi
std::vector< std::vector < OutputDivergence > >	div_phi
std::vector< std::vector < OutputShape > >	dphidxi
std::vector< std::vector < OutputShape > >	dphideta
std::vector< std::vector < OutputShape > >	dphidzeta
std::vector< std::vector < OutputShape > >	dphidx
std::vector< std::vector < OutputShape > >	dphidy
std::vector< std::vector < OutputShape > >	dphidz
std::vector< std::vector < OutputTensor > >	d2phi
std::vector< std::vector < OutputShape > >	d2phidxi2
std::vector< std::vector < OutputShape > >	d2phidxideta
std::vector< std::vector < OutputShape > >	d2phidxidzeta
std::vector< std::vector < OutputShape > >	d2phideta2
std::vector< std::vector < OutputShape > >	d2phidetadzeta
std::vector< std::vector < OutputShape > >	d2phidzeta2
std::vector< std::vector < OutputShape > >	d2phidx2
std::vector< std::vector < OutputShape > >	d2phidxdy
std::vector< std::vector < OutputShape > >	d2phidxdz
std::vector< std::vector < OutputShape > >	d2phidy2
std::vector< std::vector < OutputShape > >	d2phidydz
std::vector< std::vector < OutputShape > >	d2phidz2
std::vector< OutputGradient >	dphase
std::vector< RealGradient >	dweight
std::vector< Real >	weight
AutoPtr< FEMap >	_fe_map
const unsigned int	dim
bool	calculations_started
bool	calculate_phi
bool	calculate_dphi
bool	calculate_d2phi
bool	calculate_curl_phi
bool	calculate_div_phi
bool	calculate_dphiref
const FEType	fe_type
ElemType	elem_type
unsigned int	_p_level
QBase *	qrule
bool	shapes_on_quadrature
Static Protected Attributes
static Counts	_counts
static Threads::atomic < unsigned int >	_n_objects
static Threads::spin_mutex	_mutex
static bool	_enable_print_counter = true
Friends
class	InfFE
std::ostream &	operator<< (std::ostream &os, const FEAbstract &fe)

---

Detailed Description

template<unsigned int Dim, FEFamily T>
class libMesh::FE< Dim, T >

A specific instatiation of the FEBase class. This class is templated, and specific template instantiations will result in different Finite Element families. Full specialization of the template for specific dimensions(Dim) and families (T) provide support for specific finite element types. The use of templates allows for compile-time optimization, however it requires that the specific finite element family and dimension is also known at compile time. If this is too restricting for your application you can use the FEBase::build() member to create abstract (but still optimized) finite elements.

Author:

Benjamin S. Kirk

Date:

2002-2007

Definition at line 91 of file fe.h.

---

Member Typedef Documentation

typedef std::map<std::string, std::pair<unsigned int, unsigned int> > libMesh::ReferenceCounter::Counts [protected, inherited]

Data structure to log the information. The log is identified by the class name.

Definition at line 113 of file reference_counter.h.

typedef TensorTools::DecrementRank<OutputShape>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputDivergence [inherited]

Definition at line 151 of file fe_base.h.

typedef TensorTools::IncrementRank<OutputShape>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputGradient [inherited]

Definition at line 149 of file fe_base.h.

typedef TensorTools::MakeNumber<OutputShape>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputNumber [inherited]

Definition at line 152 of file fe_base.h.

typedef TensorTools::DecrementRank<OutputNumber>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputNumberDivergence [inherited]

Definition at line 155 of file fe_base.h.

typedef TensorTools::IncrementRank<OutputNumber>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputNumberGradient [inherited]

Definition at line 153 of file fe_base.h.

typedef TensorTools::IncrementRank<OutputNumberGradient>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputNumberTensor [inherited]

Definition at line 154 of file fe_base.h.

template<unsigned int Dim, FEFamily T>

typedef FEGenericBase<typename FEOutputType<T>::type>::OutputShape libMesh::FE< Dim, T >::OutputShape

Convenient typedefs for gradients of output, hessians of output, and potentially-complex-valued versions of same.

Reimplemented from libMesh::FEGenericBase< FEOutputType< T >::type >.

Definition at line 103 of file fe.h.

typedef TensorTools::IncrementRank<OutputGradient>::type libMesh::FEGenericBase< FEOutputType< T >::type >::OutputTensor [inherited]

Definition at line 150 of file fe_base.h.

---

Constructor & Destructor Documentation

template<unsigned int Dim, FEFamily T>

libMesh::FE< Dim, T >::FE

(

const FEType &

fet

)

[inline, explicit]

Constructor.

Definition at line 904 of file fe.h.

References libMesh::FEAbstract::get_family().

                                :
  FEGenericBase<typename FEOutputType<T>::type> (Dim,fet),
  last_side(INVALID_ELEM),
  last_edge(libMesh::invalid_uint)
{
  // Sanity check.  Make sure the
  // Family specified in the template instantiation
  // matches the one in the FEType object
  libmesh_assert_equal_to (T, this->get_family());
}

---

Member Function Documentation

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::attach_quadrature_rule

(

QBase *

q

)

[inline, virtual]

Provides the class with the quadrature rule, which provides the locations (on a reference element) where the shape functions are to be calculated.

Implements libMesh::FEAbstract.

Definition at line 44 of file fe.C.

References libMesh::FEAbstract::elem_type, libMeshEnums::INVALID_ELEM, and libMesh::FEAbstract::qrule.

{
  libmesh_assert(q);
  this->qrule = q;
  // make sure we don't cache results from a previous quadrature rule
  this->elem_type = INVALID_ELEM;
  return;
}

static AutoPtr<FEGenericBase> libMesh::FEGenericBase< FEOutputType< T >::type >::build	(	const unsigned int	dim,
		const FEType &	type
	)			[static, inherited]

Builds a specific finite element type. A AutoPtr<FEGenericBase> is returned to prevent a memory leak. This way the user need not remember to delete the object.

The build call will fail if the OutputType of this class is not compatible with the output required for the requested type

Reimplemented from libMesh::FEAbstract.

static AutoPtr<FEGenericBase> libMesh::FEGenericBase< FEOutputType< T >::type >::build_InfFE	(	const unsigned int	dim,
		const FEType &	type
	)			[static, inherited]

Builds a specific infinite element type. A AutoPtr<FEGenericBase> is returned to prevent a memory leak. This way the user need not remember to delete the object.

The build call will fail if the OutputShape of this class is not compatible with the output required for the requested type

static void libMesh::FEGenericBase< FEOutputType< T >::type >::coarsened_dof_values	(	const NumericVector< Number > &	global_vector,
		const DofMap &	dof_map,
		const Elem *	coarse_elem,
		DenseVector< Number > &	coarse_dofs,
		const unsigned int	var,
		const bool	use_old_dof_indices = false
	)			[static, inherited]

Creates a local projection on coarse_elem, based on the DoF values in global_vector for it's children.

template<>

void libMesh::FE< 3, XYZ >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 941 of file fe_xyz.C.

{}

template<>

void libMesh::FE< 2, XYZ >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 940 of file fe_xyz.C.

{}

template<>

void libMesh::FE< 3, SCALAR >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 129 of file fe_scalar.C.

  { }

template<>

void libMesh::FE< 2, SCALAR >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 122 of file fe_scalar.C.

  { }

template<>

void libMesh::FE< 3, NEDELEC_ONE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 603 of file fe_nedelec_one.C.

  { nedelec_one_compute_constraints(constraints, dof_map, variable_number, elem, /*Dim=*/3); }

template<>

void libMesh::FE< 2, NEDELEC_ONE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 596 of file fe_nedelec_one.C.

  { nedelec_one_compute_constraints(constraints, dof_map, variable_number, elem, /*Dim=*/2); }

template<>

void libMesh::FE< 1, NEDELEC_ONE >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 589 of file fe_nedelec_one.C.

  { NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

void libMesh::FE< 0, NEDELEC_ONE >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 582 of file fe_nedelec_one.C.

  { NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

void libMesh::FE< 3, MONOMIAL >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 418 of file fe_monomial.C.

{}

template<>

void libMesh::FE< 2, MONOMIAL >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 417 of file fe_monomial.C.

{}

template<>

void libMesh::FE< 3, LAGRANGE_VEC >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 958 of file fe_lagrange_vec.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { //libmesh_not_implemented();
    FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem); 
  }

template<>

void libMesh::FE< 2, LAGRANGE_VEC >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 949 of file fe_lagrange_vec.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { //libmesh_not_implemented();
    FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem); 
  }

template<>

void libMesh::FE< 3, LAGRANGE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 871 of file fe_lagrange.C.

  { lagrange_compute_constraints(constraints, dof_map, variable_number, elem, /*Dim=*/3); }

template<>

void libMesh::FE< 2, LAGRANGE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 864 of file fe_lagrange.C.

  { lagrange_compute_constraints(constraints, dof_map, variable_number, elem, /*Dim=*/2); }

template<>

void libMesh::FE< 3, L2_LAGRANGE >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 510 of file fe_l2_lagrange.C.

  { }

template<>

void libMesh::FE< 2, L2_LAGRANGE >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 503 of file fe_l2_lagrange.C.

  { }

template<>

void libMesh::FE< 3, L2_HIERARCHIC >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 209 of file fe_l2_hierarchic.C.

  { }

template<>

void libMesh::FE< 2, L2_HIERARCHIC >::compute_constraints	(	DofConstraints &	,
		DofMap &	,
		const unsigned	int,
		const Elem *
	)			[inline]

Definition at line 202 of file fe_l2_hierarchic.C.

  { }

template<>

void libMesh::FE< 3, HIERARCHIC >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 380 of file fe_hierarchic.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<>

void libMesh::FE< 2, HIERARCHIC >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 373 of file fe_hierarchic.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<>

void libMesh::FE< 3, HERMITE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 368 of file fe_hermite.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<>

void libMesh::FE< 2, HERMITE >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 361 of file fe_hermite.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<>

void libMesh::FE< 3, CLOUGH >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 359 of file fe_clough.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<>

void libMesh::FE< 2, CLOUGH >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[inline]

Definition at line 352 of file fe_clough.C.

References libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints().

  { compute_proj_constraints(constraints, dof_map, variable_number, elem); }

template<unsigned int Dim, FEFamily T>

static void libMesh::FE< Dim, T >::compute_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[static]

Computes the constraint matrix contributions (for non-conforming adapted meshes) corresponding to variable number var_number, using element-specific optimizations if possible.

void libMesh::FEAbstract::compute_node_constraints	(	NodeConstraints &	constraints,
		const Elem *	elem
	)			[static, inherited]

Computes the nodal constraint contributions (for non-conforming adapted meshes), using Lagrange geometry

Definition at line 878 of file fe_abstract.C.

References std::abs(), libMesh::Elem::build_side(), libMesh::Elem::default_order(), libMesh::Elem::dim(), libMesh::FEAbstract::fe_type, libMesh::FEInterface::inverse_map(), libMeshEnums::LAGRANGE, libMesh::Elem::level(), libMesh::FEInterface::n_dofs(), libMesh::Elem::n_sides(), libMesh::Elem::neighbor(), libMesh::Elem::parent(), libMesh::Real, libMesh::remote_elem, libMesh::FEInterface::shape(), libMesh::Threads::spin_mtx, and libMesh::Elem::subactive().

{
  libmesh_assert(elem);

  const unsigned int Dim = elem->dim();

  // Only constrain elements in 2,3D.
  if (Dim == 1)
    return;

  // Only constrain active and ancestor elements
  if (elem->subactive())
    return;

  // We currently always use LAGRANGE mappings for geometry
  const FEType fe_type(elem->default_order(), LAGRANGE);

  std::vector<const Node*> my_nodes, parent_nodes;

  // Look at the element faces.  Check to see if we need to
  // build constraints.
  for (unsigned int s=0; s<elem->n_sides(); s++)
    if (elem->neighbor(s) != NULL &&
        elem->neighbor(s) != remote_elem)
      if (elem->neighbor(s)->level() < elem->level()) // constrain dofs shared between
        {                                                     // this element and ones coarser
                                                              // than this element.
          // Get pointers to the elements of interest and its parent.
          const Elem* parent = elem->parent();

          // This can't happen...  Only level-0 elements have NULL
          // parents, and no level-0 elements can be at a higher
          // level than their neighbors!
          libmesh_assert(parent);

          const AutoPtr<Elem> my_side     (elem->build_side(s));
          const AutoPtr<Elem> parent_side (parent->build_side(s));

          const unsigned int n_side_nodes = my_side->n_nodes();

          my_nodes.clear();
          my_nodes.reserve (n_side_nodes);
          parent_nodes.clear();
          parent_nodes.reserve (n_side_nodes);

          for (unsigned int n=0; n != n_side_nodes; ++n)
            my_nodes.push_back(my_side->get_node(n));

          for (unsigned int n=0; n != n_side_nodes; ++n)
            parent_nodes.push_back(parent_side->get_node(n));

          for (unsigned int my_side_n=0;
               my_side_n < n_side_nodes;
               my_side_n++)
            {
              libmesh_assert_less (my_side_n, FEInterface::n_dofs(Dim-1, fe_type, my_side->type()));

              const Node* my_node = my_nodes[my_side_n];

              // The support point of the DOF
              const Point& support_point = *my_node;

              // Figure out where my node lies on their reference element.
              const Point mapped_point = FEInterface::inverse_map(Dim-1, fe_type,
                                                                  parent_side.get(),
                                                                  support_point);

              // Compute the parent's side shape function values.
              for (unsigned int their_side_n=0;
                   their_side_n < n_side_nodes;
                   their_side_n++)
                {
                  libmesh_assert_less (their_side_n, FEInterface::n_dofs(Dim-1, fe_type, parent_side->type()));

                  const Node* their_node = parent_nodes[their_side_n];
                  libmesh_assert(their_node);

                  const Real their_value = FEInterface::shape(Dim-1,
                                                              fe_type,
                                                              parent_side->type(),
                                                              their_side_n,
                                                              mapped_point);

                  const Real their_mag = std::abs(their_value);
#ifdef DEBUG
                  // Protect for the case u_i ~= u_j,
                  // in which case i better equal j.
                  if (their_mag > 0.999)
                    {
                      libmesh_assert_equal_to (my_node, their_node);
                      libmesh_assert_less (std::abs(their_value - 1.), 0.001);
                    }
                  else
#endif
                  // To make nodal constraints useful for constructing
                  // sparsity patterns faster, we need to get EVERY
                  // POSSIBLE constraint coupling identified, even if
                  // there is no coupling in the isoparametric
                  // Lagrange case.
                  if (their_mag < 1.e-5)
                    {
                      // since we may be running this method concurretly
                      // on multiple threads we need to acquire a lock
                      // before modifying the shared constraint_row object.
                      Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);

                      // A reference to the constraint row.
                      NodeConstraintRow& constraint_row = constraints[my_node].first;

                      constraint_row.insert(std::make_pair (their_node,
                                                            0.));
                    }
                  // To get nodal coordinate constraints right, only
                  // add non-zero and non-identity values for Lagrange
                  // basis functions.
                  else // (1.e-5 <= their_mag <= .999)
                    {
                      // since we may be running this method concurretly
                      // on multiple threads we need to acquire a lock
                      // before modifying the shared constraint_row object.
                      Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);

                      // A reference to the constraint row.
                      NodeConstraintRow& constraint_row = constraints[my_node].first;

                      constraint_row.insert(std::make_pair (their_node,
                                                            their_value));
                    }
                }
            }
        }
}

static void libMesh::FEGenericBase< FEOutputType< T >::type >::compute_periodic_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const PeriodicBoundaries &	boundaries,
		const MeshBase &	mesh,
		const PointLocatorBase *	point_locator,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[static, inherited]

Computes the constraint matrix contributions (for meshes with periodic boundary conditions) corresponding to variable number var_number, using generic projections.

void libMesh::FEAbstract::compute_periodic_node_constraints	(	NodeConstraints &	constraints,
		const PeriodicBoundaries &	boundaries,
		const MeshBase &	mesh,
		const PointLocatorBase *	point_locator,
		const Elem *	elem
	)			[static, inherited]

Computes the node position constraint equation contributions (for meshes with periodic boundary conditions)

Definition at line 1021 of file fe_abstract.C.

References libMesh::Elem::active(), libMesh::PeriodicBoundaries::boundary(), libMesh::MeshBase::boundary_info, libMesh::Elem::build_side(), libMesh::Elem::default_order(), libMesh::Elem::dim(), libMesh::FEAbstract::fe_type, libMesh::PeriodicBoundaryBase::get_corresponding_pos(), libMesh::invalid_uint, libMesh::FEInterface::inverse_map(), libMeshEnums::LAGRANGE, libMesh::Elem::level(), libMesh::FEInterface::n_dofs(), libMesh::Elem::n_sides(), libMesh::PeriodicBoundaries::neighbor(), libMesh::Elem::neighbor(), libMesh::PeriodicBoundaryBase::pairedboundary, libMesh::Real, libMesh::FEInterface::shape(), and libMesh::Threads::spin_mtx.

{
  // Only bother if we truly have periodic boundaries
  if (boundaries.empty())
    return;

  libmesh_assert(elem);

  // Only constrain active elements with this method
  if (!elem->active())
    return;

  const unsigned int Dim = elem->dim();

  // We currently always use LAGRANGE mappings for geometry
  const FEType fe_type(elem->default_order(), LAGRANGE);

  std::vector<const Node*> my_nodes, neigh_nodes;

  // Look at the element faces.  Check to see if we need to
  // build constraints.
  for (unsigned int s=0; s<elem->n_sides(); s++)
    {
      if (elem->neighbor(s))
        continue;

      const std::vector<boundary_id_type>& bc_ids = mesh.boundary_info->boundary_ids (elem, s);
      for (std::vector<boundary_id_type>::const_iterator id_it=bc_ids.begin(); id_it!=bc_ids.end(); ++id_it)
        {
          const boundary_id_type boundary_id = *id_it;
          const PeriodicBoundaryBase *periodic = boundaries.boundary(boundary_id);
          if (periodic)
            {
              libmesh_assert(point_locator);

              // Get pointers to the element's neighbor.
              const Elem* neigh = boundaries.neighbor(boundary_id, *point_locator, elem, s);

              // h refinement constraints:
              // constrain dofs shared between
              // this element and ones as coarse
              // as or coarser than this element.
              if (neigh->level() <= elem->level())
                {
                  unsigned int s_neigh =
                    mesh.boundary_info->side_with_boundary_id (neigh, periodic->pairedboundary);
                  libmesh_assert_not_equal_to (s_neigh, libMesh::invalid_uint);

#ifdef LIBMESH_ENABLE_AMR
                  libmesh_assert(neigh->active());
#endif // #ifdef LIBMESH_ENABLE_AMR

                  const AutoPtr<Elem> my_side    (elem->build_side(s));
                  const AutoPtr<Elem> neigh_side (neigh->build_side(s_neigh));

                  const unsigned int n_side_nodes = my_side->n_nodes();

                  my_nodes.clear();
                  my_nodes.reserve (n_side_nodes);
                  neigh_nodes.clear();
                  neigh_nodes.reserve (n_side_nodes);

                  for (unsigned int n=0; n != n_side_nodes; ++n)
                    my_nodes.push_back(my_side->get_node(n));

                  for (unsigned int n=0; n != n_side_nodes; ++n)
                    neigh_nodes.push_back(neigh_side->get_node(n));

                  // Make sure we're not adding recursive constraints
                  // due to the redundancy in the way we add periodic
                  // boundary constraints, or adding constraints to
                  // nodes that already have AMR constraints
                  std::vector<bool> skip_constraint(n_side_nodes, false);

                  for (unsigned int my_side_n=0;
                       my_side_n < n_side_nodes;
                       my_side_n++)
                    {
                      libmesh_assert_less (my_side_n, FEInterface::n_dofs(Dim-1, fe_type, my_side->type()));

                      const Node* my_node = my_nodes[my_side_n];

                      // Figure out where my node lies on their reference element.
                      const Point neigh_point = periodic->get_corresponding_pos(*my_node);

                      const Point mapped_point = FEInterface::inverse_map(Dim-1, fe_type,
                                                                          neigh_side.get(),
                                                                          neigh_point);

                      // If we've already got a constraint on this
                      // node, then the periodic constraint is
                      // redundant
                      {
                        Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);

                        if (constraints.count(my_node))
                          {
                            skip_constraint[my_side_n] = true;
                            continue;
                          }
                      }

                      // Compute the neighbors's side shape function values.
                      for (unsigned int their_side_n=0;
                           their_side_n < n_side_nodes;
                           their_side_n++)
                        {
                          libmesh_assert_less (their_side_n, FEInterface::n_dofs(Dim-1, fe_type, neigh_side->type()));

                          const Node* their_node = neigh_nodes[their_side_n];

                          // If there's a constraint on an opposing node,
                          // we need to see if it's constrained by
                          // *our side* making any periodic constraint
                          // on us recursive
                          {
                            Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);

                            if (!constraints.count(their_node))
                              continue;

                            const NodeConstraintRow& their_constraint_row =
                              constraints[their_node].first;

                            for (unsigned int orig_side_n=0;
                                 orig_side_n < n_side_nodes;
                                 orig_side_n++)
                              {
                                libmesh_assert_less (orig_side_n, FEInterface::n_dofs(Dim-1, fe_type, my_side->type()));

                                const Node* orig_node = my_nodes[orig_side_n];

                                if (their_constraint_row.count(orig_node))
                                  skip_constraint[orig_side_n] = true;
                              }
                          }
                        }
                    }
                  for (unsigned int my_side_n=0;
                       my_side_n < n_side_nodes;
                       my_side_n++)
                    {
                      libmesh_assert_less (my_side_n, FEInterface::n_dofs(Dim-1, fe_type, my_side->type()));

                      if (skip_constraint[my_side_n])
                        continue;

                      const Node* my_node = my_nodes[my_side_n];

                      // Figure out where my node lies on their reference element.
                      const Point neigh_point = periodic->get_corresponding_pos(*my_node);

                      // Figure out where my node lies on their reference element.
                      const Point mapped_point = FEInterface::inverse_map(Dim-1, fe_type,
                                                                          neigh_side.get(),
                                                                          neigh_point);

                      for (unsigned int their_side_n=0;
                           their_side_n < n_side_nodes;
                           their_side_n++)
                        {
                          libmesh_assert_less (their_side_n, FEInterface::n_dofs(Dim-1, fe_type, neigh_side->type()));

                          const Node* their_node = neigh_nodes[their_side_n];
                          libmesh_assert(their_node);

                          const Real their_value = FEInterface::shape(Dim-1,
                                                                      fe_type,
                                                                      neigh_side->type(),
                                                                      their_side_n,
                                                                      mapped_point);

                          // since we may be running this method concurretly
                          // on multiple threads we need to acquire a lock
                          // before modifying the shared constraint_row object.
                          {
                            Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);

                            NodeConstraintRow& constraint_row =
                              constraints[my_node].first;

                            constraint_row.insert(std::make_pair(their_node,
                                                                 their_value));
                          }
                        }
                    }
                }
            }
        }
    }
}

static void libMesh::FEGenericBase< FEOutputType< T >::type >::compute_proj_constraints	(	DofConstraints &	constraints,
		DofMap &	dof_map,
		const unsigned int	variable_number,
		const Elem *	elem
	)			[static, inherited]

Computes the constraint matrix contributions (for non-conforming adapted meshes) corresponding to variable number var_number, using generic projections.

Referenced by libMesh::FE< Dim, T >::compute_constraints().

virtual void libMesh::FEGenericBase< FEOutputType< T >::type >::compute_shape_functions	(	const Elem *	elem,
		const std::vector< Point > &	qp
	)			[protected, virtual, inherited]

After having updated the jacobian and the transformation from local to global coordinates in FEAbstract::compute_map(), the first derivatives of the shape functions are transformed to global coordinates, giving dphi, dphidx, dphidy, and dphidz. This method should rarely be re-defined in derived classes, but still should be usable for children. Therefore, keep it protected.

Implements libMesh::FEAbstract.

Reimplemented in libMesh::FEXYZ< Dim >.

Referenced by libMesh::FE< Dim, T >::reinit().

void libMesh::ReferenceCounter::disable_print_counter_info

(

)

[static, inherited]

Definition at line 106 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

{
  _enable_print_counter = false;
  return;
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::dofs_on_edge	(	const Elem *const	elem,
		const Order	o,
		unsigned int	e,
		std::vector< unsigned int > &	di
	)			[inline, static]

Fills the vector di with the local degree of freedom indices associated with edge e of element elem

On a p-refined element, o should be the base order of the element.

Definition at line 89 of file fe.C.

References libMesh::Elem::is_node_on_edge(), libMesh::FE< Dim, T >::n_dofs(), libMesh::FE< Dim, T >::n_dofs_at_node(), libMesh::Elem::n_edges(), libMesh::Elem::n_nodes(), n_nodes, libMesh::Elem::p_level(), and libMesh::Elem::type().

{
  libmesh_assert(elem);
  libmesh_assert_less (e, elem->n_edges());

  di.clear();
  unsigned int nodenum = 0;
  const unsigned int n_nodes = elem->n_nodes();
  for (unsigned int n = 0; n != n_nodes; ++n)
    {
      const unsigned int n_dofs = n_dofs_at_node(elem->type(),
                                                 static_cast<Order>(o + elem->p_level()), n);
      if (elem->is_node_on_edge(n, e))
        for (unsigned int i = 0; i != n_dofs; ++i)
          di.push_back(nodenum++);
      else
        nodenum += n_dofs;
    }
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::dofs_on_side	(	const Elem *const	elem,
		const Order	o,
		unsigned int	s,
		std::vector< unsigned int > &	di
	)			[inline, static]

Fills the vector di with the local degree of freedom indices associated with side s of element elem

On a p-refined element, o should be the base order of the element.

Definition at line 63 of file fe.C.

References libMesh::Elem::is_node_on_side(), libMesh::FE< Dim, T >::n_dofs(), libMesh::FE< Dim, T >::n_dofs_at_node(), libMesh::Elem::n_nodes(), n_nodes, libMesh::Elem::n_sides(), libMesh::Elem::p_level(), and libMesh::Elem::type().

{
  libmesh_assert(elem);
  libmesh_assert_less (s, elem->n_sides());

  di.clear();
  unsigned int nodenum = 0;
  const unsigned int n_nodes = elem->n_nodes();
  for (unsigned int n = 0; n != n_nodes; ++n)
    {
      const unsigned int n_dofs = n_dofs_at_node(elem->type(),
                                                 static_cast<Order>(o + elem->p_level()), n);
      if (elem->is_node_on_side(n, s))
        for (unsigned int i = 0; i != n_dofs; ++i)
          di.push_back(nodenum++);
      else
        nodenum += n_dofs;
    }
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::edge_reinit	(	const Elem *	elem,
		const unsigned int	edge,
		const Real	tolerance = TOLERANCE,
		const std::vector< Point > *const	pts = NULL,
		const std::vector< Real > *const	weights = NULL
	)			[inline, virtual]

Reinitializes all the physical element-dependent data based on the edge. The tolerance paremeter is passed to the involved call to inverse_map(). By default the shape functions and associated data are computed at the quadrature points specified by the quadrature rule qrule, but may be any points specified on the reference side element specified in the optional argument pts.

Implements libMesh::FEAbstract.

Definition at line 248 of file fe_boundary.C.

References libMesh::FEAbstract::_fe_map, libMesh::Elem::build_edge(), libMesh::FEAbstract::elem_type, libMesh::QBase::get_points(), libMesh::FEAbstract::get_type(), libMesh::QBase::get_weights(), libMesh::QBase::init(), libMesh::FE< Dim, T >::inverse_map(), libMesh::FE< Dim, T >::last_edge, libMesh::Elem::p_level(), libMesh::FEAbstract::qrule, libMesh::FE< Dim, T >::reinit(), libMesh::FE< Dim, T >::shapes_need_reinit(), libMesh::QBase::shapes_need_reinit(), libMesh::FEAbstract::shapes_on_quadrature, and libMesh::Elem::type().

{
  libmesh_assert(elem);
  libmesh_assert (this->qrule != NULL || pts != NULL);
  // We don't do this for 1D elements!
  libmesh_assert_not_equal_to (Dim, 1);

  // Build the side of interest
  const AutoPtr<Elem> edge(elem->build_edge(e));

  // Initialize the shape functions at the user-specified
  // points
  if (pts != NULL)
    {
      // The shape functions do not correspond to the qrule
      this->shapes_on_quadrature = false;

      // Initialize the edge shape functions
      this->_fe_map->template init_edge_shape_functions<Dim> (*pts, edge.get());

      // Compute the Jacobian*Weight on the face for integration
      if (weights != NULL)
        {
          this->_fe_map->compute_edge_map (Dim, *weights, edge.get());
        }
      else
        {
          std::vector<Real> dummy_weights (pts->size(), 1.);
          this->_fe_map->compute_edge_map (Dim, dummy_weights, edge.get());
        }
    }
  // If there are no user specified points, we use the
  // quadrature rule
  else
    {
      // initialize quadrature rule
      this->qrule->init(edge->type(), elem->p_level());

      if(this->qrule->shapes_need_reinit())
        this->shapes_on_quadrature = false;

      // We might not need to reinitialize the shape functions
      if ((this->get_type() != elem->type())                   ||
          (edge->type() != static_cast<int>(last_edge))        || // Comparison between enum and unsigned, cast the unsigned to int
          this->shapes_need_reinit()                           ||
          !this->shapes_on_quadrature)
        {
          // Set the element type
          this->elem_type = elem->type();

          // Set the last_edge
          last_edge = edge->type();

          // Initialize the edge shape functions
          this->_fe_map->template init_edge_shape_functions<Dim> (this->qrule->get_points(), edge.get());
        }

      // Compute the Jacobian*Weight on the face for integration
      this->_fe_map->compute_edge_map (Dim, this->qrule->get_weights(), edge.get());

      // The shape functions correspond to the qrule
      this->shapes_on_quadrature = true;
    }

  // make a copy of the Jacobian for integration
  const std::vector<Real> JxW_int(this->_fe_map->get_JxW());

  // Find where the integration points are located on the
  // full element.
  std::vector<Point> qp; 
  this->inverse_map (elem, this->_fe_map->get_xyz(), qp, tolerance);

  // compute the shape function and derivative values
  // at the points qp
  this->reinit  (elem, &qp);

  // copy back old data
  this->_fe_map->get_JxW() = JxW_int;
}

void libMesh::ReferenceCounter::enable_print_counter_info

(

)

[static, inherited]

Methods to enable/disable the reference counter output from print_info()

Definition at line 100 of file reference_counter.C.

References libMesh::ReferenceCounter::_enable_print_counter.

{
  _enable_print_counter = true;
  return;
}

template<>

FEContinuity libMesh::FE< 3, XYZ >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 926 of file fe_xyz.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 2, XYZ >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 925 of file fe_xyz.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 1, XYZ >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 924 of file fe_xyz.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 0, XYZ >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 923 of file fe_xyz.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 3, SZABAB >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 1261 of file fe_szabab.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 2, SZABAB >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 1260 of file fe_szabab.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 1, SZABAB >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 1259 of file fe_szabab.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 0, SZABAB >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 1258 of file fe_szabab.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 3, SCALAR >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 110 of file fe_scalar.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 2, SCALAR >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 109 of file fe_scalar.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 1, SCALAR >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 108 of file fe_scalar.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 0, SCALAR >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 107 of file fe_scalar.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 3, NEDELEC_ONE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 566 of file fe_nedelec_one.C.

References libMeshEnums::H_CURL.

{ return H_CURL; }

template<>

FEContinuity libMesh::FE< 2, NEDELEC_ONE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 565 of file fe_nedelec_one.C.

References libMeshEnums::H_CURL.

{ return H_CURL; }

template<>

FEContinuity libMesh::FE< 1, NEDELEC_ONE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 564 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

FEContinuity libMesh::FE< 0, NEDELEC_ONE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 563 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

FEContinuity libMesh::FE< 3, MONOMIAL >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 403 of file fe_monomial.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 2, MONOMIAL >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 402 of file fe_monomial.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 1, MONOMIAL >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 401 of file fe_monomial.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 0, MONOMIAL >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 400 of file fe_monomial.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 3, LAGRANGE_VEC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 929 of file fe_lagrange_vec.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 2, LAGRANGE_VEC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 928 of file fe_lagrange_vec.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 1, LAGRANGE_VEC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 927 of file fe_lagrange_vec.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 0, LAGRANGE_VEC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 926 of file fe_lagrange_vec.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 3, LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 844 of file fe_lagrange.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 2, LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 843 of file fe_lagrange.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 1, LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 842 of file fe_lagrange.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 0, LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 841 of file fe_lagrange.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 3, L2_LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 485 of file fe_l2_lagrange.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 2, L2_LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 484 of file fe_l2_lagrange.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 1, L2_LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 483 of file fe_l2_lagrange.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 0, L2_LAGRANGE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 482 of file fe_l2_lagrange.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 3, L2_HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 191 of file fe_l2_hierarchic.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 2, L2_HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 190 of file fe_l2_hierarchic.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 1, L2_HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 189 of file fe_l2_hierarchic.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 0, L2_HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 188 of file fe_l2_hierarchic.C.

References libMeshEnums::DISCONTINUOUS.

{ return DISCONTINUOUS; }

template<>

FEContinuity libMesh::FE< 3, HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 362 of file fe_hierarchic.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 2, HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 361 of file fe_hierarchic.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 1, HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 360 of file fe_hierarchic.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 0, HIERARCHIC >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 359 of file fe_hierarchic.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 3, HERMITE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 349 of file fe_hermite.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 2, HERMITE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 348 of file fe_hermite.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 1, HERMITE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 347 of file fe_hermite.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 0, HERMITE >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 346 of file fe_hermite.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 3, CLOUGH >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 341 of file fe_clough.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 2, CLOUGH >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 340 of file fe_clough.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 1, CLOUGH >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 339 of file fe_clough.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 0, CLOUGH >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 338 of file fe_clough.C.

References libMeshEnums::C_ONE.

{ return C_ONE; }

template<>

FEContinuity libMesh::FE< 3, BERNSTEIN >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 426 of file fe_bernstein.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 2, BERNSTEIN >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 425 of file fe_bernstein.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 1, BERNSTEIN >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 424 of file fe_bernstein.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<>

FEContinuity libMesh::FE< 0, BERNSTEIN >::get_continuity

(

)

const [inline, virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

Definition at line 423 of file fe_bernstein.C.

References libMeshEnums::C_ZERO.

{ return C_ZERO; }

template<unsigned int Dim, FEFamily T>

virtual FEContinuity libMesh::FE< Dim, T >::get_continuity

(

)

const [virtual]

Returns:

the continuity level of the finite element.

Implements libMesh::FEAbstract.

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_curl_phi

(

)

const [inline, inherited]

Returns:

the curl of the shape function at the quadrature points.

Definition at line 238 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_curl_phi);
    calculate_curl_phi = calculate_dphiref = true; return curl_phi; }

const std::vector<Real>& libMesh::FEAbstract::get_curvatures

(

)

const [inline, inherited]

Returns:

the curvatures for use in face integration.

Definition at line 380 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_curvatures();}

const std::vector<std::vector<OutputTensor> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phi

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 304 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phi; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phideta2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 384 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phideta2; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidetadzeta

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 392 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidetadzeta; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidx2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 312 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidx2; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidxdy

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 320 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidxdy; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidxdz

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 328 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidxdz; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidxi2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 360 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidxi2; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidxideta

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 368 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidxideta; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidxidzeta

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 376 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidxidzeta; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidy2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 336 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi =  calculate_dphiref = true; return d2phidy2; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidydz

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 344 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidydz; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidz2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points.

Definition at line 352 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidz2; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_d2phidzeta2

(

)

const [inline, inherited]

Returns:

the shape function second derivatives at the quadrature points, in reference coordinates

Definition at line 400 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_d2phi);
    calculate_d2phi = calculate_dphiref = true; return d2phidzeta2; }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdeta2

(

)

const [inline, inherited]

Returns:

the second partial derivatives in eta.

Definition at line 267 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdeta2(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdetadzeta

(

)

const [inline, inherited]

Returns:

the second partial derivatives in eta-zeta.

Definition at line 297 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdetadzeta(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdxi2

(

)

const [inline, inherited]

Returns:

the second partial derivatives in xi.

Definition at line 261 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdxi2(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdxideta

(

)

const [inline, inherited]

Returns:

the second partial derivatives in xi-eta.

Definition at line 283 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdxideta(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdxidzeta

(

)

const [inline, inherited]

Returns:

the second partial derivatives in xi-zeta.

Definition at line 291 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdxidzeta(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_d2xyzdzeta2

(

)

const [inline, inherited]

Returns:

the second partial derivatives in zeta.

Definition at line 275 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_d2xyzdzeta2(); }

const std::vector<Real>& libMesh::FEAbstract::get_detadx

(

)

const [inline, inherited]

Returns:

the deta/dx entry in the transformation matrix from physical to local coordinates.

Definition at line 327 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_detadx(); }

const std::vector<Real>& libMesh::FEAbstract::get_detady

(

)

const [inline, inherited]

Returns:

the deta/dy entry in the transformation matrix from physical to local coordinates.

Definition at line 334 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_detady(); }

const std::vector<Real>& libMesh::FEAbstract::get_detadz

(

)

const [inline, inherited]

Returns:

the deta/dz entry in the transformation matrix from physical to local coordinates.

Definition at line 341 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_detadz(); }

const std::vector<std::vector<OutputDivergence> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_div_phi

(

)

const [inline, inherited]

Returns:

the divergence of the shape function at the quadrature points.

Definition at line 246 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_div_phi);
    calculate_div_phi = calculate_dphiref = true; return div_phi; }

const std::vector<OutputGradient>& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphase

(

)

const [inline, inherited]

Returns:

the global first derivative of the phase term which is used in infinite elements, evaluated at the quadrature points.

In case of the general finite element class FE this field is initialized to all zero, so that the variational formulation for an infinite element returns correct element matrices for a mesh using both finite and infinite elements.

Definition at line 418 of file fe_base.h.

      { return dphase; }

const std::vector<std::vector<OutputGradient> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphi

(

)

const [inline, inherited]

Returns:

the shape function derivatives at the quadrature points.

Definition at line 230 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphi);
    calculate_dphi = calculate_dphiref = true; return dphi; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphideta

(

)

const [inline, inherited]

Returns:

the shape function eta-derivative at the quadrature points.

Definition at line 286 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphiref);
    calculate_dphiref = true; return dphideta; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphidx

(

)

const [inline, inherited]

Returns:

the shape function x-derivative at the quadrature points.

Definition at line 254 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphi);
    calculate_dphi = calculate_dphiref = true; return dphidx; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphidxi

(

)

const [inline, inherited]

Returns:

the shape function xi-derivative at the quadrature points.

Definition at line 278 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphiref);
    calculate_dphiref = true; return dphidxi; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphidy

(

)

const [inline, inherited]

Returns:

the shape function y-derivative at the quadrature points.

Definition at line 262 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphi);
    calculate_dphi = calculate_dphiref = true; return dphidy; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphidz

(

)

const [inline, inherited]

Returns:

the shape function z-derivative at the quadrature points.

Definition at line 270 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphi);
    calculate_dphi = calculate_dphiref = true; return dphidz; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_dphidzeta

(

)

const [inline, inherited]

Returns:

the shape function zeta-derivative at the quadrature points.

Definition at line 294 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_dphiref);
    calculate_dphiref = true; return dphidzeta; }

const std::vector<Real>& libMesh::FEAbstract::get_dxidx

(

)

const [inline, inherited]

Returns:

the dxi/dx entry in the transformation matrix from physical to local coordinates.

Definition at line 306 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dxidx(); }

const std::vector<Real>& libMesh::FEAbstract::get_dxidy

(

)

const [inline, inherited]

Returns:

the dxi/dy entry in the transformation matrix from physical to local coordinates.

Definition at line 313 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dxidy(); }

const std::vector<Real>& libMesh::FEAbstract::get_dxidz

(

)

const [inline, inherited]

Returns:

the dxi/dz entry in the transformation matrix from physical to local coordinates.

Definition at line 320 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dxidz(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_dxyzdeta

(

)

const [inline, inherited]

Returns:

the element tangents in eta-direction at the quadrature points.

Definition at line 248 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dxyzdeta(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_dxyzdxi

(

)

const [inline, inherited]

Returns:

the element tangents in xi-direction at the quadrature points.

Definition at line 241 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dxyzdxi(); }

const std::vector<RealGradient>& libMesh::FEAbstract::get_dxyzdzeta

(

)

const [inline, inherited]

Returns:

the element tangents in zeta-direction at the quadrature points.

Definition at line 255 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return _fe_map->get_dxyzdzeta(); }

const std::vector<Real>& libMesh::FEAbstract::get_dzetadx

(

)

const [inline, inherited]

Returns:

the dzeta/dx entry in the transformation matrix from physical to local coordinates.

Definition at line 348 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dzetadx(); }

const std::vector<Real>& libMesh::FEAbstract::get_dzetady

(

)

const [inline, inherited]

Returns:

the dzeta/dy entry in the transformation matrix from physical to local coordinates.

Definition at line 355 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dzetady(); }

const std::vector<Real>& libMesh::FEAbstract::get_dzetadz

(

)

const [inline, inherited]

Returns:

the dzeta/dz entry in the transformation matrix from physical to local coordinates.

Definition at line 362 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_dzetadz(); }

FEFamily libMesh::FEAbstract::get_family

(

)

const [inline, inherited]

Returns:

the finite element family of this element.

Definition at line 439 of file fe_abstract.h.

References libMesh::FEType::family, and libMesh::FEAbstract::fe_type.

Referenced by libMesh::FE< Dim, T >::FE().

{ return fe_type.family; }

const FEMap& libMesh::FEAbstract::get_fe_map

(

)

const [inline, inherited]

Returns:

the mapping object

Definition at line 444 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map, and libMesh::AutoPtr< Tp >::get().

Referenced by libMesh::InfFE< Dim, T_radial, T_map >::combine_base_radial(), libMesh::InfFE< Dim, T_radial, T_map >::init_face_shape_functions(), libMesh::HCurlFETransformation< OutputShape >::map_curl(), libMesh::H1FETransformation< OutputShape >::map_curl(), libMesh::H1FETransformation< OutputShape >::map_d2phi(), libMesh::H1FETransformation< OutputShape >::map_div(), libMesh::H1FETransformation< OutputShape >::map_dphi(), and libMesh::HCurlFETransformation< OutputShape >::map_phi().

{ return *_fe_map.get(); }

FEType libMesh::FEAbstract::get_fe_type

(

)

const [inline, inherited]

Returns:

the FE Type (approximation order and family) of the finite element.

Definition at line 418 of file fe_abstract.h.

References libMesh::FEAbstract::fe_type.

Referenced by libMesh::FEMContext::build_new_fe(), libMesh::HCurlFETransformation< OutputShape >::map_phi(), libMesh::H1FETransformation< OutputShape >::map_phi(), and libMesh::ProjectFEMSolution::operator()().

{ return fe_type; }

std::string libMesh::ReferenceCounter::get_info

(

)

[static, inherited]

Gets a string containing the reference information.

Definition at line 47 of file reference_counter.C.

References libMesh::ReferenceCounter::_counts, and libMesh::Quality::name().

Referenced by libMesh::ReferenceCounter::print_info().

{
#if defined(LIBMESH_ENABLE_REFERENCE_COUNTING) && defined(DEBUG)

  std::ostringstream oss;

  oss << '\n'
      << " ---------------------------------------------------------------------------- \n"
      << "| Reference count information                                                |\n"
      << " ---------------------------------------------------------------------------- \n";

  for (Counts::iterator it = _counts.begin();
       it != _counts.end(); ++it)
    {
      const std::string name(it->first);
      const unsigned int creations    = it->second.first;
      const unsigned int destructions = it->second.second;

      oss << "| " << name << " reference count information:\n"
          << "|  Creations:    " << creations    << '\n'
          << "|  Destructions: " << destructions << '\n';
    }

  oss << " ---------------------------------------------------------------------------- \n";

  return oss.str();

#else

  return "";

#endif
}

const std::vector<Real>& libMesh::FEAbstract::get_JxW

(

)

const [inline, inherited]

Returns:

the element Jacobian times the quadrature weight for each quadrature point.

Definition at line 234 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

Referenced by libMesh::ExactErrorEstimator::find_squared_element_error(), and libMesh::ProjectFEMSolution::operator()().

  { return this->_fe_map->get_JxW(); }

const std::vector<Point>& libMesh::FEAbstract::get_normals

(

)

const [inline, inherited]

Returns:

the normal vectors for face integration.

Definition at line 374 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_normals(); }

Order libMesh::FEAbstract::get_order

(

)

const [inline, inherited]

Returns:

the approximation order of the finite element.

Definition at line 423 of file fe_abstract.h.

References libMesh::FEAbstract::_p_level, libMesh::FEAbstract::fe_type, and libMesh::FEType::order.

Referenced by libMesh::FEXYZ< Dim >::compute_face_values(), libMesh::FEXYZ< Dim >::init_shape_functions(), and libMesh::FE< Dim, T >::init_shape_functions().

{ return static_cast<Order>(fe_type.order + _p_level); }

unsigned int libMesh::FEAbstract::get_p_level

(

)

const [inline, inherited]

Returns:

the p refinement level that the current shape functions have been calculated for.

Definition at line 413 of file fe_abstract.h.

References libMesh::FEAbstract::_p_level.

{ return _p_level; }

const std::vector<std::vector<OutputShape> >& libMesh::FEGenericBase< FEOutputType< T >::type >::get_phi

(

)

const [inline, inherited]

Returns:

the shape function values at the quadrature points on the element.

Definition at line 222 of file fe_base.h.

  { libmesh_assert(!calculations_started || calculate_phi);
    calculate_phi = true; return phi; }

void libMesh::FEAbstract::get_refspace_nodes	(	const ElemType	t,
		std::vector< Point > &	nodes
	)			[static, inherited]

returns the reference space nodes coordinates given the element type

Definition at line 415 of file fe_abstract.C.

References libMeshEnums::EDGE2, libMeshEnums::EDGE3, libMesh::err, libMeshEnums::HEX20, libMeshEnums::HEX27, libMeshEnums::HEX8, libMeshEnums::PRISM15, libMeshEnums::PRISM18, libMeshEnums::PRISM6, libMeshEnums::PYRAMID5, libMeshEnums::QUAD4, libMeshEnums::QUAD8, libMeshEnums::QUAD9, libMeshEnums::TET10, libMeshEnums::TET4, libMeshEnums::TRI3, and libMeshEnums::TRI6.

Referenced by libMesh::FE< Dim, T >::side_map().

{
  switch(itemType)
  {
    case EDGE2:
    {
       nodes.resize(2);
       nodes[0] = Point (-1.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       return;
    }
    case EDGE3:
    {
       nodes.resize(3);
       nodes[0] = Point (-1.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       nodes[2] = Point (0.,0.,0.);
       return;
    }
    case TRI3:
    {
       nodes.resize(3);
       nodes[0] = Point (0.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       nodes[2] = Point (0.,1.,0.);
       return;
    }
    case TRI6:
    {
       nodes.resize(6);
       nodes[0] = Point (0.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       nodes[2] = Point (0.,1.,0.);
       nodes[3] = Point (.5,0.,0.);
       nodes[4] = Point (.5,.5,0.);
       nodes[5] = Point (0.,.5,0.);
       return;
    }
    case QUAD4:
    {
       nodes.resize(4);
       nodes[0] = Point (-1.,-1.,0.);
       nodes[1] = Point (1.,-1.,0.);
       nodes[2] = Point (1.,1.,0.);
       nodes[3] = Point (-1.,1.,0.);
       return;
    }
    case QUAD8:
    {
       nodes.resize(8);
       nodes[0] = Point (-1.,-1.,0.);
       nodes[1] = Point (1.,-1.,0.);
       nodes[2] = Point (1.,1.,0.);
       nodes[3] = Point (-1.,1.,0.);
       nodes[4] = Point (0.,-1.,0.);
       nodes[5] = Point (1.,0.,0.);
       nodes[6] = Point (0.,1.,0.);
       nodes[7] = Point (-1.,0.,0.);
       return;
    }
    case QUAD9:
    {
       nodes.resize(9);
       nodes[0] = Point (-1.,-1.,0.);
       nodes[1] = Point (1.,-1.,0.);
       nodes[2] = Point (1.,1.,0.);
       nodes[3] = Point (-1.,1.,0.);
       nodes[4] = Point (0.,-1.,0.);
       nodes[5] = Point (1.,0.,0.);
       nodes[6] = Point (0.,1.,0.);
       nodes[7] = Point (-1.,0.,0.);
       nodes[8] = Point (0.,0.,0.);
       return;
    }
    case TET4:
    {
       nodes.resize(4);
       nodes[0] = Point (0.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       nodes[2] = Point (0.,1.,0.);
       nodes[3] = Point (0.,0.,1.);
       return;
    }
    case TET10:
    {
       nodes.resize(10);
       nodes[0] = Point (0.,0.,0.);
       nodes[1] = Point (1.,0.,0.);
       nodes[2] = Point (0.,1.,0.);
       nodes[3] = Point (0.,0.,1.);
       nodes[4] = Point (.5,0.,0.);
       nodes[5] = Point (.5,.5,0.);
       nodes[6] = Point (0.,.5,0.);
       nodes[7] = Point (0.,0.,.5);
       nodes[8] = Point (.5,0.,.5);
       nodes[9] = Point (0.,.5,.5);
       return;
    }
    case HEX8:
    {
       nodes.resize(8);
       nodes[0] = Point (-1.,-1.,-1.);
       nodes[1] = Point (1.,-1.,-1.);
       nodes[2] = Point (1.,1.,-1.);
       nodes[3] = Point (-1.,1.,-1.);
       nodes[4] = Point (-1.,-1.,1.);
       nodes[5] = Point (1.,-1.,1.);
       nodes[6] = Point (1.,1.,1.);
       nodes[7] = Point (-1.,1.,1.);
       return;
    }
    case HEX20:
    {
       nodes.resize(20);
       nodes[0] = Point (-1.,-1.,-1.);
       nodes[1] = Point (1.,-1.,-1.);
       nodes[2] = Point (1.,1.,-1.);
       nodes[3] = Point (-1.,1.,-1.);
       nodes[4] = Point (-1.,-1.,1.);
       nodes[5] = Point (1.,-1.,1.);
       nodes[6] = Point (1.,1.,1.);
       nodes[7] = Point (-1.,1.,1.);
       nodes[8] = Point (0.,-1.,-1.);
       nodes[9] = Point (1.,0.,-1.);
       nodes[10] = Point (0.,1.,-1.);
       nodes[11] = Point (-1.,0.,-1.);
       nodes[12] = Point (-1.,-1.,0.);
       nodes[13] = Point (1.,-1.,0.);
       nodes[14] = Point (1.,1.,0.);
       nodes[15] = Point (-1.,1.,0.);
       nodes[16] = Point (0.,-1.,1.);
       nodes[17] = Point (1.,0.,1.);
       nodes[18] = Point (0.,1.,1.);
       nodes[19] = Point (-1.,0.,1.);
       return;
    }
    case HEX27:
    {
       nodes.resize(27);
       nodes[0] = Point (-1.,-1.,-1.);
       nodes[1] = Point (1.,-1.,-1.);
       nodes[2] = Point (1.,1.,-1.);
       nodes[3] = Point (-1.,1.,-1.);
       nodes[4] = Point (-1.,-1.,1.);
       nodes[5] = Point (1.,-1.,1.);
       nodes[6] = Point (1.,1.,1.);
       nodes[7] = Point (-1.,1.,1.);
       nodes[8] = Point (0.,-1.,-1.);
       nodes[9] = Point (1.,0.,-1.);
       nodes[10] = Point (0.,1.,-1.);
       nodes[11] = Point (-1.,0.,-1.);
       nodes[12] = Point (-1.,-1.,0.);
       nodes[13] = Point (1.,-1.,0.);
       nodes[14] = Point (1.,1.,0.);
       nodes[15] = Point (-1.,1.,0.);
       nodes[16] = Point (0.,-1.,1.);
       nodes[17] = Point (1.,0.,1.);
       nodes[18] = Point (0.,1.,1.);
       nodes[19] = Point (-1.,0.,1.);
       nodes[20] = Point (0.,0.,-1.);
       nodes[21] = Point (0.,-1.,0.);
       nodes[22] = Point (1.,0.,0.);
       nodes[23] = Point (0.,1.,0.);
       nodes[24] = Point (-1.,0.,0.);
       nodes[25] = Point (0.,0.,1.);
       nodes[26] = Point (0.,0.,0.);
       return;
    }
    case PRISM6:
    {
       nodes.resize(6);
       nodes[0] = Point (0.,0.,-1.);
       nodes[1] = Point (1.,0.,-1.);
       nodes[2] = Point (0.,1.,-1.);
       nodes[3] = Point (0.,0.,1.);
       nodes[4] = Point (1.,0.,1.);
       nodes[5] = Point (0.,1.,1.);
       return;
    }
    case PRISM15:
    {
       nodes.resize(15);
       nodes[0] = Point (0.,0.,-1.);
       nodes[1] = Point (1.,0.,-1.);
       nodes[2] = Point (0.,1.,-1.);
       nodes[3] = Point (0.,0.,1.);
       nodes[4] = Point (1.,0.,1.);
       nodes[5] = Point (0.,1.,1.);
       nodes[6] = Point (.5,0.,-1.);
       nodes[7] = Point (.5,.5,-1.);
       nodes[8] = Point (0.,.5,-1.);
       nodes[9] = Point (0.,0.,0.);
       nodes[10] = Point (1.,0.,0.);
       nodes[11] = Point (0.,1.,0.);
       nodes[12] = Point (.5,0.,1.);
       nodes[13] = Point (.5,.5,1.);
       nodes[14] = Point (0.,.5,1.);
       return;
    }
    case PRISM18:
    {
       nodes.resize(18);
       nodes[0] = Point (0.,0.,-1.);
       nodes[1] = Point (1.,0.,-1.);
       nodes[2] = Point (0.,1.,-1.);
       nodes[3] = Point (0.,0.,1.);
       nodes[4] = Point (1.,0.,1.);
       nodes[5] = Point (0.,1.,1.);
       nodes[6] = Point (.5,0.,-1.);
       nodes[7] = Point (.5,.5,-1.);
       nodes[8] = Point (0.,.5,-1.);
       nodes[9] = Point (0.,0.,0.);
       nodes[10] = Point (1.,0.,0.);
       nodes[11] = Point (0.,1.,0.);
       nodes[12] = Point (.5,0.,1.);
       nodes[13] = Point (.5,.5,1.);
       nodes[14] = Point (0.,.5,1.);
       nodes[15] = Point (.5,0.,0.);
       nodes[16] = Point (.5,.5,0.);
       nodes[17] = Point (0.,.5,0.);
       return;
    }
    case PYRAMID5:
    {
       nodes.resize(5);
       nodes[0] = Point (-1.,-1.,0.);
       nodes[1] = Point (1.,-1.,0.);
       nodes[2] = Point (1.,1.,0.);
       nodes[3] = Point (-1.,1.,0.);
       nodes[4] = Point (-1.,-1.,1.);
       return;
    }
    default:
    {
      libMesh::err << "ERROR: Unknown element type " << itemType << std::endl;
      libmesh_error();
    }
  }
  return;
}

const std::vector<RealGradient>& libMesh::FEGenericBase< FEOutputType< T >::type >::get_Sobolev_dweight

(

)

const [inline, inherited]

Returns:

the first global derivative of the multiplicative weight at each quadrature point. See get_Sobolev_weight() for details. In case of FE initialized to all zero.

Definition at line 442 of file fe_base.h.

      { return dweight; }

const std::vector<Real>& libMesh::FEGenericBase< FEOutputType< T >::type >::get_Sobolev_weight

(

)

const [inline, inherited]

Returns:

the multiplicative weight at each quadrature point. This weight is used for certain infinite element weak formulations, so that weighted Sobolev spaces are used for the trial function space. This renders the variational form easily computable.

In case of the general finite element class FE this field is initialized to all ones, so that the variational formulation for an infinite element returns correct element matrices for a mesh using both finite and infinite elements.

Definition at line 434 of file fe_base.h.

      { return weight; }

const std::vector<std::vector<Point> >& libMesh::FEAbstract::get_tangents

(

)

const [inline, inherited]

Returns:

the tangent vectors for face integration.

Definition at line 368 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

  { return this->_fe_map->get_tangents(); }

ElemType libMesh::FEAbstract::get_type

(

)

const [inline, inherited]

Returns:

the element type that the current shape functions have been calculated for. Useful in determining when shape functions must be recomputed.

Definition at line 407 of file fe_abstract.h.

References libMesh::FEAbstract::elem_type.

Referenced by libMesh::FEXYZ< Dim >::compute_face_values(), libMesh::FE< Dim, T >::edge_reinit(), libMesh::FEXYZ< Dim >::init_shape_functions(), libMesh::FE< Dim, T >::init_shape_functions(), and libMesh::InfFE< Dim, T_radial, T_map >::reinit().

{ return elem_type; }

const std::vector<Point>& libMesh::FEAbstract::get_xyz

(

)

const [inline, inherited]

Returns:

the xyz spatial locations of the quadrature points on the element.

Definition at line 227 of file fe_abstract.h.

References libMesh::FEAbstract::_fe_map.

Referenced by libMesh::FEXYZ< Dim >::compute_shape_functions(), libMesh::ExactErrorEstimator::find_squared_element_error(), and libMesh::ProjectFEMSolution::operator()().

  { return this->_fe_map->get_xyz(); }

void libMesh::ReferenceCounter::increment_constructor_count

(

const std::string &

name

)

[inline, protected, inherited]

Increments the construction counter. Should be called in the constructor of any derived class that will be reference counted.

Definition at line 163 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::ReferenceCountedObject().

{
  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
  std::pair<unsigned int, unsigned int>& p = _counts[name];

  p.first++;
}

void libMesh::ReferenceCounter::increment_destructor_count

(

const std::string &

name

)

[inline, protected, inherited]

Increments the destruction counter. Should be called in the destructor of any derived class that will be reference counted.

Definition at line 176 of file reference_counter.h.

References libMesh::ReferenceCounter::_counts, and libMesh::Threads::spin_mtx.

Referenced by libMesh::ReferenceCountedObject< RBParametrized >::~ReferenceCountedObject().

{
  Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
  std::pair<unsigned int, unsigned int>& p = _counts[name];

  p.second++;
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::init_base_shape_functions	(	const std::vector< Point > &	qp,
		const Elem *	e
	)			[inline, protected, virtual]

Initialize the data fields for the base of an an infinite element.

Implements libMesh::FEGenericBase< FEOutputType< T >::type >.

Definition at line 530 of file fe.C.

References libMesh::FEAbstract::_fe_map, libMesh::FEAbstract::calculate_d2phi, libMesh::FEAbstract::calculate_dphi, libMesh::FEAbstract::calculate_phi, libMesh::FEAbstract::elem_type, libMesh::FE< Dim, T >::init_shape_functions(), and libMesh::Elem::type().

{
  // I don't understand infinite elements well enough to risk
  // calculating too little.  :-(  RHS
  this->calculate_phi = this->calculate_dphi = this->calculate_d2phi = true;

  this->elem_type = e->type();
  this->_fe_map->template init_reference_to_physical_map<Dim>(qp, e);
  init_shape_functions(qp, e);
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::init_shape_functions	(	const std::vector< Point > &	qp,
		const Elem *	e
	)			[inline, protected, virtual]

Update the various member data fields phi, dphidxi, dphideta, dphidzeta, etc. for the current element. These data will be computed at the points qp, which are generally (but need not be) the quadrature points.

Reimplemented in libMesh::FEXYZ< Dim >.

Definition at line 241 of file fe.C.

References libMesh::FEAbstract::calculate_curl_phi, libMesh::FEAbstract::calculate_d2phi, libMesh::FEAbstract::calculate_div_phi, libMesh::FEAbstract::calculate_dphi, libMesh::FEAbstract::calculate_dphiref, libMesh::FEAbstract::calculate_phi, libMesh::FEAbstract::calculations_started, libMesh::FEGenericBase< FEOutputType< T >::type >::curl_phi, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phi, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phideta2, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidetadzeta, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidx2, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidxdy, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidxdz, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidxi2, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidxideta, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidxidzeta, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidy2, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidydz, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidz2, libMesh::FEGenericBase< FEOutputType< T >::type >::d2phidzeta2, libMesh::FEGenericBase< FEOutputType< T >::type >::div_phi, libMesh::FEGenericBase< FEOutputType< T >::type >::dphase, libMesh::FEGenericBase< FEOutputType< T >::type >::dphi, libMesh::FEGenericBase< FEOutputType< T >::type >::dphideta, libMesh::FEGenericBase< FEOutputType< T >::type >::dphidx, libMesh::FEGenericBase< FEOutputType< T >::type >::dphidxi, libMesh::FEGenericBase< FEOutputType< T >::type >::dphidy, libMesh::FEGenericBase< FEOutputType< T >::type >::dphidz, libMesh::FEGenericBase< FEOutputType< T >::type >::dphidzeta, libMesh::FEGenericBase< FEOutputType< T >::type >::dweight, libMesh::FEAbstract::fe_type, libMesh::FEInterface::field_type(), libMesh::FEAbstract::get_order(), libMesh::FEAbstract::get_type(), libMesh::FE< Dim, T >::n_shape_functions(), libMesh::FEType::order, libMesh::FEGenericBase< FEOutputType< T >::type >::phi, libMesh::FE< Dim, T >::shape_deriv(), libMesh::FE< Dim, T >::shape_second_deriv(), libMeshEnums::TYPE_VECTOR, and libMesh::FEGenericBase< FEOutputType< T >::type >::weight.

Referenced by libMesh::FE< Dim, T >::init_base_shape_functions(), and libMesh::FE< Dim, T >::reinit().

{
  libmesh_assert(elem);
  this->calculations_started = true;

  // If the user forgot to request anything, we'll be safe and
  // calculate everything:
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
  if (!this->calculate_phi && !this->calculate_dphi && !this->calculate_d2phi 
      && !this->calculate_curl_phi && !this->calculate_div_phi)
    {
      this->calculate_phi = this->calculate_dphi = this->calculate_d2phi = this->calculate_dphiref = true;
      if( FEInterface::field_type(T) == TYPE_VECTOR )
        {
          this->calculate_curl_phi = true;
          this->calculate_div_phi  = true;
        }
    }
#else
  if (!this->calculate_phi && !this->calculate_dphi && !this->calculate_curl_phi && !this->calculate_div_phi)
    {
      this->calculate_phi = this->calculate_dphi = this->calculate_dphiref = true;
      if( FEInterface::field_type(T) == TYPE_VECTOR )
        {
          this->calculate_curl_phi = true;
          this->calculate_div_phi  = true;
        }
    }
#endif // LIBMESH_ENABLE_SECOND_DERIVATIVES

  // Start logging the shape function initialization
  START_LOG("init_shape_functions()", "FE");


  // The number of quadrature points.
  const unsigned int n_qp = libmesh_cast_int<unsigned int>(qp.size());

  // Number of shape functions in the finite element approximation
  // space.
  const unsigned int n_approx_shape_functions =
    this->n_shape_functions(this->get_type(),
                            this->get_order());

  // resize the vectors to hold current data
  // Phi are the shape functions used for the FE approximation
  // Phi_map are the shape functions used for the FE mapping
  if (this->calculate_phi)
    this->phi.resize     (n_approx_shape_functions);

  if (this->calculate_dphi)
    {
      this->dphi.resize    (n_approx_shape_functions);
      this->dphidx.resize  (n_approx_shape_functions);
      this->dphidy.resize  (n_approx_shape_functions);
      this->dphidz.resize  (n_approx_shape_functions);
    }

  if(this->calculate_dphiref)
    {
      if (Dim > 0)
        this->dphidxi.resize (n_approx_shape_functions);

      if (Dim > 1)
        this->dphideta.resize      (n_approx_shape_functions);

      if (Dim > 2)
        this->dphidzeta.resize     (n_approx_shape_functions);
    }
  
  if( this->calculate_curl_phi && (FEInterface::field_type(T) == TYPE_VECTOR) )
    this->curl_phi.resize(n_approx_shape_functions);

  if( this->calculate_div_phi && (FEInterface::field_type(T) == TYPE_VECTOR) )
    this->div_phi.resize(n_approx_shape_functions);

#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
  if (this->calculate_d2phi)
    {
      this->d2phi.resize     (n_approx_shape_functions);
      this->d2phidx2.resize  (n_approx_shape_functions);
      this->d2phidxdy.resize (n_approx_shape_functions);
      this->d2phidxdz.resize (n_approx_shape_functions);
      this->d2phidy2.resize  (n_approx_shape_functions);
      this->d2phidydz.resize (n_approx_shape_functions);
      this->d2phidz2.resize  (n_approx_shape_functions);

      if (Dim > 0)
        this->d2phidxi2.resize (n_approx_shape_functions);

      if (Dim > 1)
        {
          this->d2phidxideta.resize (n_approx_shape_functions);
          this->d2phideta2.resize   (n_approx_shape_functions);
        }
      if (Dim > 2)
        {
          this->d2phidxidzeta.resize  (n_approx_shape_functions);
          this->d2phidetadzeta.resize (n_approx_shape_functions);
          this->d2phidzeta2.resize    (n_approx_shape_functions);
        }
    }
#endif // ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES

  for (unsigned int i=0; i<n_approx_shape_functions; i++)
    {
      if (this->calculate_phi)
        this->phi[i].resize         (n_qp);
      if (this->calculate_dphi)
        {
          this->dphi[i].resize        (n_qp);
          this->dphidx[i].resize      (n_qp);
          this->dphidy[i].resize      (n_qp);
          this->dphidz[i].resize      (n_qp);
        }
      
      if(this->calculate_dphiref)
        {
          if (Dim > 0)
            this->dphidxi[i].resize(n_qp);

          if (Dim > 1)
            this->dphideta[i].resize(n_qp);

          if (Dim > 2)
            this->dphidzeta[i].resize(n_qp);
        }

      if(this->calculate_curl_phi && (FEInterface::field_type(T) == TYPE_VECTOR) )
        this->curl_phi[i].resize(n_qp);

      if(this->calculate_div_phi && (FEInterface::field_type(T) == TYPE_VECTOR) )
        this->div_phi[i].resize(n_qp);

#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
      if (this->calculate_d2phi)
        {
          this->d2phi[i].resize     (n_qp);
          this->d2phidx2[i].resize  (n_qp);
          this->d2phidxdy[i].resize (n_qp);
          this->d2phidxdz[i].resize (n_qp);
          this->d2phidy2[i].resize  (n_qp);
          this->d2phidydz[i].resize (n_qp);
          this->d2phidz2[i].resize  (n_qp);
          if (Dim > 0)
            this->d2phidxi2[i].resize (n_qp);
          if (Dim > 1)
            {
              this->d2phidxideta[i].resize (n_qp);
              this->d2phideta2[i].resize   (n_qp);
            }
          if (Dim > 2)
            {
              this->d2phidxidzeta[i].resize  (n_qp);
              this->d2phidetadzeta[i].resize (n_qp);
              this->d2phidzeta2[i].resize    (n_qp);
            }
        }
#endif // ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
    }


#ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS
  //------------------------------------------------------------
  // Initialize the data fields, which should only be used for infinite
  // elements, to some sensible values, so that using a FE with the
  // variational formulation of an InfFE, correct element matrices are
  // returned

 {
    this->weight.resize  (n_qp);
    this->dweight.resize (n_qp);
    this->dphase.resize  (n_qp);

    for (unsigned int p=0; p<n_qp; p++)
      {
        this->weight[p] = 1.;
        this->dweight[p].zero();
        this->dphase[p].zero();
      }

 }
#endif // ifdef LIBMESH_ENABLE_INFINITE_ELEMENTS

  switch (Dim)
    {

      //------------------------------------------------------------
      // 0D
    case 0:
      {
        break;
      }

      //------------------------------------------------------------
      // 1D
    case 1:
      {
        // Compute the value of the approximation shape function i at quadrature point p
        if (this->calculate_dphiref)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              this->dphidxi[i][p]  = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 0, qp[p]);
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
        if (this->calculate_d2phi)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              this->d2phidxi2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 0, qp[p]);
#endif // ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES

        break;
      }



      //------------------------------------------------------------
      // 2D
    case 2:
      {
        // Compute the value of the approximation shape function i at quadrature point p
        if (this->calculate_dphiref)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              {
                this->dphidxi[i][p]  = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 0, qp[p]);
                this->dphideta[i][p] = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 1, qp[p]);
              }
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
        if (this->calculate_d2phi)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              {
                this->d2phidxi2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 0, qp[p]);
                this->d2phidxideta[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 1, qp[p]);
                this->d2phideta2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 2, qp[p]);
              }
#endif // ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES

        
        break;
      }



      //------------------------------------------------------------
      // 3D
    case 3:
      {
        // Compute the value of the approximation shape function i at quadrature point p
        if (this->calculate_dphiref)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              {
                this->dphidxi[i][p]   = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 0, qp[p]);
                this->dphideta[i][p]  = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 1, qp[p]);
                this->dphidzeta[i][p] = FE<Dim,T>::shape_deriv (elem, this->fe_type.order, i, 2, qp[p]);
              }
#ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
        if (this->calculate_d2phi)
          for (unsigned int i=0; i<n_approx_shape_functions; i++)
            for (unsigned int p=0; p<n_qp; p++)
              {
                this->d2phidxi2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 0, qp[p]);
                this->d2phidxideta[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 1, qp[p]);
                this->d2phideta2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 2, qp[p]);
                this->d2phidxidzeta[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 3, qp[p]);
                this->d2phidetadzeta[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 4, qp[p]);
                this->d2phidzeta2[i][p] = FE<Dim,T>::shape_second_deriv (elem, this->fe_type.order, i, 5, qp[p]);
              }
#endif // ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES

        break;
      }


    default:
      libmesh_error();
    }

  // Stop logging the shape function initialization
  STOP_LOG("init_shape_functions()", "FE");
}

template<unsigned int Dim, FEFamily T>

void libMesh::FE< Dim, T >::inverse_map	(	const Elem *	elem,
		const std::vector< Point > &	physical_points,
		std::vector< Point > &	reference_points,
		const Real	tolerance = TOLERANCE,
		const bool	secure = true
	)			[inline, static]

Takes a number points in physical space (in the physical_points vector) and finds their location on the reference element for the input element elem. The values on the reference element are returned in the vector reference_points. The optional parameter tolerance defines how close is "good enough." The map inversion iteration computes the sequence , and the iteration is terminated when

Definition at line 1208 of file fe_map.C.

References libMesh::TypeVector< T >::size().

{
  // The number of points to find the
  // inverse map of
  const std::size_t n_points = physical_points.size();

  // Resize the vector to hold the points
  // on the reference element
  reference_points.resize(n_points);

  // Find the coordinates on the reference
  // element of each point in physical space
  for (std::size_t p=0; p<n_points; p++)
    reference_points[p] =
      FE<Dim,T>::inverse_map (elem, physical_points[p], tolerance, secure);
}

template<unsigned int Dim, FEFamily T>

Point libMesh::FE< Dim, T >::inverse_map	(	const Elem *	elem,
		const Point &	p,
		const Real	tolerance = TOLERANCE,
		const bool	secure = true
	)			[inline, static]

Returns:

the location (on the reference element) of the point p located in physical space. This function requires inverting the (possibly nonlinear) transformation map, so it is not trivial. The optional parameter tolerance defines how close is "good enough." The map inversion iteration computes the sequence , and the iteration is terminated when

Definition at line 818 of file fe_map.C.

References libMesh::TypeVector< T >::add(), libMesh::err, libMesh::DofObject::id(), libMesh::FE< Dim, T >::map(), libMesh::FE< Dim, T >::map_eta(), libMesh::FE< Dim, T >::map_xi(), libMesh::FE< Dim, T >::map_zeta(), libMesh::FEAbstract::on_reference_element(), libMesh::Elem::print_info(), libMesh::Real, libMesh::TypeVector< T >::size(), and libMesh::Elem::type().

Referenced by libMesh::FE< Dim, T >::edge_reinit().

{
  libmesh_assert(elem);
  libmesh_assert_greater_equal (tolerance, 0.);


  // Start logging the map inversion.
  START_LOG("inverse_map()", "FE");

  // How much did the point on the reference
  // element change by in this Newton step?
  Real inverse_map_error = 0.;

  //  The point on the reference element.  This is
  //  the "initial guess" for Newton's method.  The
  //  centroid seems like a good idea, but computing
  //  it is a little more intensive than, say taking
  //  the zero point.
  //
  //  Convergence should be insensitive of this choice
  //  for "good" elements.
  Point p; // the zero point.  No computation required

  //  The number of iterations in the map inversion process.
  unsigned int cnt = 0;

  //  The number of iterations after which we give up and declare
  //  divergence
  const unsigned int max_cnt = 10;

  //  The distance (in master element space) beyond which we give up
  //  and declare divergence.  This is no longer used...
  // Real max_step_length = 4.;



  //  Newton iteration loop.
  do
    {
      //  Where our current iterate \p p maps to.
      const Point physical_guess = FE<Dim,T>::map (elem, p);

      //  How far our current iterate is from the actual point.
      const Point delta = physical_point - physical_guess;

      //  Increment in current iterate \p p, will be computed.
      Point dp;


      //  The form of the map and how we invert it depends
      //  on the dimension that we are in.
      switch (Dim)
        {
          // ------------------------------------------------------------------
          //  0D map inversion is trivial
        case 0:
          {
            break;
          }

          // ------------------------------------------------------------------
          //  1D map inversion
          //
          //  Here we find the point on a 1D reference element that maps to
          //  the point \p physical_point in the domain.  This is a bit tricky
          //  since we do not want to assume that the point \p physical_point
          //  is also in a 1D domain.  In particular, this method might get
          //  called on the edge of a 3D element, in which case
          //  \p physical_point actually lives in 3D.
        case 1:
          {
            const Point dxi = FE<Dim,T>::map_xi (elem, p);

            //  Newton's method in this case looks like
            //
            //  {X} - {X_n} = [J]*dp
            //
            //  Where {X}, {X_n} are 3x1 vectors, [J] is a 3x1 matrix
            //  d(x,y,z)/dxi, and we seek dp, a scalar.  Since the above
            //  system is either overdetermined or rank-deficient, we will
            //  solve the normal equations for this system
            //
            //  [J]^T ({X} - {X_n}) = [J]^T [J] {dp}
            //
            //  which involves the trivial inversion of the scalar
            //  G = [J]^T [J]
            const Real G = dxi*dxi;

            if (secure)
              libmesh_assert_greater (G, 0.);

            const Real Ginv = 1./G;

            const Real  dxidelta = dxi*delta;

            dp(0) = Ginv*dxidelta;

            // No master elements have radius > 4, but sometimes we
            // can take a step that big while still converging
            // if (secure)
              // libmesh_assert_less (dp.size(), max_step_length);

            break;
          }



          // ------------------------------------------------------------------
          //  2D map inversion
          //
          //  Here we find the point on a 2D reference element that maps to
          //  the point \p physical_point in the domain.  This is a bit tricky
          //  since we do not want to assume that the point \p physical_point
          //  is also in a 2D domain.  In particular, this method might get
          //  called on the face of a 3D element, in which case
          //  \p physical_point actually lives in 3D.
        case 2:
          {
            const Point dxi  = FE<Dim,T>::map_xi  (elem, p);
            const Point deta = FE<Dim,T>::map_eta (elem, p);

            //  Newton's method in this case looks like
            //
            //  {X} - {X_n} = [J]*{dp}
            //
            //  Where {X}, {X_n} are 3x1 vectors, [J] is a 3x2 matrix
            //  d(x,y,z)/d(xi,eta), and we seek {dp}, a 2x1 vector.  Since
            //  the above system is either overdermined or rank-deficient,
            //  we will solve the normal equations for this system
            //
            //  [J]^T ({X} - {X_n}) = [J]^T [J] {dp}
            //
            //  which involves the inversion of the 2x2 matrix
            //  [G] = [J]^T [J]
            const Real
              G11 = dxi*dxi,  G12 = dxi*deta,
              G21 = dxi*deta, G22 = deta*deta;


            const Real det = (G11*G22 - G12*G21);

            if (secure)
              libmesh_assert_not_equal_to (det, 0.);

            const Real inv_det = 1./det;

            const Real
              Ginv11 =  G22*inv_det,
              Ginv12 = -G12*inv_det,

              Ginv21 = -G21*inv_det,
              Ginv22 =  G11*inv_det;


            const Real  dxidelta  = dxi*delta;
            const Real  detadelta = deta*delta;

            dp(0) = (Ginv11*dxidelta + Ginv12*detadelta);
            dp(1) = (Ginv21*dxidelta + Ginv22*detadelta);

            // No master elements have radius > 4, but sometimes we
            // can take a step that big while still converging
            // if (secure)
              // libmesh_assert_less (dp.size(), max_step_length);

            break;
          }



          // ------------------------------------------------------------------
          //  3D map inversion
          //
          //  Here we find the point in a 3D reference element that maps to
          //  the point \p physical_point in a 3D domain. Nothing special
          //  has to happen here, since (unless the map is singular because
          //  you have a BAD element) the map will be invertable and we can
          //  apply Newton's method directly.
        case 3:
          {
            const Point dxi   = FE<Dim,T>::map_xi   (elem, p);
            const Point deta  = FE<Dim,T>::map_eta  (elem, p);
            const Point dzeta = FE<Dim,T>::map_zeta (elem, p);

            //  Newton's method in this case looks like
            //
            //  {X} = {X_n} + [J]*{dp}
            //
            //  Where {X}, {X_n} are 3x1 vectors, [J] is a 3x3 matrix
            //  d(x,y,z)/d(xi,eta,zeta), and we seek {dp}, a 3x1 vector.
            //  Since the above system is nonsingular for invertable maps
            //  we will solve
            //
            //  {dp} = [J]^-1 ({X} - {X_n})
            //
            //  which involves the inversion of the 3x3 matrix [J]
            const Real
              J11 = dxi(0), J12 = deta(0), J13 = dzeta(0),
              J21 = dxi(1), J22 = deta(1), J23 = dzeta(1),
              J31 = dxi(2), J32 = deta(2), J33 = dzeta(2);

            const Real det = (J11*(J22*J33 - J23*J32) +
                              J12*(J23*J31 - J21*J33) +
                              J13*(J21*J32 - J22*J31));

            if (secure)
              libmesh_assert_not_equal_to (det, 0.);

            const Real inv_det = 1./det;

            const Real
              Jinv11 =  (J22*J33 - J23*J32)*inv_det,
              Jinv12 = -(J12*J33 - J13*J32)*inv_det,
              Jinv13 =  (J12*J23 - J13*J22)*inv_det,

              Jinv21 = -(J21*J33 - J23*J31)*inv_det,
              Jinv22 =  (J11*J33 - J13*J31)*inv_det,
              Jinv23 = -(J11*J23 - J13*J21)*inv_det,

              Jinv31 =  (J21*J32 - J22*J31)*inv_det,
              Jinv32 = -(J11*J32 - J12*J31)*inv_det,
              Jinv33 =  (J11*J22 - J12*J21)*inv_det;

            dp(0) = (Jinv11*delta(0) +
                     Jinv12*delta(1) +
                     Jinv13*delta(2));

            dp(1) = (Jinv21*delta(0) +
                     Jinv22*delta(1) +
                     Jinv23*delta(2));

            dp(2) = (Jinv31*delta(0) +
                     Jinv32*delta(1) +
                     Jinv33*delta(2));

            // No master elements have radius > 4, but sometimes we
            // can take a step that big while still converging
            // if (secure)
              // libmesh_assert_less (dp.size(), max_step_length);

            break;
          }


          //  Some other dimension?
        default:
          libmesh_error();
        } // end switch(Dim), dp now computed



      //  ||P_n+1 - P_n||
      inverse_map_error = dp.size();

      //  P_n+1 = P_n + dp
      p.add (dp);

      //  Increment the iteration count.
      cnt++;

      //  Watch for divergence of Newton's
      //  method.  Here's how it goes:
      //  (1) For good elements, we expect convergence in 10
      //      iterations, with no too-large steps.
      //      - If called with (secure == true) and we have not yet converged
      //        print out a warning message.
      //      - If called with (secure == true) and we have not converged in
      //        20 iterations abort
      //  (2) This method may be called in cases when the target point is not
      //      inside the element and we have no business expecting convergence.
      //      For these cases if we have not converged in 10 iterations forget
      //      about it.
      if (cnt > max_cnt)
        {
          //  Warn about divergence when secure is true - this
          //  shouldn't happen
          if (secure)
            {
              // Print every time in devel/dbg modes
#ifndef NDEBUG
              libmesh_here();
              libMesh::err << "WARNING: Newton scheme has not converged in "
                            << cnt << " iterations:" << std::endl
                            << "   physical_point="
                            << physical_point
                            << "   physical_guess="
                            << physical_guess
                            << "   dp="
                            << dp
                            << "   p="
                            << p
                            << "   error=" << inverse_map_error
                            << "   in element " << elem->id()
                            << std::endl;

              elem->print_info(libMesh::err);
#else
              // In optimized mode, just print once that an inverse_map() call
              // had trouble converging its Newton iteration.
              libmesh_do_once(libMesh::err << "WARNING: At least one element took more than "
                              << max_cnt
                              << " iterations to converge in inverse_map()...\n"
                              << "Rerun in devel/dbg mode for more details."
                              << std::endl;);

#endif // NDEBUG

              if (cnt > 2*max_cnt)
                {
                  libMesh::err << "ERROR: Newton scheme FAILED to converge in "
                                << cnt
                                << " iterations!"
                                << " in element " << elem->id()
                                << std::endl;

                  elem->print_info(libMesh::err);

                  libmesh_error();
                }
            }
          //  Return a far off point when secure is false - this
          //  should only happen when we're trying to map a point
          //  that's outside the element
          else
            {
              for (unsigned int i=0; i != Dim; ++i)
                p(i) = 1e6;

              STOP_LOG("inverse_map()", "FE");
              return p;
            }
        }
    }
  while (inverse_map_error > tolerance);



  //  If we are in debug mode do two sanity checks.
#ifdef DEBUG

  if (secure)
    {
      // Make sure the point \p p on the reference element actually
      // does map to the point \p physical_point within a tolerance.

      const Point check = FE<Dim,T>::map (elem, p);
      const Point diff  = physical_point - check;

      if (diff.size() > tolerance)
        {
          libmesh_here();
          libMesh::err << "WARNING:  diff is "
                        << diff.size()
                        << std::endl
                        << " point="
                        << physical_point;
          libMesh::err << " local=" << check;
          libMesh::err << " lref= " << p;

          elem->print_info(libMesh::err);
        }

      // Make sure the point \p p on the reference element actually
      // is

      if (!FEAbstract::on_reference_element(p, elem->type(), 2*tolerance))
        {
          libmesh_here();
          libMesh::err << "WARNING:  inverse_map of physical point "
                       << physical_point
                       << "is not on element." << '\n';
          elem->print_info(libMesh::err);
        }
    }

#endif



  //  Stop logging the map inversion.
  STOP_LOG("inverse_map()", "FE");

  return p;
}

template<>

bool libMesh::FE< 3, XYZ >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 933 of file fe_xyz.C.

{ return true; }

template<>

bool libMesh::FE< 2, XYZ >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 932 of file fe_xyz.C.

{ return true; }

template<>

bool libMesh::FE< 1, XYZ >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 931 of file fe_xyz.C.

{ return true; }

template<>

bool libMesh::FE< 0, XYZ >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 930 of file fe_xyz.C.

{ return true; }

template<>

bool libMesh::FE< 3, SZABAB >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 1267 of file fe_szabab.C.

{ return false; }

template<>

bool libMesh::FE< 2, SZABAB >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 1266 of file fe_szabab.C.

{ return false; }

template<>

bool libMesh::FE< 1, SZABAB >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 1265 of file fe_szabab.C.

{ return false; }

template<>

bool libMesh::FE< 0, SZABAB >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 1264 of file fe_szabab.C.

{ return false; }

template<>

bool libMesh::FE< 3, SCALAR >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 116 of file fe_scalar.C.

{ return false; }

template<>

bool libMesh::FE< 2, SCALAR >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 115 of file fe_scalar.C.

{ return false; }

template<>

bool libMesh::FE< 1, SCALAR >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 114 of file fe_scalar.C.

{ return false; }

template<>

bool libMesh::FE< 0, SCALAR >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 113 of file fe_scalar.C.

{ return false; }

template<>

bool libMesh::FE< 3, NEDELEC_ONE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 572 of file fe_nedelec_one.C.

{ return false; }

template<>

bool libMesh::FE< 2, NEDELEC_ONE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 571 of file fe_nedelec_one.C.

{ return false; }

template<>

bool libMesh::FE< 1, NEDELEC_ONE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 570 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

bool libMesh::FE< 0, NEDELEC_ONE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 569 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

bool libMesh::FE< 3, MONOMIAL >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 410 of file fe_monomial.C.

{ return true; }

template<>

bool libMesh::FE< 2, MONOMIAL >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 409 of file fe_monomial.C.

{ return true; }

template<>

bool libMesh::FE< 1, MONOMIAL >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 408 of file fe_monomial.C.

{ return true; }

template<>

bool libMesh::FE< 0, MONOMIAL >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 407 of file fe_monomial.C.

{ return true; }

template<>

bool libMesh::FE< 3, LAGRANGE_VEC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 935 of file fe_lagrange_vec.C.

{ return false; }

template<>

bool libMesh::FE< 2, LAGRANGE_VEC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 934 of file fe_lagrange_vec.C.

{ return false; }

template<>

bool libMesh::FE< 1, LAGRANGE_VEC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 933 of file fe_lagrange_vec.C.

{ return false; }

template<>

bool libMesh::FE< 0, LAGRANGE_VEC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 932 of file fe_lagrange_vec.C.

{ return false; }

template<>

bool libMesh::FE< 3, LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 850 of file fe_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 2, LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 849 of file fe_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 1, LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 848 of file fe_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 0, LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 847 of file fe_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 3, L2_LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 491 of file fe_l2_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 2, L2_LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 490 of file fe_l2_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 1, L2_LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 489 of file fe_l2_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 0, L2_LAGRANGE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 488 of file fe_l2_lagrange.C.

{ return false; }

template<>

bool libMesh::FE< 3, L2_HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 197 of file fe_l2_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 2, L2_HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 196 of file fe_l2_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 1, L2_HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 195 of file fe_l2_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 0, L2_HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 194 of file fe_l2_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 3, HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 368 of file fe_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 2, HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 367 of file fe_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 1, HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 366 of file fe_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 0, HIERARCHIC >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 365 of file fe_hierarchic.C.

{ return true; }

template<>

bool libMesh::FE< 3, HERMITE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 355 of file fe_hermite.C.

{ return true; }

template<>

bool libMesh::FE< 2, HERMITE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 354 of file fe_hermite.C.

{ return true; }

template<>

bool libMesh::FE< 1, HERMITE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 353 of file fe_hermite.C.

{ return true; }

template<>

bool libMesh::FE< 0, HERMITE >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 352 of file fe_hermite.C.

{ return true; }

template<>

bool libMesh::FE< 3, CLOUGH >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 347 of file fe_clough.C.

{ return false; } // FIXME - this will be changed

template<>

bool libMesh::FE< 2, CLOUGH >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 346 of file fe_clough.C.

{ return false; } // FIXME - this will be changed

template<>

bool libMesh::FE< 1, CLOUGH >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 345 of file fe_clough.C.

{ return false; } // FIXME - this will be changed

template<>

bool libMesh::FE< 0, CLOUGH >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 344 of file fe_clough.C.

{ return false; } // FIXME - this will be changed

template<>

bool libMesh::FE< 3, BERNSTEIN >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 432 of file fe_bernstein.C.

{ return false; }

template<>

bool libMesh::FE< 2, BERNSTEIN >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 431 of file fe_bernstein.C.

{ return false; }

template<>

bool libMesh::FE< 1, BERNSTEIN >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 430 of file fe_bernstein.C.

{ return false; }

template<>

bool libMesh::FE< 0, BERNSTEIN >::is_hierarchic

(

)

const [inline, virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

Definition at line 429 of file fe_bernstein.C.

{ return false; }

template<unsigned int Dim, FEFamily T>

virtual bool libMesh::FE< Dim, T >::is_hierarchic

(

)

const [virtual]

Returns:

true if the finite element's higher order shape functions are hierarchic

Implements libMesh::FEAbstract.

template<unsigned int Dim, FEFamily T>

Point libMesh::FE< Dim, T >::map	(	const Elem *	elem,
		const Point &	reference_point
	)			[inline, static]

Returns:

the location (in physical space) of the point p located on the reference element.

Definition at line 1231 of file fe_map.C.

References libMesh::TypeVector< T >::add_scaled(), libMesh::Elem::default_order(), libMesh::FE< Dim, T >::n_shape_functions(), libMesh::Elem::point(), libMesh::FE< Dim, T >::shape(), and libMesh::Elem::type().

Referenced by libMesh::FE< Dim, T >::inverse_map().

{
  libmesh_assert(elem);

  Point p;

  const ElemType type     = elem->type();
  const Order order       = elem->default_order();
  const unsigned int n_sf = FE<Dim,LAGRANGE>::n_shape_functions(type, order);

  // Lagrange basis functions are used for mapping
  for (unsigned int i=0; i<n_sf; i++)
    p.add_scaled (elem->point(i),
                  FE<Dim,LAGRANGE>::shape(type,
                                          order,
                                          i,
                                          reference_point)
                  );

  return p;
}

template<unsigned int Dim, FEFamily T>

Point libMesh::FE< Dim, T >::map_eta	(	const Elem *	elem,
		const Point &	reference_point
	)			[inline, static]

Returns:

d(xyz)/deta (in physical space) of the point p located on the reference element.

Definition at line 1283 of file fe_map.C.

References libMesh::TypeVector< T >::add_scaled(), libMesh::Elem::default_order(), libMesh::FE< Dim, T >::n_shape_functions(), libMesh::Elem::point(), libMesh::FE< Dim, T >::shape_deriv(), and libMesh::Elem::type().

Referenced by libMesh::FE< Dim, T >::inverse_map().

{
  libmesh_assert(elem);

  Point p;

  const ElemType type     = elem->type();
  const Order order       = elem->default_order();
  const unsigned int n_sf = FE<Dim,LAGRANGE>::n_shape_functions(type, order);

  // Lagrange basis functions are used for mapping
  for (unsigned int i=0; i<n_sf; i++)
    p.add_scaled (elem->point(i),
                  FE<Dim,LAGRANGE>::shape_deriv(type,
                                                order,
                                                i,
                                                1,
                                                reference_point)
                  );

  return p;
}

template<unsigned int Dim, FEFamily T>

Point libMesh::FE< Dim, T >::map_xi	(	const Elem *	elem,
		const Point &	reference_point
	)			[inline, static]

Returns:

d(xyz)/dxi (in physical space) of the point p located on the reference element.

Definition at line 1257 of file fe_map.C.

References libMesh::TypeVector< T >::add_scaled(), libMesh::Elem::default_order(), libMesh::FE< Dim, T >::n_shape_functions(), libMesh::Elem::point(), libMesh::FE< Dim, T >::shape_deriv(), and libMesh::Elem::type().

Referenced by libMesh::FE< Dim, T >::inverse_map().

{
  libmesh_assert(elem);

  Point p;

  const ElemType type     = elem->type();
  const Order order       = elem->default_order();
  const unsigned int n_sf = FE<Dim,LAGRANGE>::n_shape_functions(type, order);

  // Lagrange basis functions are used for mapping
  for (unsigned int i=0; i<n_sf; i++)
    p.add_scaled (elem->point(i),
                  FE<Dim,LAGRANGE>::shape_deriv(type,
                                                order,
                                                i,
                                                0,
                                                reference_point)
                  );

  return p;
}

template<unsigned int Dim, FEFamily T>

Point libMesh::FE< Dim, T >::map_zeta	(	const Elem *	elem,
		const Point &	reference_point
	)			[inline, static]

Returns:

d(xyz)/dzeta (in physical space) of the point p located on the reference element.

Definition at line 1309 of file fe_map.C.

References libMesh::TypeVector< T >::add_scaled(), libMesh::Elem::default_order(), libMesh::FE< Dim, T >::n_shape_functions(), libMesh::Elem::point(), libMesh::FE< Dim, T >::shape_deriv(), and libMesh::Elem::type().

Referenced by libMesh::FE< Dim, T >::inverse_map().

{
  libmesh_assert(elem);

  Point p;

  const ElemType type     = elem->type();
  const Order order       = elem->default_order();
  const unsigned int n_sf = FE<Dim,LAGRANGE>::n_shape_functions(type, order);

  // Lagrange basis functions are used for mapping
  for (unsigned int i=0; i<n_sf; i++)
    p.add_scaled (elem->point(i),
                  FE<Dim,LAGRANGE>::shape_deriv(type,
                                                order,
                                                i,
                                                2,
                                                reference_point)
                  );

  return p;
}

template<>

unsigned int libMesh::FE< 3, XYZ >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 907 of file fe_xyz.C.

{ return xyz_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, XYZ >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 906 of file fe_xyz.C.

{ return xyz_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, XYZ >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 905 of file fe_xyz.C.

{ return xyz_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, XYZ >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 904 of file fe_xyz.C.

{ return xyz_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, SZABAB >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 1243 of file fe_szabab.C.

{ return szabab_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, SZABAB >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 1242 of file fe_szabab.C.

{ return szabab_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, SZABAB >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 1241 of file fe_szabab.C.

{ return szabab_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, SZABAB >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 1240 of file fe_szabab.C.

{ return szabab_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, SCALAR >::n_dofs	(	const	ElemType,
		const Order	o
	)			[inline]

Definition at line 90 of file fe_scalar.C.

{ return o; }

template<>

unsigned int libMesh::FE< 2, SCALAR >::n_dofs	(	const	ElemType,
		const Order	o
	)			[inline]

Definition at line 89 of file fe_scalar.C.

{ return o; }

template<>

unsigned int libMesh::FE< 1, SCALAR >::n_dofs	(	const	ElemType,
		const Order	o
	)			[inline]

Definition at line 88 of file fe_scalar.C.

{ return o; }

template<>

unsigned int libMesh::FE< 0, SCALAR >::n_dofs	(	const	ElemType,
		const Order	o
	)			[inline]

Definition at line 87 of file fe_scalar.C.

{ return o; }

template<>

unsigned int libMesh::FE< 3, NEDELEC_ONE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 544 of file fe_nedelec_one.C.

{ return nedelec_one_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, NEDELEC_ONE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 543 of file fe_nedelec_one.C.

{ return nedelec_one_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, NEDELEC_ONE >::n_dofs	(	const	ElemType,
		const	Order
	)			[inline]

Definition at line 542 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

unsigned int libMesh::FE< 0, NEDELEC_ONE >::n_dofs	(	const	ElemType,
		const	Order
	)			[inline]

Definition at line 541 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

unsigned int libMesh::FE< 3, MONOMIAL >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 383 of file fe_monomial.C.

{ return monomial_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, MONOMIAL >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 382 of file fe_monomial.C.

{ return monomial_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, MONOMIAL >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 381 of file fe_monomial.C.

{ return monomial_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, MONOMIAL >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 380 of file fe_monomial.C.

{ return monomial_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, LAGRANGE_VEC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 907 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs().

{ return 3*FE<3,LAGRANGE>::n_dofs(t,o); }

template<>

unsigned int libMesh::FE< 2, LAGRANGE_VEC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 906 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs().

{ return 2*FE<2,LAGRANGE>::n_dofs(t,o); }

template<>

unsigned int libMesh::FE< 1, LAGRANGE_VEC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 905 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs().

{ return FE<1,LAGRANGE>::n_dofs(t,o); }

template<>

unsigned int libMesh::FE< 0, LAGRANGE_VEC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 904 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs().

{ return FE<0,LAGRANGE>::n_dofs(t,o); }

template<>

unsigned int libMesh::FE< 3, LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 822 of file fe_lagrange.C.

{ return lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 821 of file fe_lagrange.C.

{ return lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 820 of file fe_lagrange.C.

{ return lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 819 of file fe_lagrange.C.

{ return lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, L2_LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 465 of file fe_l2_lagrange.C.

{ return l2_lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, L2_LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 464 of file fe_l2_lagrange.C.

{ return l2_lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, L2_LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 463 of file fe_l2_lagrange.C.

{ return l2_lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, L2_LAGRANGE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 462 of file fe_l2_lagrange.C.

{ return l2_lagrange_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, L2_HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 172 of file fe_l2_hierarchic.C.

{ return l2_hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, L2_HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 171 of file fe_l2_hierarchic.C.

{ return l2_hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, L2_HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 170 of file fe_l2_hierarchic.C.

{ return l2_hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, L2_HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 169 of file fe_l2_hierarchic.C.

{ return l2_hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 344 of file fe_hierarchic.C.

{ return hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 343 of file fe_hierarchic.C.

{ return hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 342 of file fe_hierarchic.C.

{ return hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, HIERARCHIC >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 341 of file fe_hierarchic.C.

{ return hierarchic_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, HERMITE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 330 of file fe_hermite.C.

{ return hermite_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, HERMITE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 329 of file fe_hermite.C.

{ return hermite_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, HERMITE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 328 of file fe_hermite.C.

{ return hermite_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, HERMITE >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 327 of file fe_hermite.C.

{ return hermite_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, CLOUGH >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 322 of file fe_clough.C.

{ return clough_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, CLOUGH >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 321 of file fe_clough.C.

{ return clough_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, CLOUGH >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 320 of file fe_clough.C.

{ return clough_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, CLOUGH >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 319 of file fe_clough.C.

{ return clough_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 3, BERNSTEIN >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 408 of file fe_bernstein.C.

{ return bernstein_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 2, BERNSTEIN >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 407 of file fe_bernstein.C.

{ return bernstein_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 1, BERNSTEIN >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 406 of file fe_bernstein.C.

{ return bernstein_n_dofs(t, o); }

template<>

unsigned int libMesh::FE< 0, BERNSTEIN >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[inline]

Definition at line 405 of file fe_bernstein.C.

{ return bernstein_n_dofs(t, o); }

template<unsigned int Dim, FEFamily T>

static unsigned int libMesh::FE< Dim, T >::n_dofs	(	const ElemType	t,
		const Order	o
	)			[static]

Returns:

the number of shape functions associated with this finite element.

On a p-refined element, o should be the total order of the element.

Referenced by libMesh::FE< Dim, T >::dofs_on_edge(), libMesh::FE< Dim, T >::dofs_on_side(), libMesh::FE< Dim, T >::n_dofs(), libMesh::FE< Dim, LAGRANGE_VEC >::n_shape_functions(), and libMesh::FE< Dim, T >::n_shape_functions().

template<>

unsigned int libMesh::FE< 3, XYZ >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 914 of file fe_xyz.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 2, XYZ >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 913 of file fe_xyz.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 1, XYZ >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 912 of file fe_xyz.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 0, XYZ >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 911 of file fe_xyz.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 3, SZABAB >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 1249 of file fe_szabab.C.

{ return szabab_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, SZABAB >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 1248 of file fe_szabab.C.

{ return szabab_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, SZABAB >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 1247 of file fe_szabab.C.

{ return szabab_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 0, SZABAB >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 1246 of file fe_szabab.C.

{ return szabab_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 3, SCALAR >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 97 of file fe_scalar.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 2, SCALAR >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 96 of file fe_scalar.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 1, SCALAR >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 95 of file fe_scalar.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 0, SCALAR >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 94 of file fe_scalar.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 3, NEDELEC_ONE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 552 of file fe_nedelec_one.C.

{ return nedelec_one_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, NEDELEC_ONE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 551 of file fe_nedelec_one.C.

{ return nedelec_one_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, NEDELEC_ONE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 550 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

unsigned int libMesh::FE< 0, NEDELEC_ONE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 549 of file fe_nedelec_one.C.

{ NEDELEC_LOW_D_ERROR_MESSAGE }

template<>

unsigned int libMesh::FE< 3, MONOMIAL >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 390 of file fe_monomial.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 2, MONOMIAL >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 389 of file fe_monomial.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 1, MONOMIAL >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 388 of file fe_monomial.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 0, MONOMIAL >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 387 of file fe_monomial.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 3, LAGRANGE_VEC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 915 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs_at_node().

{ return 3*FE<2,LAGRANGE>::n_dofs_at_node(t,o,n); }

template<>

unsigned int libMesh::FE< 2, LAGRANGE_VEC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 914 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs_at_node().

{ return 2*FE<2,LAGRANGE>::n_dofs_at_node(t,o,n); }

template<>

unsigned int libMesh::FE< 1, LAGRANGE_VEC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 913 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs_at_node().

{ return FE<1,LAGRANGE>::n_dofs_at_node(t,o,n); }

template<>

unsigned int libMesh::FE< 0, LAGRANGE_VEC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 912 of file fe_lagrange_vec.C.

References libMesh::FE< Dim, T >::n_dofs_at_node().

{ return FE<0,LAGRANGE>::n_dofs_at_node(t,o,n); }

template<>

unsigned int libMesh::FE< 3, LAGRANGE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 830 of file fe_lagrange.C.

{ return lagrange_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, LAGRANGE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 829 of file fe_lagrange.C.

{ return lagrange_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, LAGRANGE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 828 of file fe_lagrange.C.

{ return lagrange_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 0, LAGRANGE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 827 of file fe_lagrange.C.

{ return lagrange_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 3, L2_LAGRANGE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 472 of file fe_l2_lagrange.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 2, L2_LAGRANGE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 471 of file fe_l2_lagrange.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 1, L2_LAGRANGE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 470 of file fe_l2_lagrange.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 0, L2_LAGRANGE >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 469 of file fe_l2_lagrange.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 3, L2_HIERARCHIC >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 179 of file fe_l2_hierarchic.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 2, L2_HIERARCHIC >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 178 of file fe_l2_hierarchic.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 1, L2_HIERARCHIC >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 177 of file fe_l2_hierarchic.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 0, L2_HIERARCHIC >::n_dofs_at_node	(	const	ElemType,
		const	Order,
		const unsigned	int
	)			[inline]

Definition at line 176 of file fe_l2_hierarchic.C.

{ return 0; }

template<>

unsigned int libMesh::FE< 3, HIERARCHIC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 350 of file fe_hierarchic.C.

{ return hierarchic_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, HIERARCHIC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 349 of file fe_hierarchic.C.

{ return hierarchic_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, HIERARCHIC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 348 of file fe_hierarchic.C.

{ return hierarchic_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 0, HIERARCHIC >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 347 of file fe_hierarchic.C.

{ return hierarchic_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 3, HERMITE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 337 of file fe_hermite.C.

{ return hermite_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, HERMITE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 336 of file fe_hermite.C.

{ return hermite_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, HERMITE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 335 of file fe_hermite.C.

{ return hermite_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 0, HERMITE >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 334 of file fe_hermite.C.

{ return hermite_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 3, CLOUGH >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 329 of file fe_clough.C.

{ return clough_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, CLOUGH >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 328 of file fe_clough.C.

{ return clough_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, CLOUGH >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 327 of file fe_clough.C.

{ return clough_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 0, CLOUGH >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 326 of file fe_clough.C.

{ return clough_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 3, BERNSTEIN >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 414 of file fe_bernstein.C.

{ return bernstein_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 2, BERNSTEIN >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)			[inline]

Definition at line 413 of file fe_bernstein.C.

{ return bernstein_n_dofs_at_node(t, o, n); }

template<>

unsigned int libMesh::FE< 1, BERNSTEIN >::n_dofs_at_node	(	const ElemType	t,
		const Order	o,
		const unsigned int	n
	)