Intrepid
Intrepid_HCURL_TET_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>
53  {
54  this -> basisCardinality_ = 6;
55  this -> basisDegree_ = 1;
56  this -> basisCellTopology_ = shards::CellTopology(shards::getCellTopologyData<shards::Tetrahedron<4> >() );
57  this -> basisType_ = BASIS_FEM_DEFAULT;
58  this -> basisCoordinates_ = COORDINATES_CARTESIAN;
59  this -> basisTagsAreSet_ = false;
60  }
61 
62 template<class Scalar, class ArrayScalar>
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  1, 3, 0, 1,
77  1, 4, 0, 1,
78  1, 5, 0, 1 };
79 
80  // Basis-independent function sets tag and enum data in tagToOrdinal_ and ordinalToTag_ arrays:
81  Intrepid::setOrdinalTagData(this -> tagToOrdinal_,
82  this -> ordinalToTag_,
83  tags,
84  this -> basisCardinality_,
85  tagSize,
86  posScDim,
87  posScOrd,
88  posDfOrd);
89 }
90 
91 
92 
93 template<class Scalar, class ArrayScalar>
95  const ArrayScalar & inputPoints,
96  const EOperator operatorType) const {
97 
98 // Verify arguments
99 #ifdef HAVE_INTREPID_DEBUG
100  Intrepid::getValues_HCURL_Args<Scalar, ArrayScalar>(outputValues,
101  inputPoints,
102  operatorType,
103  this -> getBaseCellTopology(),
104  this -> getCardinality() );
105 #endif
106 
107  // Number of evaluation points = dim 0 of inputPoints
108  int dim0 = inputPoints.dimension(0);
109 
110  // Temporaries: (x,y,z) coordinates of the evaluation point
111  Scalar x = 0.0;
112  Scalar y = 0.0;
113  Scalar z = 0.0;
114 
115  switch (operatorType) {
116  case OPERATOR_VALUE:
117  for (int i0 = 0; i0 < dim0; i0++) {
118  x = inputPoints(i0, 0);
119  y = inputPoints(i0, 1);
120  z = inputPoints(i0, 2);
121 
122  // outputValues is a rank-3 array with dimensions (basisCardinality_, dim0, spaceDim)
123  outputValues(0, i0, 0) = 1.0 - y - z;
124  outputValues(0, i0, 1) = x;
125  outputValues(0, i0, 2) = x;
126 
127  outputValues(1, i0, 0) =-y;
128  outputValues(1, i0, 1) = x;
129  outputValues(1, i0, 2) = 0.0;
130 
131  outputValues(2, i0, 0) = -y;
132  outputValues(2, i0, 1) = -1.0 + x + z;
133  outputValues(2, i0, 2) = -y;
134 
135  outputValues(3, i0, 0) = z;
136  outputValues(3, i0, 1) = z;
137  outputValues(3, i0, 2) = 1.0 - x - y;
138 
139  outputValues(4, i0, 0) =-z;
140  outputValues(4, i0, 1) = 0.0;
141  outputValues(4, i0, 2) = x;
142 
143  outputValues(5, i0, 0) = 0.0;
144  outputValues(5, i0, 1) =-z;
145  outputValues(5, i0, 2) = y;
146  }
147  break;
148 
149  case OPERATOR_CURL:
150  for (int i0 = 0; i0 < dim0; i0++) {
151  x = inputPoints(i0, 0);
152  y = inputPoints(i0, 1);
153  z = inputPoints(i0, 2);
154 
155  // outputValues is a rank-3 array with dimensions (basisCardinality_, dim0, spaceDim)
156  outputValues(0, i0, 0) = 0.0;
157  outputValues(0, i0, 1) =-2.0;
158  outputValues(0, i0, 2) = 2.0;
159 
160  outputValues(1, i0, 0) = 0.0;
161  outputValues(1, i0, 1) = 0.0;
162  outputValues(1, i0, 2) = 2.0;
163 
164  outputValues(2, i0, 0) =-2.0;
165  outputValues(2, i0, 1) = 0.0;
166  outputValues(2, i0, 2) = 2.0;
167 
168  outputValues(3, i0, 0) =-2.0;
169  outputValues(3, i0, 1) = 2.0;
170  outputValues(3, i0, 2) = 0.0;
171 
172  outputValues(4, i0, 0) = 0.0;
173  outputValues(4, i0, 1) =-2.0;
174  outputValues(4, i0, 2) = 0.0;
175 
176  outputValues(5, i0, 0) = 2.0;
177  outputValues(5, i0, 1) = 0.0;
178  outputValues(5, i0, 2) = 0.0;
179  }
180  break;
181 
182  case OPERATOR_DIV:
183  TEUCHOS_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_DIV), std::invalid_argument,
184  ">>> ERROR (Basis_HCURL_TET_I1_FEM): DIV is invalid operator for HCURL Basis Functions");
185  break;
186 
187  case OPERATOR_GRAD:
188  TEUCHOS_TEST_FOR_EXCEPTION( (operatorType == OPERATOR_GRAD), std::invalid_argument,
189  ">>> ERROR (Basis_HCURL_TET_I1_FEM): GRAD is invalid operator for HCURL Basis Functions");
190  break;
191 
192  case OPERATOR_D1:
193  case OPERATOR_D2:
194  case OPERATOR_D3:
195  case OPERATOR_D4:
196  case OPERATOR_D5:
197  case OPERATOR_D6:
198  case OPERATOR_D7:
199  case OPERATOR_D8:
200  case OPERATOR_D9:
201  case OPERATOR_D10:
202  TEUCHOS_TEST_FOR_EXCEPTION( ( (operatorType == OPERATOR_D1) ||
203  (operatorType == OPERATOR_D2) ||
204  (operatorType == OPERATOR_D3) ||
205  (operatorType == OPERATOR_D4) ||
206  (operatorType == OPERATOR_D5) ||
207  (operatorType == OPERATOR_D6) ||
208  (operatorType == OPERATOR_D7) ||
209  (operatorType == OPERATOR_D8) ||
210  (operatorType == OPERATOR_D9) ||
211  (operatorType == OPERATOR_D10) ),
212  std::invalid_argument,
213  ">>> ERROR (Basis_HCURL_TET_I1_FEM): Invalid operator type");
214  break;
215 
216  default:
217  TEUCHOS_TEST_FOR_EXCEPTION( ( (operatorType != OPERATOR_VALUE) &&
218  (operatorType != OPERATOR_GRAD) &&
219  (operatorType != OPERATOR_CURL) &&
220  (operatorType != OPERATOR_DIV) &&
221  (operatorType != OPERATOR_D1) &&
222  (operatorType != OPERATOR_D2) &&
223  (operatorType != OPERATOR_D3) &&
224  (operatorType != OPERATOR_D4) &&
225  (operatorType != OPERATOR_D5) &&
226  (operatorType != OPERATOR_D6) &&
227  (operatorType != OPERATOR_D7) &&
228  (operatorType != OPERATOR_D8) &&
229  (operatorType != OPERATOR_D9) &&
230  (operatorType != OPERATOR_D10) ),
231  std::invalid_argument,
232  ">>> ERROR (Basis_HCURL_TET_I1_FEM): Invalid operator type");
233  }
234 }
235 
236 
237 
238 template<class Scalar, class ArrayScalar>
240  const ArrayScalar & inputPoints,
241  const ArrayScalar & cellVertices,
242  const EOperator operatorType) const {
243  TEUCHOS_TEST_FOR_EXCEPTION( (true), std::logic_error,
244  ">>> ERROR (Basis_HCURL_TET_I1_FEM): FEM Basis calling an FVD member function");
245  }
246 
247 template<class Scalar, class ArrayScalar>
249 #ifdef HAVE_INTREPID_DEBUG
250  // Verify rank of output array.
251  TEUCHOS_TEST_FOR_EXCEPTION( !(DofCoords.rank() == 2), std::invalid_argument,
252  ">>> ERROR: (Intrepid::Basis_HCURL_TET_I1_FEM::getDofCoords) rank = 2 required for DofCoords array");
253  // Verify 0th dimension of output array.
254  TEUCHOS_TEST_FOR_EXCEPTION( !( DofCoords.dimension(0) == this -> basisCardinality_ ), std::invalid_argument,
255  ">>> ERROR: (Intrepid::Basis_HCURL_TET_I1_FEM::getDofCoords) mismatch in number of DoF and 0th dimension of DofCoords array");
256  // Verify 1st dimension of output array.
257  TEUCHOS_TEST_FOR_EXCEPTION( !( DofCoords.dimension(1) == (int)(this -> basisCellTopology_.getDimension()) ), std::invalid_argument,
258  ">>> ERROR: (Intrepid::Basis_HCURL_TET_I1_FEM::getDofCoords) incorrect reference cell (1st) dimension in DofCoords array");
259 #endif
260 
261  DofCoords(0,0) = 0.5; DofCoords(0,1) = 0.0; DofCoords(0,2) = 0.0;
262  DofCoords(1,0) = 0.5; DofCoords(1,1) = 0.5; DofCoords(1,2) = 0.0;
263  DofCoords(2,0) = 0.0; DofCoords(2,1) = 0.5; DofCoords(2,2) = 0.0;
264  DofCoords(3,0) = 0.0; DofCoords(3,1) = 0.0; DofCoords(3,2) = 0.5;
265  DofCoords(4,0) = 0.5; DofCoords(4,1) = 0.0; DofCoords(4,2) = 0.5;
266  DofCoords(5,0) = 0.0; DofCoords(5,1) = 0.5; DofCoords(5,2) = 0.5;
267 }
268 
269 }// namespace Intrepid
270 
271 #if defined(Intrepid_SHOW_DEPRECATED_WARNINGS)
272 #ifdef __GNUC__
273 #warning "The Intrepid package is deprecated"
274 #endif
275 #endif
276 
void getDofCoords(ArrayScalar &DofCoords) const
Returns spatial locations (coordinates) of degrees of freedom on a reference Triangle.
void getValues(ArrayScalar &outputValues, const ArrayScalar &inputPoints, const EOperator operatorType) const
Evaluation of a FEM basis on a reference Tetrahedron cell.
void initializeTags()
Initializes tagToOrdinal_ and ordinalToTag_ lookup arrays.