Intrepid
Intrepid_HGRAD_TET_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_TET_Cn_FEM_HPP
50#define INTREPID_HGRAD_TET_Cn_FEM_HPP
51
52#include "Intrepid_Types.hpp"
53#include "Intrepid_Basis.hpp"
55#include "Intrepid_HGRAD_TET_Cn_FEM_ORTH.hpp"
56#include "Teuchos_SerialDenseMatrix.hpp"
57#include "Teuchos_SerialDenseSolver.hpp"
58
59
60namespace Intrepid {
61
87template<class Scalar, class ArrayScalar>
88class Basis_HGRAD_TET_Cn_FEM: public Basis<Scalar, ArrayScalar> {
89private:
90
93 virtual void initializeTags();
94
97 Basis_HGRAD_TET_Cn_FEM_ORTH<Scalar,FieldContainer<Scalar> > Phis;
100 FieldContainer<Scalar> V;
104 FieldContainer<Scalar> Vinv;
107 FieldContainer<Scalar> latticePts;
108
109public:
110
113 Basis_HGRAD_TET_Cn_FEM( const int n , const EPointType pointType );
114
125 void getValues(ArrayScalar & outputValues,
126 const ArrayScalar & inputPoints,
127 const EOperator operatorType) const;
128
129
132 void getValues(ArrayScalar & outputValues,
133 const ArrayScalar & inputPoints,
134 const ArrayScalar & cellVertices,
135 const EOperator operatorType = OPERATOR_VALUE) const;
136};
137
138}// namespace Intrepid
139
140#include "Intrepid_HGRAD_TET_Cn_FEMDef.hpp"
141
142#endif
Header file for the abstract base class Intrepid::Basis.
Header file for utility class to provide point tools, such as barycentric coordinates,...
Contains definitions of custom data types in Intrepid.
Implementation of the default H(grad)-compatible Lagrange basis of arbitrary degree on Tetrahedron ce...
Basis_HGRAD_TET_Cn_FEM(const int n, const EPointType pointType)
Constructor.
FieldContainer< Scalar > latticePts
stores the points at which degrees of freedom are located.
FieldContainer< Scalar > Vinv
The inverse of V. The columns of Vinv express the Lagrange basis in terms of the orthogonal basis.
Basis_HGRAD_TET_Cn_FEM_ORTH< Scalar, FieldContainer< Scalar > > Phis
The orthogonal basis on triangles, out of which the nodal basis is constructed.
virtual void initializeTags()
Initializes tagToOrdinal_ and ordinalToTag_ lookup arrays.
void getValues(ArrayScalar &outputValues, const ArrayScalar &inputPoints, const EOperator operatorType) const
Evaluation of a FEM basis on a reference Triangle cell.
FieldContainer< Scalar > V
The Vandermonde matrix with V_{ij} = phi_i(x_j), where x_j is the j_th point in the lattice.
An abstract base class that defines interface for concrete basis implementations for Finite Element (...