Why do I have a scalar base element when initializing FE_Q<dim> with dim=3?
I thought that I then have the finite element as 3d-element?
Am Donnerstag, 3. August 2017 16:32:03 UTC+2 schrieb Wolfgang Bangerth:
>
>
> >
> > |
> > template<intdim>
> >
> >
> std::vector<constFiniteElement<dim>*>Diffusion<dim>::create_fe_list(constunsignedintpolynomial_degree)
>
>
> > {
> > std::vector<constFiniteElement<dim>*>fe_list;
> > fe_list.push_back(newFE_Q<dim>(polynomial_degree));
> > returnfe_list;
>
> So you have one scalar base element...
>
>
> > template<intdim>
> > std::vector<unsignedint>Diffusion<dim>::create_fe_multiplicities()
> > {
> > std::vector<unsignedint>multiplicities;
> > multiplicities.push_back(3);
>
> ...with multiplicity three, for a total of three components. As a
> consequence...
>
> > constFEValuesExtractors::VectorE_temperature(dim);
> > constFEValuesExtractors::VectorL_temperature(2*dim);
>
> ...using these fill have to fail if dim>=2. Indeed, this happens:
>
>
> > --------------------------------------------------------
> > Anerror occurred inline <3996>of file
> > </opt/dealII/include/deal.II/fe/fe_values.h>infunction
> >
> constdealii::FEValuesViews::Vector<dim,spacedim>&dealii::FEValuesBase<dim,spacedim>::operator[](constdealii::FEValuesExtractors::Vector&)const[withintdim
>
>
> > =3;intspacedim =3]
> > Theviolated condition was:
> > vector.first_vector_component <fe_values_views_cache.vectors.size()
> > Additionalinformation:
> > Index3isnotinthe half-open range [0,1).
>
> Best
> W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> <javascript:>
> www: http://www.math.colostate.edu/~bangerth/
>
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
