Why do you need pointwise? The topological should be fine, especially if you are a version of DOLFIN later than this changeset:
revno: 4483 committer: Anders Logg <[email protected]> branch nick: dolfin-dev timestamp: Thu 2010-02-11 19:24:46 +0100 message: Bug fix in SubDomain::mark, fixes bug in DirichletBC based on SubDomain::inside This fixed a bug that we had when setting boundary conditions for a channel flow problem. -- Anders On Fri, Mar 26, 2010 at 11:33:07AM +0100, Patrick Riesen wrote: > hello, > > i'm under way of creating in-/outflow bc conditions by transferring > a discrete 1D flow profile function onto a boundary subdomain of a > 2D mesh. I marked the subdomain by a meshfunction and use an > expression in which i evaluate the bcs values by extracting them > from the 1D discrete solution as > > BCInFlowProfile::eval (Array<double>& values, const Array<double>& x) const > { > [....] > velocity1Dprofile.eval(values, x_eval); > [....] > > } > i want to use "pointwise" method for applying the bcs, so that the > velocities evaluated and specified on the vertices, but i get > > dolfin/fem/DirichletBC.cpp:635: void > dolfin::DirichletBC::compute_bc_pointwise(std::map<unsigned int, > double, std::less<unsigned int>, std::allocator<std::pair<const > unsigned int, double> > >&, dolfin::BoundaryCondition::LocalData&) > const: Assertion `user_sub_domain' failed. > > > is using DirichletBC "pointwise" with a meshfunction and a subdomain > id not possible? > > regards, > patrick > > > _______________________________________________ > Mailing list: https://launchpad.net/~dolfin > Post to : [email protected] > Unsubscribe : https://launchpad.net/~dolfin > More help : https://help.launchpad.net/ListHelp
signature.asc
Description: Digital signature
_______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
