Anders Logg wrote:
> On Mon, Mar 31, 2008 at 03:52:57PM +0200, Murtazo Nazarov wrote:
>
>> 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);
>> ...}
>>
>
> Do you remember to set everything to 1 (number of subdomains) here?
>
>
Yes I do:
// Mark everything as sub domain
1
(*sub_domains) = 1;
// Mark the sub domain as sub domain
0
sub_domain.mark(*sub_domains, 0);
/murtazo
_______________________________________________
DOLFIN-dev mailing list
[email protected]
http://www.fenics.org/mailman/listinfo/dolfin-dev
