[firedrake] Matrix multiplication in bilinear form

Dear all, Can I use an assembled matrix into a bilinear form? E.g. I want to write something like mu1 = Function(V) mu = TrialFunction(V) v = TestFunction(V) Q = assemble(v*dx) lhs = inner(grad(mu),grad(v))*dx + Q*mu*dx rhs = .... solve(lhs, mu1, rhs) Shall/can I just assemble each term in the lhs independently, and then add/multiply the resulting matrices? Best, Anna.