¿Se sigue trabajando activamente en QUADPACK para incorporar nuevos avances en integración numérica, o la biblioteca sigue basándose únicamente en los algoritmos que existían en el momento en que se escribió inicialmente?
Respuesta
¿Demasiados anuncios?
No se mantiene, la última actualización que puedo ver es de 1987 . Dicho esto, la integración numérica 1D no ha cambiado increíblemente rápido en los últimos 30 años. (La cuadratura de alta dimensión estaría irremediablemente desfasada.) El único algoritmo que obviamente falta es el doble exponencial cuadratura para integradas situadas en espacios de Hardy.
Curtis-Clenshaw y Gauss-Kronrod son caballos de batalla que probablemente serán útiles durante siglos. Sin embargo, el principal problema de QUADPACK es que está escrito en Fortran 77. Esto no sería un problema si la gente siguiera leyendo y aprendiendo Fortran 77, pero no es así. g77 ni siquiera se mantiene, por lo que no puede generar nuevos conjuntos de instrucciones. Así que tendrás que trabajar mucho para conseguir proyectos que funcionen con QUADPACK en bruto. scipy mantiene un QUADPACK envoltura que es muy fácil de usar, por lo que puede evitar la molestia de aprender Fortran 77.
