Intrepid
http://trilinos.sandia.gov/packages/docs/r10.12/packages/intrepid/src/Discretization/Basis/Intrepid_HGRAD_POLY_C1_FEM.hpp
Go to the documentation of this file.
00001 // @HEADER
00002 // ************************************************************************
00003 //
00004 //                           Intrepid Package
00005 //                 Copyright (2007) Sandia Corporation
00006 //
00007 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
00008 // license for use of this work by or on behalf of the U.S. Government.
00009 //
00010 // This library is free software; you can redistribute it and/or modify
00011 // it under the terms of the GNU Lesser General Public License as
00012 // published by the Free Software Foundation; either version 2.1 of the
00013 // License, or (at your option) any later version.
00014 //
00015 // This library is distributed in the hope that it will be useful, but
00016 // WITHOUT ANY WARRANTY; without even the implied warranty of
00017 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00018 // Lesser General Public License for more details.
00019 //
00020 // You should have received a copy of the GNU Lesser General Public
00021 // License along with this library; if not, write to the Free Software
00022 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
00023 // USA
00024 // Questions? Contact Pavel Bochev (pbboche@sandia.gov) or
00025 //                    Denis Ridzal (dridzal@sandia.gov).
00026 //
00027 // ************************************************************************
00028 // @HEADER
00029 
00035 #ifndef INTREPID_HGRAD_POLY_C1_FEM_HPP
00036 #define INTREPID_HGRAD_POLY_C1_FEM_HPP
00037 
00038 #include "Intrepid_Basis.hpp"
00039 #include "Shards_CellTopology.hpp"
00040 
00041 namespace Intrepid{
00048   template<class Scalar, class ArrayScalar> 
00049   class Basis_HGRAD_POLY_C1_FEM : public Basis<Scalar, ArrayScalar> {
00050   public:
00055     Basis_HGRAD_POLY_C1_FEM(const shards::CellTopology& cellTopology);
00056     
00059     void getValues(ArrayScalar& outputValues,
00060                    const ArrayScalar& inputPoints,
00061                    const EOperator operatorType) const;
00062 
00070     void getValues(ArrayScalar& outputValues,
00071                    const ArrayScalar& inputPoints,
00072                    const ArrayScalar& cellVertices,
00073                    const EOperator operatorType = OPERATOR_VALUE) const;
00074     
00075   private:
00076     
00079     void initializeTags();
00080 
00081     
00084     template<class Scalar1, class ArrayScalar1>
00085     Scalar1 computeArea(const ArrayScalar1& p1,
00086                        const ArrayScalar1& p2,
00087                        const ArrayScalar1& p3) const;
00088 
00091     template<class Scalar1, class ArrayScalar1>
00092     void evaluateWeightFunctions(ArrayScalar1& outputValues,
00093                                  const ArrayScalar1& inputValues,
00094                                  const ArrayScalar1& cellVertices) const;
00095     
00096     
00097     
00100     template<class Scalar1, class ArrayScalar1>
00101     void shapeFunctions(ArrayScalar1& outputValues,
00102                         const ArrayScalar1& inputValues,
00103                         const ArrayScalar1& cellVertices)const;
00104   }; // end class Basis_HGRAD_POLY_C1_FEM
00105 
00106 } // namespace Intrepid
00107 
00108 #include "Intrepid_HGRAD_POLY_C1_FEMDef.hpp"
00109 #endif