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.

Reply via email to