O estudo de fluxo de fluidos em grande escala em meios porosos exige cálculos intensos e ocorre nas mais diversas aplicações, incluindo o fluxo de águas subterrâneas e a recuperação de petróleo. Este trabalho apresenta novas estratégias computacionais aplicadas à geomecânica de reservatórios. São propostos avanços para a montagem eficiente de matrizes de elementos finitos e a solução de sistemas lineares usando código altamente vetorizado no MATLAB. Na versão para CPU, a montagem da matriz de elementos é realizada usando procedimentosconvencionais de vetorização, com base em duas estratégias: as matrizes explícitas e os produtos multidimensionais. A montagem adicional da matriz esparsa global é obtida usando a função esparsa nativa. Para a versão GPU, o cálculo do conjunto completo de matrizes deelementos é realizado com as mesmas estratégias da abordagem da CPU, usando estruturas gpuArraye o suporte nativo CUDA fornecido pelo MATLAB Parallel Computing Toolbox. A solução do sistema linear resultante nas versões de CPU e GPU é obtida com duas estratégias usando uma abordagem unidirecional: o solucionador de gradiente conjugado nativo (pcg) e o fornecido pela biblioteca Eigen. Uma ampla discussão é apresentada em um benchmark dedicado, em que as diferentes estratégias usando CPU e GPU são comparadas em termos de tempo de processamento e requisitos de memória. Essas análises apresentam aumentos de velocidades significativos em relação aos códigos seriais. Por fim, são fornecidas recomendações.