29 typedef basis_rep<T> base;
33 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
38 typedef basis_rep<T> base;
42 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
47 typedef basis_rep<T> base;
51 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
56 typedef basis_rep<T> base;
60 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
65 typedef basis_rep<T> base;
69 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
74 typedef basis_rep<T> base;
78 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
83 typedef basis_rep<T> base;
87 static void hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>&);
93 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
100 basis_P0_p<T>::grad_evaluate (
108 basis_P0_p<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
117 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
124 basis_P0_e<T>::grad_evaluate (
133 basis_P0_e<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
142 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
149 basis_P0_t<T>::grad_evaluate (
159 basis_P0_t<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
168 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
175 basis_P0_q<T>::grad_evaluate (
185 basis_P0_q<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
194 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
201 basis_P0_T<T>::grad_evaluate (
212 basis_P0_T<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
215 x[0] =
point_basic<T>(2.5000000000000000e-01, 2.5000000000000000e-01, 2.5000000000000000e-01);
221 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
228 basis_P0_P<T>::grad_evaluate (
239 basis_P0_P<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
248 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
255 basis_P0_H<T>::grad_evaluate (
266 basis_P0_H<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
295 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
const
320 error_macro (
"evaluate: unsupported `" << hat_K.
name() <<
"' element type");
333 return basis_P0_p<T>::grad_evaluate (hat_x, values);
336 return basis_P0_e<T>::grad_evaluate (hat_x, values);
339 return basis_P0_t<T>::grad_evaluate (hat_x, values);
342 return basis_P0_q<T>::grad_evaluate (hat_x, values);
345 return basis_P0_T<T>::grad_evaluate (hat_x, values);
348 return basis_P0_P<T>::grad_evaluate (hat_x, values);
351 return basis_P0_H<T>::grad_evaluate (hat_x, values);
354 error_macro (
"grad_evaluate: unsupported `" << hat_K.
name() <<
"' element type");
362 const Eigen::Matrix<T,Eigen::Dynamic,1>& f_xnod,
363 Eigen::Matrix<T,Eigen::Dynamic,1>& dof)
const
373 base::is_continuous(),
374 base::_ndof_on_subgeo,
375 base::_nnod_on_subgeo,
376 base::_first_idof_by_dimension,
377 base::_first_inod_by_dimension);
386 return basis_P0_p<T>::hat_node (_hat_node[hat_K.
variant()]);
389 return basis_P0_e<T>::hat_node (_hat_node[hat_K.
variant()]);
392 return basis_P0_t<T>::hat_node (_hat_node[hat_K.
variant()]);
395 return basis_P0_q<T>::hat_node (_hat_node[hat_K.
variant()]);
398 return basis_P0_T<T>::hat_node (_hat_node[hat_K.
variant()]);
401 return basis_P0_P<T>::hat_node (_hat_node[hat_K.
variant()]);
404 return basis_P0_H<T>::hat_node (_hat_node[hat_K.
variant()]);
407 error_macro (
"hat_node: unsupported `" << hat_K.
name() <<
"' element type");
field::size_type size_type
base::size_type size_type
void _initialize_cstor_sizes() 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 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 _initialize_data(reference_element hat_K) const
static void initialize_local_first(size_type k, bool is_continuous, std::array< std::array< size_type, reference_element::max_variant >, 4 > &ndof_on_subgeo, std::array< std::array< size_type, reference_element::max_variant >, 4 > &nnod_on_subgeo, std::array< std::array< size_type, 5 >, reference_element::max_variant > &first_idof_by_dimension, std::array< std::array< size_type, 5 >, reference_element::max_variant > &first_inod_by_dimension)
see the basis_option page for the full documentation
void set_continuous(bool c=true)
piola_fem< T > _piola_fem
see the reference_element page for the full documentation
static const variant_type H
static const variant_type q
static const variant_type e
static const variant_type p
variant_type variant() const
static const variant_type T
static const variant_type P
static const variant_type t
#define error_macro(message)
This file is part of Rheolef.
void evaluate(const geo_basic< float_type, M > &omega_K, const geo_element &K, Eigen::Matrix< Result, Eigen::Dynamic, 1 > &value) const