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 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? 

On Friday, January 6, 2023 at 7:29:17 PM UTC+1 Wolfgang Bangerth wrote:

>
> > 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 subset of vector components). I will see if I can at least add 
> assertions for some simple cases.
>
> Best
> W.
>
> -- 
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bang...@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/582227f6-425d-4a06-acbc-5db23dcea79bn%40googlegroups.com.

Reply via email to