Um dos principais desafios em algoritmos de visualização volumétrica é calcular a integral volumétrica de maneira eficiente, mantendo uma precisão mínima adequada. Geralmente, métodos de integração numérica utilizam passos de tamanho constante, não incluindo nenhuma estratégia de controle numérico. Como uma possível solução, métodos numéricos adaptativos podem ser utilizados, pois conseguem adaptar o tamanho do passo de integração dada uma tolerância de erro pré-definida. Em CPU, os algoritmos adaptativos de integração numérica são, normalmente, implementados recursivamente. Já em GPU, é desejável eliminar implementações recursivas. O presente trabalho propõe um algoritmo adaptativo e iterativo para a avaliação da integral volumétrica em malhas regulares, apresentando soluções para manter o controle do passo da integral interna e externa. Os resultados do trabalho buscaram comparar a precisão e eficiência do método proposto com o modelo de integração com passo de tamanho constante, utilizando a soma de Riemann. Verificou-se que o algoritmo proposto gerou resultados precisos, com desempenho competitivo. As comparações foram feitas em CPU e GPU.