[GeoMechanicsApplication] Performance issue in D-GeoFlow

[avdg81]

It appears that the performance of D-GeoFlow has decreased significantly. Walter has prepared a set of input and output files that should demonstrate the issue. I haven't attempted to reproduce it yet. But to avoid that the files will get lost (they were stored in my postbox), I have created this issue.

Below you'll find Walter's original message (in Dutch).

Things that need to be done:

• Ask Walter the date (or commit hash) of the "old" Kratos installation.
• Build the "old" Kratos package (using release build).
• Profile the code by running the simulation (using Intel's VTune).
• Build the current Kratos package (again using release build).
• Again profile the code by running the simulation (using Intel's VTune).
• Compare the performance reports to see where the hot spots have changed.
• Create a new issue to address the performance bottlenecks.
• Document how to profile code using VTune (try to keep it practical).

---

Hoi Anne,

Volgens mij had Ray je hier ook al eens kort over bij gepraat. Zouden jullie eens willen kijken naar de volgende files:
N:\Deltabox\Postbox\Graaf,van de Anne\vanWalter\DGeoFlowCriticalHead

Deze D-Geo Flow / critical head berekening, duurt in de laatste release ongeveer 40 seconden. Wanneer ik draai met de meest recente Kratos file, duurt deze ongeveer 1 minuut langer. Aangezien er verder geen functionele verandering zijn geweest, blijven we voor de komende release werken met de oude kernel. Het zou wel fijn zijn als jullie eens kunnen uitzoeken waar dit snelheidsverschil door komt en of daar mogelijk iets aan gedaan moet worden. Echter (misschien ten overvloede), werk aan D-Setttlement heeft wat ons betreft altijd prioriteit boven dit issue.

Alvast bedankt voor uitzoeken en mochten er nog vragen over zijn, dan hoor ik dat graag.

Groeten,
Walter

---

We have been able to improve the performance of the supplied benchmark by:

• Calling reserve to set the required capacity of vectors (to avoid reallocating memory).
• Using noalias for some assignments.
• Using fixed-size matrices rather than dynamic matrices.
• Not copying smart pointers when there is no need to participate in the ownership of the guarded object.

[avdg81]

DGeoFlowCriticalHead.zip