Cuadratura Adaptativa

La cuadratura adaptativa subdivide automáticamente el intervalo de integración donde el integrando es difícil, utilizando estimaciones de error local para alcanzar la precisión solicitada con el menor número posible de evaluaciones de función.

Definition

La cuadratura adaptativa es cualquier estrategia de integración numérica que utiliza estimaciones del error de aproximación local para decidir dónde y con qué precisión subdividir el dominio de integración, de modo que se logre de manera eficiente una tolerancia de error global prescrita.

Scope

Este tema abarca la estimación de errores locales mediante la comparación de reglas de diferentes órdenes o niveles de refinamiento, la bisección recursiva de intervalos (Simpson adaptativo y Gauss-Kronrod adaptativo), los presupuestos de error global y los criterios de parada, el manejo de singularidades y características pronunciadas, y el diseño de integradores automáticos de producción como los de la biblioteca QUADPACK.

Core questions

• ¿Cómo se calcula el error local de una estimación de cuadratura sin conocer la integral exacta?
• ¿Cómo concentra la subdivisión recursiva el esfuerzo donde el integrando varía más?
• ¿Qué criterios de parada logran de manera fiable la tolerancia solicitada evitando el trabajo innecesario?
• ¿Cómo se detectan y manejan de forma robusta las singularidades y discontinuidades integrables?

Key theories

Estimación de error local y subdivisión

La comparación de una estimación gruesa y una más fina (o de orden superior) en un subintervalo produce una estimación del error local; si excede la parte de la tolerancia asignada a ese subintervalo, el subintervalo se divide y el procedimiento se repite, de lo contrario se acepta su contribución.

Estrategias globalmente adaptativas

En lugar de tratar los subintervalos de forma independiente, los integradores globalmente adaptativos mantienen una cola de subintervalos ordenados por error estimado y siempre refinan el peor, lo que maneja las singularidades localizadas de manera eficiente y es la base de las rutinas QUADPACK.

Mechanisms

En cada subintervalo, el integrador evalúa un par de reglas incrustadas —por ejemplo, un par de Gauss-Kronrod o dos estimaciones de Simpson con diferentes refinamientos— cuya diferencia estima el error local. Un método localmente adaptativo recurre bisecando cualquier subintervalo cuyo error estimado sea demasiado grande. Un método globalmente adaptativo mantiene una cola de prioridad de subintervalos clasificados por error estimado y subdivide repetidamente el peor infractor actual hasta que la estimación de error sumada cumple con la tolerancia. Se añaden extrapolación y manejo especializado de pesos para hacer frente a singularidades en los puntos finales e integrandos oscilatorios.

Clinical relevance

La cuadratura adaptativa es en lo que se basan las rutinas de integración de propósito general en el software científico para entregar un resultado con una precisión especificada por el usuario sin que este tenga que analizar el integrando; es esencial para integrandos con picos, comportamiento de capa límite o singularidades integrables que anularían una regla fija, y es la base de los integradores automáticos en paquetes numéricos y estadísticos ampliamente utilizados.

History

La integración automática y controlada por errores maduró en las décadas de 1970 y principios de 1980, culminando en el paquete QUADPACK (1983), cuyas rutinas adaptativas de Gauss-Kronrod con extrapolación se convirtieron en el estándar de facto y fueron posteriormente adoptadas, portadas o reimplementadas en muchos sistemas de software numérico y estadístico.

Key figures

• Robert Piessens
• Philip J. Davis
• Philip Rabinowitz

Seminal works

• davis1984
• piessens1983

Frequently asked questions

¿Cómo sabe un integrador adaptativo el error si no conoce la respuesta?

Estima el error local comparando dos aproximaciones de diferente precisión en el mismo subintervalo —por ejemplo, una regla de orden superior y otra de orden inferior. Su diferencia aproxima el error y guía dónde refinar, aunque la integral verdadera sea desconocida.

¿Cuándo tiene dificultades la cuadratura adaptativa?

Puede ser engañada por integrandos que son suaves en los puntos de muestreo pero tienen características ocultas entre ellos, por integrandos fuertemente oscilatorios o por singularidades no integrables. En tales casos, se necesitan reglas especializadas, transformaciones o métodos de integración oscilatoria.