### Kronecker product of matrices

109

views

0

Hello,

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!

Community: FEniCS Project

Please login to add an answer/comment or follow this question.