dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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 <logg@xxxxxxxxx>
 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     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help   : https://help.launchpad.net/ListHelp

Attachment: signature.asc
Description: Digital signature