Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-03-13 Thread Abbas
Thank you Prof. Bangerth. I am always greatfull for your input. On Friday, March 10, 2023 at 10:27:49 PM UTC+1 Wolfgang Bangerth wrote: > On 3/8/23 07:49, Abbas wrote: > > > > I am solving a nonlinear problem with DG and my BCs are in the weak > form. It > > would be helpful if I start with a

Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-03-10 Thread Wolfgang Bangerth
On 3/8/23 07:49, Abbas wrote: I am solving a nonlinear problem with DG and my BCs are in the weak form. It would be helpful if I start with a solution vector that is closer to the prescribed solution so I won't have to do iterations otherwise would have been unnecessary. For a purely linear p

Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-03-08 Thread Abbas
I would still want to initialise my solution vector to the values which are prescribed at the boundary before solving. Context: I am solving a nonlinear problem with DG and my BCs are in the weak form. It would be helpful if I start with a solution vector that is closer to the prescribed sol

Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-01-06 Thread Wolfgang Bangerth
It's cool how dealii is consitent with the maths. Maybe @Wolfgang would agree that this should have thrown an exception in Debug? In principle yes, but I don't know how to write such an assertion given the generality of the function (for the hp and non-hp case, and when selecting only a sub

Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-01-04 Thread Abbas
Thank you Daniel. I was using penalties before taking this route, because I had artifacts at the boundary faces present at a corner between a Dirichlet and a Neuman BC edge. A problem I didn't have with CG. My first intuition was to apply BCs strongly to fix this, but I don't think it will ev

Re: [deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-01-04 Thread Daniel Arndt
Abbas, Yes, FE_DGQ only defines dofs in the interior of cells. Thus, trying to interpolate the boundary has no effect. You would normally use penalty terms to prescribe boundary values for DG methods. Best, Daniel On Wed, Jan 4, 2023 at 6:04 PM Abbas wrote: > Hello, > > The VectorTools::interp

[deal.II] VectorTools::interpolate_boundary_values doesn't work for DG

2023-01-04 Thread Abbas
Hello, The VectorTools::interpolate_boundary_values function works for interpolating fe_Q functions but not fe_DGQ . Was this meant to be the case? I attached a minimal code that illustrates this. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options,