NumLib::IntegrationGaussLegendreRegular< N_DIM > Class Template Reference

Detailed Description

template<unsigned N_DIM>
class NumLib::IntegrationGaussLegendreRegular< N_DIM >

Gauss-Legendre quadrature rule for regular shape elements: line, quad and hex.

Template Parameters

N_DIM

Spatial dimension

Definition at line 28 of file IntegrationGaussLegendreRegular.h.

#include <IntegrationGaussLegendreRegular.h>

Public Types
using	WeightedPoint = typename MathLib::TemplateWeightedPoint< double, double, N_DIM >

Public Member Functions
	IntegrationGaussLegendreRegular (unsigned order=2)
void	setIntegrationOrder (unsigned order)
	Change the integration order. More...
unsigned	getIntegrationOrder () const
	return current integration order. More...
unsigned	getNumberOfPoints () const
	return the number of sampling points More...
WeightedPoint	getWeightedPoint (unsigned igp) const
std::array< unsigned, 1 >	getPositionIndices (unsigned, unsigned igp)
std::array< unsigned, 2 >	getPositionIndices (unsigned order, unsigned igp)
std::array< unsigned, 3 >	getPositionIndices (unsigned order, unsigned igp)
template<typename Method >
MathLib::TemplateWeightedPoint< double, double, N_DIM >	getWeightedPoint (std::array< unsigned, N_DIM > const &pos)

Static Public Member Functions
static std::array< unsigned, N_DIM >	getPositionIndices (unsigned order, unsigned igp)
static WeightedPoint	getWeightedPoint (unsigned order, unsigned igp)

Static Private Member Functions
template<typename Method >
static WeightedPoint	getWeightedPoint (std::array< unsigned, N_DIM > const &pos)

Private Attributes
unsigned	_order
unsigned	_n_sampl_pt {0}

Member Typedef Documentation

WeightedPoint

template<unsigned N_DIM>

using NumLib::IntegrationGaussLegendreRegular< N_DIM >::WeightedPoint = typename MathLib::TemplateWeightedPoint<double, double, N_DIM>

Definition at line 31 of file IntegrationGaussLegendreRegular.h.

Constructor & Destructor Documentation

IntegrationGaussLegendreRegular()

template<unsigned N_DIM>

NumLib::IntegrationGaussLegendreRegular< N_DIM >::IntegrationGaussLegendreRegular ( unsigned  order = 2 )

inlineexplicit

Create IntegrationGaussLegendreRegular of the given Gauss-Legendre integration order.

Parameters

order

integration order (default 2)

Definition at line 38 of file IntegrationGaussLegendreRegular.h.

                                                                 : _order(order)
    {
        this->setIntegrationOrder(order);
    }

References NumLib::IntegrationGaussLegendreRegular< N_DIM >::setIntegrationOrder().

Member Function Documentation

getIntegrationOrder()

template<unsigned N_DIM>

unsigned NumLib::IntegrationGaussLegendreRegular< N_DIM >::getIntegrationOrder ( ) const

inline

return current integration order.

Definition at line 51 of file IntegrationGaussLegendreRegular.h.

{ return _order; }

References NumLib::IntegrationGaussLegendreRegular< N_DIM >::_order.

Referenced by NumLib::IntegrationGaussLegendreRegular< N_DIM >::getWeightedPoint().

getNumberOfPoints()

template<unsigned N_DIM>

unsigned NumLib::IntegrationGaussLegendreRegular< N_DIM >::getNumberOfPoints ( ) const

inline

return the number of sampling points

Definition at line 54 of file IntegrationGaussLegendreRegular.h.

{ return _n_sampl_pt; }

References NumLib::IntegrationGaussLegendreRegular< N_DIM >::_n_sampl_pt.

getPositionIndices() [1/4]

std::array< unsigned, 2 > NumLib::IntegrationGaussLegendreRegular< 2 >::getPositionIndices ( unsigned  order, unsigned  igp  )

inline

Definition at line 29 of file IntegrationGaussLegendreRegular-impl.h.

{
    assert(igp < order * order);
    std::array<unsigned, 2> result;
    result[0] = igp / order;
    result[1] = igp % order;
    return result;
}

getPositionIndices() [2/4]

std::array< unsigned, 3 > NumLib::IntegrationGaussLegendreRegular< 3 >::getPositionIndices ( unsigned  order, unsigned  igp  )

inline

Definition at line 41 of file IntegrationGaussLegendreRegular-impl.h.

{
    assert(igp < order * order * order);
    unsigned const gp_r = igp / (order * order);
    unsigned const gp_s = igp % (order * order);
    std::array<unsigned, 3> result;
    result[0] = gp_r;
    result[1] = gp_s / order;
    result[2] = gp_s % order;
    return result;
}

getPositionIndices() [3/4]

template<unsigned N_DIM>

static std::array<unsigned, N_DIM> NumLib::IntegrationGaussLegendreRegular< N_DIM >::getPositionIndices ( unsigned  order, unsigned  igp  )

static

Get position indexes in r-s-t axis.

Parameters

order	The number of integration points
igp	The integration point index

Returns

