Currently have order, which corresponds to the order of Lebedev quadrature, as a template argument to all of the different functions which calculate the singular potential. This was so that the Lebedev coordinates could be a static attribute, but that's unnecessary because we only use one object to calculate the singular potential anyways. Would make it much easier to specify parameters from the command line, and would also cut compilation time considerably.
Currently have
order
, which corresponds to the order of Lebedev quadrature, as a template argument to all of the different functions which calculate the singular potential. This was so that the Lebedev coordinates could be a static attribute, but that's unnecessary because we only use one object to calculate the singular potential anyways. Would make it much easier to specify parameters from the command line, and would also cut compilation time considerably.