Flow123d
JB_transport-9331eee
|
Go to the documentation of this file.
40 template<
unsigned int spacedim>
50 template<
unsigned int spacedim>
53 unsigned int first_component_idx,
57 for (
unsigned int ip=0; ip<
n_points; ip++)
58 for (
unsigned int id=0;
id<dof_indices.size();
id++)
67 template<
unsigned int spacedim>
68 template<
unsigned int DIM>
88 ASSERT(fe_sys !=
nullptr).error(
"Mixed system must be represented by FESystem.");
93 unsigned int comp_offset = 0;
94 for (
auto fe : fe_sys->
fe())
110 ASSERT(
false).error(
"Not implemented.");
122 template<
unsigned int spacedim>
130 template<
unsigned int spacedim>
135 template<
unsigned int spacedim>
136 template<
unsigned int DIM>
145 allocate( q.
size(), _fe, _flags);
146 elm_values = std::make_shared<ElementValues<spacedim> >(q, update_flags, DIM);
152 ASSERT(fe !=
nullptr).error(
"Mixed system must be represented by FESystem.");
154 fe_values_vec.resize(fe->
fe().size());
155 for (
unsigned int f=0; f<fe->
fe().size(); f++)
156 fe_values_vec[f].initialize(q, *fe->
fe()[f], update_flags);
160 if ( q.dim() == DIM )
162 fe_data = init_fe_data(_fe, q);
164 else if ( q.dim() + 1 == DIM )
167 for (
unsigned int sid = 0; sid < RefElement<DIM>::n_sides; sid++)
169 side_fe_data[sid] = init_fe_data(_fe, q.make_from_side<DIM>(sid));
173 ASSERT(
false)(q.dim())(DIM).error(
"Dimension mismatch in FEValues::initialize().");
178 template<
unsigned int spacedim>
179 template<
unsigned int DIM>
181 unsigned int n_points,
188 ASSERT(_fe.
n_components() == spacedim).error(
"FEVector must have spacedim components.");
190 ASSERT(_fe.
n_components() == spacedim*spacedim).error(
"FETensor must have spacedim*spacedim components.");
193 fe_sys_dofs_.clear();
194 fe_sys_n_components_.clear();
195 fe_sys_n_space_components_.clear();
198 n_points_ = n_points;
201 fe_type_ = _fe.
type_;
203 if (fe_sys !=
nullptr)
205 for (
unsigned int f=0; f<fe_sys->
fe().size(); f++)
207 fe_sys_dofs_.push_back(fe_sys->
fe_dofs(f));
208 fe_sys_n_components_.push_back(fe_sys->
fe()[f]->n_components());
209 fe_sys_n_space_components_.push_back(fe_sys->
fe()[f]->n_space_components(spacedim));
217 shape_values.resize(n_points_,
vector<double>(n_dofs_*n_components_));
220 shape_gradients.resize(n_points_,
vector<arma::vec::fixed<spacedim> >(n_dofs_*n_components_));
222 views_cache_.initialize(*
this, _fe);
227 template<
unsigned int spacedim>
228 template<
unsigned int DIM>
233 std::shared_ptr<FEInternalData> data = std::make_shared<FEInternalData>(q.
size(), n_dofs_);
236 for (
unsigned int i=0; i<q.
size(); i++)
238 for (
unsigned int j=0; j<n_dofs_; j++)
243 data->ref_shape_values[i][j] = trans(shape_values.row(j));
248 for (
unsigned int i=0; i<q.
size(); i++)
250 for (
unsigned int j=0; j<n_dofs_; j++)
256 data->ref_shape_grads[i][j] = grad;
276 template<
unsigned int spacedim>
278 const unsigned int point_no,
279 const unsigned int comp)
const
284 return shape_gradients[point_no][function_no*n_components_+comp];
493 template<
unsigned int spacedim>
498 this->fill_data_specialized<MapScalar<spacedim>>(elm_values, fe_data);
501 this->fill_data_specialized<MapVector<spacedim>>(elm_values, fe_data);
504 this->fill_data_specialized<MapContravariant<spacedim>>(elm_values, fe_data);
507 this->fill_data_specialized<MapPiola<spacedim>>(elm_values, fe_data);
510 this->fill_data_specialized<MapTensor<spacedim>>(elm_values, fe_data);
513 this->fill_data_specialized<MapSystem<spacedim>>(elm_values, fe_data);
522 template<
unsigned int spacedim>
523 template<
class MapType>
526 map_type.fill_values_vec(*
this, elm_values, fe_data);
528 map_type.update_values(*
this, elm_values, fe_data);
530 map_type.update_gradients(*
this, elm_values, fe_data);
536 template<
unsigned int spacedim>
541 if (!elm_values->cell().is_valid() ||
542 elm_values->cell() != cell)
544 elm_values->reinit(cell);
547 fill_data(*elm_values, *fe_data);
551 template<
unsigned int spacedim>
556 if (!elm_values->side().is_valid() ||
557 elm_values->side() != cell_side)
559 elm_values->reinit(cell_side);
565 fill_data(*elm_values, *(side_fe_data[sid]) );
576 fv[0].initialize(quadrature[0], *fe[0_d], flags);
577 fv[1].initialize(quadrature[1], *fe[1_d], flags);
578 fv[2].initialize(quadrature[2], *fe[2_d], flags);
579 fv[3].initialize(quadrature[3], *fe[3_d], flags);
~FEValues()
Correct deallocation of objects created by 'initialize' methods.
@ update_gradients
Shape function gradients.
Helper class allows update values and gradients of FEValues of FEVectorPiola type.
unsigned int n_points() const
Returns the number of quadrature points.
arma::vec::fixed< dim > shape_grad(const unsigned int i, const arma::vec::fixed< dim > &p, const unsigned int comp=0) const
Calculates the comp-th component of the gradient of the i-th shape function at the point p on the ref...
FEValues()
Default constructor with postponed initialization.
static UpdateFlags update_each(UpdateFlags flags)
Determines which additional quantities have to be computed.
Class FEValues calculates finite element data on the actual cells such as shape function values,...
@ update_values
Shape function values.
unsigned int n_dofs_
Number of finite element dofs.
unsigned int size() const
Returns number of quadrature points.
void initialize(Quadrature &_quadrature, FiniteElement< DIM > &_fe, UpdateFlags _flags)
Initialize structures and calculates cell-independent data.
unsigned int dim() const
Returns dimension of the side, that is dimension of the element minus one.
void initialize(const FEValues &fv, const FiniteElement< DIM > &fe)
Class FESystem for compound finite elements.
#define ASSERT_LT(a, b)
Definition of comparative assert macro (Less Than) only for debug mode.
Helper class allows update values and gradients of FEValues of FEVector type.
Calculates finite element data on the actual cell.
unsigned int dim_
Dimension of reference space.
unsigned int n_components() const
Returns numer of components of the basis function.
std::vector< FEValues< 3 > > mixed_fe_values(QGauss::array &quadrature, MixedPtr< FiniteElement > fe, UpdateFlags flags)
void fill_data(const ElementValues< spacedim > &elm_values, const FEInternalData &fe_data)
Computes the shape function values and gradients on the actual cell and fills the FEValues structure.
const std::vector< std::shared_ptr< FiniteElement< dim > > > & fe() const
virtual UpdateFlags update_each(UpdateFlags flags)
Decides which additional quantities have to be computed for each cell.
Structure for storing the precomputed finite element data.
ArmaMat< double, N, M > mat
#define ASSERT_PERMANENT(expr)
Allow use shorter versions of macro names if these names is not used with external library.
FEType type_
Type of FiniteElement.
std::vector< std::vector< arma::mat > > ref_shape_grads
Precomputed gradients of basis functions at the quadrature points.
Abstract class for description of finite elements.
Helper class allows update values and gradients of FEValues of FETensor type.
std::array< QGauss, 4 > array
#define ASSERT_EQ(a, b)
Definition of comparative assert macro (EQual) only for debug mode.
unsigned int side_idx() const
Returns local index of the side on the element.
Helper class allows update values and gradients of FEValues of FEMixedSystem type.
Abstract class for the description of a general finite element on a reference simplex in dim dimensio...
unsigned int n_dofs() const
Returns the number of degrees of freedom needed by the finite element.
Class for computation of data on cell and side.
std::vector< std::vector< arma::vec > > ref_shape_values
Precomputed values of basis functions at the quadrature points.
unsigned int n_points_
Number of integration points.
int LongIdx
Define type that represents indices of large arrays (elements, nodes, dofs etc.)
Helper class allows update values and gradients of FEValues of FEScalar type.
Armor::ArmaVec< double, point_dim > point(unsigned int i) const
Returns the ith quadrature point.
std::vector< unsigned int > fe_dofs(unsigned int fe_index)
Return dof indices belonging to given sub-FE.
void allocate(unsigned int n_points, FiniteElement< DIM > &_fe, UpdateFlags flags)
Allocates space for computed data.
Basic definitions of numerical quadrature rules.
void fill_data_specialized(const ElementValues< spacedim > &elm_values, const FEInternalData &fe_data)
Computes the shape function values and gradients on the actual cell and fills the FEValues structure....
arma::vec::fixed< spacedim > shape_grad_component(const unsigned int function_no, const unsigned int point_no, const unsigned int comp) const
Return the gradient of the function_no-th shape function at the point_no-th quadrature point.
Class ElementValues calculates data related to transformation of reference cell to actual cell (Jacob...
UpdateFlags
Enum type UpdateFlags indicates which quantities are to be recomputed on each finite element cell.
Class MappingP1 implements the affine transformation of the unit cell onto the actual cell.
void reinit(const ElementAccessor< spacedim > &cell)
Update cell-dependent data (gradients, Jacobians etc.)
Base class for quadrature rules on simplices in arbitrary dimensions.
Compound finite element on dim dimensional simplex.
std::shared_ptr< FEInternalData > init_fe_data(const FiniteElement< DIM > &fe, const Quadrature &q)
Precompute finite element data on reference element.
Helper class allows update values and gradients of FEValues of FEVectorContravariant type.
unsigned int n_space_components(unsigned int spacedim)
Number of components of FE in a mapped space with dimension spacedim.
FEInternalData(unsigned int np, unsigned int nd)
unsigned int n_points
Number of quadrature points.
double shape_value(const unsigned int i, const arma::vec::fixed< dim > &p, const unsigned int comp=0) const
Calculates the value of the comp-th component of the i-th shape function at the point p on the refere...