BEM++
2.0
|
Default family of quadrature rules over pairs of elements. More...
#include </home/wojtek/Projects/BEM/bempp-sven/bempp/lib/fiber/default_double_quadrature_rule_family.hpp>
Public Member Functions | |
virtual void | fillQuadraturePointsAndWeights (const DoubleQuadratureDescriptor &desc, arma::Mat< CoordinateType > &testPoints, arma::Mat< CoordinateType > &trialPoints, std::vector< CoordinateType > &testWeights, std::vector< CoordinateType > &trialWeights, bool &isTensor) const |
Fill arrays of quadrature points and weights. More... | |
![]() | |
virtual | ~DoubleQuadratureRuleFamily () |
Destructor. | |
Default family of quadrature rules over pairs of elements.
Regular integrals are treated with tensor-product quadrature rules, singular integrals with non-tensor-product rules taken from Sauter and Schwab, "Boundary Element Methods", Springer 2011.
|
virtual |
Fill arrays of quadrature points and weights.
[in] | desc | Quadrature descriptor. |
[out] | testPoints | 2D array whose (i, j)th element will contain the ith (local) coordinate of the jth quadrature point on the test element. |
[out] | trialPoints | 2D array whose (i, j)th element will contain the ith (local) coordinate of the jth quadrature point on the trial element. |
[out] | testWeights | If isTensor is set to true , this parameter will be a vector of test point weights of a separable quadrature rule; otherwise it will be a vector of point weights of a nonseparable quadarture rule. |
[out] | trialWeights | If isTensor is set to true , this parameter will be a vector of trial point weights of a separable quadrature rule; otherwise it will be empty. |
[out] | isTensor | If set to true , the returned quadrature rule is a tensor-product rule, otherwise it is not. |
Implements Fiber::DoubleQuadratureRuleFamily< CoordinateType >.
References Fiber::DoubleQuadratureDescriptor::testOrder, Fiber::ElementPairTopology::testVertexCount, Fiber::DoubleQuadratureDescriptor::topology, Fiber::DoubleQuadratureDescriptor::trialOrder, Fiber::ElementPairTopology::trialVertexCount, and Fiber::ElementPairTopology::type.