Dear all,
I want a way to set a NumercVector<Number> on a sideset but I'm not sure how I
can access the DOF indices of a given sideset.
Here is a snippet of what that part of my code looks like:
for ( ; el1 != end_el; ++el1)// loop on all elements
{
const Elem* elem = *el1;
dof_map.dof_indices (elem, dof_indices, var);
const unsigned int n_dofs = dof_indices.size();
local_vec.resize (n_dofs); // local vector I'm using to set the Global
vector
local_vec.zero();
for (unsigned int s=0; s<elem->n_sides(); s++)// loop on all sides
if (elem->neighbor(s) == NULL)// my side set is always on a
boundary
{
boundary_id_type bc_id =
mesh.boundary_info->boundary_id (elem,s);
if(bc_id == given_id )// see if the side "s" is on
sideset "given_id"
{
//here is what I miss ?
}
.
.
.
}
dof_map.constrain_element_vector (local_vec, dof_indices);
Global_vec.insert(local_vec, dof_indices_side); // dof_indices_side
?????
}
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"?
As always thanks for the help,
Ata
------------------------------------------------------------------------------
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
