|
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]
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.
