[deal.II] Re: Use of 'VectorTools::interpolate_boundary_values' using numerical solution.

Hi Jaekwang One way I can think of is to construct your ConstraintMatrix for Part B manually. Iterate over your triangulation, find all the dofs on the target boundaries, then use ConstraintMatrix::add_line and ConstraintMatrix::set_inhomogeinity to set the BC. Another method might be putting

[deal.II] Different time steps lead to different asymptotic values for the residual in a heat equation

This question might be related to earlier of my questions, but finally I could not find an explanation for the behaviour yet, thus this question. I intended to calculate the time-dependent heat equation with constant factors, using the newton method (such that time-dependent factors can be