Public Member Functions | Static Public Member Functions | List of all members
FEGenericBase< T > Class Template Reference

#include <exact_error_estimator.h>

Public Member Functions

template<>
AutoPtr< FEGenericBase
< RealGradient > > 		build (const unsigned int dim, const FEType &fet)
 
template<>
AutoPtr< FEGenericBase< Real > > build_InfFE (const unsigned int dim, const FEType &fet)
 
template<>
AutoPtr< FEGenericBase
< RealGradient > > 		build_InfFE (const unsigned int, const FEType &)
 

Static Public Member Functions

template<>
AutoPtr< FEGenericBase< Real > > build (const unsigned int dim, const FEType &fet)
 

Detailed Description

template<typename T>
class FEGenericBase< T >

Definition at line 36 of file exact_error_estimator.h.

Member Function Documentation

template<>
AutoPtr< FEGenericBase< Real > > FEGenericBase< Real >::build ( const unsigned int  dim,
const FEType &  fet 
)
static

Definition at line 184 of file fe_base.C.

References libMeshEnums::BERNSTEIN, libMeshEnums::CLOUGH, libMesh::FEType::family, libMeshEnums::HERMITE, libMeshEnums::HIERARCHIC, libMeshEnums::L2_HIERARCHIC, libMeshEnums::L2_LAGRANGE, libMeshEnums::LAGRANGE, libMeshEnums::MONOMIAL, libMesh::out, libMeshEnums::SCALAR, libMeshEnums::SZABAB, and libMeshEnums::XYZ.

