00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #include "SundanceBasisDOFTopologyBase.hpp"
00032 #include "PlayaExceptions.hpp"
00033
00034
00035 using namespace Sundance;
00036 using namespace Teuchos;
00037
00038
00039 int BasisDOFTopologyBase::nReferenceDOFsWithFacets(
00040 const CellType& maximalCellType,
00041 const CellType& cellType
00042 ) const
00043 {
00044 switch(cellType)
00045 {
00046 case TetCell:
00047 return nReferenceDOFsWithoutFacets(maximalCellType, TetCell)
00048 + 4*nReferenceDOFsWithoutFacets(maximalCellType, TriangleCell)
00049 + 6*nReferenceDOFsWithoutFacets(maximalCellType, LineCell)
00050 + 4*nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00051 case BrickCell:
00052 return nReferenceDOFsWithoutFacets(maximalCellType, BrickCell)
00053 + 6*nReferenceDOFsWithoutFacets(maximalCellType, QuadCell)
00054 + 12*nReferenceDOFsWithoutFacets(maximalCellType, LineCell)
00055 + 8*nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00056 case TriangleCell:
00057 return nReferenceDOFsWithoutFacets(maximalCellType, TriangleCell)
00058 + 3*nReferenceDOFsWithoutFacets(maximalCellType, LineCell)
00059 + 3*nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00060 case QuadCell:
00061 return nReferenceDOFsWithoutFacets(maximalCellType, QuadCell)
00062 + 4*nReferenceDOFsWithoutFacets(maximalCellType, LineCell)
00063 + 4*nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00064 case LineCell:
00065 return nReferenceDOFsWithoutFacets(maximalCellType, LineCell)
00066 + 2*nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00067 case PointCell:
00068 return nReferenceDOFsWithoutFacets(maximalCellType, PointCell);
00069 default:
00070 TEUCHOS_TEST_FOR_EXCEPTION(true, std::runtime_error,
00071 "case cellType=" << cellType << " not defined");
00072 }
00073 return -1;
00074 }