|
Intrepid
|
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., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 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
1.7.6.1