The IntPuff model does not actually integrate the Gaussian Puff: It provides an interface to calculate the finite sum of puffs or uses an approximate integral.
An addition to this could be a keyword argument, which defaults to the approximate integral to maintain backwards compatibility, which allows the user to specify that actual numerical integration should happen.
This could take the form of pre-generating a set of points for Gaussian quadrature, and using that to calculate the integral on demand or just using quadgk with adaptive quadrature.
The IntPuff model does not actually integrate the Gaussian Puff: It provides an interface to calculate the finite sum of puffs or uses an approximate integral.
An addition to this could be a keyword argument, which defaults to the approximate integral to maintain backwards compatibility, which allows the user to specify that actual numerical integration should happen.
This could take the form of pre-generating a set of points for Gaussian quadrature, and using that to calculate the integral on demand or just using quadgk with adaptive quadrature.