dolfin team mailing list archive

[Patrick Riesen <priesen@xxxxxxxxxxxxxxx>]

similarly, i'm solving channel flow problems and the like.
no, i don't need pointwise, you are right, it's fine now also with 
topological, i found the self-inflicted problem ;-)

thank you for the help,

best regards,
patrick


Anders Logg wrote:
> 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