On Jan 22, 2013, at 1:24 PM, "Kirk, Benjamin (JSC-EG311)" <[email protected]> wrote:
> On Jan 22, 2013, at 9:22 AM, Ataollah Mesgarnejad <[email protected]> > wrote: > >> To reiterate my question is: is there a way to get dof_indices of the >> sideset to use for setting the global vector using my "local_vec"? > > Not sure I understand 100% why you need to do this, but at any rate you can > get what you want by building a side from the interior element. This will > construct a valid element in its own right that lives in (dim-1). You can > then feed *this side element* to the DofMap instead, and yo'll get the > indices you desire. > Works like a charm, can't thank you enough. PS: for future reference in case anybody was searching the mailing list: The reason I was looking to do this was that I'm using PETSc dm nonlinear solver and I needed to set the upper & lower bounds for my variable to something specific based on the side set id of the boundary i.e. I'm imposing a Dirichlet boundary on a boxed constrained nonlinear system. Best, Ata > -Ben ------------------------------------------------------------------------------ Master Visual Studio, SharePoint, SQL, ASP.NET, C# 2012, HTML5, CSS, MVC, Windows 8 Apps, JavaScript and much more. Keep your skills current with LearnDevNow - 3,200 step-by-step video tutorials by Microsoft MVPs and experts. ON SALE this month only -- learn more at: http://p.sf.net/sfu/learnnow-d2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
