En cuanto a las EDP (deterministas), tenemos la noción de "solución fuerte" (que resuelve directamente el operador diferencial en la formulación fuerte del problema) y la "solución débil", que se ocupa de una formulación débil del problema.
Para la formulación fuerte, las diferencias finitas son el camino a seguir ya que son la discretización natural del operador diferencial.
Para la formulación débil, los métodos de elementos finitos son el camino a seguir, ya que abordan directamente la formulación débil del problema restringiéndola a un espacio dimensional finito (dependiendo del tipo de "funciones-elemento" o funciones-base que se elija para este espacio).
El método de las diferencias finitas tiene problemas con las geometrías complejas y las mallas adaptativas - la geometría no será un problema en la valoración de opciones ya que siempre se considera el rectángulo $[0,T]\times[S_\text{min}, S_\text{max}]$ . Perfeccionamiento local puede ser un problema - pero depende de la ecuación y de la condición inicial/de límite. Además, en el caso de algunas ecuaciones diferenciales (no lineales), se producen problemas de disconticidad y se pueden producir efectos de oscilación. Hay esquemas que evitan este tipo de problemas, pero hay que invertir en el algoritmo.
El método de elementos finitos es una noción más general. Dado que las SPDE se definen por su formulación integral (Ito-), un enfoque que se aproxime a una formulación integral resultará más natural. Esto se debe a que el resultado de una llamada europea (variable S) y las trayectorias del movimiento browniano (variable t) no son diferenciables. El uso de métodos de diferencias finitas para SPDEs discretizaciones más naturales para los "operadores diferenciales" no le da el esquema correcto por lo que yo sé. Esa sería otra pista que apunta en la dirección de los elementos finitos un poco.
Un problema al comparar el rendimiento de ambos es, sin duda, que hay tantos esquemas diferentes (implícitos/explícitos de diferentes ordas para las diferencias finitas y elecciones de funciones base y técnicas de refinamiento local para los elementos finitos) que es imposible decirlo. Quizás para alguna elección de modelo y condición inicial/limítrofe un método supere al otro, pero creo que es difícil de generalizar.