[deal.II] Neumann bc in step-20

Hi, I’m trying to implement homogeneous Neumann bc on part of the boundary in tutorial example step-20. The Neumann bc should be imposed strongly (as for Dirichlet bc in normal, non-mixed, formulation). Is this correct? In any case, I’m trying to impose the Neumann bc strongly using the

[deal.II] Re: Vector-valued gradient of solution vector

Joel, I did try it, but I cant just change the local_velocity_values from a > std::vector> to std::vector. Then the > fe_values[velocities].get_function_values function doesnt compile. > That's true, but you should use std::vector and not std::vector, see also [1]. >

[deal.II] Re: Vector-valued gradient of solution vector

Dear Daniel, I did try it, but I cant just change the local_velocity_values from a std::vector> to std::vector. Then the fe_values[velocities].get_function_values function doesnt compile. So I commented that part out, the results you see is from the code as it is, I used the for

[deal.II] Re: More power outages

Everything should be available again. Am Mittwoch, 28. September 2016 15:19:15 UTC+2 schrieb sotelo...@gmail.com: > > it is already 28th, the server still down. is it going to take longer > than announced? > Edith > > On Monday, September 12, 2016 at 11:21:13 AM UTC-5, Guido Kanschat wrote: >>

[deal.II] Re: More power outages

it is already 28th, the server still down. is it going to take longer than announced? Edith On Monday, September 12, 2016 at 11:21:13 AM UTC-5, Guido Kanschat wrote: > > Dear all, > > We are expecting more power outages on September 26 to 28. They are trying > to fix the emergency supply such

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

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

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

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] Unique (Continuous) Gradient at Edges/Vertices

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

Re: [deal.II] Finite Difference with deal.ii and Parallel Distributed Triangulation

Dear Rajat, You could do finite differences using the p4est wrapper in deal.II in principle, but it would probably not be very efficient. You want to know the indices which are not available in our wrapper at least, not sure about p4est. To get them, I would build a DoFHandler based on some