• Francesco Montanari's avatar
    Fix Simpson's rule · 4b8d3549
    Francesco Montanari authored
    * cosmo.el (cosmo--trapzd, cosmo-qsimp): Check Simpson's rule for
      convergence and stop if the desired accuracy is reached. The
      previous algorithm did't perform any convergence check. At low
      redshift z~1 numerical integrals are now about one order of
      magnitude faster, at large redshift z~1 one order of magnitude
      slower (but it is guaranteed to reach the desired accuracy).
    * cosmo.el (cosmo, cosmo-int-prec, cosmo-int-maxsteps): Define
      customization group to allow users to set numerical integrals
      precision variable through the `customize' command.
    * README.org: Explain how to customize variables.
cosmo.el 16.1 KB