Varying all the parameters

Just for completeness, Here is a run where all the parameters are varied. This was 40k points, run in two separate chains. I found the bug. The minimum value of $\sigma_8$ that I was using was wrong, so every once in a while the chain would run out of range and the code would crash. I can now, in principle, run arbitrarily long chains (which is more accurate than running a short chain 15 times).



So, what we see here is that this looks esentially like the plot where $\Omega_b h^2$ is fixed, which is consistent with what we were seeing before, $\Omega_b h^2$ is less important than the other parameters in terms of it's effect on the $\Omega_m h^2 - \sigma_8$ plane. I think this is because the major driver for the boundaries of the degeneracy range are the cutoff values of k, whereas $\Omega_b h^2$ is controlling the relative values of the peak heights. It is probably true that the points which first exceed the $1\sigma$ threshold are either at one end or the other of the allowed range in k. I think this is also the reason that the degenerate region has straight edges, but we will have to see.