On 7 June 2011 13:55, Chong Luo <luo.ch...@yahoo.com> wrote: > Suppose we have a square domain. Is it possible to impose Dirichlet boundary > conditions only on the 4 corners?
Have a look at dolfin/demo/undocumented/dg-advection-diffusion/cpp/main.cpp DirichletBC bc(V, g, boundary, "geometric"); replace "geometric" with "pointwise" and make sure that your SubDomain is defined in terms of geometry only, i.e., it shouldn't use the 'on_boundary' variable in the inside() function. Kristian > Thanks! > > _______________________________________________ > Mailing list: https://launchpad.net/~dolfin > Post to : dolfin@lists.launchpad.net > Unsubscribe : https://launchpad.net/~dolfin > More help : https://help.launchpad.net/ListHelp > _______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : dolfin@lists.launchpad.net Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp