geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
217 stars 232 forks source link

Make Peierls approximation more efficient #5945

Open bobmyhill opened 5 days ago

bobmyhill commented 5 days ago

Most of the calculations in the approximate Peierls rheology are operating on constants, e.g.:

c_cutoff = std::pow(p.stress_cutoff/p.peierls_stress, p.glide_parameter_p);

Doing these calculations millions of times is inefficient and we should just store the derived parameters instead.