Hi,
I want to apply a boundary condition to the vertices on a boundary. For that I
define:
// Sub domain for MyBC
class MyBC_Boundary2D : public SubDomain
{
public:
bool inside(const real* p, bool on_boundary) const
{
return on_boundary && (p[0] < xmax - bmarg) && (p[0] > xmin + bmarg);
}
};
Then I initialize:
//-----------------------------------------------------------------------------
void MyBC::init(SubDomain& sub_domain)
{ ...
mesh.init(0);
sub_domains = new MeshFunction<uint>(mesh, 0);
...}
and skip interior nodes:
void MyBC::apply(...)
{...
for (VertexIterator vertex(mesh); !vertex.end();++vertex)
{
// Skip facets not inside the sub domain
if ( (*sub_domains)(*vertex) != sub_domain )
{
continue;
}
cout << "vertex = " << (*vertex).index() << endl;
}
...}
But, it still does not skip interior nodes. Do I do something wrong or there
should be some changes in subdomain.mark?
/murtazo
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
