Hello Bruno,
Here is the error message that I am getting:
--------------------------------------------------------
An error occurred in line <388> of file
</g/a-relliott/COMMON/senxx051/dealii_source/dealii_library/include/deal.II/base/tensor.h>
in function
dealii::Tensor<rank, dim, Number>::value_type& dealii::Tensor<rank,
dim, Number>::operator[](unsigned int) [with int rank_ = 2; int dim = 3;
Number = double; dealii::Tensor<rank, dim, Number>::value_type =
dealii::Tensor<1, 3, double>]
The violated condition was:
i<dim
The name and call sequence of the exception was:
ExcIndexRange(i, 0, dim)
Additional Information:
Index 3 is not in the half-open range [0,3).
Stacktrace:
-----------
#0 ./nodal_growth_dyn_full:
Step8::ElasticProblem<1>::assemble_system_matrix()
#1 ./nodal_growth_dyn_full: Step8::ElasticProblem<1>::run()
#2 ./nodal_growth_dyn_full: main
--------------------------------------------------------
CMakeFiles/run.dir/build.make:57: recipe for target 'CMakeFiles/run' failed
make[3]: *** [CMakeFiles/run] Aborted
CMakeFiles/Makefile2:232: recipe for target 'CMakeFiles/run.dir/all' failed
make[2]: *** [CMakeFiles/run.dir/all] Error 2
CMakeFiles/Makefile2:239: recipe for target 'CMakeFiles/run.dir/rule' failed
make[1]: *** [CMakeFiles/run.dir/rule] Error 2
Makefile:183: recipe for target 'run' failed
make: *** [run] Error 2
I am using deal.ii for a nonlinear elasticity problem. And the change in
the code that is causing this error message is:
template <int dim>
ElasticProblem<dim>::ElasticProblem ()
:
dof_handler (triangulation),
fe (FE_Q<dim>(3), dim)
{}
The code runs fine when I am using FE_Q<dim>(2) or FE_Q<dim>(1) instead.
Thanks,
Krishanu
On Tuesday, 16 July 2019 14:22:13 UTC-5, Bruno Turcksin wrote:
>
> Krishanu,
>
> It's impossible to help you with so little information. Somewhere you set
> something to the wrong size. Can you give us a backtrace? We need to see
> the error message and the code that produces it.
>
> Best,
>
> Bruno
>
> On Tuesday, July 16, 2019 at 3:08:42 PM UTC-4, Krishanu Sen wrote:
>>
>> I am trying to run a code using higher order elements. The code runs fine
>> when I am using linear [FE_Q<dim>(1)] or quadratic [FE_Q<dim>(2)] elements.
>> But when I tried to use cubic elements [FE_Q<dim>(3)], the code stops with
>> an error "ExcIndexRange", which says that the code is trying to access an
>> index of a tensor that exceeds the limit of the tensor indices. I am not
>> sure why would that happen just because of changing the polynomial degree
>> of the element. I would appreciate any help or insight regarding the issue.
>>
>> Thanks,
>> Krishanu
>>
>
--
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/0884e920-3ac7-4d63-84b5-d9993223e568%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
