libMesh
fe_hierarchic_vec.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2025 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner
3 
4 // This library is free software; you can redistribute it and/or
5 // modify it under the terms of the GNU Lesser General Public
6 // License as published by the Free Software Foundation; either
7 // version 2.1 of the License, or (at your option) any later version.
8 
9 // This library is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // Lesser General Public License for more details.
13 
14 // You should have received a copy of the GNU Lesser General Public
15 // License along with this library; if not, write to the Free Software
16 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 
18 
19 
20 // Local includes
21 #include "libmesh/dof_map.h"
22 #include "libmesh/elem.h"
23 #include "libmesh/enum_to_string.h"
24 #include "libmesh/fe.h"
25 #include "libmesh/fe_interface.h"
26 #include "libmesh/fe_macro.h"
27 #include "libmesh/tensor_value.h"
28 
29 
30 namespace libMesh
31 {
32 
33 
34 LIBMESH_DEFAULT_VECTORIZED_FE(0,HIERARCHIC_VEC)
35 LIBMESH_DEFAULT_VECTORIZED_FE(1,HIERARCHIC_VEC)
36 LIBMESH_DEFAULT_VECTORIZED_FE(2,HIERARCHIC_VEC)
37 LIBMESH_DEFAULT_VECTORIZED_FE(3,HIERARCHIC_VEC)
38 LIBMESH_DEFAULT_VECTORIZED_FE(0,L2_HIERARCHIC_VEC)
39 LIBMESH_DEFAULT_VECTORIZED_FE(1,L2_HIERARCHIC_VEC)
40 LIBMESH_DEFAULT_VECTORIZED_FE(2,L2_HIERARCHIC_VEC)
41 LIBMESH_DEFAULT_VECTORIZED_FE(3,L2_HIERARCHIC_VEC)
42 
43 
44 // ------------------------------------------------------------
45 // Hierarchic-specific implementations
46 
47 
48 // Anonymous namespace for local helper functions
49 namespace {
50 void hierarchic_vec_nodal_soln(const Elem * elem,
51  const Order order,
52  const std::vector<Number> & elem_soln,
53  const int dim,
54  std::vector<Number> & nodal_soln,
55  const bool add_p_level)
56 {
57  const unsigned int n_nodes = elem->n_nodes();
58 
59  // Constant shape functions can't be supported, even for
60  // L2_HIERARCHIC*, without breaking the "HIERARCHIC is
61  // hierarchic" guarantee
62  libmesh_assert(order != CONSTANT);
63 
64  nodal_soln.resize(dim*n_nodes);
65 
66  // Do interpolation at the nodes explicitly.
67  FEType fe_type(order, HIERARCHIC);
68 
69  const unsigned int n_sf =
70  FEInterface::n_shape_functions(fe_type, elem, add_p_level);
71 
72  std::vector<Point> refspace_nodes;
73  FEBase::get_refspace_nodes(elem->type(),refspace_nodes);
74  libmesh_assert_equal_to (refspace_nodes.size(), n_nodes);
75  libmesh_assert_equal_to (elem_soln.size(), n_sf*dim);
76 
77  // Zero before summation
78  std::fill(nodal_soln.begin(), nodal_soln.end(), 0);
79 
80  for (unsigned int n=0; n<n_nodes; n++)
81  for (int d=0; d != dim; ++d)
82  // u_i = Sum (alpha_i phi_i); we're here only looking
83  // at vector components in direction d
84  for (unsigned int i=0; i<n_sf; i++)
85  nodal_soln[n*dim+d] += elem_soln[i*dim+d] *
86  FEInterface::shape(fe_type, elem, i, refspace_nodes[n]);
87 
88 }// void hierarchic_vec_nodal_soln
89 
90 } // anonymous namespace
91 
92 
93  // Do full-specialization for every dimension, instead
94  // of explicit instantiation at the end of this file.
95  // This could be macro-ified so that it fits on one line...
96 LIBMESH_FE_NODAL_SOLN_DIM(HIERARCHIC_VEC, (FE<0, HIERARCHIC>::nodal_soln), 0)
97 LIBMESH_FE_NODAL_SOLN_DIM(HIERARCHIC_VEC, (FE<1, HIERARCHIC>::nodal_soln), 1)
98 
99 template <>
100 void FE<2,HIERARCHIC_VEC>::nodal_soln(const Elem * elem,
101  const Order order,
102  const std::vector<Number> & elem_soln,
103  std::vector<Number> & nodal_soln,
104  const bool add_p_level,
105  const unsigned)
106 { hierarchic_vec_nodal_soln(elem, order, elem_soln, 2 /*dim*/, nodal_soln, add_p_level); }
107 
108 template <>
109 void FE<3,HIERARCHIC_VEC>::nodal_soln(const Elem * elem,
110  const Order order,
111  const std::vector<Number> & elem_soln,
112  std::vector<Number> & nodal_soln,
113  const bool add_p_level,
114  const unsigned)
115 { hierarchic_vec_nodal_soln(elem, order, elem_soln, 3 /*dim*/, nodal_soln, add_p_level); }
116 
117 LIBMESH_FE_SIDE_NODAL_SOLN(HIERARCHIC_VEC)
118 
119 LIBMESH_FE_NODAL_SOLN_DIM(L2_HIERARCHIC_VEC, (FE<0, HIERARCHIC_VEC>::nodal_soln), 0)
120 LIBMESH_FE_NODAL_SOLN_DIM(L2_HIERARCHIC_VEC, (FE<1, HIERARCHIC_VEC>::nodal_soln), 1)
121 LIBMESH_FE_NODAL_SOLN_DIM(L2_HIERARCHIC_VEC, (FE<2, HIERARCHIC_VEC>::nodal_soln), 2)
122 LIBMESH_FE_NODAL_SOLN_DIM(L2_HIERARCHIC_VEC, (FE<3, HIERARCHIC_VEC>::nodal_soln), 3)
123 
124 LIBMESH_FE_SIDE_NODAL_SOLN(L2_HIERARCHIC_VEC)
125 
126 
127 // Specialize for shape function routines by leveraging scalar HIERARCHIC elements
128 
129 // 0-D
130 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
131  const unsigned int i, const Point & p)
132 {
133  Real value = FE<0,HIERARCHIC>::shape( type, order, i, p );
134  return libMesh::RealGradient( value );
135 }
136 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
137  const unsigned int i, const unsigned int j,
138  const Point & p)
139 {
140  Real value = FE<0,HIERARCHIC>::shape_deriv( type, order, i, j, p );
141  return libMesh::RealGradient( value );
142 }
143 
144 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
145 
146 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
147  const unsigned int i, const unsigned int j,
148  const Point & p)
149 {
150  Real value = FE<0,HIERARCHIC>::shape_second_deriv( type, order, i, j, p );
151  return libMesh::RealGradient( value );
152 }
153 #endif
154 
155 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
156  const unsigned int i, const Point & p)
157 {
158  return FE<0,HIERARCHIC_VEC>::shape(type, order, i, p);
159 }
160 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
161  const unsigned int i, const unsigned int j,
162  const Point & p)
163 {
164  return FE<0,HIERARCHIC_VEC>::shape_deriv(type, order, i, j, p);
165 }
166 
167 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
168 
169 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
170  const unsigned int i, const unsigned int j,
171  const Point & p)
172 {
173  return FE<0,HIERARCHIC_VEC>::shape_second_deriv(type, order, i, j, p);
174 }
175 #endif
176 
177 // 1-D
178 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
179  const unsigned int i, const Point & p)
180 {
181  Real value = FE<1,HIERARCHIC>::shape( type, order, i, p );
182  return libMesh::RealGradient( value );
183 }
184 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
185  const unsigned int i, const unsigned int j,
186  const Point & p)
187 {
188  Real value = FE<1,HIERARCHIC>::shape_deriv( type, order, i, j, p );
189  return libMesh::RealGradient( value );
190 }
191 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
192 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
193  const unsigned int i, const unsigned int j,
194  const Point & p)
195 {
196  Real value = FE<1,HIERARCHIC>::shape_second_deriv( type, order, i, j, p );
197  return libMesh::RealGradient( value );
198 }
199 
200 #endif
201 
202 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
203  const unsigned int i, const Point & p)
204 {
205  return FE<1,HIERARCHIC_VEC>::shape(type, order, i, p);
206 }
207 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
208  const unsigned int i, const unsigned int j,
209  const Point & p)
210 {
211  return FE<1,HIERARCHIC_VEC>::shape_deriv(type, order, i, j, p);
212 }
213 
214 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
215 
216 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
217  const unsigned int i, const unsigned int j,
218  const Point & p)
219 {
220  return FE<1,HIERARCHIC_VEC>::shape_second_deriv(type, order, i, j, p);
221 }
222 #endif
223 
224 // 2-D
225 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
226  const unsigned int i, const Point & p)
227 {
228  Real value = FE<2,HIERARCHIC>::shape( type, order, i/2, p );
229 
230  switch( i%2 )
231  {
232  case 0:
233  return libMesh::RealGradient( value );
234 
235  case 1:
236  return libMesh::RealGradient( Real(0), value );
237 
238  default:
239  libmesh_error_msg("i%2 must be either 0 or 1!");
240  }
241 
242  //dummy
243  return libMesh::RealGradient();
244 }
245 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
246  const unsigned int i, const unsigned int j,
247  const Point & p)
248 {
249  Real value = FE<2,HIERARCHIC>::shape_deriv( type, order, i/2, j, p );
250 
251  switch( i%2 )
252  {
253  case 0:
254  return libMesh::RealGradient( value );
255 
256  case 1:
257  return libMesh::RealGradient( Real(0), value );
258 
259  default:
260  libmesh_error_msg("i%2 must be either 0 or 1!");
261  }
262 
263  //dummy
264  return libMesh::RealGradient();
265 }
266 
267 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
268 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
269  const unsigned int i, const unsigned int j,
270  const Point & p)
271 {
272  Real value = FE<2,HIERARCHIC>::shape_second_deriv( type, order, i/2, j, p );
273 
274  switch( i%2 )
275  {
276  case 0:
277  return libMesh::RealGradient( value );
278 
279  case 1:
280  return libMesh::RealGradient( Real(0), value );
281 
282  default:
283  libmesh_error_msg("i%2 must be either 0 or 1!");
284  }
285 
286  //dummy
287  return libMesh::RealGradient();
288 }
289 
290 #endif
291 
292 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
293  const unsigned int i, const Point & p)
294 {
295  return FE<2,HIERARCHIC_VEC>::shape(type, order, i, p);
296 }
297 
298 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
299  const unsigned int i, const unsigned int j,
300  const Point & p)
301 {
302  return FE<2,HIERARCHIC_VEC>::shape_deriv(type, order, i, j, p);
303 }
304 
305 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
306 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
307  const unsigned int i, const unsigned int j,
308  const Point & p)
309 {
310  return FE<2,HIERARCHIC_VEC>::shape_second_deriv(type, order, i, j, p);
311 }
312 
313 #endif
314 
315 // 3-D
316 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
317  const unsigned int i, const Point & p)
318 {
319  Real value = FE<3,HIERARCHIC>::shape( type, order, i/3, p );
320 
321  switch( i%3 )
322  {
323  case 0:
324  return libMesh::RealGradient( value );
325 
326  case 1:
327  return libMesh::RealGradient( Real(0), value );
328 
329  case 2:
330  return libMesh::RealGradient( Real(0), Real(0), value );
331 
332  default:
333  libmesh_error_msg("i%3 must be 0, 1, or 2!");
334  }
335 
336  //dummy
337  return libMesh::RealGradient();
338 }
339 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
340  const unsigned int i, const unsigned int j,
341  const Point & p)
342 {
343  Real value = FE<3,HIERARCHIC>::shape_deriv( type, order, i/3, j, p );
344 
345  switch( i%3 )
346  {
347  case 0:
348  return libMesh::RealGradient( value );
349 
350  case 1:
351  return libMesh::RealGradient( Real(0), value );
352 
353  case 2:
354  return libMesh::RealGradient( Real(0), Real(0), value );
355 
356  default:
357  libmesh_error_msg("i%3 must be 0, 1, or 2!");
358  }
359 
360  //dummy
361  return libMesh::RealGradient();
362 }
363 
364 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
365 
366 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
367  const unsigned int i, const unsigned int j,
368  const Point & p)
369 {
370  Real value = FE<3,HIERARCHIC>::shape_second_deriv( type, order, i/3, j, p );
371 
372  switch( i%3 )
373  {
374  case 0:
375  return libMesh::RealGradient( value );
376 
377  case 1:
378  return libMesh::RealGradient( Real(0), value );
379 
380  case 2:
381  return libMesh::RealGradient( Real(0), Real(0), value );
382 
383  default:
384  libmesh_error_msg("i%3 must be 0, 1, or 2!");
385  }
386 
387  //dummy
388  return libMesh::RealGradient();
389 }
390 
391 #endif
392 
393 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape(const ElemType type, const Order order,
394  const unsigned int i, const Point & p)
395 {
396  return FE<3,HIERARCHIC_VEC>::shape(type, order, i, p);
397 }
398 
399 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape_deriv(const ElemType type, const Order order,
400  const unsigned int i, const unsigned int j,
401  const Point & p)
402 {
403  return FE<3,HIERARCHIC_VEC>::shape_deriv(type, order, i, j, p);
404 }
405 
406 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
407 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape_second_deriv(const ElemType type, const Order order,
408  const unsigned int i, const unsigned int j,
409  const Point & p)
410 {
411  return FE<3,HIERARCHIC_VEC>::shape_second_deriv(type, order, i, j, p);
412 }
413 
414 #endif
415 
416 
417 // 0-D
418 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
419  const unsigned int i, const Point & p,
420  const bool add_p_level)
421 {
422  const Real value = FE<0,HIERARCHIC>::shape(elem, order, i, p, add_p_level);
423  return libMesh::RealGradient( value );
424 }
425 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
426  const unsigned int i, const unsigned int j,
427  const Point & p,
428  const bool add_p_level)
429 {
430  const Real value = FE<0,HIERARCHIC>::shape_deriv(elem, order, i, j, p, add_p_level);
431  return libMesh::RealGradient( value );
432 }
433 
434 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
435 
436 template <> RealGradient FE<0,HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
437  const unsigned int i, const unsigned int j,
438  const Point & p,
439  const bool add_p_level)
440 {
441  Real value = FE<0,HIERARCHIC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
442  return libMesh::RealGradient( value );
443 }
444 
445 #endif
446 
447 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
448  const unsigned int i, const Point & p,
449  const bool add_p_level)
450 {
451  return FE<0,HIERARCHIC_VEC>::shape(elem, order, i, p, add_p_level);
452 }
453 
454 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
455  const unsigned int i, const unsigned int j,
456  const Point & p,
457  const bool add_p_level)
458 {
459  return FE<0,HIERARCHIC_VEC>::shape_deriv(elem, order, i, j, p, add_p_level);
460 }
461 
462 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
463 
464 template <> RealGradient FE<0,L2_HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
465  const unsigned int i, const unsigned int j,
466  const Point & p,
467  const bool add_p_level)
468 {
469  return FE<0,HIERARCHIC_VEC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
470 }
471 
472 #endif
473 
474 // 1-D
475 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
476  const unsigned int i, const Point & p,
477  const bool add_p_level)
478 {
479  Real value = FE<1,HIERARCHIC>::shape(elem, order, i, p, add_p_level);
480  return libMesh::RealGradient( value );
481 }
482 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
483  const unsigned int i, const unsigned int j,
484  const Point & p,
485  const bool add_p_level)
486 {
487  Real value = FE<1,HIERARCHIC>::shape_deriv(elem, order, i, j, p, add_p_level);
488  return libMesh::RealGradient( value );
489 }
490 
491 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
492 template <> RealGradient FE<1,HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
493  const unsigned int i, const unsigned int j,
494  const Point & p,
495  const bool add_p_level)
496 {
497  Real value = FE<1,HIERARCHIC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
498  return libMesh::RealGradient( value );
499 }
500 
501 #endif
502 
503 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
504  const unsigned int i, const Point & p,
505  const bool add_p_level)
506 {
507  return FE<1,HIERARCHIC_VEC>::shape(elem, order, i, p, add_p_level);
508 }
509 
510 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
511  const unsigned int i, const unsigned int j,
512  const Point & p,
513  const bool add_p_level)
514 {
515  return FE<1,HIERARCHIC_VEC>::shape_deriv(elem, order, i, j, p, add_p_level);
516 }
517 
518 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
519 
520 template <> RealGradient FE<1,L2_HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
521  const unsigned int i, const unsigned int j,
522  const Point & p,
523  const bool add_p_level)
524 {
525  return FE<1,HIERARCHIC_VEC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
526 }
527 
528 #endif
529 
530 // 2-D
531 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
532  const unsigned int i, const Point & p,
533  const bool add_p_level)
534 {
535  const Real value = FE<2,HIERARCHIC>::shape(elem, order, i/2, p, add_p_level);
536 
537  switch( i%2 )
538  {
539  case 0:
540  return libMesh::RealGradient( value );
541 
542  case 1:
543  return libMesh::RealGradient( Real(0), value );
544 
545  default:
546  libmesh_error_msg("i%2 must be either 0 or 1!");
547  }
548 
549  //dummy
550  return libMesh::RealGradient();
551 }
552 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
553  const unsigned int i, const unsigned int j,
554  const Point & p,
555  const bool add_p_level)
556 {
557  const Real value = FE<2,HIERARCHIC>::shape_deriv(elem, order, i/2, j, p, add_p_level);
558 
559  switch( i%2 )
560  {
561  case 0:
562  return libMesh::RealGradient( value );
563 
564  case 1:
565  return libMesh::RealGradient( Real(0), value );
566 
567  default:
568  libmesh_error_msg("i%2 must be either 0 or 1!");
569  }
570 
571  //dummy
572  return libMesh::RealGradient();
573 }
574 
575 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
576 template <> RealGradient FE<2,HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
577  const unsigned int i, const unsigned int j,
578  const Point & p,
579  const bool add_p_level)
580 {
581  const Real value = FE<2,HIERARCHIC>::shape_second_deriv(elem, order, i/2, j, p, add_p_level);
582 
583  switch( i%2 )
584  {
585  case 0:
586  return libMesh::RealGradient( value );
587 
588  case 1:
589  return libMesh::RealGradient( Real(0), value );
590 
591  default:
592  libmesh_error_msg("i%2 must be either 0 or 1!");
593  }
594 
595  //dummy
596  return libMesh::RealGradient();
597 }
598 
599 #endif
600 
601 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
602  const unsigned int i, const Point & p,
603  const bool add_p_level)
604 {
605  return FE<2,HIERARCHIC_VEC>::shape(elem, order, i, p, add_p_level);
606 }
607 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
608  const unsigned int i, const unsigned int j,
609  const Point & p,
610  const bool add_p_level)
611 {
612  return FE<2,HIERARCHIC_VEC>::shape_deriv(elem, order, i, j, p, add_p_level);
613 }
614 
615 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
616 
617 template <> RealGradient FE<2,L2_HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
618  const unsigned int i, const unsigned int j,
619  const Point & p,
620  const bool add_p_level)
621 {
622  return FE<2,HIERARCHIC_VEC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
623 }
624 
625 #endif
626 
627 // 3-D
628 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
629  const unsigned int i, const Point & p,
630  const bool add_p_level)
631 {
632  const Real value = FE<3,HIERARCHIC>::shape(elem, order, i/3, p, add_p_level);
633 
634  switch( i%3 )
635  {
636  case 0:
637  return libMesh::RealGradient( value );
638 
639  case 1:
640  return libMesh::RealGradient( Real(0), value );
641 
642  case 2:
643  return libMesh::RealGradient( Real(0), Real(0), value );
644 
645  default:
646  libmesh_error_msg("i%3 must be 0, 1, or 2!");
647  }
648 
649  //dummy
650  return libMesh::RealGradient();
651 }
652 
653 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
654  const unsigned int i, const unsigned int j,
655  const Point & p,
656  const bool add_p_level)
657 {
658  const Real value = FE<3,HIERARCHIC>::shape_deriv(elem, order, i/3, j, p, add_p_level);
659 
660  switch( i%3 )
661  {
662  case 0:
663  return libMesh::RealGradient( value );
664 
665  case 1:
666  return libMesh::RealGradient( Real(0), value );
667 
668  case 2:
669  return libMesh::RealGradient( Real(0), Real(0), value );
670 
671  default:
672  libmesh_error_msg("i%3 must be 0, 1, or 2!");
673  }
674 
675  //dummy
676  return libMesh::RealGradient();
677 }
678 
679 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
680 
681 template <> RealGradient FE<3,HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
682  const unsigned int i, const unsigned int j,
683  const Point & p,
684  const bool add_p_level)
685 {
686  const Real value = FE<3,HIERARCHIC>::shape_second_deriv(elem, order, i/3, j, p, add_p_level);
687 
688  switch( i%3 )
689  {
690  case 0:
691  return libMesh::RealGradient( value );
692 
693  case 1:
694  return libMesh::RealGradient( Real(0), value );
695 
696  case 2:
697  return libMesh::RealGradient( Real(0), Real(0), value );
698 
699  default:
700  libmesh_error_msg("i%3 must be 0, 1, or 2!");
701  }
702 
703  //dummy
704  return libMesh::RealGradient();
705 }
706 
707 #endif
708 
709 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape(const Elem * elem, const Order order,
710  const unsigned int i, const Point & p,
711  const bool add_p_level)
712 {
713  return FE<3,HIERARCHIC_VEC>::shape(elem, order, i, p, add_p_level);
714 }
715 
716 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape_deriv(const Elem * elem, const Order order,
717  const unsigned int i, const unsigned int j,
718  const Point & p,
719  const bool add_p_level)
720 {
721  return FE<3,HIERARCHIC_VEC>::shape_deriv(elem, order, i, j, p, add_p_level);
722 }
723 
724 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
725 
726 template <> RealGradient FE<3,L2_HIERARCHIC_VEC>::shape_second_deriv(const Elem * elem, const Order order,
727  const unsigned int i, const unsigned int j,
728  const Point & p,
729  const bool add_p_level)
730 {
731  return FE<3,HIERARCHIC_VEC>::shape_second_deriv(elem, order, i, j, p, add_p_level);
732 }
733 
734 #endif
735 
736 
737 // Instantiate n_dofs*() functions for every dimension
738 LIBMESH_DEFAULT_VEC_NDOFS(HIERARCHIC)
739 
740 // L2 Hierarchic has a simpler but non-default implementation, with
741 // all their dofs on the element
742 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs(const ElemType t, const Order o) { return FE<0,L2_HIERARCHIC>::n_dofs(t,o); }
743 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs(const ElemType t, const Order o) { return FE<1,L2_HIERARCHIC>::n_dofs(t,o); }
744 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs(const ElemType t, const Order o) { return 2*FE<2,L2_HIERARCHIC>::n_dofs(t,o); }
745 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs(const ElemType t, const Order o) { return 3*FE<3,L2_HIERARCHIC>::n_dofs(t,o); }
746 
747 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs(const Elem * e, const Order o) { return FE<0,L2_HIERARCHIC>::n_dofs(e,o); }
748 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs(const Elem * e, const Order o) { return FE<1,L2_HIERARCHIC>::n_dofs(e,o); }
749 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs(const Elem * e, const Order o) { return 2*FE<2,L2_HIERARCHIC>::n_dofs(e,o); }
750 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs(const Elem * e, const Order o) { return 3*FE<3,L2_HIERARCHIC>::n_dofs(e,o); }
751 
752 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs_at_node(const ElemType, const Order, const unsigned int) { return 0; }
753 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs_at_node(const ElemType, const Order, const unsigned int) { return 0; }
754 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs_at_node(const ElemType, const Order, const unsigned int) { return 0; }
755 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs_at_node(const ElemType, const Order, const unsigned int) { return 0; }
756 
757 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs_at_node(const Elem &, const Order, const unsigned int) { return 0; }
758 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs_at_node(const Elem &, const Order, const unsigned int) { return 0; }
759 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs_at_node(const Elem &, const Order, const unsigned int) { return 0; }
760 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs_at_node(const Elem &, const Order, const unsigned int) { return 0; }
761 
762 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const ElemType t, const Order o) { return FE<0,L2_HIERARCHIC_VEC>::n_dofs(t, o); }
763 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const ElemType t, const Order o) { return FE<1,L2_HIERARCHIC_VEC>::n_dofs(t, o); }
764 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const ElemType t, const Order o) { return FE<2,L2_HIERARCHIC_VEC>::n_dofs(t, o); }
765 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const ElemType t, const Order o) { return FE<3,L2_HIERARCHIC_VEC>::n_dofs(t, o); }
766 
767 template <> unsigned int FE<0,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const Elem & e, const Order o) { return FE<0,L2_HIERARCHIC_VEC>::n_dofs(&e, o); }
768 template <> unsigned int FE<1,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const Elem & e, const Order o) { return FE<1,L2_HIERARCHIC_VEC>::n_dofs(&e, o); }
769 template <> unsigned int FE<2,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const Elem & e, const Order o) { return FE<2,L2_HIERARCHIC_VEC>::n_dofs(&e, o); }
770 template <> unsigned int FE<3,L2_HIERARCHIC_VEC>::n_dofs_per_elem(const Elem & e, const Order o) { return FE<3,L2_HIERARCHIC_VEC>::n_dofs(&e, o); }
771 
772 // Hierarchic FEMs are always C^0 continuous
773 template <> FEContinuity FE<0,HIERARCHIC_VEC>::get_continuity() const { return C_ZERO; }
774 template <> FEContinuity FE<1,HIERARCHIC_VEC>::get_continuity() const { return C_ZERO; }
775 template <> FEContinuity FE<2,HIERARCHIC_VEC>::get_continuity() const { return C_ZERO; }
776 template <> FEContinuity FE<3,HIERARCHIC_VEC>::get_continuity() const { return C_ZERO; }
777 
778 // L2 Hierarchic FEMs are always discontinuous
779 template <> FEContinuity FE<0,L2_HIERARCHIC_VEC>::get_continuity() const { return DISCONTINUOUS; }
780 template <> FEContinuity FE<1,L2_HIERARCHIC_VEC>::get_continuity() const { return DISCONTINUOUS; }
781 template <> FEContinuity FE<2,L2_HIERARCHIC_VEC>::get_continuity() const { return DISCONTINUOUS; }
782 template <> FEContinuity FE<3,L2_HIERARCHIC_VEC>::get_continuity() const { return DISCONTINUOUS; }
783 
784 // Hierarchic FEMs are hierarchic
785 template <> bool FE<0,HIERARCHIC_VEC>::is_hierarchic() const { return true; }
786 template <> bool FE<1,HIERARCHIC_VEC>::is_hierarchic() const { return true; }
787 template <> bool FE<2,HIERARCHIC_VEC>::is_hierarchic() const { return true; }
788 template <> bool FE<3,HIERARCHIC_VEC>::is_hierarchic() const { return true; }
789 template <> bool FE<0,L2_HIERARCHIC_VEC>::is_hierarchic() const { return true; }
790 template <> bool FE<1,L2_HIERARCHIC_VEC>::is_hierarchic() const { return true; }
791 template <> bool FE<2,L2_HIERARCHIC_VEC>::is_hierarchic() const { return true; }
792 template <> bool FE<3,L2_HIERARCHIC_VEC>::is_hierarchic() const { return true; }
793 
794 // Hierarchic FEM shapes need reinit
795 template <> bool FE<0,HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
796 template <> bool FE<1,HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
797 template <> bool FE<2,HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
798 template <> bool FE<3,HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
799 template <> bool FE<0,L2_HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
800 template <> bool FE<1,L2_HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
801 template <> bool FE<2,L2_HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
802 template <> bool FE<3,L2_HIERARCHIC_VEC>::shapes_need_reinit() const { return true; }
803 
804 // Methods for computing Hierarchic constraints. Note: we pass the
805 // dimension as the last argument to the anonymous helper function.
806 // Also note: we only need instantiations of this function for
807 // Dim==2 and 3.
808 #ifdef LIBMESH_ENABLE_AMR
809 template <>
810 void FE<2,HIERARCHIC_VEC>::compute_constraints (DofConstraints & constraints,
811  DofMap & dof_map,
812  const unsigned int variable_number,
813  const Elem * elem)
814 { //libmesh_not_implemented();
815  FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem);
816 }
817 
818 template <>
819 void FE<3,HIERARCHIC_VEC>::compute_constraints (DofConstraints & constraints,
820  DofMap & dof_map,
821  const unsigned int variable_number,
822  const Elem * elem)
823 { //libmesh_not_implemented();
824  FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem);
825 }
826 
827 template <>
828 void FE<2,L2_HIERARCHIC_VEC>::compute_constraints (DofConstraints & constraints,
829  DofMap & dof_map,
830  const unsigned int variable_number,
831  const Elem * elem)
832 { //libmesh_not_implemented();
833  FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem);
834 }
835 
836 template <>
837 void FE<3,L2_HIERARCHIC_VEC>::compute_constraints (DofConstraints & constraints,
838  DofMap & dof_map,
839  const unsigned int variable_number,
840  const Elem * elem)
841 { //libmesh_not_implemented();
842  FEVectorBase::compute_proj_constraints(constraints, dof_map, variable_number, elem);
843 }
844 #endif // LIBMESH_ENABLE_AMR
845 
846 } // namespace libMesh
libMesh::FEType
class FEType hides (possibly multiple) FEFamily and approximation orders, thereby enabling specialize...
Definition: fe_type.h:196
libMesh::FE::n_dofs
static unsigned int n_dofs(const ElemType t, const Order o)
libMesh::ElemType
ElemType
Defines an enum for geometric element types.
Definition: enum_elem_type.h:33
libMesh::RealGradient
RealVectorValue RealGradient
Definition: hp_coarsentest.h:49
libMesh::Order
Order
defines an enum for polynomial orders.
Definition: enum_order.h:40
libMesh::LIBMESH_FE_NODAL_SOLN_DIM
LIBMESH_FE_NODAL_SOLN_DIM(LIBMESH_FE_NODAL_SOLN_DIM(HIERARCHIC_VEC,(FE< 0, HIERARCHIC >::nodal_soln), 0)
Definition: fe_hierarchic_vec.C:96
libMesh::FE::shape
static OutputShape shape(const ElemType t, const Order o, const unsigned int i, const Point &p)
dim
unsigned int dim
Definition: adaptivity_ex3.C:113
libMesh::Elem
This is the base class from which all geometric element types are derived.
Definition: elem.h:94
libMesh::FE::n_dofs_at_node
static unsigned int n_dofs_at_node(const ElemType t, const Order o, const unsigned int n)
libMesh::FE::shape_deriv
static OutputShape shape_deriv(const ElemType t, const Order o, const unsigned int i, const unsigned int j, const Point &p)
libMesh::FE::shapes_need_reinit
virtual bool shapes_need_reinit() const override
libMesh
The libMesh namespace provides an interface to certain functionality in the library.
libMesh::FE::is_hierarchic
virtual bool is_hierarchic() const override
libMesh::DofMap
This class handles the numbering of degrees of freedom on a mesh.
Definition: dof_map.h:179
libMesh::LIBMESH_DEFAULT_VECTORIZED_FE
LIBMESH_DEFAULT_VECTORIZED_FE(template<>Real FE< 0, BERNSTEIN)
Definition: fe_bernstein_shape_0D.C:30
libMesh::FE
A specific instantiation of the FEBase class.
Definition: fe.h:127
libMesh::LIBMESH_FE_SIDE_NODAL_SOLN
LIBMESH_FE_SIDE_NODAL_SOLN(HIERARCHIC_VEC)
Definition: fe_hierarchic_vec.C:117
n_nodes
const dof_id_type n_nodes
Definition: tecplot_io.C:67
libMesh::FEInterface::n_shape_functions
static unsigned int n_shape_functions(const unsigned int dim, const FEType &fe_t, const ElemType t)
Definition: fe_interface.C:272
libMesh::Elem::n_nodes
virtual unsigned int n_nodes() const =0
libMesh::FEInterface::shape
static Real shape(const unsigned int dim, const FEType &fe_t, const ElemType t, const unsigned int i, const Point &p)
Definition: fe_interface.C:760
libMesh::FEAbstract::get_refspace_nodes
static void get_refspace_nodes(const ElemType t, std::vector< Point > &nodes)
Definition: fe_abstract.C:400
libMesh::libmesh_assert
libmesh_assert(ctx)
libMesh::FEGenericBase::compute_proj_constraints
static void compute_proj_constraints(DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
Computes the constraint matrix contributions (for non-conforming adapted meshes) corresponding to var...
Definition: fe_base.C:1570
libMesh::FE::n_dofs_per_elem
static unsigned int n_dofs_per_elem(const ElemType t, const Order o)
libMesh::FE::get_continuity
virtual FEContinuity get_continuity() const override
libMesh::FE::compute_constraints
static void compute_constraints(DofConstraints &constraints, DofMap &dof_map, const unsigned int variable_number, const Elem *elem)
Computes the constraint matrix contributions (for non-conforming adapted meshes) corresponding to var...
libMesh::FE::nodal_soln
static void nodal_soln(const Elem *elem, const Order o, const std::vector< Number > &elem_soln, std::vector< Number > &nodal_soln, bool add_p_level=true, const unsigned vdim=1)
Build the nodal soln from the element soln.
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Definition: libmesh_common.h:144
libMesh::LIBMESH_DEFAULT_VEC_NDOFS
LIBMESH_DEFAULT_VEC_NDOFS(HIERARCHIC)
Definition: fe_hierarchic_vec.C:738
libMesh::FEContinuity
FEContinuity
defines an enum for finite element types to libmesh_assert a certain level (or type? Hcurl?) of continuity.
Definition: enum_fe_family.h:84
value
static const bool value
Definition: xdr_io.C:54
libMesh::Elem::type
virtual ElemType type() const =0
libMesh::Point
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:39
libMesh::DofConstraints
The constraint matrix storage format.
Definition: dof_map.h:108
libMesh::FE::shape_second_deriv
static OutputShape shape_second_deriv(const ElemType t, const Order o, const unsigned int i, const unsigned int j, const Point &p)
Generated on Thu Jul 17 2025 01:29:06 for libMesh by   doxygen 1.8.14