Sundance::CurveQuadratureIntegral Class Reference

Inheritance diagram for Sundance::CurveQuadratureIntegral:

List of all members.

Public Member Functions
	CurveQuadratureIntegral (const CellType &maxCellType, const bool isConstantIntegral, const QuadratureFamily &quad, const ParametrizedCurve &globalCurve, const Mesh &mesh, int verb)
	CurveQuadratureIntegral (const CellType &maxCellType, const bool isConstantIntegral, const BasisFamily &testBasis, int alpha, int testDerivOrder, const QuadratureFamily &quad, const ParametrizedCurve &globalCurve, const Mesh &mesh, int verb)
	CurveQuadratureIntegral (const CellType &maxCellType, const bool isConstantIntegral, const BasisFamily &testBasis, int alpha, int testDerivOrder, const BasisFamily &unkBasis, int beta, int unkDerivOrder, const QuadratureFamily &quad, const ParametrizedCurve &globalCurve, const Mesh &mesh, int verb)
virtual	~CurveQuadratureIntegral ()
virtual void	transform (const CellJacobianBatch &JTrans, const CellJacobianBatch &JVol, const Array< int > &isLocalFlag, const Array< int > &facetNum, const RCP< Array< int > > &cellLIDs, const double constCoeff, const double *const coeff, RCP< Array< double > > &A) const
virtual void	transformZeroForm (const CellJacobianBatch &JTrans, const CellJacobianBatch &JVol, const Array< int > &isLocalFlag, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, const double constCoeff, const double *const coeff, RCP< Array< double > > &A) const
virtual void	transformTwoForm (const CellJacobianBatch &JTrans, const CellJacobianBatch &JVol, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, const double constCoeff, const double *const coeff, RCP< Array< double > > &A) const
void	transformOneForm (const CellJacobianBatch &JTrans, const CellJacobianBatch &JVol, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, const double constCoeff, const double *const coeff, RCP< Array< double > > &A) const
Private Member Functions
void	updateRefCellInformation (int maxCellLID, const ParametrizedCurve &curve) const
void	updateRefCellIntegralOneForm (int maxCellLID, int cellInBatch) const
void	updateRefCellIntegralTwoForm (int maxCellLID, int cellInBatch) const
void	transformSummingFirst (int nCells, const CellJacobianBatch &JVol, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, const double constCoeff, const double *const GPtr, const double *const coeff, RCP< Array< double > > &A) const
void	transformSummingLast (int nCells, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, const double *const GPtr, const double *const coeff, RCP< Array< double > > &A) const
bool	useSumFirstMethod () const
double &	wValue (int q, int testDerivDir, int testNode, int unkDerivDir, int unkNode) const
const double &	wValue (int facetCase, int q, int testDerivDir, int testNode, int unkDerivDir, int unkNode) const
double &	wValue (int q, int testDerivDir, int testNode) const
Private Attributes
QuadratureFamily	quad_
Array< Point >	quadPts_
Array< double >	quadWeights_
Array< double >	W_
bool	useSumFirstMethod_
bool	useConstCoeff_
Array< Point >	quadCurveDerivs_
Array< Point >	quadCurveNormals_

---

Detailed Description

Integration by quadrature over a maximal cell

Definition at line 56 of file SundanceCurveQuadratureIntegral.hpp.

---

Constructor & Destructor Documentation

CurveQuadratureIntegral::CurveQuadratureIntegral	(	const CellType &	maxCellType,
		const bool	isConstantIntegral,
		const QuadratureFamily &	quad,
		const ParametrizedCurve &	globalCurve,
		const Mesh &	mesh,
		int	verb
	)

Construct a zero-form to be computed by quadrature

Definition at line 77 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::Mesh::cellType(), Sundance::ElementIntegral::describe(), Sundance::QuadratureFamily::getPoints(), Playa::Out::os(), quadCurveDerivs_, quadCurveNormals_, quadPts_, quadWeights_, Sundance::ElementIntegral::setupVerb(), Sundance::Mesh::spatialDim(), SUNDANCE_MSG1, and W_.

