Is it possible to calculate the derivative of a solution u with respect to a piecewise constant parameter c? My goal is to calculate the sensitivity of u by doing something like:
u
c
sL = diff(L, c) - action(diff(a,c), u)