[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread A.Z Ihsan
Great, it works now. THanks Simon. BR, Ihsan On Wednesday, June 10, 2020 at 12:50:30 PM UTC+2, Simon Sticko wrote: > > Hi, > Yes. Here, you add a number of scalar elements of different degree: > > > for (unsigned int deg=1; deg <= max_fe_degree; ++deg) > > { > >

[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread Simon Sticko
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

[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread A.Z Ihsan
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));

[deal.II] Re: Error calculating shape gradient

2020-06-10 Thread Simon Sticko
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