[deal.II] Finite Element Software deal.II Version 9.2.0 released

2020-06-10 Thread Matthias Maier
Version 9.2.0 of deal.II, the object-oriented finite element library awarded the J. H. Wilkinson Prize for Numerical Software, has been released. It is available for free under an Open Source license from the deal.II homepage at https://www.dealii.org/ The major changes of this

Re: [deal.II] Particle Class Implementation

2020-06-10 Thread Wolfgang Bangerth
On 6/9/20 9:05 PM, Bruno Blais wrote: In additional to what Jean-Paul suggested, you can look at the preliminary version of step-68 which does exactly what you would like to achieve with particles. The code is available on the following pull request. Rene and I have put some work into it and

Re: [deal.II] Particle Class Implementation

2020-06-10 Thread Victoria W.
Dear Bruno and Jean-Paul, The new syntax worked great, and thank you for the extra resources Bruno. Best, Victoria On Tuesday, June 9, 2020 at 11:05:50 PM UTC-4, Bruno Blais wrote: > > Dear Victoria, > > In additional to what Jean-Paul suggested, you can look at the preliminary > version of

[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) > > { > > fe_collection.push_back(FE_Q(d

[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)); face_quadrature_collection.

[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 shou