I'm sorry, please let me ask you one more question.... > Why is that expected? It should be the same. The new sparsity_pattern > should be the same as the old one > > > According to your comment, I remove system_matrix.reinit(sparsity_pattern) only (system_rhs.reinit and solution.reinit are not removed)
The result is as follows... Number of active cells: 5120 Number of degrees of freedom: 20609 iteration 0 ||u_{k}-u_{k-1}||42.5271 iteration 1 ||u_{k}-u_{k-1}||30.1063 iteration 2 ||u_{k}-u_{k-1}||11.1065 iteration 3 ||u_{k}-u_{k-1}||0.159375 iteration 4 ||u_{k}-u_{k-1}||0.0176957 iteration 5 ||u_{k}-u_{k-1}||0.0149181 iteration 6 ||u_{k}-u_{k-1}||0.0129168 iteration 7 ||u_{k}-u_{k-1}||0.0113656 iteration 8 ||u_{k}-u_{k-1}||0.0101467 iteration 9 ||u_{k}-u_{k-1}||0.00916819 It seems that the error is greatly reduced... So, It seems that I can't ignore this difference. Could you please explain to me what makes this big difference?? Thank you Kyusik. -- 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.