A 3D front propagation algorithm for the simulation of geological folding is presented. The algorithm is implemented on multicore and GPU architectures. Numerical experiments based on OpenMP show an almost perfect scaling. Using the automatic code translator Mint, Cuda code is created from seqential C code. On a GPU the Cuda code runs 67-95 times faster than sequential code, and allows 3D grids to be computed at interactive times.