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)
	MaxQpsThread (MaxQpsThread &x, Threads::split split)
void	operator() (const libMesh::ConstElemRange &range)
void	join (const MaxQpsThread &y)
unsigned int	max () const

Protected Attributes
FEProblemBase &	_fe_problem
THREAD_ID	_tid
unsigned int	_max
	Maximum number of qps 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 27 of file MaxQpsThread.h.

Constructor & Destructor Documentation

MaxQpsThread() [1/2]

LIBMESH_DEFINE_HASH_POINTERS MaxQpsThread::MaxQpsThread

(

FEProblemBase &

fe_problem

)

Definition at line 20 of file MaxQpsThread.C.

: _fe_problem(fe_problem), _max(0) {}

MaxQpsThread() [2/2]

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

Definition at line 23 of file MaxQpsThread.C.

  : _fe_problem(x._fe_problem), _max(x._max)
{
}

Member Function Documentation

join()

void MaxQpsThread::join

(

const MaxQpsThread &

y

)

Definition at line 82 of file MaxQpsThread.C.

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

max()

unsigned int MaxQpsThread::max ( ) const

inline

Definition at line 39 of file MaxQpsThread.h.

Referenced by FEProblemBase::updateMaxQps().

{ return _max; }

operator()()

void MaxQpsThread::operator()

(

const libMesh::ConstElemRange &

range

)

Definition at line 29 of file MaxQpsThread.C.

{
  ParallelUniqueId puid;
  _tid = puid.id;

  // Not actually using any pre-existing data so it shouldn't matter which assembly we use
  auto & assem = _fe_problem.assembly(_tid, 0);

  // For short circuiting reinit.  With potential block-specific qrules we
  // need to track "seen" element types by their subdomains as well.
  std::set<std::pair<ElemType, SubdomainID>> seen_it;

  for (const auto & elem : range)
  {
    // Only reinit if the element type has not previously been seen
    if (!seen_it.insert(std::make_pair(elem->type(), elem->subdomain_id())).second)
      continue;

    // This ensures we can access the correct qrules if any block-specific
    // qrules have been created.
    assem.setCurrentSubdomainID(elem->subdomain_id());
    assem.setVolumeQRule(elem);

    auto & qrule = assem.writeableQRule();
    qrule->init(*elem);
    if (qrule->n_points() > _max)
      _max = qrule->n_points();

    // We used to check side quadrature rules too - badly (assuming
    // side 0 doesn't have a smaller rule than the others, which is
    // often untrue on prisms) and generally unnecessarily (a side
    // rule will always use fewer points than a corresponding interior
    // rule).
    //
    // We should handle the possibility of users manually specifying a
    // higher order for side than for interior integration.  Doing
    // this efficiently will need to wait for a new libMesh API,
    // though.

    // 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();
  }

  // Clear the cached quadrature rules because we may add FE objects in between now and simulation
  // start and we only ensure we set all the FE quadrature rules if a quadrature rule is different
  // from the cached quadrature rule
  assem.clearCachedQRules();
}

Member Data Documentation

_fe_problem

FEProblemBase& MaxQpsThread::_fe_problem

protected

Definition at line 42 of file MaxQpsThread.h.

Referenced by operator()().

_max

unsigned int MaxQpsThread::_max

protected

Maximum number of qps encountered.

Definition at line 47 of file MaxQpsThread.h.

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

_tid

THREAD_ID MaxQpsThread::_tid

protected

Definition at line 44 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