Find an adaptive refinement algorithm

[duetosymmetry]

We should be able to subdivide intervals in a where we think the functions \omega(a) or A(a) are varying on smaller scales than currently sampled. See e.g. this discussion where graph refinement is discussed. Here's another discussion I found on SE, which points to this implementation. There is another implementation in axch/adaptive-plot.