00001 /* @HEADER@ */ 00002 // ************************************************************************ 00003 // 00004 // Sundance 00005 // Copyright (2005) Sandia Corporation 00006 // 00007 // Copyright (year first published) Sandia Corporation. Under the terms 00008 // of Contract DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government 00009 // retains certain rights in this software. 00010 // 00011 // This library is free software; you can redistribute it and/or modify 00012 // it under the terms of the GNU Lesser General Public License as 00013 // published by the Free Software Foundation; either version 2.1 of the 00014 // License, or (at your option) any later version. 00015 // 00016 // This library is distributed in the hope that it will be useful, but 00017 // WITHOUT ANY WARRANTY; without even the implied warranty of 00018 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00019 // Lesser General Public License for more details. 00020 // 00021 // You should have received a copy of the GNU Lesser General Public 00022 // License along with this library; if not, write to the Free Software 00023 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00024 // USA 00025 // Questions? Contact Kevin Long (krlong@sandia.gov), 00026 // Sandia National Laboratories, Livermore, California, USA 00027 // 00028 // ************************************************************************ 00029 /* @HEADER@ */ 00030 00031 #ifndef SUNDANCE_PARTIALELEMENTDOFMAP_H 00032 #define SUNDANCE_PARTIALELEMENTDOFMAP_H 00033 00034 #include "SundanceDefs.hpp" 00035 #include "SundanceDOFMapBase.hpp" 00036 #include "SundanceCellSet.hpp" 00037 #include "SundanceCellFilter.hpp" 00038 #include "SundanceBasisFamily.hpp" 00039 #include "SundanceObjectWithVerbosity.hpp" 00040 00041 00042 namespace Sundance 00043 { 00044 using namespace Teuchos; 00045 00046 /** 00047 * PartialElementDOFMap is a DOF map specialized to the case of element-based 00048 * DOFs on a subset of cells in the domain. All elements must have the same 00049 * set of functions. 00050 */ 00051 class PartialElementDOFMap : public DOFMapBase 00052 { 00053 public: 00054 /** */ 00055 PartialElementDOFMap(const Mesh& mesh, 00056 const CellFilter& subdomain, 00057 int nFuncs, 00058 int setupVerb); 00059 00060 /** */ 00061 virtual ~PartialElementDOFMap(){;} 00062 00063 /** */ 00064 RCP<const MapStructure> 00065 getDOFsForCellBatch(int cellDim, 00066 const Array<int>& cellLID, 00067 const Set<int>& requestedFuncSet, 00068 Array<Array<int> >& dofs, 00069 Array<int>& nNodes, 00070 int verbosity) const ; 00071 00072 /** */ 00073 RCP<const Set<int> > 00074 allowedFuncsOnCellBatch(int cellDim, 00075 const Array<int>& cellLID) const ; 00076 00077 /** */ 00078 const Array<CellFilter>& funcDomains() const {return funcDomains_;} 00079 00080 /** */ 00081 virtual void print(std::ostream& os) const ; 00082 00083 00084 protected: 00085 00086 void init(); 00087 00088 void computeOffsets(int localCount) ; 00089 00090 void shareRemoteDOFs(const Array<Array<int> >& remoteElems); 00091 00092 int dim_; 00093 int nFuncs_; 00094 int nElems_; 00095 CellFilter subdomain_; 00096 Array<CellFilter> funcDomains_; 00097 Array<int> elemDofs_; 00098 RCP<MapStructure> structure_; 00099 RCP<const Set<int> > allFuncs_; 00100 }; 00101 00102 } 00103 00104 00105 #endif