Hi,
Yes. Here, you add a number of scalar elements of different degree:
> for (unsigned int deg=1; deg <= max_fe_degree; ++deg)
> {
> fe_collection.push_back(FE_Q(deg));
> quadrature_collection.push_back(QGauss(deg+1));
> face_quadrature_collection.push_back(QGauss(deg+1));
> }
which makes sense
Hi Simon,
i followed the tutorial step-27 about hp-fem, using the similar constructor
with dim = 3
for (unsigned int deg=1; deg <= max_fe_degree; ++deg)
{
fe_collection.push_back(FE_Q(deg));
quadrature_collection.push_back(QGauss(deg+1));
Hi,
from the error message you can see that the element you are using only has
1 component. You get an error because you are trying to access component 1,
which doesn't exist. Since your element should have dim components, there
is likely something wrong with how you create your element. It
