ReferenceElementSets.h
Go to the documentation of this file.
/*****************************************************************************
*
* Copyright (c) 2003-2020 by The University of Queensland
* http://www.uq.edu.au
*
* Primary Business: Queensland, Australia
* Licensed under the Apache License, version 2.0
* http://www.apache.org/licenses/LICENSE-2.0
*
* Development until 2012 by Earth Systems Science Computational Center (ESSCC)
* Development 2012-2013 by School of Earth Sciences
* Development from 2014-2017 by Centre for Geoscience Computing (GeoComp)
* Development from 2019 by School of Earth and Environmental Sciences
**
*****************************************************************************/
#ifndef __FINLEY_REFERENCEELEMENTSETS_H__
#define __FINLEY_REFERENCEELEMENTSETS_H__
#include "ReferenceElements.h"
namespace finley {
struct ReferenceElementSet {
ReferenceElementSet(ElementTypeId id, int order, int reduced_order)
{
const ReferenceElementInfo* id_info = ReferenceElement::getInfo(id);
const ShapeFunctionInfo* bf_info = ShapeFunction::getInfo(
id_info->BasisFunctions);
if (order<0)
order=std::max(2*bf_info->numOrder, 0);
referenceElement.reset(new ReferenceElement(id, order));
if (reduced_order<0)
reduced_order=std::max(2*(bf_info->numOrder-1), 0);
referenceElementReducedQuadrature.reset(new ReferenceElement(id,
reduced_order));
if (referenceElement->getNumNodes() != referenceElementReducedQuadrature->getNumNodes()) {
throw escript::ValueError("ReferenceElementSet: numNodes in referenceElement and referenceElementReducedQuadrature don't match.");
}
}
const_ShapeFunction_ptr borrowBasisFunctions(bool reducedShapefunction,
bool reducedIntegrationOrder) const
{
if (reducedShapefunction) {
return (reducedIntegrationOrder ?
referenceElementReducedQuadrature->LinearBasisFunctions :
referenceElement->LinearBasisFunctions);
}
return (reducedIntegrationOrder ?
referenceElementReducedQuadrature->BasisFunctions :
referenceElement->BasisFunctions);
}
const_ShapeFunction_ptr borrowParametrization(bool reducedIntegrationOrder) const
{
return (reducedIntegrationOrder ?
referenceElementReducedQuadrature->Parametrization :
referenceElement->Parametrization);
}
const_ReferenceElement_ptr borrowReferenceElement(bool reducedIntOrder) const
{
return (reducedIntOrder ? referenceElementReducedQuadrature :
referenceElement);
}
inline int getNumNodes() const { return referenceElement->getNumNodes(); }
ReferenceElement_ptr referenceElementReducedQuadrature;
ReferenceElement_ptr referenceElement;
};
typedef boost::shared_ptr<const ReferenceElementSet> const_ReferenceElementSet_ptr;
} // namespace finley
#endif // __FINLEY_REFERENCEELEMENTSETS_H__
ReferenceElements.h
escript::ValueError
An exception class that signals an invalid argument value.
Definition: EsysException.h:100
finley
A suite of factory methods for creating various finley domains.
Definition: finley/src/Assemble.h:32
finley::const_ReferenceElement_ptr
boost::shared_ptr< const ReferenceElement > const_ReferenceElement_ptr
Definition: ReferenceElements.h:216
finley::const_ShapeFunction_ptr
boost::shared_ptr< const ShapeFunction > const_ShapeFunction_ptr
Definition: ShapeFunctions.h:102
finley::ReferenceElement_ptr
boost::shared_ptr< ReferenceElement > ReferenceElement_ptr
Definition: ReferenceElements.h:215
finley::ElementTypeId
ElementTypeId
Definition: ReferenceElements.h:41
finley::const_ReferenceElementSet_ptr
boost::shared_ptr< const ReferenceElementSet > const_ReferenceElementSet_ptr
Definition: ReferenceElementSets.h:80
finley::ReferenceElementInfo
this struct holds the definition of the reference element
Definition: ReferenceElements.h:123
finley::ReferenceElementInfo::BasisFunctions
ShapeFunctionTypeId BasisFunctions
shape function for the basis functions
Definition: ReferenceElements.h:148
finley::ReferenceElementSet
Definition: ReferenceElementSets.h:27
finley::ReferenceElementSet::borrowParametrization
const_ShapeFunction_ptr borrowParametrization(bool reducedIntegrationOrder) const
Definition: ReferenceElementSets.h:60
finley::ReferenceElementSet::getNumNodes
int getNumNodes() const
Definition: ReferenceElementSets.h:73
finley::ReferenceElementSet::referenceElementReducedQuadrature
ReferenceElement_ptr referenceElementReducedQuadrature
Definition: ReferenceElementSets.h:75
finley::ReferenceElementSet::referenceElement
ReferenceElement_ptr referenceElement
Definition: ReferenceElementSets.h:76
finley::ReferenceElementSet::borrowBasisFunctions
const_ShapeFunction_ptr borrowBasisFunctions(bool reducedShapefunction, bool reducedIntegrationOrder) const
Definition: ReferenceElementSets.h:47
finley::ReferenceElementSet::ReferenceElementSet
ReferenceElementSet(ElementTypeId id, int order, int reduced_order)
Definition: ReferenceElementSets.h:28
finley::ReferenceElementSet::borrowReferenceElement
const_ReferenceElement_ptr borrowReferenceElement(bool reducedIntOrder) const
Definition: ReferenceElementSets.h:67
finley::ReferenceElement
this struct holds the realization of a reference element
Definition: ReferenceElements.h:179
finley::ReferenceElement::getInfo
static const ReferenceElementInfo * getInfo(ElementTypeId id)
returns the element information structure for the given type id
Definition: ReferenceElements.cpp:664
finley::ShapeFunctionInfo
this struct holds the definition of the shape functions on an element
Definition: ShapeFunctions.h:60
finley::ShapeFunctionInfo::numOrder
int numOrder
order of the shape functions
Definition: ShapeFunctions.h:70
finley::ShapeFunction::getInfo
static const ShapeFunctionInfo * getInfo(ShapeFunctionTypeId id)
Definition: ShapeFunctions.cpp:105