Make Peierls approximation more efficient

[bobmyhill]

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.