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.