On Sun, 2 Mar 2014, Francesco Ballarin wrote:
> I have a NumericVector v that stores a discretized function and I would
> like to use this vector v as a Dirichlet BC for a system s via the
> DirichletBoundary class.
> I guess I could use something like
> DirichletBoundary( MeshFunction(v) )
> but, at least if the vector v and the system s are defined on the same mesh
> and with the same FE discretization, is there any more efficient way to do
> this?
At the moment, no. Ideally we ought to have a DirichletBoundary
capable of optionally taking an FEMFunctionBase rather than just a
FunctionBase, but that would be a significant patch to write; I
wouldn't suggest you bother with it unless profiling tells you you
have a bottleneck in the constraint creation. The inefficient code
path will be O(N^(2/3) log N) instead of O(N^{2/3}), but that will
still be asymptotically cheaper than your solve.
---
Roy
------------------------------------------------------------------------------
Subversion Kills Productivity. Get off Subversion & Make the Move to Perforce.
With Perforce, you get hassle-free workflows. Merge that actually works.
Faster operations. Version large binaries. Built-in WAN optimization and the
freedom to use Git, Perforce or both. Make the move to Perforce.
http://pubads.g.doubleclick.net/gampad/clk?id=122218951&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
