El enfoque lattice/FD sugerido por @LocalVolatility funcionará bien.
Sin embargo, dado que usted menciona específicamente " basado en las trayectorias simuladas de un movimiento browniano geométrico ", se podría considerar alternativamente el Monte-Carlo por mínimos cuadrados.
Más concretamente, trabajando hacia atrás en el tiempo desde el vencimiento hasta el inicio del contrato, el algoritmo de Longstaff-Schwartz le permitirá calcular los valores de continuación a partir de las trayectorias simuladas -- bajo el capó, esta expectativa condicional se evalúa en el sentido de mínimos cuadrados utilizando técnicas de regresión estándar (de ahí el nombre original del método), pero también podría utilizar métodos de aprendizaje supervisado más elaborados en este punto, supongo. La intersección del valor de continuación y el valor intrínseco para cada tiempo $t$ entonces, por definición, constituye el límite óptimo de ejercicio.
La parte complicada es que el algoritmo estándar de Longstaff-Schwartz en realidad identifica una estrategia de ejercicio subóptima (los precios resultantes tendrán un sesgo bajo). El lado positivo es que las estructuras generales de dividendos pueden manejarse casi sin problemas en Monte Carlo.