Re: [firedrake] retrieve dof and entities of assembled matrix entries

On 25 Aug 2020, at 14:31, Karin&NiKo <niko.karin@gmail.com> wrote:

I feel sorry to bother you again but something is going wrong. Here is my script :

from firedrake import *

nbx = 1 nby = 1 lx = 0.5 ly = 0.5 mesh = RectangleMesh(nbx, nby, lx, ly, quadrilateral=False)

Vu = VectorFunctionSpace(mesh, "Lagrange", 2) Vp = FunctionSpace(mesh, "Lagrange", 1) Vt = FunctionSpace(mesh, "Lagrange", 1) Z = Vu * Vp * Vt

dof_coord = Function(Z) for fn in dof_coord.split(): fn.interpolate(mesh.coordinates)

The mesh.coordinates function is in a vector space, and Vp and Vt are both scalar spaces, so you can't interpolate from one into the other. Lawrence