That’s exactly what I was doing. I am playing with the example https://github.com/libMesh/libmesh/blob/master/examples/systems_of_equations/systems_of_equations_ex7/systems_of_equations_ex7.C where the mesh is built with HEX27, but the finite element is first order lagrange. Why is there this mismatch? I didn’t know that meshes could have an order, I thought it was up to the finite element space.
Miguel From: David Knezevic <david.kneze...@akselos.com<mailto:david.kneze...@akselos.com>> Date: Thursday, November 10, 2016 at 12:02 PM To: Miguel Salazar <salazardet...@llnl.gov<mailto:salazardet...@llnl.gov>> Cc: "libmesh-users@lists.sourceforge.net<mailto:libmesh-users@lists.sourceforge.net>" <libmesh-users@lists.sourceforge.net<mailto:libmesh-users@lists.sourceforge.net>> Subject: Re: [Libmesh-users] Implementing rigid body modes On Thu, Nov 10, 2016 at 3:01 PM, David Knezevic <david.kneze...@akselos.com<mailto:david.kneze...@akselos.com>> wrote: On Thu, Nov 10, 2016 at 2:57 PM, Salazar De Troya, Miguel <salazardet...@llnl.gov<mailto: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
