GalSim-developers/GalSim — Issue #2: Replacing Simpson's rule integration from SBProfile Simpson.h
Author: rmandelb
Created: 2012-02-24 16:28:10+00:00
State: open

Currently Simpson.h has a NR-based routine for integration by Simpson's rule, which Gary modified so it's written as a template with the integrand as the template argument. We should replace this with some alternative, such as Mike's integration package. Gary points out the decision to be made: whether to adopt the templating interface, or to change sbprofile to do more traditional function calls. The advantage of the template is that the integrand can be a class instead of a function, which allows better control of internal variables.