Public Member Functions | Private Member Functions | Private Attributes
Sundance::MaximalQuadratureIntegral Class Reference
Inheritance diagram for Sundance::MaximalQuadratureIntegral:
Sundance::ElementIntegral

List of all members.

Public Member Functions

 MaximalQuadratureIntegral (const CellType &maxCellType, const QuadratureFamily &quad, const ParametrizedCurve &globalCurve, const Mesh &mesh, int verb)
 MaximalQuadratureIntegral (const CellType &maxCellType, const BasisFamily &testBasis, int alpha, int testDerivOrder, const QuadratureFamily &quad, const ParametrizedCurve &globalCurve, const Mesh &mesh, int verb)
 MaximalQuadratureIntegral (const CellType &maxCellType, 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 ~MaximalQuadratureIntegral ()
virtual void transform (const CellJacobianBatch &JTrans, const CellJacobianBatch &JVol, const Array< int > &isLocalFlag, const Array< int > &facetNum, const RCP< Array< int > > &cellLIDs, 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 *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 *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 *const coeff, RCP< Array< double > > &A) const

Private Member Functions

void transformSummingFirst (int nCells, const Array< int > &facetIndex, const RCP< Array< int > > &cellLIDs, 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 double & wValue (int facetCase, int q, int testDerivDir, int testNode, int unkDerivDir, int unkNode) const
double & wValue (int q, int testDerivDir, int testNode)
const double & wValue (int q, int testDerivDir, int testNode) const

Private Attributes

Array< double > W_
bool useSumFirstMethod_
Array< Array< Array< double > > > W_ACI_F1_
Array< Array< Array< Array
< Array< double > > > > > 		W_ACI_F2_
QuadratureFamily quad_
Array< Point > quadPts_
Array< double > quadWeights_

Detailed Description

Integration by quadrature over a maximal cell

Definition at line 56 of file SundanceMaximalQuadratureIntegral.hpp.

Constructor & Destructor Documentation

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

Construct a zero-form to be computed by quadrature

Definition at line 89 of file SundanceMaximalQuadratureIntegral.cpp.

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

MaximalQuadratureIntegral::MaximalQuadratureIntegral ( const CellType &  maxCellType,
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 123 of file SundanceMaximalQuadratureIntegral.cpp.

References Sundance::ElementIntegral::addFlops(), Sundance::ElementIntegral::assertLinearForm(), Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::describe(), Sundance::BasisFamily::dim(), Sundance::ElementIntegral::evalCellType(), Sundance::QuadratureFamily::getPoints(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nRefDerivTest(), Playa::Out::os(), quadPts_, quadWeights_, Sundance::BasisFamily::refEval(), Sundance::ElementIntegral::setupVerb(), SUNDANCE_MSG1, W_, W_ACI_F1_, and wValue().

MaximalQuadratureIntegral::MaximalQuadratureIntegral ( const CellType &  maxCellType,
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 195 of file SundanceMaximalQuadratureIntegral.cpp.

References Sundance::ElementIntegral::addFlops(), Sundance::ElementIntegral::assertBilinearForm(), Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::describe(), Sundance::BasisFamily::dim(), Sundance::ElementIntegral::evalCellType(), Sundance::QuadratureFamily::getPoints(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::ElementIntegral::nRefDerivTest(), Sundance::ElementIntegral::nRefDerivUnk(), Playa::Out::os(), quadPts_, quadWeights_, Sundance::BasisFamily::refEval(), Sundance::ElementIntegral::setupVerb(), SUNDANCE_MSG1, W_, W_ACI_F2_, and wValue().

virtual Sundance::MaximalQuadratureIntegral::~MaximalQuadratureIntegral ( ) [inline, virtual]

virtual dtor

Definition at line 93 of file SundanceMaximalQuadratureIntegral.hpp.

Member Function Documentation

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

Definition at line 96 of file SundanceMaximalQuadratureIntegral.hpp.

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

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

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

Definition at line 387 of file SundanceMaximalQuadratureIntegral.cpp.

References Sundance::ElementIntegral::addFlops(), Sundance::ElementIntegral::alpha(), Sundance::ElementIntegral::cellType(), Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::createOneFormTransformationMatrix(), Sundance::CellJacobianBatch::detJ(), Sundance::ElementIntegral::dim(), Sundance::ElementIntegral::G(), Sundance::QuadratureFamily::getAdaptedWeights(), Sundance::ElementIntegral::globalCurve(), Sundance::ElementIntegral::integrationVerb(), maxCellQuadrature1Timer(), Sundance::ElementIntegral::mesh(), Sundance::ElementIntegral::nNodes(), Sundance::CellJacobianBatch::numCells(), Sundance::ElementIntegral::order(), Playa::Out::os(), quad_, quadPts_, quadWeights_, SUNDANCE_MSG2, SUNDANCE_MSG4, SUNDANCE_MSG5, Sundance::ElementIntegral::testDerivOrder(), transformSummingFirst(), Sundance::ElementIntegral::transformVerb(), W_, W_ACI_F1_, and Sundance::writeTable().

Referenced by transform().

void MaximalQuadratureIntegral::transformSummingFirst ( 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 summing reference quantities over quadrature points and then transforming the sum to physical quantities.

Definition at line 673 of file SundanceMaximalQuadratureIntegral.cpp.

Referenced by transformOneForm(), and transformTwoForm().

void Sundance::MaximalQuadratureIntegral::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 MaximalQuadratureIntegral::transformTwoForm ( const CellJacobianBatch &  JTrans,
const CellJacobianBatch &  JVol,
const Array< int > &  facetIndex,
const RCP< Array< int > > &  cellLIDs,
const double *const  coeff,
RCP< Array< double > > &  A 
) const [virtual]

Definition at line 540 of file SundanceMaximalQuadratureIntegral.cpp.

References Sundance::ElementIntegral::addFlops(), Sundance::ElementIntegral::alpha(), Sundance::ElementIntegral::beta(), Sundance::ElementIntegral::cellType(), Sundance::ElementIntegral::chop(), Sundance::ElementIntegral::createTwoFormTransformationMatrix(), Sundance::CellJacobianBatch::detJ(), Sundance::ElementIntegral::dim(), Sundance::ElementIntegral::G(), Sundance::QuadratureFamily::getAdaptedWeights(), Sundance::ElementIntegral::globalCurve(), Sundance::ElementIntegral::integrationVerb(), maxCellQuadrature2Timer(), Sundance::ElementIntegral::mesh(), Sundance::ElementIntegral::nNodes(), Sundance::ElementIntegral::nNodesTest(), Sundance::ElementIntegral::nNodesUnk(), Sundance::CellJacobianBatch::numCells(), Sundance::ElementIntegral::order(), quad_, quadPts_, quadWeights_, SUNDANCE_MSG2, Sundance::ElementIntegral::testDerivOrder(), transformSummingFirst(), Sundance::ElementIntegral::transformVerb(), Sundance::ElementIntegral::unkDerivOrder(), W_, and W_ACI_F2_.

Referenced by transform().

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

Definition at line 294 of file SundanceMaximalQuadratureIntegral.cpp.

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

Referenced by transform().

bool Sundance::MaximalQuadratureIntegral::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 156 of file SundanceMaximalQuadratureIntegral.hpp.

References useSumFirstMethod_.

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

Definition at line 159 of file SundanceMaximalQuadratureIntegral.hpp.

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

Referenced by MaximalQuadratureIntegral().

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

Definition at line 170 of file SundanceMaximalQuadratureIntegral.hpp.

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

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

Definition at line 183 of file SundanceMaximalQuadratureIntegral.hpp.

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

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

Definition at line 188 of file SundanceMaximalQuadratureIntegral.hpp.

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

Member Data Documentation

QuadratureFamily Sundance::MaximalQuadratureIntegral::quad_ [private]

The quadrature family needed for special integration (ACI)

Definition at line 206 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by transformOneForm(), and transformTwoForm().

Array<Point> Sundance::MaximalQuadratureIntegral::quadPts_ [private]

The quadrature points

Definition at line 209 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by MaximalQuadratureIntegral(), transformOneForm(), and transformTwoForm().

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

The standard weights (in case of ACI these might change)

Definition at line 212 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by MaximalQuadratureIntegral(), transformOneForm(), and transformTwoForm().

bool Sundance::MaximalQuadratureIntegral::useSumFirstMethod_ [private]

Definition at line 195 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by useSumFirstMethod().

Array<double> Sundance::MaximalQuadratureIntegral::W_ [private]

Definition at line 192 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by MaximalQuadratureIntegral(), transformOneForm(), transformTwoForm(), and wValue().

Array<Array<Array<double> > > Sundance::MaximalQuadratureIntegral::W_ACI_F1_ [private]

For ACI (ACI = Adaptive Cell Integration), store the reference integral values for one form The indexes facet, quadPoints, nRefDerivTest , nNodesTest

Definition at line 199 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by MaximalQuadratureIntegral(), and transformOneForm().

Array<Array<Array<Array<Array<double> > > > > Sundance::MaximalQuadratureIntegral::W_ACI_F2_ [private]

For ACI (ACI = Adaptive Cell Integration), store the reference integral values for two form The indexes quadPoints, nRefDerivTest , nNodesTest , nRefDerivUnk , nNodesUnk

Definition at line 203 of file SundanceMaximalQuadratureIntegral.hpp.

Referenced by MaximalQuadratureIntegral(), and transformTwoForm().

Site Contact