51template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
56template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
59 CubatureTensor<Scalar,ArrayPoint,ArrayWeight>(cubature1,cubature2) {}
62template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
66 CubatureTensor<Scalar,ArrayPoint,ArrayWeight>(cubature1,cubature2, cubature3) {}
69template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
74template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
76 ArrayWeight & cubWeights)
const {
79 int numCubPoints = getNumPoints();
81 for(
int i=0; i<numCubPoints; i++)
83 Scalar zTerm = 0.5*(1.0 - cubPoints(i,2));
84 cubPoints(i,0) = zTerm*cubPoints(i,0);
85 cubPoints(i,1) = zTerm*cubPoints(i,1);
86 cubPoints(i,2) = 1.0 - zTerm;
92template<
class Scalar,
class ArrayPo
int,
class ArrayWeight>
94 ArrayWeight& cubWeights,
95 ArrayPoint& cellCoords)
const
97 TEUCHOS_TEST_FOR_EXCEPTION( (
true), std::logic_error,
98 ">>> ERROR (CubatureTensorPyr): Cubature defined in reference space calling method for physical space cubature.");
101template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
107template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
113template <
class Scalar,
class ArrayPo
int,
class ArrayWeight>
Defines direct cubature (integration) rules in Intrepid.
virtual void getAccuracy(std::vector< int > °ree) const
Returns max. degree of polynomials that are integrated exactly. The return vector has the size of the...
virtual int getDimension() const
Returns dimension of integration domain.
void getCubature(ArrayPoint &cubPoints, ArrayWeight &cubWeights) const
Returns cubature points and weights (return arrays must be pre-sized/pre-allocated).
virtual int getNumPoints() const
Returns the number of cubature points.
CubatureTensorPyr(std::vector< Teuchos::RCP< Cubature< Scalar, ArrayPoint, ArrayWeight > > > cubatures)
Constructor.
Defines tensor-product cubature (integration) rules in Intrepid.
virtual int getNumPoints() const
Returns the number of cubature points.
virtual void getAccuracy(std::vector< int > °ree) const
Returns max. degree of polynomials that are integrated exactly. The return vector has the size of the...
virtual int getDimension() const
Returns dimension of integration domain.
virtual void getCubature(ArrayPoint &cubPoints, ArrayWeight &cubWeights) const
Returns cubature points and weights (return arrays must be pre-sized/pre-allocated).
Defines the base class for cubature (integration) rules in Intrepid.