SundanceLocalDOFMap.hpp
Go to the documentation of this file.
00001 /* @HEADER@ */
00002 // ************************************************************************
00003 // 
00004 //                             Sundance
00005 //                 Copyright 2011 Sandia Corporation
00006 // 
00007 // Under the terms of Contract DE-AC04-94AL85000 with Sandia Corporation,
00008 // the U.S. Government retains certain rights in this software.
00009 //
00010 // Redistribution and use in source and binary forms, with or without
00011 // modification, are permitted provided that the following conditions are
00012 // met:
00013 //
00014 // 1. Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 //
00017 // 2. Redistributions in binary form must reproduce the above copyright
00018 // notice, this list of conditions and the following disclaimer in the
00019 // documentation and/or other materials provided with the distribution.
00020 //
00021 // 3. Neither the name of the Corporation nor the names of the
00022 // contributors may be used to endorse or promote products derived from
00023 // this software without specific prior written permission.
00024 //
00025 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
00026 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00027 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
00028 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
00029 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00030 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00031 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00032 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00033 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00034 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00035 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00036 //
00037 // Questions? Contact Kevin Long (kevin.long@ttu.edu)
00038 // 
00039 
00040 /* @HEADER@ */
00041 
00042 
00043 
00044 #ifndef SUNDANCE_LOCALDOFMAP_H
00045 #define SUNDANCE_LOCALDOFMAP_H
00046 
00047 #include "SundanceDefs.hpp"
00048 #include "SundanceDOFMapBase.hpp"
00049 
00050 
00051 namespace Sundance
00052 {
00053 
00054 /** 
00055  * LocalDOFMap bundles several compact tables used for fast lookup of
00056  * local DOFs. 
00057  */
00058 class LocalDOFMap
00059 {
00060 public:
00061   /** */
00062   LocalDOFMap(int numBlocks, int verb);
00063 
00064   /** */
00065   void setCells(int cellDim, int maxCellDim,
00066     const RCP<const Array<int> >& cellLID);
00067 
00068   /** */
00069   int nCells() const ;
00070 
00071   /** */
00072   bool isUsed(int b) const {return isUsed_[b];}
00073 
00074   /** */
00075   bool isUnused(int b) const {return !isUsed(b);}
00076 
00077   /** */
00078   bool isUnused() const ;
00079 
00080   /** */
00081   void markAsUnused() ;
00082 
00083   /** */
00084   bool hasCells() const {return hasCells_;}
00085 
00086   /** */
00087   const RCP<const Array<int> >& cellLIDs() const {return cellLID_;}
00088 
00089   /** */
00090   void markAsUsed(int b) {isUsed_[b] = true;}
00091 
00092   /** */
00093   int numBlocks() const {return mapStruct_->size();}
00094 
00095   /** */
00096   const Array<int>& nLocalNodesPerChunk(int b) const 
00097     {return (*nLocalNodesPerChunk_)[b];}
00098 
00099   /** */
00100   const RCP<const MapStructure>& mapStruct(int b) const 
00101     {return (*mapStruct_)[b];}
00102 
00103   /** */
00104   const Array<Array<int> >& localDOFs(int b) const 
00105     {return (*localDOFs_)[b];}
00106 
00107   /** */
00108   std::ostream& print(std::ostream& os) const ;
00109 
00110   /** */
00111   void fillBlock(int b, const RCP<DOFMapBase>& globalMap,
00112     const Array<Set<int> >& requiredFunc);
00113 
00114   /** */
00115   void setVerb(int v) const {verb_ = v;}
00116 
00117 private:
00118 
00119 
00120   /** */
00121   Array<int>& nLocalNodesPerChunk(int b)
00122     {return (*nLocalNodesPerChunk_)[b];}
00123 
00124   /** */
00125   RCP<const MapStructure>& mapStruct(int b) 
00126     {return (*mapStruct_)[b];}
00127 
00128   /** */
00129   Array<Array<int> >& localDOFs(int b) 
00130     {return (*localDOFs_)[b];}
00131 
00132   /** */
00133   void verifyValidBlock(int b) const ;
00134 
00135 
00136   mutable int verb_;
00137   Array<int> isUsed_;
00138   bool hasCells_;
00139   RCP<Array<Array<int> > > nLocalNodesPerChunk_;
00140   RCP<Array<RCP<const MapStructure> > > mapStruct_;  
00141   RCP<Array<Array<Array<int> > > > localDOFs_;
00142   RCP<const Array<int> > cellLID_;
00143   int activeCellDim_;
00144   int maxCellDim_;
00145 };
00146 
00147 
00148 /** \relates LocalDOFMap */
00149 inline std::ostream& operator<<(std::ostream& os,
00150   const LocalDOFMap& m)
00151 {
00152   return m.print(os);
00153 }
00154 
00155 }
00156 
00157 
00158 
00159 #endif

Site Contact