High order smoothing schemes for inhomogeneous parabolic problems with applications in option pricing