1 // The libMesh Finite Element Library.
2 // Copyright (C) 2002-2020 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 namespace libMesh
27 {
28
29
30 LIBMESH_DEFAULT_VECTORIZED_FE(0,SZABAB)
31
32
33 template <>
shape(const ElemType,const Order,const unsigned int libmesh_dbg_var (i),const Point &)34 Real FE<0,SZABAB>::shape(const ElemType,
35 const Order,
36 const unsigned int libmesh_dbg_var(i),
37 const Point &)
38 {
39 libmesh_assert_less (i, 1);
40 return 1.;
41 }
42
43
44
45 template <>
shape(const Elem *,const Order,const unsigned int libmesh_dbg_var (i),const Point &,const bool)46 Real FE<0,SZABAB>::shape(const Elem *,
47 const Order,
48 const unsigned int libmesh_dbg_var(i),
49 const Point &,
50 const bool)
51 {
52 libmesh_assert_less (i, 1);
53 return 1.;
54 }
55
56
57 template <>
shape(const FEType,const Elem *,const unsigned int libmesh_dbg_var (i),const Point &,const bool)58 Real FE<0,SZABAB>::shape(const FEType,
59 const Elem *,
60 const unsigned int libmesh_dbg_var(i),
61 const Point &,
62 const bool)
63 {
64 libmesh_assert_less (i, 1);
65 return 1.;
66 }
67
68
69
70
71 template <>
shape_deriv(const ElemType,const Order,const unsigned int,const unsigned int,const Point &)72 Real FE<0,SZABAB>::shape_deriv(const ElemType,
73 const Order,
74 const unsigned int,
75 const unsigned int,
76 const Point &)
77 {
78 libmesh_error_msg("No spatial derivatives in 0D!");
79 return 0.;
80 }
81
82
83
84 template <>
shape_deriv(const Elem *,const Order,const unsigned int,const unsigned int,const Point &,const bool)85 Real FE<0,SZABAB>::shape_deriv(const Elem *,
86 const Order,
87 const unsigned int,
88 const unsigned int,
89 const Point &,
90 const bool)
91 {
92 libmesh_error_msg("No spatial derivatives in 0D!");
93 return 0.;
94 }
95
96
97 template <>
shape_deriv(const FEType,const Elem *,const unsigned int,const unsigned int,const Point &,const bool)98 Real FE<0,SZABAB>::shape_deriv(const FEType,
99 const Elem *,
100 const unsigned int,
101 const unsigned int,
102 const Point &,
103 const bool)
104 {
105 libmesh_error_msg("No spatial derivatives in 0D!");
106 return 0.;
107 }
108
109
110
111 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES
112
113 template <>
shape_second_deriv(const ElemType,const Order,const unsigned int,const unsigned int,const Point &)114 Real FE<0,SZABAB>::shape_second_deriv(const ElemType,
115 const Order,
116 const unsigned int,
117 const unsigned int,
118 const Point &)
119 {
120 libmesh_error_msg("No spatial derivatives in 0D!");
121 return 0.;
122 }
123
124
125
126 template <>
shape_second_deriv(const Elem *,const Order,const unsigned int,const unsigned int,const Point &,const bool)127 Real FE<0,SZABAB>::shape_second_deriv(const Elem *,
128 const Order,
129 const unsigned int,
130 const unsigned int,
131 const Point &,
132 const bool)
133 {
134 libmesh_error_msg("No spatial derivatives in 0D!");
135 return 0.;
136 }
137
138
139
140 template <>
shape_second_deriv(const FEType,const Elem *,const unsigned int,const unsigned int,const Point &,const bool)141 Real FE<0,SZABAB>::shape_second_deriv(const FEType,
142 const Elem *,
143 const unsigned int,
144 const unsigned int,
145 const Point &,
146 const bool)
147 {
148 libmesh_error_msg("No spatial derivatives in 0D!");
149 return 0.;
150 }
151
152 #endif
153
154 } // namespace libMesh
155