On Thu, 16 Dec 2010, R P wrote:
> I am interested if it is posible to imposse periodicity on the attached mesh > (created in Gmsh) in circular (torioidal) direction by this code: > > DofMap & dof_map = system.get_dof_map(); > PeriodicBoundary vert(RealVectorValue(0.0, 0.000001, 0.0)); > > vert.myboundary = 5; > vert.pairedboundary = 6; > > dof_map.add_periodic_boundary(vert); > > If yes, could you please advise me, how to set correct boundary_info IDs > (e.g. 5, 6) on the apropriate sides? MeshBase::boundary_info->remove_side() and MeshBase::boundary_info->add_side() are the methods for changing boundary IDs. The best thing you can usually do is loop over all elements, find those which match some criteria (location, face normal, previous boundary ID, etc), and set their new boundary ID accordingly. I've put one of my own little utilities up as an example: http://users.ices.utexas.edu/~roystgnr/meshbcid.C If anyone is interested I could add this utility to the libMesh src/apps/, or accept a patch adding its functionality to src/apps/meshtool.C --- Roy ------------------------------------------------------------------------------ Gaining the trust of online customers is vital for the success of any company that requires sensitive data to be transmitted over the Web. Learn how to best implement a security strategy that keeps consumers' information secure and instills the confidence they need to proceed with transactions. http://p.sf.net/sfu/oracle-sfdevnl _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
