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 problem, I need to do two non_linear iterations to
converge instead of 1.
If I can't use VectorTools::interpolate_boundary_values what are my options?
Does VectorTools::project_boundary_values work in that case?
No, it also needs to know which DoFs are located on the boundary, but for DG
elements no DoFs are logically on the boundary.
I don't see a good solution to your problem that uses what's already in the
library. I imagine that using code similar to what's in
interpolate_boundary_values() but using a different way of finding out which
DoF's support point is physically (though not logically) located on the
boundary, this all could be implemented in 200 lines of code. But there is
nothing simple the obvious already provides you for this case.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/75e595c3-deba-c1e6-5cb4-9b6fb964d7b5%40colostate.edu.
