MaxQpsThread Class Reference

This class determines the maximum number of Quadrature Points and Shape Functions used for a given simulation based on the variable discretizations, and quadrature rules used for all variables in the system. More...

#include <MaxQpsThread.h>

Public Member Functions
	MaxQpsThread (FEProblemBase &fe_problem, QuadratureType type, Order order, Order face_order)
	MaxQpsThread (MaxQpsThread &x, Threads::split split)
void	operator() (const ConstElemRange &range)
void	join (const MaxQpsThread &y)
unsigned int	max () const
unsigned int	max_shape_funcs () const

Protected Attributes
FEProblemBase &	_fe_problem
QuadratureType	_qtype
Order	_order
Order	_face_order
THREAD_ID	_tid
unsigned int	_max
	Maximum number of qps encountered. More...
unsigned int	_max_shape_funcs
	Maximum number of shape functions encountered. More...

Detailed Description

This class determines the maximum number of Quadrature Points and Shape Functions used for a given simulation based on the variable discretizations, and quadrature rules used for all variables in the system.

Definition at line 28 of file MaxQpsThread.h.

Constructor & Destructor Documentation

MaxQpsThread() [1/2]

LIBMESH_DEFINE_HASH_POINTERS MaxQpsThread::MaxQpsThread	(	FEProblemBase &	fe_problem,
		QuadratureType	type,
		Order	order,
		Order	face_order
	)

Definition at line 19 of file MaxQpsThread.C.

  : _fe_problem(fe_problem),
    _qtype(qtype),
    _order(order),
    _face_order(face_order),
    _max(0),
    _max_shape_funcs(0)
{
}

MaxQpsThread() [2/2]

MaxQpsThread::MaxQpsThread	(	MaxQpsThread &	x,
		Threads::split	split
	)

Definition at line 33 of file MaxQpsThread.C.

  : _fe_problem(x._fe_problem),
    _qtype(x._qtype),
    _order(x._order),
    _face_order(x._face_order),
    _max(x._max),
    _max_shape_funcs(x._max_shape_funcs)
{
}

Member Function Documentation

join()

void MaxQpsThread::join

(

const MaxQpsThread &

y

)

Definition at line 97 of file MaxQpsThread.C.

{
  if (y._max > _max)
    _max = y._max;

  if (y._max_shape_funcs > _max_shape_funcs)
    _max_shape_funcs = y._max_shape_funcs;
}

max()

unsigned int MaxQpsThread::max ( ) const

inline

Definition at line 40 of file MaxQpsThread.h.

Referenced by FEProblemBase::createQRules().

{ return _max; }

max_shape_funcs()

unsigned int MaxQpsThread::max_shape_funcs ( ) const

inline

Definition at line 42 of file MaxQpsThread.h.

Referenced by FEProblemBase::createQRules().

{ return _max_shape_funcs; }

operator()()

void MaxQpsThread::operator()

(

const ConstElemRange &

range

)

Definition at line 44 of file MaxQpsThread.C.

{
  ParallelUniqueId puid;
  _tid = puid.id;

  // For short circuiting reinit
  std::set<ElemType> seen_it;
  for (const auto & elem : range)
  {
    // Only reinit if the element type has not previously been seen
    if (seen_it.insert(elem->type()).second)
    {
      FEType fe_type(FIRST, LAGRANGE);
      unsigned int dim = elem->dim();
      unsigned int side = 0; // we assume that any element will have at least one side ;)

      // We cannot mess with the FE objects in Assembly, because we might need to request second
      // derivatives
      // later on. If we used them, we'd call reinit on them, thus making the call to request second
      // derivatives harmful (i.e. leading to segfaults/asserts). Thus, we have to use a locally
      // allocated object here.
      std::unique_ptr<FEBase> fe(FEBase::build(dim, fe_type));

      // figure out the number of qps for volume
      std::unique_ptr<QBase> qrule(QBase::build(_qtype, dim, _order));
      fe->attach_quadrature_rule(qrule.get());
      fe->reinit(elem);
      if (qrule->n_points() > _max)
        _max = qrule->n_points();

      unsigned int n_shape_funcs = fe->n_shape_functions();
      if (n_shape_funcs > _max_shape_funcs)
        _max_shape_funcs = n_shape_funcs;

      // figure out the number of qps for the face
      // NOTE: user might specify higher order rule for faces, thus possibly ending up with more qps
      // than in the volume
      std::unique_ptr<QBase> qrule_face(QBase::build(_qtype, dim - 1, _face_order));
      fe->attach_quadrature_rule(qrule_face.get());
      fe->reinit(elem, side);
      if (qrule_face->n_points() > _max)
        _max = qrule_face->n_points();

      // In initial conditions nodes are enumerated as pretend quadrature points
      // using the _qp index to access coupled variables. In order to be able to
      // use _zero (resized according to _max_qps) with _qp, we need to count nodes.
      if (elem->n_nodes() > _max)
        _max = elem->n_nodes();
    }
  }
}

Member Data Documentation

_face_order

Order MaxQpsThread::_face_order

protected

Definition at line 49 of file MaxQpsThread.h.

Referenced by operator()().

_fe_problem

FEProblemBase& MaxQpsThread::_fe_problem

protected

Definition at line 45 of file MaxQpsThread.h.

_max

unsigned int MaxQpsThread::_max

protected

Maximum number of qps encountered.

Definition at line 54 of file MaxQpsThread.h.

Referenced by join(), max(), and operator()().

_max_shape_funcs

unsigned int MaxQpsThread::_max_shape_funcs

protected

Maximum number of shape functions encountered.

Definition at line 57 of file MaxQpsThread.h.

Referenced by join(), max_shape_funcs(), and operator()().

_order

Order MaxQpsThread::_order

protected

Definition at line 48 of file MaxQpsThread.h.

Referenced by operator()().

_qtype

QuadratureType MaxQpsThread::_qtype

protected

Definition at line 47 of file MaxQpsThread.h.

Referenced by operator()().

_tid

THREAD_ID MaxQpsThread::_tid

protected

Definition at line 51 of file MaxQpsThread.h.

Referenced by operator()().

---

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

• include/loops/MaxQpsThread.h
• src/loops/MaxQpsThread.C