Hi all,<div><br></div><div>Perhaps this may be a simple question, but say I have this bilinear and linear form:</div><div><br></div><div>a = grad(u)*D*grad(v)*dx</div>L = F*v*dx<div><br></div><div>where u,v is trial/test function on CG1 space, and D and F are coefficients that live in DG0 space.</div><div><br></div><div>Say I have a velocity vector function q (CG1) on the same mesh and want D to be the element-wise divergence of said velocity.</div><div><br></div>How would I formulate D in Firedrake? <div><br></div><div>Thanks,</div><div>Justin</div>