SheafSystem  0.0.0.0
uniform_3d.h
1 
2 //
3 // Copyright (c) 2014 Limit Point Systems, Inc.
4 //
5 // Licensed under the Apache License, Version 2.0 (the "License");
6 // you may not use this file except in compliance with the License.
7 // You may obtain a copy of the License at
8 //
9 // http://www.apache.org/licenses/LICENSE-2.0
10 //
11 // Unless required by applicable law or agreed to in writing, software
12 // distributed under the License is distributed on an "AS IS" BASIS,
13 // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 // See the License for the specific language governing permissions and
15 // limitations under the License.
16 //
17 
18 // Interface for class uniform_3d
19 
20 #ifndef UNIFORM_3D_H
21 #define UNIFORM_3D_H
22 
23 #ifndef SHEAF_DLL_SPEC_H
24 #include "SheafSystem/sheaf_dll_spec.h"
25 #endif
26 
27 #ifndef LINEAR_FCN_SPACE_H
28 #include "SheafSystem/linear_fcn_space.h"
29 #endif
30 
31 namespace fiber_bundle
32 {
33 
39 class SHEAF_DLL_SPEC uniform_3d : public linear_fcn_space
40 {
41 
42  // ===========================================================
44  // ===========================================================
46 
47 public:
48 
52  uniform_3d();
53 
57  uniform_3d(const uniform_3d& xother);
58 
62  virtual ~uniform_3d();
63 
64 protected:
65 
66 private:
67 
73  enum static_const_int {DL = 8};
74 
78  value_type _basis_value_buffer[DL];
79 
83  value_type _basis_deriv_value_buffer[3*DL];
84 
88  value_type _jacobian_value_buffer[128];
89 
91 
92  // ===========================================================
94  // ===========================================================
96 
97 public:
98 
102  virtual int dl() const;
103 
108  virtual void basis_at_coord(const dof_type xlocal_coord[],
109  size_type xlocal_coord_ub);
110 
115  virtual void basis_derivs_at_coord(const dof_type xlocal_coords[],
116  size_type xlocal_coords_ub);
117 
118 protected:
119 
120 private:
121 
123 
124  // ===========================================================
126  // ===========================================================
128 
129 public:
130 
135  virtual value_type jacobian_determinant(const dof_type xcoord_dofs[],
136  size_type xcoord_dofs_ub,
137  size_type xdf,
138  const coord_type xlocal_coords[],
139  size_type xlocal_coords_ub);
144  virtual value_type volume(const dof_type xcoord_dofs[],
145  size_type xcoord_dofs_ub,
146  size_type xdf);
147 
151  virtual void integrate(const dof_type xcoord_dofs[],
152  size_type xcoord_dofs_ub,
153  size_type xdf,
154  const dof_type xintegrands[],
155  size_type xintegrands_ub,
156  value_type xresult_integrals[],
157  size_type xresult_integrals_ub);
158 
162  virtual void integrate(const dof_type xcoord_dofs[],
163  size_type xcoord_dofs_ub,
164  size_type xdf,
165  const dof_type xintegrand,
166  value_type xresult_integrals[],
167  size_type xresult_integrals_ub);
168 
172  virtual void gauss_point(pod_index_type xindex,
173  coord_type xresult[],
174  size_type xresult_ub);
175 
176 protected:
177 
178 private:
179 
181 
182  // ===========================================================
184  // ===========================================================
186 
187 public:
188 
193  void dxi_local(size_type xlocal_coord_index,
194  const dof_type xsource_dofs[],
195  size_type xsource_dofs_ub,
196  dof_type xresult_dofs[],
197  size_type xresult_dofs_ub) const;
198 
203  virtual void jacobian(const dof_type xcoord_dofs[],
204  size_type xcoord_dofs_ub,
205  size_type xdf,
206  const dof_type xlocal_coords[],
207  size_type xlocal_coords_ub);
208 
209 protected:
210 
211 private:
212 
214 
215  // ===========================================================
217  // ===========================================================
219 
220 public:
221 
226  virtual int db() const;
227 
231  virtual void local_coordinates(pod_index_type xindex,
232  coord_type xresult[],
233  size_type xresult_ub) const;
238  virtual bool in_standard_domain(const dof_type xlocal_coords[], size_type xlocal_coords_ub) const;
239 
240 protected:
241 
242 private:
243 
245 
246  // ===========================================================
248  // ===========================================================
250 
251 public:
252 
258  virtual void coord_at_value(const dof_type xdofs[],
259  size_type xdofs_ub,
260  const dof_type xglobal_coords[],
261  size_type xglobal_coord_ub,
262  dof_type xlocal_coords[],
263  size_type xlocal_coords_ub) const;
264 
265 protected:
266 
267 private:
268 
270 
271  // ===========================================================
273  // ===========================================================
275 
276 public:
277 
281  virtual uniform_3d* clone() const;
282 
286  virtual uniform_3d& operator=(const section_evaluator& xother);
287 
291  uniform_3d& operator=(const uniform_3d& xother);
292 
296  virtual bool invariant() const;
297 
301  virtual bool is_ancestor_of(const any* xother) const;
302 
303 protected:
304 
305 private:
306 
308 
309 };
310 
311 } // namespace fiber_bundle
312 
313 #endif // ifndef UNIFORM_3D_H
A section evaluator using trilinear interpolation over a cubic 3D domain. Intended for use with unifo...
Definition: uniform_3d.h:39
sec_vd_dof_type dof_type
The type of degree of freedom.
Abstract base class with useful features for all objects.
Definition: any.h:39
An abstract integrable section evaluator which is a member of a linear function space.
unsigned long size_type
An unsigned integral type used to represent sizes and capacities.
Definition: sheaf.h:52
chart_point_coord_type coord_type
The type of local coordinate; the scalar type for the local coordinate vector space.
An abstract local section evaluator; a map from {local coordinates x dofs} to section value...
int_type pod_index_type
The plain old data index type.
Definition: pod_types.h:49
Namespace for the fiber_bundles component of the sheaf system.