[Libmesh-users] DirichletBoundary question

[David K. <dkn...@se...>]

I'm passing a subclass of FunctionBase<Number> to a DirichletBoundary 
object, and the FunctionBase<Number> uses a MeshFunction from a separate 
(dim-1)-dimensional mesh to return a value. In this context, the 
lower-dimensional mesh "matches" the relevant boundary of the 
dim-dimensional mesh.

I sometimes have some tolerance issues, so that I have to find the 
"nearest point" in the (dim-1)-dimensional mesh. At present, I just find 
the "nearest node" on the (dim-1)-dimensional mesh (rather than the 
nearest interior point or whatever), and it works fine.

I gather that this is because (at least in my case) the 
DirichletBoundary only evaluates the FunctionBase at the nodes on the 
boundary of the dim-dimensional mesh. I'm wondering if this will be true 
in general, or if this is just due to the fact that I'm using nodal 
(i.e. LAGRANGE) basis functions in this case?

Thanks!
David