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)
100basis_P0_p<T>::grad_evaluate (
108basis_P0_p<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
117 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
124basis_P0_e<T>::grad_evaluate (
133basis_P0_e<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
142 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
149basis_P0_t<T>::grad_evaluate (
159basis_P0_t<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
168 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
175basis_P0_q<T>::grad_evaluate (
185basis_P0_q<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
194 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
201basis_P0_T<T>::grad_evaluate (
212basis_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)
228basis_P0_P<T>::grad_evaluate (
239basis_P0_P<T>::hat_node (Eigen::Matrix<
point_basic<T>,Eigen::Dynamic,1>& x)
248 Eigen::Matrix<T,Eigen::Dynamic,1>& values)
255basis_P0_H<T>::grad_evaluate (
266basis_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_internal,
375 base::_ndof_on_subgeo,
376 base::_nnod_on_subgeo_internal,
377 base::_nnod_on_subgeo,
378 base::_first_idof_by_dimension_internal,
379 base::_first_idof_by_dimension,
380 base::_first_inod_by_dimension_internal,
381 base::_first_inod_by_dimension);
390 return basis_P0_p<T>::hat_node (_hat_node[hat_K.
variant()]);
393 return basis_P0_e<T>::hat_node (_hat_node[hat_K.
variant()]);
396 return basis_P0_t<T>::hat_node (_hat_node[hat_K.
variant()]);
399 return basis_P0_q<T>::hat_node (_hat_node[hat_K.
variant()]);
402 return basis_P0_T<T>::hat_node (_hat_node[hat_K.
variant()]);
405 return basis_P0_P<T>::hat_node (_hat_node[hat_K.
variant()]);
408 return basis_P0_H<T>::hat_node (_hat_node[hat_K.
variant()]);
411 error_macro (
"hat_node: unsupported `" << hat_K.
name() <<
"' element 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_internal, 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_internal, 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_internal, 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_internal, 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)
reference_element::size_type size_type
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