[deal.II] Quadrature on Lagrange Polynomials

2016-11-24 Thread Andreas Krämer
Hi everybody, Does deal.II provide a way to use Gauss-Lobatto-Chebychev quadratures when using FE_DGQArbitraryNodes<2>? Specifically, I want to integrate functions f(x) using Lagrange polynomials on the GLC support points. The quadrature there is usually defined with a weighting function, as \in

[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] Parallel hp::DoFHandler

2016-09-30 Thread Andreas Krämer
Dear developers, can you already foresee when there will be a version of hp::DoFHandler that works with parallel::distributed::Triangulation? Best, Andreas -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?

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

2016-09-28 Thread Andreas Krämer
eldFunction. > > I hope this helps, > J-P > > On Wednesday, September 28, 2016 at 10:18:30 AM UTC+2, Andreas Krämer > wrote: >> >> Hi everybody, >> >> I am currently implementing a semi-Lagrangian advection solver in deal.II. >> >> Everything w

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

2016-09-28 Thread Andreas Krämer
Hi everybody, I am currently implementing a semi-Lagrangian advection solver in deal.II. Everything works fine in periodic domains, but now I want to define boundary conditions that depend on the gradient of the solution at the boundary. When I use Lagrangian FE_Q elements, the gradients are di

[deal.II] Mapping and round-off errors

2016-05-19 Thread Andreas Krämer
Hi, I have just encountered an issue with projecting points into a cell. Say, I have a point x which is very close to a cell (e.g. distance = 10^-16) and I want this point to be inside the cell. My first attempt was to map the point to the unit cell, round the coordinates to 0 or 1, and map