Compute averages in some directions.

Author

Michel Mehrenberger

Generic function which can be used for computing charge density. By default we suppose a uniform mesh node based. Trapezoid quadrature formula is used. We can also integrate by providing a function whose signature is sll_integration_discrete_1d.