{
// The stupid AutoPtr<FEBase> ap(); return ap;
// construct is required to satisfy IBM's xlC
switch (dim)
{
// 0D
case 0:
{
switch (fet.family)
{
case CLOUGH:
{
AutoPtr<FEBase> ap(new FE<0,CLOUGH>(fet));
return ap;
}
case HERMITE:
{
AutoPtr<FEBase> ap(new FE<0,HERMITE>(fet));
return ap;
}
case LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<0,LAGRANGE>(fet));
return ap;
}
case L2_LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<0,L2_LAGRANGE>(fet));
return ap;
}
case HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<0,HIERARCHIC>(fet));
return ap;
}
case L2_HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<0,L2_HIERARCHIC>(fet));
return ap;
}
case MONOMIAL:
{
AutoPtr<FEBase> ap(new FE<0,MONOMIAL>(fet));
return ap;
}
#ifdef LIBMESH_ENABLE_HIGHER_ORDER_SHAPES
case SZABAB:
{
AutoPtr<FEBase> ap(new FE<0,SZABAB>(fet));
return ap;
}
case BERNSTEIN:
{
AutoPtr<FEBase> ap(new FE<0,BERNSTEIN>(fet));
return ap;
}
#endif
case XYZ:
{
AutoPtr<FEBase> ap(new FEXYZ<0>(fet));
return ap;
}
case SCALAR:
{
AutoPtr<FEBase> ap(new FEScalar<0>(fet));
return ap;
}
default:
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
// 1D
case 1:
{
switch (fet.family)
{
case CLOUGH:
{
AutoPtr<FEBase> ap(new FE<1,CLOUGH>(fet));
return ap;
}
case HERMITE:
{
AutoPtr<FEBase> ap(new FE<1,HERMITE>(fet));
return ap;
}
case LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<1,LAGRANGE>(fet));
return ap;
}
case L2_LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<1,L2_LAGRANGE>(fet));
return ap;
}
case HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<1,HIERARCHIC>(fet));
return ap;
}
case L2_HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<1,L2_HIERARCHIC>(fet));
return ap;
}
case MONOMIAL:
{
AutoPtr<FEBase> ap(new FE<1,MONOMIAL>(fet));
return ap;
}
#ifdef LIBMESH_ENABLE_HIGHER_ORDER_SHAPES
case SZABAB:
{
AutoPtr<FEBase> ap(new FE<1,SZABAB>(fet));
return ap;
}
case BERNSTEIN:
{
AutoPtr<FEBase> ap(new FE<1,BERNSTEIN>(fet));
return ap;
}
#endif
case XYZ:
{
AutoPtr<FEBase> ap(new FEXYZ<1>(fet));
return ap;
}
case SCALAR:
{
AutoPtr<FEBase> ap(new FEScalar<1>(fet));
return ap;
}
default:
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
// 2D
case 2:
{
switch (fet.family)
{
case CLOUGH:
{
AutoPtr<FEBase> ap(new FE<2,CLOUGH>(fet));
return ap;
}
case HERMITE:
{
AutoPtr<FEBase> ap(new FE<2,HERMITE>(fet));
return ap;
}
case LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<2,LAGRANGE>(fet));
return ap;
}
case L2_LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<2,L2_LAGRANGE>(fet));
return ap;
}
case HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<2,HIERARCHIC>(fet));
return ap;
}
case L2_HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<2,L2_HIERARCHIC>(fet));
return ap;
}
case MONOMIAL:
{
AutoPtr<FEBase> ap(new FE<2,MONOMIAL>(fet));
return ap;
}
#ifdef LIBMESH_ENABLE_HIGHER_ORDER_SHAPES
case SZABAB:
{
AutoPtr<FEBase> ap(new FE<2,SZABAB>(fet));
return ap;
}
case BERNSTEIN:
{
AutoPtr<FEBase> ap(new FE<2,BERNSTEIN>(fet));
return ap;
}
#endif
case XYZ:
{
AutoPtr<FEBase> ap(new FEXYZ<2>(fet));
return ap;
}
case SCALAR:
{
AutoPtr<FEBase> ap(new FEScalar<2>(fet));
return ap;
}
default:
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
// 3D
case 3:
{
switch (fet.family)
{
case CLOUGH:
{
libMesh::out << "ERROR: Clough-Tocher elements currently only support 1D and 2D"
<< std::endl;
libmesh_error();
}
case HERMITE:
{
AutoPtr<FEBase> ap(new FE<3,HERMITE>(fet));
return ap;
}
case LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<3,LAGRANGE>(fet));
return ap;
}
case L2_LAGRANGE:
{
AutoPtr<FEBase> ap(new FE<3,L2_LAGRANGE>(fet));
return ap;
}
case HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<3,HIERARCHIC>(fet));
return ap;
}
case L2_HIERARCHIC:
{
AutoPtr<FEBase> ap(new FE<3,L2_HIERARCHIC>(fet));
return ap;
}
case MONOMIAL:
{
AutoPtr<FEBase> ap(new FE<3,MONOMIAL>(fet));
return ap;
}
#ifdef LIBMESH_ENABLE_HIGHER_ORDER_SHAPES
case SZABAB:
{
AutoPtr<FEBase> ap(new FE<3,SZABAB>(fet));
return ap;
}
case BERNSTEIN:
{
AutoPtr<FEBase> ap(new FE<3,BERNSTEIN>(fet));
return ap;
}
#endif
case XYZ:
{
AutoPtr<FEBase> ap(new FEXYZ<3>(fet));
return ap;
}
case SCALAR:
{
AutoPtr<FEBase> ap(new FEScalar<3>(fet));
return ap;
}
default:
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
default:
libmesh_error();
}
libmesh_error();
AutoPtr<FEBase> ap(NULL);
return ap;
}
template<>
AutoPtr< FEGenericBase< RealGradient > > FEGenericBase< RealGradient >::build ( const unsigned int  dim,
const FEType &  fet 
)

Definition at line 522 of file fe_base.C.

References libMesh::FEType::family, libMeshEnums::LAGRANGE_VEC, libMeshEnums::NEDELEC_ONE, and libMesh::out.

