On Thu, Nov 10, 2016 at 3:01 PM, David Knezevic <david.kneze...@akselos.com> wrote:
> 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? > "P1 basis functions on a second order mesh" I mean 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