a tuple of position indexes

getPositionIndices() [4/4]

std::array< unsigned, 1 > NumLib::IntegrationGaussLegendreRegular< 1 >::getPositionIndices ( unsigned  , unsigned  igp  )

inline

Definition at line 19 of file IntegrationGaussLegendreRegular-impl.h.

{
    std::array<unsigned, 1> result;
    result[0] = igp;
    return result;
}

getWeightedPoint() [1/4]

template<unsigned N_DIM>

template<typename Method >

MathLib::TemplateWeightedPoint<double, double, N_DIM> NumLib::IntegrationGaussLegendreRegular< N_DIM >::getWeightedPoint ( std::array< unsigned, N_DIM > const &  pos )

inline

Definition at line 81 of file IntegrationGaussLegendreRegular-impl.h.

{
    std::array<double, N_DIM> coords;
    double weight = 1;
    for (unsigned d = 0; d < N_DIM; d++)
    {
        coords[d] = Method::X[pos[d]];
        weight *= Method::W[pos[d]];
    }
 
    return MathLib::TemplateWeightedPoint<double, double, N_DIM>(coords,
                                                                 weight);
}

getWeightedPoint() [2/4]

template<unsigned N_DIM>

template<typename Method >

static WeightedPoint NumLib::IntegrationGaussLegendreRegular< N_DIM >::getWeightedPoint ( std::array< unsigned, N_DIM > const &  pos )

staticprivate

Computes weighted point using given integration method.

Template Parameters

Method

Integration method to use.

Parameters

pos	Point indices computed by getPositionIndices.

getWeightedPoint() [3/4]

template<unsigned N_DIM>

WeightedPoint NumLib::IntegrationGaussLegendreRegular< N_DIM >::getWeightedPoint ( unsigned  igp ) const

inline

Get coordinates of the integration point.

Parameters

igp	The integration point index

Returns

a weighted point

Definition at line 60 of file IntegrationGaussLegendreRegular.h.

    {
        return getWeightedPoint(getIntegrationOrder(), igp);
    }

References NumLib::IntegrationGaussLegendreRegular< N_DIM >::getIntegrationOrder().

getWeightedPoint() [4/4]

template<unsigned N_DIM>

MathLib::TemplateWeightedPoint< double, double, N_DIM > NumLib::IntegrationGaussLegendreRegular< N_DIM >::getWeightedPoint ( unsigned  order, unsigned  igp  )

inlinestatic

Get coordinates of the integration point.

Parameters

order	The number of integration points
igp	The integration point index

Returns

a weighted point

Definition at line 56 of file IntegrationGaussLegendreRegular-impl.h.

{
    assert(igp < std::pow(order, N_DIM));
    std::array<unsigned, N_DIM> const pos = getPositionIndices(order, igp);
 
    switch (order)
    {
        case 1:
            return getWeightedPoint<MathLib::GaussLegendre<1>>(pos);
        case 2:
            return getWeightedPoint<MathLib::GaussLegendre<2>>(pos);
        case 3:
            return getWeightedPoint<MathLib::GaussLegendre<3>>(pos);
        case 4:
            return getWeightedPoint<MathLib::GaussLegendre<4>>(pos);
    }
 
    return MathLib::TemplateWeightedPoint<double, double, N_DIM>(
        std::array<double, N_DIM>(), 0);
}

setIntegrationOrder()

template<unsigned N_DIM>

void NumLib::IntegrationGaussLegendreRegular< N_DIM >::setIntegrationOrder ( unsigned  order )

inline

Change the integration order.

Definition at line 44 of file IntegrationGaussLegendreRegular.h.

    {
        this->_n_sampl_pt = static_cast<unsigned>(std::pow(order, N_DIM));
        this->_order = order;
    }

References NumLib::IntegrationGaussLegendreRegular< N_DIM >::_n_sampl_pt, and NumLib::IntegrationGaussLegendreRegular< N_DIM >::_order.

Referenced by NumLib::IntegrationGaussLegendreRegular< N_DIM >::IntegrationGaussLegendreRegular().

Member Data Documentation

_n_sampl_pt

template<unsigned N_DIM>

unsigned NumLib::IntegrationGaussLegendreRegular< N_DIM >::_n_sampl_pt {0}

private

Definition at line 91 of file IntegrationGaussLegendreRegular.h.

Referenced by NumLib::IntegrationGaussLegendreRegular< N_DIM >::getNumberOfPoints(), and NumLib::IntegrationGaussLegendreRegular< N_DIM >::setIntegrationOrder().

_order

template<unsigned N_DIM>

unsigned NumLib::IntegrationGaussLegendreRegular< N_DIM >::_order

private

Definition at line 90 of file IntegrationGaussLegendreRegular.h.

Referenced by NumLib::IntegrationGaussLegendreRegular< N_DIM >::getIntegrationOrder(), and NumLib::IntegrationGaussLegendreRegular< N_DIM >::setIntegrationOrder().

---

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

• NumLib/Fem/Integration/IntegrationGaussLegendreRegular.h
• NumLib/Fem/Integration/IntegrationGaussLegendreRegular-impl.h