Intrepid
Intrepid_HGRAD_TRI_Cn_FEM.hpp
Go to the documentation of this file.
1 // @HEADER
2 // ************************************************************************
3 //
4 // Intrepid Package
5 // Copyright (2007) Sandia Corporation
6 //
7 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive
8 // license for use of this work by or on behalf of the U.S. Government.
9 //
10 // Redistribution and use in source and binary forms, with or without
11 // modification, are permitted provided that the following conditions are
12 // met:
13 //
14 // 1. Redistributions of source code must retain the above copyright
15 // notice, this list of conditions and the following disclaimer.
16 //
17 // 2. Redistributions in binary form must reproduce the above copyright
18 // notice, this list of conditions and the following disclaimer in the
19 // documentation and/or other materials provided with the distribution.
20 //
21 // 3. Neither the name of the Corporation nor the names of the
22 // contributors may be used to endorse or promote products derived from
23 // this software without specific prior written permission.
24 //
25 // THIS SOFTWARE IS PROVIDED BY SANDIA CORPORATION "AS IS" AND ANY
26 // EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
27 // IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR
28 // PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL SANDIA CORPORATION OR THE
29 // CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
30 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
31 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
32 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
33 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
34 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
35 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
36 //
37 // Questions? Contact Pavel Bochev (pbboche@sandia.gov)
38 // Denis Ridzal (dridzal@sandia.gov), or
39 // Kara Peterson (kjpeter@sandia.gov)
40 //
41 // ************************************************************************
42 // @HEADER
43 
49 #ifndef INTREPID_HGRAD_TRI_Cn_FEM_HPP
50 #define INTREPID_HGRAD_TRI_Cn_FEM_HPP
51 
52 #include "Intrepid_Types.hpp"
53 #include "Intrepid_Basis.hpp"
54 #include "Intrepid_PointTools.hpp"
55 #include "Intrepid_HGRAD_TRI_Cn_FEM_ORTH.hpp"
56 #include "Teuchos_SerialDenseMatrix.hpp"
57 #include "Teuchos_SerialDenseSolver.hpp"
58 
59 
60 namespace Intrepid {
61 
83 template<class Scalar, class ArrayScalar>
84 class Basis_HGRAD_TRI_Cn_FEM: public Basis<Scalar, ArrayScalar> {
85 private:
86 
89  virtual void initializeTags();
90 
104 
105 public:
106 
109  Basis_HGRAD_TRI_Cn_FEM( const int n , const EPointType pointType );
110 
121  void getValues(ArrayScalar & outputValues,
122  const ArrayScalar & inputPoints,
123  const EOperator operatorType) const;
124 
125 
128  void getValues(ArrayScalar & outputValues,
129  const ArrayScalar & inputPoints,
130  const ArrayScalar & cellVertices,
131  const EOperator operatorType = OPERATOR_VALUE) const;
132 };
133 
134 }// namespace Intrepid
135 
136 #include "Intrepid_HGRAD_TRI_Cn_FEMDef.hpp"
137 
138 #endif
Implementation of the default H(grad)-compatible orthogonal basis (Dubiner) of arbitrary degree on tr...
FieldContainer< Scalar > latticePts
stores the points at which degrees of freedom are located.
Header file for utility class to provide point tools, such as barycentric coordinates, equispaced lattices, and warp-blend point distrubtions.
Basis_HGRAD_TRI_Cn_FEM(const int n, const EPointType pointType)
Constructor.
Contains definitions of custom data types in Intrepid.
FieldContainer< Scalar > V
The Vandermonde matrix with V_{ij} = phi_i(x_j), where x_j is the j_th point in the lattice...
Header file for the abstract base class Intrepid::Basis.
FieldContainer< Scalar > Vinv
The inverse of V. The columns of Vinv express the Lagrange basis in terms of the orthogonal basis...
An abstract base class that defines interface for concrete basis implementations for Finite Element (...
void getValues(ArrayScalar &outputValues, const ArrayScalar &inputPoints, const EOperator operatorType) const
Evaluation of a FEM basis on a reference Triangle cell.
Implementation of the default H(grad)-compatible Lagrange basis of arbitrary degree on Triangle cell...
virtual void initializeTags()
Initializes tagToOrdinal_ and ordinalToTag_ lookup arrays.
Basis_HGRAD_TRI_Cn_FEM_ORTH< Scalar, FieldContainer< Scalar > > Phis
The orthogonal basis on triangles, out of which the nodal basis is constructed.