Intrepid
/usr/src/RPM/BUILD/trilinos10-10.6.4/packages/intrepid/src/Discretization/Basis/Intrepid_HGRAD_LINE_Cn_FEM.hpp
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_LINE_Cn_FEM_HPP
00036 #define INTREPID_HGRAD_LINE_Cn_FEM_HPP
00037 
00038 #include "Intrepid_Basis.hpp"
00039 #include "Intrepid_PointTools.hpp"
00040 #include "Intrepid_Polylib.hpp"
00041 #include "Intrepid_HGRAD_LINE_Cn_FEM_JACOBI.hpp"
00042 #include "Teuchos_SerialDenseMatrix.hpp"
00043 #include "Teuchos_SerialDenseSolver.hpp"
00044 
00045 
00046 namespace Intrepid {
00047   
00063 template<class Scalar, class ArrayScalar> 
00064 class Basis_HGRAD_LINE_Cn_FEM: public Basis<Scalar, ArrayScalar> {
00065 private:
00067   FieldContainer<Scalar> latticePts_;
00068 
00070   Basis_HGRAD_LINE_Cn_FEM_JACOBI<Scalar, FieldContainer<Scalar> > Phis_;
00071   
00073   FieldContainer<Scalar> V_;
00074 
00077   FieldContainer<Scalar> Vinv_;
00078 
00081   void initializeTags();
00082   
00083 public:
00084   
00087   Basis_HGRAD_LINE_Cn_FEM(int order , const ArrayScalar &pts );  
00088 
00092   Basis_HGRAD_LINE_Cn_FEM(int order , const EPointType &pointType );  
00093   
00094   
00105   void getValues(ArrayScalar &          outputValues,
00106                  const ArrayScalar &    inputPoints,
00107                  const EOperator        operatorType) const;
00108   
00109   
00112   void getValues(ArrayScalar &          outputValues,
00113                  const ArrayScalar &    inputPoints,
00114                  const ArrayScalar &    cellVertices,
00115                  const EOperator        operatorType = OPERATOR_VALUE) const;
00116 
00117 };
00118 
00119 }// namespace Intrepid
00120 
00121 #include "Intrepid_HGRAD_LINE_Cn_FEMDef.hpp"
00122 
00123 #endif