CuboidMeshConstructor.cpp
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 #include "CuboidMeshConstructor.hpp"
00031 #include "UblasCustomFunctions.hpp"
00032
00033
00035
00037
00038
00039 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00040 void CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::ConstructHyperCube(AbstractTetrahedralMesh<1,SPACE_DIM>& rMesh, unsigned width)
00041 {
00042 rMesh.ConstructLinearMesh(width);
00043 }
00044
00045 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00046 void CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::ConstructHyperCube(AbstractTetrahedralMesh<2,2>& rMesh, unsigned width)
00047 {
00048 rMesh.ConstructRectangularMesh(width, width);
00049 }
00050
00051 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00052 void CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::ConstructHyperCube(AbstractTetrahedralMesh<3,3>& rMesh, unsigned width)
00053 {
00054 rMesh.ConstructCuboid(width, width, width);
00055 }
00056
00057 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00058 void CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::Construct(AbstractTetrahedralMesh<ELEMENT_DIM,SPACE_DIM>& rMesh, unsigned meshRefinementNum, double meshWidth)
00059 {
00060
00061 assert(ELEMENT_DIM == SPACE_DIM || (ELEMENT_DIM == 1 && SPACE_DIM == 3));
00062
00063 mMeshWidth = meshWidth;
00064 assert(meshRefinementNum < 30);
00065
00066
00067 unsigned mesh_size = (unsigned) SmallPow(2, meshRefinementNum+2);
00068 double scaling = mMeshWidth/(double) mesh_size;
00069 ConstructHyperCube(rMesh, mesh_size);
00070 rMesh.Scale(scaling, scaling, scaling);
00071 mNumElements = rMesh.GetNumElements();
00072 mNumNodes = rMesh.GetNumNodes();
00073 }
00074
00075 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00076 double CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::GetWidth()
00077 {
00078 return mMeshWidth;
00079 }
00080
00081 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00082 unsigned CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::GetNumElements()
00083 {
00084 return mNumElements;
00085 }
00086
00087 template<unsigned ELEMENT_DIM, unsigned SPACE_DIM>
00088 unsigned CuboidMeshConstructor<ELEMENT_DIM, SPACE_DIM>::GetNumNodes()
00089 {
00090 return mNumNodes;
00091 }
00092
00093
00094
00096
00098
00099 template class CuboidMeshConstructor<1>;
00100 template class CuboidMeshConstructor<2>;
00101 template class CuboidMeshConstructor<3>;
00102 template class CuboidMeshConstructor<1,3>;