Hello everyone. I am trying to solve elasticity and laplace equations one after the other. The elasticity solution feeds as input to the laplace equation (stiffness matrix of laplace equation depends on elasticity solution) and vice versa. I am using 65536 elements. For step-8 the assembly takes very less time (around 0.15second) while for my assemble_elastic, it takes around 5 seconds. The only difference between my assemble_elastic function and the assemble function of step-8 is that for each Gauss point, I additionally call a function(H_plus) that takes the laplace solution, the current cell and Gauss point as input and evaluates some quantity using this information. The H_plus function is called 4*65536 times but the function is very simple. My doubt is whether such a huge increase in cost (from 0.15 sec to 5 sec) is expected for this problem or is there something that I am doing that is increasing the cost so much.
I am running the code in release mode. Thanks and regards Wasim -- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/7ff1d083-ef44-43f0-b74f-99634a9ab851n%40googlegroups.com.
