Kyusik, On Sunday, March 5, 2017 at 10:53:29 PM UTC-5, hanks0...@gmail.com wrote: > > I have a question in the lecture 31.65 that is about Picard iteration. > > In the lecture, the professor doesn't include system_matrix.reinit > (sparsity_pattern), system_rhs.reinit (dof_handler.n_dofs()), and > solution.reinit (dof_handler.n_dofs()) in assemble_system(). > > However, I was worried that these three values might not be substituted > with new values(after every iteration). > > So I added these three reinit in assemble_system(). > > The two results are a little bit different as expected. > Why is that expected? It should be the same. The new sparsity_pattern should be the same as the old one and you still have the same number of degrees of freedom. That's why there is no need to do the reinit. The values of the matrix and the vectors should be reset but not their size.
Best, Bruno -- 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. For more options, visit https://groups.google.com/d/optout.
