> Elems already have ids... elem->id(). If you are careful when adding your > elems from your mesh file you can make them match what's in the mesh file.
Yes, this is true. But I wasn't sure how the elements get renumbered eventually and I didnt want to mess with that. Moreover, my application does not rely on element id's for anything and so it did not make sense to me to rely on this. > As for surfaces.... those are probably used for boundary conditions. In > that case you should add those surface ids to the BoundaryInfo class > associated with the Mesh. You can look in the Exodus reader / writer for an > idea of how this can work. This will allow you to get lists of sides back > from the BoundaryInfo object when you want to apply a boundary condition. The surface/volume id's are specifically to group different elements, of the same physical material (same subdomain_id) that are in a certain region of the domain. For example, you could have a 3x4 lattice with patches of same material as 1 2 3 1 2 3 1 2 1 3 2 2 This is very useful in transients I've come across and makes the process of code verification so much easier rather than relying on physical coordinate points themselves. I've not used this for any kind of boundary conditions yet but I'll look in to Exodus IO to see how you've used it there. On Wed, Oct 28, 2009 at 11:52 AM, Derek Gaston <[email protected]> wrote: > On Oct 28, 2009, at 10:38 AM, Vijay S. Mahadevan wrote: > >> Arvind, If you are interested in storing the surface/volume ids for >> each element, you need to create additional data in Elem class. This >> would be very similar to the subdomain_id() method. > > Elems already have ids... elem->id(). If you are careful when adding your > elems from your mesh file you can make them match what's in the mesh file. > > As for surfaces.... those are probably used for boundary conditions. In > that case you should add those surface ids to the BoundaryInfo class > associated with the Mesh. You can look in the Exodus reader / writer for an > idea of how this can work. This will allow you to get lists of sides back > from the BoundaryInfo object when you want to apply a boundary condition. > > Derek > ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