CurveQuadratureIntegral::CurveQuadratureIntegral	(	const CellType &	maxCellType,
		const bool	isConstantIntegral,
		const BasisFamily &	testBasis,
		int	alpha,
		int	testDerivOrder,
		const QuadratureFamily &	quad,
		const ParametrizedCurve &	globalCurve,
		const Mesh &	mesh,
		int	verb
	)

Construct a one form to be computed by quadrature

Definition at line 117 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::assertLinearForm(), Sundance::Mesh::cellType(), Sundance::ElementIntegral::describe(), Sundance::QuadratureFamily::getPoints(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nRefDerivTest(), Playa::Out::os(), quadCurveDerivs_, quadCurveNormals_, quadPts_, quadWeights_, Sundance::ElementIntegral::setupVerb(), Sundance::Mesh::spatialDim(), SUNDANCE_MSG1, and W_.

CurveQuadratureIntegral::CurveQuadratureIntegral	(	const CellType &	maxCellType,
		const bool	isConstantIntegral,
		const BasisFamily &	testBasis,
		int	alpha,
		int	testDerivOrder,
		const BasisFamily &	unkBasis,
		int	beta,
		int	unkDerivOrder,
		const QuadratureFamily &	quad,
		const ParametrizedCurve &	globalCurve,
		const Mesh &	mesh,
		int	verb
	)

Construct a two-form to be computed by quadrature

Definition at line 161 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::assertBilinearForm(), Sundance::Mesh::cellType(), Sundance::ElementIntegral::describe(), Sundance::QuadratureFamily::getPoints(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::ElementIntegral::nRefDerivTest(), Sundance::ElementIntegral::nRefDerivUnk(), Playa::Out::os(), quadCurveDerivs_, quadCurveNormals_, quadPts_, quadWeights_, Sundance::ElementIntegral::setupVerb(), Sundance::Mesh::spatialDim(), SUNDANCE_MSG1, and W_.

virtual Sundance::CurveQuadratureIntegral::~CurveQuadratureIntegral

(

)

[inline, virtual]

virtual dtor

Definition at line 96 of file SundanceCurveQuadratureIntegral.hpp.

---

Member Function Documentation

virtual void Sundance::CurveQuadratureIntegral::transform	(	const CellJacobianBatch &	JTrans,
		const CellJacobianBatch &	JVol,
		const Array< int > &	isLocalFlag,
		const Array< int > &	facetNum,
		const RCP< Array< int > > &	cellLIDs,
		const double	constCoeff,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const [inline, virtual]

Definition at line 99 of file SundanceCurveQuadratureIntegral.hpp.

References Sundance::ElementIntegral::order(), transformOneForm(), transformTwoForm(), and transformZeroForm().

Referenced by Sundance::IntegralGroup::evaluate().

void CurveQuadratureIntegral::transformOneForm	(	const CellJacobianBatch &	JTrans,
		const CellJacobianBatch &	JVol,
		const Array< int > &	facetIndex,
		const RCP< Array< int > > &	cellLIDs,
		const double	constCoeff,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const

Definition at line 435 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::addFlops(), Sundance::ElementIntegral::alpha(), Sundance::ElementIntegral::createOneFormTransformationMatrix(), Sundance::CellJacobianBatch::detJ(), Sundance::ElementIntegral::G(), Sundance::ElementIntegral::integrationVerb(), maxCellQuadratureTimer(), Sundance::ElementIntegral::nNodes(), Sundance::CellJacobianBatch::numCells(), Sundance::ElementIntegral::order(), Playa::Out::os(), quadCurveDerivs_, quadWeights_, SUNDANCE_MSG2, SUNDANCE_MSG4, SUNDANCE_MSG5, Sundance::ElementIntegral::testDerivOrder(), transformSummingFirst(), Sundance::ElementIntegral::transformVerb(), updateRefCellIntegralOneForm(), useConstCoeff_, W_, and Sundance::writeTable().

Referenced by transform().

void CurveQuadratureIntegral::transformSummingFirst	(	int	nCells,
		const CellJacobianBatch &	JVol,
		const Array< int > &	facetIndex,
		const RCP< Array< int > > &	cellLIDs,
		const double	constCoeff,
		const double *const	GPtr,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const [private]

Do the integration by summing reference quantities over quadrature points and then transforming the sum to physical quantities.

Definition at line 619 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::CellJacobianBatch::detJ(), and Sundance::CellJacobianBatch::numCells().

Referenced by transformOneForm(), and transformTwoForm().

void Sundance::CurveQuadratureIntegral::transformSummingLast	(	int	nCells,
		const Array< int > &	facetIndex,
		const RCP< Array< int > > &	cellLIDs,
		const double *const	GPtr,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const [private]

Do the integration by transforming to physical coordinates at each quadrature point, and then summing

void CurveQuadratureIntegral::transformTwoForm	(	const CellJacobianBatch &	JTrans,
		const CellJacobianBatch &	JVol,
		const Array< int > &	facetIndex,
		const RCP< Array< int > > &	cellLIDs,
		const double	constCoeff,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const [virtual]

Definition at line 532 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::alpha(), Sundance::ElementIntegral::beta(), Sundance::ElementIntegral::createTwoFormTransformationMatrix(), Sundance::ElementIntegral::G(), Sundance::ElementIntegral::integrationVerb(), maxCellQuadratureTimer(), Sundance::ElementIntegral::nNodes(), Sundance::CellJacobianBatch::numCells(), Sundance::ElementIntegral::order(), quadCurveDerivs_, quadWeights_, SUNDANCE_MSG2, Sundance::ElementIntegral::testDerivOrder(), transformSummingFirst(), Sundance::ElementIntegral::transformVerb(), Sundance::ElementIntegral::unkDerivOrder(), updateRefCellIntegralTwoForm(), useConstCoeff_, and W_.

Referenced by transform().

void CurveQuadratureIntegral::transformZeroForm	(	const CellJacobianBatch &	JTrans,
		const CellJacobianBatch &	JVol,
		const Array< int > &	isLocalFlag,
		const Array< int > &	facetIndex,
		const RCP< Array< int > > &	cellLIDs,
		const double	constCoeff,
		const double *const	coeff,
		RCP< Array< double > > &	A
	)		const [virtual]

Definition at line 341 of file SundanceCurveQuadratureIntegral.cpp.

References maxCellQuadratureTimer(), Sundance::CellJacobianBatch::numCells(), Playa::Out::os(), SUNDANCE_MSG1, SUNDANCE_MSG5, and Sundance::writeTable().

Referenced by transform().

void CurveQuadratureIntegral::updateRefCellInformation	(	int	maxCellLID,
		const ParametrizedCurve &	curve
	)		const [private]

updates the reference cell information (quadPoints,derivatives,normals)

Definition at line 207 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::cellType(), Sundance::Mesh::getCurvePoints(), Sundance::CurveIntegralCalc::getCurveQuadPoints(), Sundance::Mesh::hasCurvePoints(), Sundance::ElementIntegral::integrationVerb(), Sundance::ElementIntegral::mesh(), quad_, quadCurveDerivs_, quadCurveNormals_, quadPts_, Sundance::Mesh::setCurvePoints(), Sundance::ElementIntegral::spatialDim(), and SUNDANCE_MSG3.

Referenced by updateRefCellIntegralOneForm(), and updateRefCellIntegralTwoForm().

void CurveQuadratureIntegral::updateRefCellIntegralOneForm	(	int	maxCellLID,
		int	cellInBatch
	)		const [private]

updates the W_ for the given cell for one form Integral

Definition at line 232 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::dim(), Sundance::ElementIntegral::evalCellType(), Sundance::ElementIntegral::globalCurve(), Sundance::ElementIntegral::integrationVerb(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nRefDerivTest(), quadPts_, quadWeights_, Sundance::BasisFamily::refEval(), Sundance::ElementIntegral::setupVerb(), SUNDANCE_MSG1, SUNDANCE_MSG2, Sundance::ElementIntegral::testBasis(), Sundance::ElementIntegral::testDerivOrder(), updateRefCellInformation(), W_, and wValue().

Referenced by transformOneForm().

void CurveQuadratureIntegral::updateRefCellIntegralTwoForm	(	int	maxCellLID,
		int	cellInBatch
	)		const [private]

updates the W_ for the given cell for two form Integral

Definition at line 276 of file SundanceCurveQuadratureIntegral.cpp.

References Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::dim(), Sundance::ElementIntegral::evalCellType(), Sundance::ElementIntegral::globalCurve(), Sundance::ElementIntegral::integrationVerb(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::ElementIntegral::nRefDerivTest(), Sundance::ElementIntegral::nRefDerivUnk(), quadPts_, quadWeights_, Sundance::BasisFamily::refEval(), Sundance::ElementIntegral::setupVerb(), SUNDANCE_MSG1, SUNDANCE_MSG2, Sundance::ElementIntegral::testBasis(), Sundance::ElementIntegral::testDerivOrder(), Sundance::ElementIntegral::unkBasis(), Sundance::ElementIntegral::unkDerivOrder(), updateRefCellInformation(), W_, and wValue().

Referenced by transformTwoForm().

bool Sundance::CurveQuadratureIntegral::useSumFirstMethod

(

)

const [inline, private]

Determine whether to do this batch of integrals using the sum-first method or the sum-last method

Definition at line 174 of file SundanceCurveQuadratureIntegral.hpp.

References useSumFirstMethod_.

double& Sundance::CurveQuadratureIntegral::wValue	(	int	q,
		int	testDerivDir,
		int	testNode,
		int	unkDerivDir,
		int	unkNode
	)		const [inline, private]

Definition at line 177 of file SundanceCurveQuadratureIntegral.hpp.

References Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::ElementIntegral::nRefDerivTest(), Sundance::ElementIntegral::nRefDerivUnk(), and W_.

Referenced by updateRefCellIntegralOneForm(), and updateRefCellIntegralTwoForm().

const double& Sundance::CurveQuadratureIntegral::wValue	(	int	facetCase,
		int	q,
		int	testDerivDir,
		int	testNode,
		int	unkDerivDir,
		int	unkNode
	)		const [inline, private]

Definition at line 188 of file SundanceCurveQuadratureIntegral.hpp.

References Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::ElementIntegral::nRefDerivTest(), Sundance::ElementIntegral::nRefDerivUnk(), and W_.

double& Sundance::CurveQuadratureIntegral::wValue	(	int	q,
		int	testDerivDir,
		int	testNode
	)		const [inline, private]

Definition at line 201 of file SundanceCurveQuadratureIntegral.hpp.

References Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nRefDerivTest(), and W_.

---

Member Data Documentation

QuadratureFamily Sundance::CurveQuadratureIntegral::quad_ [private]

The quadrature family needed for curve integration

Definition at line 205 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by updateRefCellInformation().

Array<Point> Sundance::CurveQuadratureIntegral::quadCurveDerivs_ [mutable, private]

The derivative of the curve at the curve

Definition at line 224 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by CurveQuadratureIntegral(), transformOneForm(), transformTwoForm(), and updateRefCellInformation().

Array<Point> Sundance::CurveQuadratureIntegral::quadCurveNormals_ [mutable, private]

The derivative of the curve at the curve

Definition at line 227 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by CurveQuadratureIntegral(), and updateRefCellInformation().

Array<Point> Sundance::CurveQuadratureIntegral::quadPts_ [mutable, private]

The quadrature points

Definition at line 208 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by CurveQuadratureIntegral(), updateRefCellInformation(), updateRefCellIntegralOneForm(), and updateRefCellIntegralTwoForm().

Array<double> Sundance::CurveQuadratureIntegral::quadWeights_ [private]

The standard weights

Definition at line 211 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by CurveQuadratureIntegral(), transformOneForm(), transformTwoForm(), updateRefCellIntegralOneForm(), and updateRefCellIntegralTwoForm().

bool Sundance::CurveQuadratureIntegral::useConstCoeff_ [private]

weather to use constant coefficients

Definition at line 221 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by transformOneForm(), and transformTwoForm().

bool Sundance::CurveQuadratureIntegral::useSumFirstMethod_ [private]

Definition at line 218 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by useSumFirstMethod().

Array<double> Sundance::CurveQuadratureIntegral::W_ [mutable, private]

Definition at line 215 of file SundanceCurveQuadratureIntegral.hpp.

Referenced by CurveQuadratureIntegral(), transformOneForm(), transformTwoForm(), updateRefCellIntegralOneForm(), updateRefCellIntegralTwoForm(), and wValue().