We evaluate an adaptive gaussian quadrature integration scheme that will be suitable for the numerical evaluation of generalized redistribution in frequency functions. The latter are indispensable ingredients for "full non-LTE" radiation transfer computations i.e., assuming potential deviations of the velocity distribution of massive particles from the usual Maxwell-Boltzmann distribution. A first validation is made with computations of the usual Voigt profile.