RESOLVING THIN BOUNDARY LAYERS IN NUMERICAL QUADRATURE ê
S. Adam, A. Tifrea, A. Neacsu · 2006
The identification of a library code able to solve numerically parametric integrals critically depends on the correctness of the diagnostic concerning the integrand behaviour at the integration domain boundaries. A second degree polynomial least squares fit of data over a minimal integrand sampling inside a mesoscopic boundary layer yields the necessary diagnostics provided solutions to the difficulties following from the lack of hardware and/or software conformity to the IEEE 754 standard. Relevant numerical examples let us infer that the developed code is robust, reliable, efficient, and provides precise outputs.