[deal.II] Re: Unique (Continuous) Gradient at Edges/Vertices

2016-09-30 Thread Andreas Krämer
Hi, I have started the implementation of such a class. The only problem is that currently there seems to be no parallel::distributed version of hp::DoFHandler. Calculating smooth gradients over the whole domain with your method would be way too costly, I guess. So I am going to have to wait un

[deal.II] Re: Unique (Continuous) Gradient at Edges/Vertices

2016-09-28 Thread Jean-Paul Pelteret
Hi Andreas, No problem, I happened to reimplement this for myself yesterday, so its very fresh in my mind :-) Best, J-P On Wednesday, September 28, 2016 at 11:08:28 AM UTC+2, Andreas Krämer wrote: > > Hi Jean-Paul, > > wow, thanks for the quick answer. > Your approach is really elegant. It sho

[deal.II] Re: Unique (Continuous) Gradient at Edges/Vertices

2016-09-28 Thread Andreas Krämer
Hi Jean-Paul, wow, thanks for the quick answer. Your approach is really elegant. It should work in principle. I will try out to get around the FEFieldFunction and get back to this thread afterwards. Best, Andreas Am Mittwoch, 28. September 2016 10:40:19 UTC+2 schrieb Jean-Paul Pelteret: > >

[deal.II] Re: Unique (Continuous) Gradient at Edges/Vertices

2016-09-28 Thread Jean-Paul Pelteret
Hi Andreas, I've had to do something similar in a different context. You could compute the L2-projection of the solution gradient using VectorTools::project. You'd need to define a new DoFHandler for the gradient components and define a function that returns the value of the solution gradient a