Hello To implement the rigid body modes I need the node coordinates and the dofs they are associated with, or vice versa, the dofs and a coordinate they are associated with. I can iterate over the nodes to obtain their coordinates and to obtain their dofs with DofObject::dof_number(), but this is triggering the assertion libmesh_assert (comp < this->n_comp(s,var)). Are there tools in libMesh to do this more elegantly or to just avoid that error? I think that this way of calculating rigid body modes would be restricted to nodal value finite elements. I don’t know how I could do it with non-lagrange finite elements, but that is not of my concern.
Thanks Miguel ------------------------------------------------------------------------------ Developer Access Program for Intel Xeon Phi Processors Access to Intel Xeon Phi processor-based developer platforms. With one year of Intel Parallel Studio XE. Training and support from Colfax. Order your platform today. http://sdm.link/xeonphi _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users