On Mon, 4 Mar 2013, Ataollah Mesgarnejad wrote: > On the other hand if I wanted to use the API how would you suggest > updating the ID? Should I do a loop over all the elements and then > faces and add faces that are in contact zone to an ID?
Right. > If so can a face have two IDs? In libMesh itself, definitely yes. I'm not sure which of our I/O formats handles that correctly, though. > And is there a place in the source code you do something similar so > I can look and get an idea? src/apps/meshbcid.C loops over elements, tests them, and changes their side bc ids. But you'll want to use the other remove_side overload to only remove a single id rather than all ids, and you'll be removing and adding ids on different tests rather than always doing a remove and an add at once. Ooh, but this actually reminds me of another concern: right now we only store id information on the coarse mesh. If you want to combine this with adaptivity and have Dirichlet boundaries which contain only portions of coarse element sides, you're out of luck unless you want to rework boundary_info.C first. --- Roy ------------------------------------------------------------------------------ Everyone hates slow websites. So do we. Make your web apps faster with AppDynamics Download AppDynamics Lite for free today: http://p.sf.net/sfu/appdyn_d2d_feb _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
