Dynamic control of step size in a

[duetosymmetry]

Using curvature of omega(a) and A(a) can give estimates of delta_a which will potentially save computation time in smooth regions, and help to resolve fine features

[duetosymmetry]

Try using d^2(ln(omega))/da^2 to estimate curvature instead of d^2(omega)/da^2