Kyusik, [...] > 3. I calculated l2 norm > Difference = solution; > Difference -= exact_solution; > norm = Difference.l2norm(); > > But, according to what you mentioned, I calculated l2norm of U_i-U_i_exact > not l2norm of u_h(x)-u_h_exact(x) > Yes, this gives you the l2-norm of U_i-U_i_exact and not the L2-norm of u_h(x)-u_h_exact(x). VectorTools::integrate_difference[1] allows you to compute the L2-norm. This function returns a Vector containing the error on each cell. The l2-norm of that Vector is the L2-norm you are looking for. Have also a look at step-7.
Best, Daniel [1] https://www.dealii.org/8.4.1/doxygen/deal.II/namespaceVectorTools.html -- 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.
