1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
// file automatically generated by "../../../rheolef/fem/lib/basis_symbolic_cxx.cc"
///
/// This file is part of Rheolef.
///
/// Copyright (C) 2000-2009 Pierre Saramito <Pierre.Saramito@imag.fr>
///
/// Rheolef is free software; you can redistribute it and/or modify
/// it under the terms of the GNU General Public License as published by
/// the Free Software Foundation; either version 2 of the License, or
/// (at your option) any later version.
///
/// Rheolef is distributed in the hope that it will be useful,
/// but WITHOUT ANY WARRANTY; without even the implied warranty of
/// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
/// GNU General Public License for more details.
///
/// You should have received a copy of the GNU General Public License
/// along with Rheolef; if not, write to the Free Software
/// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
///
/// =========================================================================
#ifndef _RHEOLEF_bubble_H
#define _RHEOLEF_bubble_H
#include "rheolef/basis.h"
#include "rheolef/tensor.h"
#include "Pk_get_local_idof_on_side.icc"
#include "basis_fem_Pk_lagrange.h"
namespace rheolef {
template<class T>
class basis_bubble: public basis_rep<T> {
public:
typedef basis_rep<T> base;
typedef typename base::size_type size_type;
basis_bubble();
~basis_bubble();
std::string name() const { return family_name(); }
bool have_index_parameter() const { return false; }
std::string family_name() const { return "bubble"; }
bool have_compact_support_inside_element() const { return true; }
size_type degree() const;
bool is_nodal() const { return true; }
void evaluate (
reference_element hat_K,
const point_basic<T>& hat_x,
Eigen::Matrix<T,Eigen::Dynamic,1>& values) const;
void grad_evaluate (
reference_element hat_K,
const point_basic<T>& hat_x,
Eigen::Matrix<point_basic<T>,Eigen::Dynamic,1>& values) const;
void _compute_dofs (
reference_element hat_K,
const Eigen::Matrix<T,Eigen::Dynamic,1>& f_xnod,
Eigen::Matrix<T,Eigen::Dynamic,1>& dof) const;
void _initialize_cstor_sizes() const;
void _initialize_data (reference_element hat_K) const;
const Eigen::Matrix<point_basic<T>,Eigen::Dynamic,1>& hat_node(reference_element hat_K) const {
base::_initialize_data_guard(hat_K);
return _hat_node [hat_K.variant()]; }
mutable std::array<
Eigen::Matrix<point_basic<T>,Eigen::Dynamic,1>,
reference_element::max_variant> _hat_node;
};
} // namespace rheolef
#endif // _RHEOLEF_bubble_H
|
