Re: [deal.II] How to use FE_FaceQ with update_gradients.

2019-08-26 Thread Lars Corbijn
Wolfgang, Thank you for the detailed information. That sounds like a good solution to the problem given the current situation. Now that I have a working checkout (but waiting for a recompile), I will see if I can find some spare time to program the solution. Lars On Tue, 20 Aug 2019 at 00:54, W

[deal.II] How to use FE_FaceQ with update_gradients.

2019-08-14 Thread Lars Corbijn
of FE_FaceQ is no problem, as its gradients are not need and (I expect) not well defined, but I do need the gradients of FE_DGQ. Is there a way around this without eliminating the degrees of freedom for the elements (like in step 51)? Regards, Lars Corbijn [1] https://epubs.siam.org/doi/abs

Re: [deal.II] Evaluate solution of face elements

2017-03-22 Thread Lars Corbijn van Willenswaard
Thank you for the answer. After some more testing I found out the code was correct, but my reasoning/math was not for cases with hanging nodes. Best, Lars On Fri, Mar 17, 2017 at 9:42 AM, Wolfgang Bangerth wrote: > On 03/14/2017 10:54 AM, l.j.corbijnvanwillenswa...@student.utwente.nl > wrote: >

Re: [deal.II] InternalError in VectorTools::compute_nonzero_normal_flux_constraints

2017-02-26 Thread Lars Corbijn van Willenswaard
On Tue, Feb 21, 2017 at 11:28 PM, Wolfgang Bangerth wrote: > On 02/21/2017 07:54 PM, Lars Corbijn van Willenswaard wrote: > >> I will try to do so somewhere this week, but it depends a bit on how much >> free >> evening time I have this week. Furthermore, I don'

Re: [deal.II] InternalError in VectorTools::compute_nonzero_normal_flux_constraints

2017-02-21 Thread Lars Corbijn van Willenswaard
I will also make a patch for it. Lars On Mon, Feb 20, 2017 at 10:56 PM, Wolfgang Bangerth wrote: > On 02/20/2017 07:22 PM, Lars Corbijn van Willenswaard wrote: > >> Thanks for the reply. I tested the change and it did not work, as >> get_generalized_face_support_points() is a

Re: [deal.II] InternalError in VectorTools::compute_nonzero_normal_flux_constraints

2017-02-20 Thread Lars Corbijn van Willenswaard
Thanks for the reply. I tested the change and it did not work, as get_generalized_face_support_points() is also empty. The cause seems that FESystem doesn't initialize the generalized_face_support_points. To test that I added some code to initialize it when initialize_unit_face_support_points() 'fa