This is a read only copy of the old FEniCS QA forum. Please visit the new QA forum to ask questions

vertex_to_dofmap() producing nonsensical indices

0 votes

Hello,

I'm getting nonsensical indices (huge numbers) from the vertex_to_dofmap() function.

from dolfin import *
mesh = UnitCubeMesh(10, 10, 10)
VV = VectorFunctionSpace(mesh, 'CG', 1)
len(mesh.coordinates()) # returns: 1331
ver2dof = vertex_to_dof_map(VV)
len(ver2dof)/3 # returns: 1331
ver2dof # returns: array([4964982195331, 6519760356485, 6528350291439, ..., 46047632, 215733567192, 0])
ver2dof.max() #returns: 17141214482326

Given the output of line 4 I would expect the indices in ver2dof in the range [0, ..., 3992]. I'm getting instead really large numbers. This is happening in fenics/1.6.0 installed on a cluster in my University. On my laptop everything works as it should, under the same version. I would appreciate any ideas what may be going wrong.

Kind regards,
Miguel A. Valdez G.

asked Nov 4, 2016 by miguelvaldez FEniCS Novice (140 points)
...