On Thu, Mar 17, 2016 at 3:31 PM, Roy Stogner <[email protected]>
wrote:

>
> On Thu, 17 Mar 2016, David Knezevic wrote:
>
> 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.
>>
>
> Would it make sense to just add that to initialize_condensed_dofs()?
>
> Just before
>   local_non_condensed_dofs_set.insert(i)
> add
>   if (!this->get_dof_map().is_constrained_dof(i))
> ?
>

This would add all constrained dofs to the list of condensed dofs. That
should do what we want, I agree.

Would you want to get rid of the argument "std::set<unsigned int> &
global_dirichlet_dofs_set"? I can imagine using that to set "extra"
condensed dofs, but I don't think that's a very common use-case, and anyway
users could just add DofConstraintRows to handle that case if needed.

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
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to