Two models were compared: Fast-Marching Method, Fast-Iterative Method, with the analytical solution. All the models were tested using a 3D block with a single initial stimulus on the upper corner It seems that when you account for anisotropy both models have similar accuracy. However, the Fast-Iterative method has more potential to be run in parallel. For that reason, I have chosen this method.