Hi Roy, I've got a branch node_and_edge_dirichlet in my fork on github (github.com/dknez/libmesh) which is working for me.
The changes to ConstrainDirichlet were indeed simple. I also added some methods in BoundaryInfo to allow us to add edge boundary ids. I've tested this on a modified version of systems_of_equations_ex6 (also in the branch) and it works well. The branch is based on your adjoint_dirichlet branch, so I'll wait for you to merge that before I make a pull request. David > As to the algorithm... cross your fingers, but I suspect that with the > current refactoring all you'd have to do is upgrade the > "do_this_side", "is_boundary_node[n]" and "is_boundary_edge[e]" > calculations in ConstrainDirichlet::apply_dirichlet_impl and you'll be > done. ------------------------------------------------------------------------------ October Webinars: Code for Performance Free Intel webinars can help you accelerate application performance. Explore tips for MPI, OpenMP, advanced profiling, and more. Get the most from the latest Intel processors and coprocessors. See abstracts and register > http://pubads.g.doubleclick.net/gampad/clk?id=60134071&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
