SundanceEdgeLocalizedBasis.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 #ifndef SUNDANCE_EDGELOCALIZEDBASIS_H
00043 #define SUNDANCE_EDGELOCALIZEDBASIS_H
00044 
00045 #include "SundanceDefs.hpp"
00046 #include "Teuchos_RefCountPtr.hpp"
00047 #include "SundanceBasisFamilyBase.hpp"
00048 
00049 namespace Sundance 
00050 {
00051 /** 
00052  * EdgeLocalizedBasis is a basis function defined only on edges. 
00053  */
00054 class EdgeLocalizedBasis : public ScalarBasis
00055 {
00056 public:
00057   /** */
00058   EdgeLocalizedBasis();
00059 
00060   /**   
00061    * \brief Inform caller as to whether a given cell type is supported 
00062    */
00063   bool supportsCellTypePair(
00064     const CellType& maximalCellType,
00065     const CellType& cellType
00066     ) const ;
00067 
00068   /** */
00069   void print(std::ostream& os) const ;
00070 
00071   /** */
00072   int order() const {return 0;}
00073 
00074   /** return the number of nodes for this basis on the given cell type */
00075   int nReferenceDOFsWithoutFacets(
00076     const CellType& maximalCellType,
00077     const CellType& cellType
00078     ) const ;
00079 
00080   /** */
00081   void getReferenceDOFs(
00082     const CellType& maximalCellType,
00083     const CellType& cellType,
00084     Array<Array<Array<int> > >& dofs) const ;
00085 
00086   /** */
00087   void refEval(
00088     const CellType& cellType,
00089     const Array<Point>& pts,
00090     const SpatialDerivSpecifier& deriv,
00091     Array<Array<Array<double> > >& result,
00092     int verbosity=0) const ;
00093 
00094 
00095   /* Handleable boilerplate */
00096   GET_RCP(BasisFamilyBase);
00097 
00098 private:
00099 
00100   /** evaluate on a line cell  */
00101   void evalOnLine(const Point& pt,
00102     const MultiIndex& deriv,
00103     Array<double>& result) const ;
00104     
00105   /** evaluate on a triangle cell  */
00106   void evalOnTriangle(const Point& pt,
00107     const MultiIndex& deriv,
00108     Array<double>& result) const ;
00109   /** evaluate on a tet cell  */
00110   void evalOnTet(const Point& pt,
00111     const MultiIndex& deriv,
00112     Array<double>& result) const ;
00113     
00114 };
00115 }
00116 
00117 #endif

Site Contact