{
// The stupid AutoPtr<FEBase> ap(); return ap;
// construct is required to satisfy IBM's xlC
switch (dim)
{
// 0D
case 0:
{
switch (fet.family)
{
case LAGRANGE_VEC:
{
AutoPtr<FEVectorBase> ap( new FELagrangeVec<0>(fet) );
return ap;
}
default:
{
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
}
case 1:
{
switch (fet.family)
{
case LAGRANGE_VEC:
{
AutoPtr<FEVectorBase> ap( new FELagrangeVec<1>(fet) );
return ap;
}
default:
{
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
}
case 2:
{
switch (fet.family)
{
case LAGRANGE_VEC:
{
AutoPtr<FEVectorBase> ap( new FELagrangeVec<2>(fet) );
return ap;
}
case NEDELEC_ONE:
{
AutoPtr<FEVectorBase> ap( new FENedelecOne<2>(fet) );
return ap;
}
default:
{
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
}
case 3:
{
switch (fet.family)
{
case LAGRANGE_VEC:
{
AutoPtr<FEVectorBase> ap( new FELagrangeVec<3>(fet) );
return ap;
}
case NEDELEC_ONE:
{
AutoPtr<FEVectorBase> ap( new FENedelecOne<3>(fet) );
return ap;
}
default:
{
libMesh::out << "ERROR: Bad FEType.family= " << fet.family << std::endl;
libmesh_error();
}
}
}
default:
libmesh_error();
} // switch(dim)
libmesh_error();
AutoPtr<FEVectorBase> ap(NULL);
return ap;
}
template<>
AutoPtr< FEGenericBase< Real > > FEGenericBase< Real >::build_InfFE ( const unsigned int  dim,
const FEType &  fet 
)

Definition at line 627 of file fe_base.C.

References libMeshEnums::CARTESIAN, libMesh::err, libMesh::FEType::inf_map, libMeshEnums::INFINITE_MAP, libMeshEnums::JACOBI_20_00, libMeshEnums::JACOBI_30_00, libMeshEnums::LAGRANGE, libMeshEnums::LEGENDRE, and libMesh::FEType::radial_family.

{
// The stupid AutoPtr<FEBase> ap(); return ap;
// construct is required to satisfy IBM's xlC
switch (dim)
{
// 1D
case 1:
{
switch (fet.radial_family)
{
case INFINITE_MAP:
{
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with FEFamily = " << fet.radial_family << std::endl;
libmesh_error();
}
case JACOBI_20_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<1,JACOBI_20_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case JACOBI_30_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<1,JACOBI_30_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LEGENDRE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<1,LEGENDRE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LAGRANGE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<1,LAGRANGE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
default:
libMesh::err << "ERROR: Bad FEType.radial_family= " << fet.radial_family << std::endl;
libmesh_error();
}
}
// 2D
case 2:
{
switch (fet.radial_family)
{
case INFINITE_MAP:
{
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with FEFamily = " << fet.radial_family << std::endl;
libmesh_error();
}
case JACOBI_20_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<2,JACOBI_20_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case JACOBI_30_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<2,JACOBI_30_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LEGENDRE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<2,LEGENDRE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LAGRANGE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<2,LAGRANGE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
default:
libMesh::err << "ERROR: Bad FEType.radial_family= " << fet.radial_family << std::endl;
libmesh_error();
}
}
// 3D
case 3:
{
switch (fet.radial_family)
{
case INFINITE_MAP:
{
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with FEFamily = " << fet.radial_family << std::endl;
libmesh_error();
}
case JACOBI_20_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<3,JACOBI_20_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case JACOBI_30_00:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<3,JACOBI_30_00,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LEGENDRE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<3,LEGENDRE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
case LAGRANGE:
{
switch (fet.inf_map)
{
case CARTESIAN:
{
AutoPtr<FEBase> ap(new InfFE<3,LAGRANGE,CARTESIAN>(fet));
return ap;
}
default:
libMesh::err << "ERROR: Don't build an infinite element " << std::endl
<< " with InfMapType = " << fet.inf_map << std::endl;
libmesh_error();
}
}
default:
libMesh::err << "ERROR: Bad FEType.radial_family= " << fet.radial_family << std::endl;
libmesh_error();
}
}
default:
libmesh_error();
}
libmesh_error();
AutoPtr<FEBase> ap(NULL);
return ap;
}
template<>
AutoPtr< FEGenericBase< RealGradient > > FEGenericBase< RealGradient >::build_InfFE ( const unsigned  int,
const FEType &   
)

Definition at line 909 of file fe_base.C.

{
// No vector types defined... YET.
libmesh_error();
AutoPtr<FEVectorBase> ap(NULL);
return ap;
}
The documentation for this class was generated from the following file:
Site Created By: libMesh Developers
Last modified: January 12 2013 15:00:43 UTC

Hosted By:
SourceForge.net Logo