25#include "rheolef/rheostream.h"
42 _scalar_basis(scalar_basis),
50 "vector(basis): basis.option.map_dimension should be initialized for component number");
51 _n_comp = base::option().dimension();
59 const size_type unset = std::numeric_limits<basis_option::size_type>::max();
60 return (base::option().
dimension() == unset) ? map_d : base::option().dimension();
67 for (
size_type map_d = 0; map_d < 4; ++map_d) {
69 base::_ndof_on_subgeo [map_d][subgeo_variant] = _n_comp*_scalar_basis.ndof_on_subgeo (map_d, subgeo_variant);
70 base::_nnod_on_subgeo [map_d][subgeo_variant] = _scalar_basis.nnod_on_subgeo (map_d, subgeo_variant);
75 for (
size_type subgeo_d = 0; subgeo_d < 5; ++subgeo_d) {
76 base::_first_idof_by_dimension [variant][subgeo_d] = _n_comp*_scalar_basis.first_idof_by_dimension (hat_K, subgeo_d);
77 base::_first_inod_by_dimension [variant][subgeo_d] = _scalar_basis.first_inod_by_dimension (hat_K, subgeo_d);
87const Eigen::Matrix<point_basic<T>,Eigen::Dynamic,1>&
90 return _scalar_basis.hat_node (hat_K);
99 base::_initialize_data_guard (hat_K);
100 _scalar_basis.evaluate (hat_K, hat_x, _scalar_value);
101 size_type loc_comp_ndof = _scalar_value.size();
102 size_type loc_ndof = _n_comp*loc_comp_ndof;
103 value.resize (loc_ndof);
105 for (
size_type loc_comp_idof = 0; loc_comp_idof < loc_comp_ndof; ++loc_comp_idof) {
106 for (
size_type i_comp = 0; i_comp < _n_comp; ++i_comp) {
107 size_type loc_idof = _n_comp*loc_comp_idof + i_comp;
108 value[loc_idof][i_comp] = _scalar_value[loc_comp_idof];
119 base::_initialize_data_guard (hat_K);
120 _scalar_basis.grad_evaluate (hat_K, hat_x, _vector_value);
121 size_type loc_comp_ndof = _vector_value.size();
122 size_type loc_ndof = _n_comp*loc_comp_ndof;
123 value.resize (loc_ndof);
125 for (
size_type loc_comp_idof = 0; loc_comp_idof < loc_comp_ndof; ++loc_comp_idof) {
126 for (
size_type i_comp = 0; i_comp < _n_comp; ++i_comp) {
127 size_type loc_idof = _n_comp*loc_comp_idof + i_comp;
128 for (
size_type j_comp = 0; j_comp < _n_comp; ++j_comp) {
129 value[loc_idof] (i_comp,j_comp) = _vector_value[loc_comp_idof][j_comp];
139 Eigen::Matrix<T,Eigen::Dynamic,1>& dof)
const
141 base::_initialize_data_guard (hat_K);
142 size_type loc_comp_ndof = _scalar_basis.ndof (hat_K);
143 size_type loc_comp_nnod = _scalar_basis.nnod (hat_K);
144 size_type loc_ndof = _n_comp*loc_comp_ndof;
145 Eigen::Matrix<T,Eigen::Dynamic,1> f_comp_xnod (loc_comp_nnod);
146 Eigen::Matrix<T,Eigen::Dynamic,1> comp_dof (loc_comp_ndof);
147 dof.resize (loc_ndof);
148 for (
size_type i_comp = 0; i_comp < _n_comp; ++i_comp) {
149 for (
size_type loc_comp_inod = 0; loc_comp_inod < loc_comp_nnod; ++loc_comp_inod) {
150 f_comp_xnod [loc_comp_inod] = f_xnod [loc_comp_inod] [i_comp];
152 _scalar_basis.compute_dofs (hat_K, f_comp_xnod, comp_dof);
153 for (
size_type loc_comp_idof = 0; loc_comp_idof < loc_comp_ndof; ++loc_comp_idof) {
154 size_type loc_idof = _n_comp*loc_comp_idof + i_comp;
155 dof [loc_idof] = comp_dof [loc_comp_idof];
162#define _RHEOLEF_instanciation(T) \
163template class basis_fem_vector<T>;
#define _RHEOLEF_instanciation(T, M, A)
field::size_type size_type
see the Float page for the full documentation
basis_fem_vector(const basis_basic< T > &scalar_basis, const basis_option &sopt)
void grad_evaluate(reference_element hat_K, const point_basic< T > &hat_x, Eigen::Matrix< tensor_basic< T >, Eigen::Dynamic, 1 > &value) const
base::size_type size_type
void _initialize_cstor_sizes() const
basis_basic< T > _scalar_basis
void evaluate(reference_element hat_K, const point_basic< T > &hat_x, Eigen::Matrix< point_basic< T >, Eigen::Dynamic, 1 > &value) const
void _compute_dofs(reference_element hat_K, const Eigen::Matrix< point_basic< T >, Eigen::Dynamic, 1 > &f_xnod, Eigen::Matrix< T, Eigen::Dynamic, 1 > &dof) const
size_type family_index() const
const Eigen::Matrix< point_basic< T >, Eigen::Dynamic, 1 > & hat_node(reference_element hat_K) const
std::string family_name() const
void _initialize_data(reference_element hat_K) const
see the basis_option page for the full documentation
see the reference_element page for the full documentation
static const variant_type max_variant
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
This file is part of Rheolef.