Integrate Struct Reference

Public Member Functions
	Integrate (const System &sys, const FEMFunctionBase< Number > &f)
	Integrate (Integrate &other, Threads::split)
void	operator() (const ConstElemRange &range)
void	join (const Integrate &other)
Number	integral () const

Private Attributes
const System &	_sys
std::unique_ptr< FEMFunctionBase< Number > >	_f
Number	_integral

Detailed Description

Definition at line 104 of file calculator.C.

Constructor & Destructor Documentation

Integrate() [1/2]

Integrate::Integrate ( const System &  sys, const FEMFunctionBase< Number > &  f  )

inline

Definition at line 105 of file calculator.C.

                                                :
    _sys(sys), _f(f.clone()), _integral(0) {}

Integrate() [2/2]

Integrate::Integrate ( Integrate &  other, Threads::split    )

inline

Definition at line 109 of file calculator.C.

                                              :
    _sys(other._sys), _f(other._f->clone()), _integral(0) {}

Member Function Documentation

integral()

Number Integrate::integral ( ) const

inline

Definition at line 152 of file calculator.C.

{ return _integral; }

join()

void Integrate::join ( const Integrate &  other )

inline

Definition at line 147 of file calculator.C.

References _integral.

  {
    _integral += other._integral;
  }

operator()()

void Integrate::operator() ( const ConstElemRange &  range )

inline

Definition at line 112 of file calculator.C.

References libMesh::FEMContext::elem_fe_reinit(), libMesh::FEMContext::get_element_fe(), libMesh::FEAbstract::get_JxW(), libMesh::FEAbstract::get_nothing(), libMesh::FEAbstract::get_xyz(), libMesh::index_range(), and libMesh::FEMContext::pre_fe_reinit().

                                                 {
    FEMContext context(_sys);

    FEBase * elem_fe = nullptr;
    context.get_element_fe(0, elem_fe);

    const std::vector<Real> & JxW = elem_fe->get_JxW();
    const std::vector<Point> & xyz = elem_fe->get_xyz();

    _f->init_context(context);

    // If _f didn't need anything, that wasn't a fluke
    elem_fe->get_nothing();

    // We only integrate on the highest dimensional elements in the
    // mesh, lest we be adding apples to oranges
    const unsigned int mesh_dim = _sys.get_mesh().mesh_dimension();

    for (const auto & elem : range)
      {
        if (elem->dim() < mesh_dim)
          continue;

        context.pre_fe_reinit(_sys, elem);
        context.elem_fe_reinit();

        for (auto qp : index_range(JxW))
          {
            const Number output = (*_f)(context, xyz[qp]);

            _integral += output * JxW[qp];
          }
      }
  }

Member Data Documentation

_f

std::unique_ptr<FEMFunctionBase<Number> > Integrate::_f

private

Definition at line 157 of file calculator.C.

_integral

Number Integrate::_integral

private

Definition at line 159 of file calculator.C.

Referenced by join().

_sys

const System& Integrate::_sys

private

Definition at line 155 of file calculator.C.

---

The documentation for this struct was generated from the following file:

• src/apps/calculator.C