libMesh
fe_clough_shape_3D.C
Go to the documentation of this file.
1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2019 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 // C++ includes
20 
21 // Local includes
22 #include "libmesh/fe.h"
23 #include "libmesh/elem.h"
24 
25 
26 // FIXME: 3D C1 finite elements are still a work in progress
27 
28 namespace libMesh
29 {
30 
31 
32 template <>
34  const Order,
35  const unsigned int,
36  const Point &)
37 {
38  libmesh_error_msg("Clough-Tocher elements require the real element \nto construct gradient-based degrees of freedom.");
39  return 0.;
40 }
41 
42 
43 
44 template <>
45 Real FE<3,CLOUGH>::shape(const Elem * libmesh_dbg_var(elem),
46  const Order,
47  const unsigned int,
48  const Point &,
49  const bool)
50 {
51  libmesh_assert(elem);
52 
53  libmesh_not_implemented();
54  return 0.;
55 }
56 
57 
58 
59 template <>
61  const Order,
62  const unsigned int,
63  const unsigned int,
64  const Point &)
65 {
66  libmesh_error_msg("Clough-Tocher elements require the real element \nto construct gradient-based degrees of freedom.");
67  return 0.;
68 }
69 
70 
71 
72 template <>
73 Real FE<3,CLOUGH>::shape_deriv(const Elem * libmesh_dbg_var(elem),
74  const Order,
75  const unsigned int,
76  const unsigned int,
77  const Point &,
78  const bool)
79 {
80  libmesh_assert(elem);
81  libmesh_not_implemented();
82  return 0.;
83 }
84 
85 
86 
87 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
88 
89 template <>
91  const Order,
92  const unsigned int,
93  const unsigned int,
94  const Point &)
95 {
96  libmesh_error_msg("Clough-Tocher elements require the real element \nto construct gradient-based degrees of freedom.");
97  return 0.;
98 }
99 
100 
101 template <>
102 Real FE<3,CLOUGH>::shape_second_deriv(const Elem * libmesh_dbg_var(elem),
103  const Order,
104  const unsigned int,
105  const unsigned int,
106  const Point &,
107  const bool)
108 {
109  libmesh_assert(elem);
110  libmesh_not_implemented();
111  return 0.;
112 }
113 
114 #endif
115 
116 } // namespace libMesh
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)
libMesh
The libMesh namespace provides an interface to certain functionality in the library.
Definition: factoryfunction.C:55
libMesh::Order
Order
Definition: enum_order.h:40
libMesh::libmesh_assert
libmesh_assert(ctx)
libMesh::FE::shape
static OutputShape shape(const ElemType t, const Order o, const unsigned int i, const Point &p)
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::Point
A Point defines a location in LIBMESH_DIM dimensional Real space.
Definition: point.h:38
libMesh::Elem
This is the base class from which all geometric element types are derived.
Definition: elem.h:100
libMesh::Real
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real
Definition: libmesh_common.h:121
libMesh::ElemType
ElemType
Defines an enum for geometric element types.
Definition: enum_elem_type.h:33