FIx a bug (for non-standard dark energy) and optimize hubbleFunction.
Implement angularDiameterDistance for non-flat geometries (not well tested yet).