Switched order of axes in probability surfaces to (EBV, DM). This
should result in better caching, since probability surfaces are
scanned primarily along the DM axis.
int casting in l.o.s. integrals done only once
Re-arranging loop order in l.o.s. integral to improve caching.
Removing matrix bounds checks in l.o.s. integrals. Performing
checks beforehand in ln(p) functions.
Strengthening prior on E(B-V) in piecewise-linear model
Merge branch 'improved-caching' into develop