Re: [deal.II] Re: Issue regarding higher order element
Thanks a lot Bruno. I did find the issue following your suggestion! On Tuesday, 16 July 2019 15:10:40 UTC-5, Bruno Turcksin wrote: > > Krishanu > > Le mar. 16 juil. 2019 à 16:00, Krishanu Sen > a écrit : > > > > I am using deal.ii for a nonlinear elasticity problem. And the change in > the code that is causing this error message is: > > > > > > template > > ElasticProblem::ElasticProblem () > > : > > dof_handler (triangulation), > > fe (FE_Q(3), dim) > > {} > No it's not where the error comes from. It is from somewhere in > Step8::ElasticProblem<1>::assemble_system_matrix() You need to find > which line produces the error. > > Best, > > Bruno > -- 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/4bd717a4-738b-4c8b-a1da-b95d2042eb96%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Re: Issue regarding higher order element
Thanks Prof. Bangerth for the suggestion. I will definitely learn about the debugger as it seems like a great tool. Thanks, Krishanu On Thursday, 18 July 2019 03:42:54 UTC-5, Wolfga ng Bangerth wrote: > > On 7/16/19 2:10 PM, Bruno Turcksin wrote: > > No it's not where the error comes from. It is from somewhere in > > Step8::ElasticProblem<1>::assemble_system_matrix() You need to find > > which line produces the error. > > Krishanu, > indeed, Bruno knows this because he is looking at the backtrace shown in > the > error message you posted: > > #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 > > The easiest way to find these errors is to run the program in a debugger. > It > will stop in exactly the place where the problem is, and allow you to > inspect > all local variables. It's a tool so worth learning how to use! (There are > also > video lectures showing how to do that in eclipse.) > > Best > W. > > -- > > Wolfgang Bangerth email: bang...@colostate.edu > > www: http://www.math.colostate.edu/~bangerth/ > > -- 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/102a6b2e-50cc-4d09-86b4-a3195fbf6c06%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Re: Issue regarding higher order element
On 7/16/19 2:10 PM, Bruno Turcksin wrote: > No it's not where the error comes from. It is from somewhere in > Step8::ElasticProblem<1>::assemble_system_matrix() You need to find > which line produces the error. Krishanu, indeed, Bruno knows this because he is looking at the backtrace shown in the error message you posted: #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 The easiest way to find these errors is to run the program in a debugger. It will stop in exactly the place where the problem is, and allow you to inspect all local variables. It's a tool so worth learning how to use! (There are also video lectures showing how to do that in eclipse.) Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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/165fa819-3baf-0faf-cabc-5206b54fb1bf%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
Re: [deal.II] Re: Issue regarding higher order element
Krishanu Le mar. 16 juil. 2019 à 16:00, Krishanu Sen a écrit : > > I am using deal.ii for a nonlinear elasticity problem. And the change in the > code that is causing this error message is: > > > template > ElasticProblem::ElasticProblem () > : > dof_handler (triangulation), > fe (FE_Q(3), dim) > {} No it's not where the error comes from. It is from somewhere in Step8::ElasticProblem<1>::assemble_system_matrix() You need to find which line produces the error. Best, Bruno -- 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/CAGVt9eOaE9g2sehRfaT-9B8JVq9pE79jdOV5uqMBF3LFzCQuaw%40mail.gmail.com. For more options, visit https://groups.google.com/d/optout.
[deal.II] Re: Issue regarding higher order element
Hello Bruno, Here is the error message that I am getting: An error occurred in line <388> of file in function dealii::Tensor::value_type& dealii::Tensor::operator[](unsigned int) [with int rank_ = 2; int dim = 3; Number = double; dealii::Tensor::value_type = dealii::Tensor<1, 3, double>] The violated condition was: i::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 ElasticProblem::ElasticProblem () : dof_handler (triangulation), fe (FE_Q(3), dim) {} The code runs fine when I am using FE_Q(2) or FE_Q(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(1)] or quadratic [FE_Q(2)] elements. >> But when I tried to use cubic elements [FE_Q(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.
[deal.II] Re: Issue regarding higher order element
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(1)] or quadratic [FE_Q(2)] elements. > But when I tried to use cubic elements [FE_Q(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/4f277d7a-396b-43f7-b887-32c850ae6944%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.