Kronecker product of matrices

11 weeks ago by


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:

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.

Similar posts:
Search »