Kronecker product of matrices
I am trying to understand how one can compute the Kronecker product. Let us say, we have the following code:
mesh = UnitSquareMesh(8,8)
V = FunctionSpace(mesh, 'CG', 1)
u = TrialFunction(V)
v = TestFunction(V)
K = assemble( dot(grad(u),grad(v))*dx )
M = assemble(u*v*dx)
Now, I would like to compute the Kronecker product: M \otimes K.
I had a look at the documentation here: http://fenics.readthedocs.io/projects/ufl/en/latest/manual/form_language.html
However, it does not work as I do not clearly understand how the form conversion.
If somebody has an idea how do we implement this, please let me know. Thanks!