Intrepid
Intrepid_HDIV_TRI_I1_FEMDef.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 namespace Intrepid {
50 
51 template<class Scalar, class ArrayScalar>
52 Basis_HDIV_TRI_I1_FEM<Scalar,ArrayScalar>::Basis_HDIV_TRI_I1_FEM()
53  {
54  this -> basisCardinality_ = 3;
55  this -> basisDegree_ = 1;
56  this -> basisCellTopology_ = shards::CellTopology(shards::getCellTopologyData<shards::Triangle<3> >() );
57  this -> basisType_ = BASIS_FEM_DEFAULT;
58  this -> basisCoordinates_ = COORDINATES_CARTESIAN;
59  this -> basisTagsAreSet_ = false;
60  }
61 
62 template<class Scalar, class ArrayScalar>
63 void Basis_HDIV_TRI_I1_FEM<Scalar, ArrayScalar>::initializeTags() {
64 
65  // Basis-dependent intializations
66  int tagSize = 4; // size of DoF tag
67  int posScDim = 0; // position in the tag, counting from 0, of the subcell dim
68  int posScOrd = 1; // position in the tag, counting from 0, of the subcell ordinal
69  int posDfOrd = 2; // position in the tag, counting from 0, of DoF ordinal relative to the subcell
70 
71  // An array with local DoF tags assigned to basis functions, in the order of their local enumeration
72  int tags[] = {
73  1, 0, 0, 1,
74  1, 1, 0, 1,
75  1, 2, 0, 1
76  // 2, 0, 0, 1,
77  // 2, 1, 0, 1,
78  // 2, 2, 0, 1
79  };
80 
81  // Basis-independent function sets tag and enum data in tagToOrdinal_ and ordinalToTag_ arrays:
82  Intrepid::setOrdinalTagData(this -> tagToOrdinal_,
83  this -> ordinalToTag_,
84  tags,
85  this -> basisCardinality_,
86  tagSize,
87  posScDim,
88  posScOrd,
89  posDfOrd);
90 }
91 
92 
93 
94 template<class Scalar, class ArrayScalar>
95 void Basis_HDIV_TRI_I1_FEM<Scalar, ArrayScalar>::getValues(ArrayScalar & outputValues,
96  const ArrayScalar & inputPoints,
97  const EOperator operatorType) const {
98 
99 
100 // Verify arguments
101 #ifdef HAVE_INTREPID_DEBUG
102  Intrepid::getValues_HDIV_Args<Scalar, ArrayScalar>(outputValues,
103  inputPoints,
104  operatorType,
105  this -> getBaseCellTopology(),
106  this -> getCardinality() );
107 #endif
108 
109 
110  // Number of evaluation points = dim 0 of inputPoints
111  int dim0 = inputPoints.dimension(0);
112 
113  // Temporaries: (x,y) coordinates of the evaluation point
114  Scalar x = 0.0;
115  Scalar y = 0.0;
116 
117  switch (operatorType) {
118  case OPERATOR_VALUE:
119  for (int i0 = 0; i0 < dim0; i0++) {
120  x = inputPoints(i0, 0);
121  y = inputPoints(i0, 1);
122 
123  // outputValues is a rank-3 array with dimensions (basisCardinality_, dim0, spaceDim)
124  outputValues(0, i0, 0) = x;
125  outputValues(0, i0, 1) = y - 1.0;
126 
127  outputValues(1, i0, 0) = x;
128  outputValues(1, i0, 1) = y;
129 
130  outputValues(2, i0, 0) = x - 1.0;
131  outputValues(2, i0, 1) = y;
132  }
133  break;
134 
135  case OPERATOR_DIV:
136  // outputValues is a rank-2 array with dimensions (basisCardinality_, dim0)
137  for (int i0 = 0; i0 < dim0; i0++) {
138  outputValues(0, i0) = 2.0;
139  outputValues(1, i0) = 2.0;
140  outputValues(2, i0) = 2.0;
141  }
142  break;
143 
144  case OPERATOR_CURL:
145  TEUCHOS_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_CURL), std::invalid_argument,
146  ">>> ERROR (Basis_HDIV_TRI_I1_FEM): CURL is invalid operator for HDIV Basis Functions");
147  break;
148 
149  case OPERATOR_GRAD:
150  TEUCHOS_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_GRAD), std::invalid_argument,
151  ">>> ERROR (Basis_HDIV_TRI_I1_FEM): GRAD is invalid operator for HDIV Basis Functions");
152  break;
153 
154  case OPERATOR_D1:
155  case OPERATOR_D2:
156  case OPERATOR_D3:
157  case OPERATOR_D4:
158  case OPERATOR_D5:
159  case OPERATOR_D6:
160  case OPERATOR_D7:
161  case OPERATOR_D8:
162  case OPERATOR_D9:
163  case OPERATOR_D10:
164  TEUCHOS_TEST_FOR_EXCEPTION( ( (operatorType == OPERATOR_D1) &&
165  (operatorType == OPERATOR_D2) &&
166  (operatorType == OPERATOR_D3) &&
167  (operatorType == OPERATOR_D4) &&
168  (operatorType == OPERATOR_D5) &&
169  (operatorType == OPERATOR_D6) &&
170  (operatorType == OPERATOR_D7) &&
171  (operatorType == OPERATOR_D8) &&
172  (operatorType == OPERATOR_D9) &&
173  (operatorType == OPERATOR_D10) ),
174  std::invalid_argument,
175  ">>> ERROR (Basis_HDIV_TRI_I1_FEM): Invalid operator type");
176  break;
177 
178  default:
179  TEUCHOS_TEST_FOR_EXCEPTION( ( (operatorType != OPERATOR_VALUE) &&
180  (operatorType != OPERATOR_GRAD) &&
181  (operatorType != OPERATOR_CURL) &&
182  (operatorType != OPERATOR_DIV) &&
183  (operatorType != OPERATOR_D1) &&
184  (operatorType != OPERATOR_D2) &&
185  (operatorType != OPERATOR_D3) &&
186  (operatorType != OPERATOR_D4) &&
187  (operatorType != OPERATOR_D5) &&
188  (operatorType != OPERATOR_D6) &&
189  (operatorType != OPERATOR_D7) &&
190  (operatorType != OPERATOR_D8) &&
191  (operatorType != OPERATOR_D9) &&
192  (operatorType != OPERATOR_D10) ),
193  std::invalid_argument,
194  ">>> ERROR (Basis_HDIV_TRI_I1_FEM): Invalid operator type");
195  }
196 }
197 
198 
199 
200 template<class Scalar, class ArrayScalar>
201 void Basis_HDIV_TRI_I1_FEM<Scalar, ArrayScalar>::getValues(ArrayScalar& outputValues,
202  const ArrayScalar & inputPoints,
203  const ArrayScalar & cellVertices,
204  const EOperator operatorType) const {
205  TEUCHOS_TEST_FOR_EXCEPTION( (true), std::logic_error,
206  ">>> ERROR (Basis_HDIV_TRI_I1_FEM): FEM Basis calling an FVD member function");
207  }
208 
209 }// namespace Intrepid
Intrepid::Basis_HDIV_TRI_I1_FEM::getValues
void getValues(ArrayScalar &outputValues, const ArrayScalar &inputPoints, const EOperator operatorType) const
Evaluation of a FEM basis on a reference Triangle cell.
Definition: Intrepid_HDIV_TRI_I1_FEMDef.hpp:95
Intrepid::Basis_HDIV_TRI_I1_FEM::Basis_HDIV_TRI_I1_FEM
Basis_HDIV_TRI_I1_FEM()
Constructor.
Definition: Intrepid_HDIV_TRI_I1_FEMDef.hpp:52
Intrepid::Basis_HDIV_TRI_I1_FEM::initializeTags
void initializeTags()
Initializes tagToOrdinal_ and ordinalToTag_ lookup arrays.
Definition: Intrepid_HDIV_TRI_I1_FEMDef.hpp:63
Generated by   doxygen 1.8.5