Order assembly 1D matrix

mesh
1

Hi guys,

I’m a little confused when I assemble my mass (or stiffness) matrix.
I wanted to double check with FEniCS if my assembly was correct, but it is kind of reverted.
For example, what my code put in (0,0), FEniCS put in (n-1,n-1) and viceversa. So it seems like the dofmap is somehow reversed.
But I checked it and it is the normal enumeration for 1D cases, from left to right.

Is there something I’m missing, like that fenics starts from the right node for some reason ?

Alos using mesh.coordinates() gives me in return the nodes coordinates from left to right.

Any suggestion ?

2

Hi,
no Dolfin reorders dofs even in this simple case:

from dolfin import *
mesh = UnitIntervalMesh(5)
V = FunctionSpace(mesh, "CG", 1)
V.tabulate_dof_coordinates()

returns

array([[1. ],
       [0.8],
       [0.6],
       [0.4],
       [0.2],
       [0. ]])
3

Thank you! This eraseevery doubts