There should be an option to compute the first n kinetic coefficients simultaneously, especially since the calculations share a lot of work. Setting a tolerance may also become more subtle in this case. This may also benefit auxiliary integration, since A_0 seems to be a good auxiliary integrand for A_1, although some investigation of the adaptive grids needed for each one may be useful.
There should be an option to compute the first n kinetic coefficients simultaneously, especially since the calculations share a lot of work. Setting a tolerance may also become more subtle in this case. This may also benefit auxiliary integration, since A_0 seems to be a good auxiliary integrand for A_1, although some investigation of the adaptive grids needed for each one may be useful.