[Libmesh-users] A node does not have dof attached to it?

[Yuxiang W. <yw...@vi...>]

Dear all,

Sorry for the naive question. In short, I have a system with three
variables (first order Lagrangian elements). I thought that for each of the
node in the mesh, it would have a DOF attached to it with n_comp==1.
However, when I used the following code snippet to check, the n_comp is
actually 0 (not 1).

for (auto node : mesh.node_ptr_range())
  std::cout << node->n_comp(0, 0) << std::endl;

What I was trying to achieve is to iterate through all the nodes, and
assign solution values to each variable DOF attached to the node from a
Eigen::Vector3d. Then I got the error message about:

Assertion `comp < this->n_comp(s,var)' failed.
comp = 0
this->n_comp(s,var) = 0

Any help would be appreciated. Thank you!

Best,
Shawn

-- 
Yuxiang "Shawn" Wang, PhD
yw...@vi...
+1 (434) 284-0836