Intrepid
Intrepid_HGRAD_LINE_Cn_FEM.hpp
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_LINE_Cn_FEM_HPP
50#define INTREPID_HGRAD_LINE_Cn_FEM_HPP
51
52#include "Intrepid_Basis.hpp"
54#include "Intrepid_Polylib.hpp"
56#include "Teuchos_SerialDenseMatrix.hpp"
57#include "Teuchos_SerialDenseSolver.hpp"
58
59
60namespace Intrepid {
61
77template<class Scalar, class ArrayScalar>
79 public Basis<Scalar, ArrayScalar>,
80 public DofCoordsInterface<ArrayScalar> {
81private:
83 FieldContainer<Scalar> latticePts_;
84
86 Basis_HGRAD_LINE_Cn_FEM_JACOBI<Scalar, FieldContainer<Scalar> > Phis_;
87
89 FieldContainer<Scalar> V_;
90
93 FieldContainer<Scalar> Vinv_;
94
97 void initializeTags();
98
99public:
103
106 Basis_HGRAD_LINE_Cn_FEM(int order , const ArrayScalar &pts );
107
111 Basis_HGRAD_LINE_Cn_FEM(int order , const EPointType &pointType );
112
113
124 void getValues(ArrayScalar & outputValues,
125 const ArrayScalar & inputPoints,
126 const EOperator operatorType) const;
127
128
131 void getValues(ArrayScalar & outputValues,
132 const ArrayScalar & inputPoints,
133 const ArrayScalar & cellVertices,
134 const EOperator operatorType = OPERATOR_VALUE) const;
135
137 virtual void getDofCoords( ArrayScalar & DofCoords ) const;
138
139};
140
141}// namespace Intrepid
142
143#include "Intrepid_HGRAD_LINE_Cn_FEMDef.hpp"
144
145#endif
Header file for the abstract base class Intrepid::Basis.
Header file for the Intrepid::HGRAD_LINE_Cn_FEM class.
Header file for utility class to provide point tools, such as barycentric coordinates,...
Header file for a set of functions providing orthogonal polynomial polynomial calculus and interpolat...
Implementation of the locally H(grad)-compatible FEM basis of variable order on the [-1,...
void initializeTags()
Initializes tagToOrdinal_ and ordinalToTag_ lookup arrays.
Basis_HGRAD_LINE_Cn_FEM(int order, const ArrayScalar &pts)
Constructor.
FieldContainer< Scalar > latticePts_
Holds the points defining the Lagrange basis.
FieldContainer< Scalar > Vinv_
inverse of Generalized Vandermonde matrix, whose columns store the expansion coefficients of the noda...
void getValues(ArrayScalar &outputValues, const ArrayScalar &inputPoints, const EOperator operatorType) const
Evaluation of a FEM basis on a reference Line cell.
FieldContainer< Scalar > V_
Generalized Vandermonde matrix V_{ij} = phis_i(x_j)
virtual void getDofCoords(ArrayScalar &DofCoords) const
implements the dofcoords interface
Basis_HGRAD_LINE_Cn_FEM_JACOBI< Scalar, FieldContainer< Scalar > > Phis_
orthogonal basis
An abstract base class that defines interface for concrete basis implementations for Finite Element (...
This is an interface class for bases whose degrees of freedom can be associated with spatial location...