anonymous_namespace{IntegrationGaussLegendrePrism.cpp} Namespace Reference

Functions
template<int OrderGaussLegendreTri, int OrderGaussLegendre>
constexpr unsigned	getNumberOfPointsConcrete ()
template<int OrderGaussLegendreTri, int OrderGaussLegendre>
MathLib::WeightedPoint	getWeightedPointConcrete (unsigned const igp)

Function Documentation

getNumberOfPointsConcrete()

template<int OrderGaussLegendreTri, int OrderGaussLegendre>

unsigned anonymous_namespace{IntegrationGaussLegendrePrism.cpp}::getNumberOfPointsConcrete ( )

constexpr

Definition at line 21 of file IntegrationGaussLegendrePrism.cpp.

{
    return MathLib::GaussLegendreTri<OrderGaussLegendreTri>::NPoints *
           OrderGaussLegendre;
}

References MathLib::GaussLegendreTri< ORDER >::NPoints.

Referenced by getWeightedPointConcrete().

getWeightedPointConcrete()

template<int OrderGaussLegendreTri, int OrderGaussLegendre>

MathLib::WeightedPoint anonymous_namespace{IntegrationGaussLegendrePrism.cpp}::getWeightedPointConcrete

(

unsigned const

igp

)

Definition at line 28 of file IntegrationGaussLegendrePrism.cpp.

{
    using GL = MathLib::GaussLegendre<OrderGaussLegendre>;
    using GLT = MathLib::GaussLegendreTri<OrderGaussLegendreTri>;
 
    assert(igp < (getNumberOfPointsConcrete<OrderGaussLegendreTri,
                                            OrderGaussLegendre>()));
 
    const unsigned gp_r = igp % GLT::NPoints;
    const unsigned gp_t = igp / GLT::NPoints;
 
    std::array<double, 3> rst{GLT::X[gp_r][0], GLT::X[gp_r][1], GL::X[gp_t]};
 
    double const weight = GLT::W[gp_r] * 0.5 * GL::W[gp_t];
 
    return MathLib::WeightedPoint(rst, weight);
}

References getNumberOfPointsConcrete().