On Thu, Nov 10, 2016 at 2:57 PM, Salazar De Troya, Miguel < salazardet...@llnl.gov> wrote:
> 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. > Sounds like some of your nodes don't have any dofs on them. Are you using P1 elements on a second order mesh? I suggest you use something like: for(unsigned int var=0; var<node->n_dofs(sys_number); var++) David ------------------------------------------------------------------------------ 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
