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 should probably be similar as in step-8.
Best, Simon On Wednesday, June 10, 2020 at 8:29:03 AM UTC+2, A.Z Ihsan wrote: > > Hi all, > > i am getting this error while calculating the local cell matrix for my hp > fem application > > dealii::Tensor<1, spacedim> dealii::FEValuesBase<dim, >>> spacedim>::shape_grad_component(unsigned int, unsigned int, unsigned int) >>> const [with int dim = 3; int spacedim = 3] >> >> The violated condition was: >> >> component < fe->n_components() >> >> Additional information: >> >> Index 1 is not in the half-open range [0,1). >> >> Stacktrace: >> >> ----------- >> >> #0 ../build/local: dealii::FEValuesBase<3, >>> 3>::shape_grad_component(unsigned int, unsigned int, unsigned int) const >> >> #1 ../build/local: dealii::SymmetricTensor<2, 3, double> >>> LinearElasticity::get_strain<3>(dealii::FEValues<3, 3> const&, unsigned >>> int, unsigned int) >> >> #2 ../build/local: LinearElasticity::HPSolver<3, dealii::Vector<double> >>> >::assemble_cell_matrix(dealii::TriaActiveIterator<dealii::DoFCellAccessor<dealii::hp::DoFHandler<3, >>> > >>> 3>, false> > const&, dealii::FullMatrix<double>&, dealii::hp::FEValues<3, >>> 3>&) >> >> #3 ../build/local: >>> LinearElasticity::HPSerialSolver<3>::assemble_linear_system(LinearElasticity::HPSerialSolver<3>::LinearSystem&) >> >> #4 ../build/local: LinearElasticity::HPSerialSolver<3>::solve_problem() >> >> >> >> > The idea is similar to the tutorial step-27, but here i use the symmetric > tensor for returning strain. Can someone explain to me what the error is > and how to resolve it? > > BR, > Ihsan > -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/9e6af2f7-695c-4a4f-bcd6-17c9a4f008cao%40googlegroups.com.