how do you best extend black-scholes or numerical methods based on lognormal processes to handle discrete dividends?