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/83d6f2cc-b746-4cb0-bf43-7b0e1d86f3bbo%40googlegroups.com.