ParILUT - Un ILU de umbral paralelo para GPUs

En este documento, presentamos el primer algoritmo para calcular factorizaciones ILU de umbral en arquitecturas de GPU. El algoritmo propuesto ParILUT-GPU se basa en intercalar iteraciones de punto fijo paralelas que aproximan los factores incompletos para un patrón no nulo existente con una estrategia que adapta dinámicamente el patrón no nulo a las características del problema. Esto requiere la selección eficiente de umbrales que separen los valores a ser eliminados de los factores incompletos, y diseñamos un nuevo algoritmo de selección adaptado a GPUs. Todos los componentes del algoritmo ParILUT-GPU hacen un uso intensivo de las características disponibles en las últimas generaciones de GPU NVIDIA y superan las implementaciones multihilo existentes en CPU.