Resolución de Sistemas Hepta-Diagonales en GPU

Abstract

Este artículo presenta la implementación del método Strongly Implicit Procedure (SIP) para la resolución de sistemas lineales hepta-diagonales en Graphics Processing Unit (GPU). En particular, se realizaron tres implementaciones distintas del método SIP que buscan aprovechar diferentes características de las GPUs (accesos eficientes a memoria, uso de memoria compartida, transferencias de datos optimizadas, etc.). Los resultados obtenidos trabajando en simple y doble precisión muestran una importante mejora en el tiempo de ejecución alcanzado, con aceleraciones de hasta 22×, comparado con la versión paralela del método en CPU. Además, se incluyó el solver desarrollado en un modelo numérico
de mediano porte, el caffa3d.MB, obteniendo una reducción de hasta 19% en el tiempo total de ejecución del modelo.