Hi all! How compute the global stiffness matrix derivatives with respect to cell coeeficients？

Consider a bilinear form：

a(u, v)=\int_{\Omega}d_h^p\mathbb{\sigma}:\mathbb{\epsilon}d\Omega

where d_h is a cell-wise function, h\in\mathcal{T} is the cell of the fem domain \mathcal{T}, p is a constant, \sigma is the stress tensor and \epsilon is the strain tensor.

Assemble the a(u,v), the stiffness matrix K is obtained.

Now my question is how to compute the derivative of K w.r.t. d_h using dolfinx and ufl.