00001 #include "fedealwrapperq1.h" 00002 00003 00004 FEDealWrapperQ1::FEDealWrapperQ1() 00005 :m_fe(1) 00006 { 00007 wrap(m_fe); 00008 m_global_map.resize(m_fe.n_dofs_per_cell()); 00009 } 00010 00011 00012 00013 00014 00015 FEDealWrapperQ1::~FEDealWrapperQ1() 00016 { 00017 00018 } 00019 00020 00021 const VecIndex& FEDealWrapperQ1::get_global_map(OrthoMesh::Cell_It &cell) 00022 { 00023 for (unsigned i=0;i<m_global_map.size();i++) 00024 m_global_map[i] = cell->vertex_index(static_cast<VertexDirection3D>(i)); 00025 return m_global_map; 00026 }