Implement -c(t,x,a)v term in the differential operator

[anriseth]

Extend the code to handle equations of the form

v_t+\sup_{a\in A}\{b*v_x + 1/2*\sigma^2*v_{xx} - c*v + f\} = 0