On Tue, 8 Oct 2013, David Knezevic wrote: > However, I still get compilation errors with --enable-complex (on eb4d88f). > Have you tested this?
No; I just replaced the "complex won't implicitly convert to bool" problem with an "int won't implicitly convert to complex" problem. I'm replacing that with an explicit conversion and testing the --enable-complex build now. > DirichletBoundary on nodes. Do you have any suggestions on how I should get > started? I think the key point from your original email was this: > > "For parallel consistency we actually do: > 1. Interpolate at nodes > 2. L2 project any edge degrees of freedom left over when holding node > dofs fixed > 3. L2 project any face dofs left over when holding node and edge dofs > fixed " That key point relates to the constraint equation algorithm, though... but I guess that's the only hard part? I think that the first thing I'd upgrade for dirichlet nodes and edges would be the API (very straightforward "add_foo(ForwardDeclaredFoo&)" declarations); next the data structures (no modifications necessary if we go the "use boundary ids corresponding to nodes and edges" route), finally the algorithms (not too hard once the data structures are done). No, scratch that. The very first thing I'd do is add "boundary ids that apply to edges" support. 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. --- Roy ------------------------------------------------------------------------------ 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
