On Thu, Mar 17, 2016 at 2:57 PM, Roy Stogner <royst...@ices.utexas.edu> wrote:
> > On Thu, 17 Mar 2016, David Knezevic wrote: > > I'd be interested to see the code, if you can send it through. In >> particular, I'd be interested to see how you're imposing the >> boundary conditions. Are you using CondensedEigenSystem? If so, >> perhaps there is an issue with how the "indices to condense" are >> being communicated in parallel. >> > > There's one obvious bug I can see: > > If a Dirichlet node is owned by a processor which doesn't own any > elements with a Dirichlet side containing that node, then that dof > won't be condensed by get_dirichlet_dofs, because the processors that > try to add it to the condensed set don't own it and the processor that > owns it won't try to add it. > I think this is the right explanation here, thanks Roy! > IIRC I had to fix a similar bug when putting together the > DirichletBoundary support. We really ought to be adding > DirichletBoundary support to CondensedEigenSystem and not forcing the > user to try to do this themselves. I agree that it would be great to update CondensedEigenSystem so that it uses DirichletBoundary objects instead of requiring the user to set the condensed dofs manually. Actually, I just looked at one of my apps that uses CondensedEigenSystem and DirichletBoundary, and there I build up the list of dofs to condense using code based on ConstrainDirichlet::apply_dirichlet_impl. It would make sense to include something like that in CondensedEigenSystem, I think. > Hmmm... if constrained dofs were > treated properly in CondensedEigenSystem, both DirichletBoundary and > AMR would work automatically, right? Do you mean just condensing out all constrained dofs? I use that already with DirichletBoundary and that works fine. And I think it would work automatically with AMR too. David ------------------------------------------------------------------------------ Transform Data into Opportunity. Accelerate data analysis in your applications with Intel Data Analytics Acceleration Library. Click to learn more. http://pubads.g.doubleclick.net/gampad/clk?id=278785231&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
