I added an exact solution object in adaptivity_ex2 with ExactSolution exact_sol(equation_systems); exact_sol.attach_exact_value(exact_solution);
Then I computed errors with exact_sol.compute_error("Convection-Diffusion", "u") and printed out L2 error in each refinement step. The output for the first time step is: Time = 0.025, refinement step = 0, elements = 10, l2_error = 0.443873 Time = 0.025, refinement step = 1, elements = 40, l2_error = 0.045196 Time = 0.025, refinement step = 2, elements = 160, l2_error = 0.131169 Time = 0.025, refinement step = 3, elements = 640, l2_error = 0.116789 Time = 0.025, refinement step = 4, elements = 2560, l2_error = 0.118175 I am curious why sometimes L2 error gets bigger, e.g., from r_step 1 to r_step 2. Don't more refinements give smaller errors? BTW, I used uniformly_refine() in each refinement step. The input mesh has 8 triangles and 2 quadrilaterals. --Junchao Zhang ------------------------------------------------------------------------------ Site24x7 APM Insight: Get Deep Visibility into Application Performance APM + Mobile APM + RUM: Monitor 3 App instances at just $35/Month Monitor end-to-end web transactions and take corrective actions now Troubleshoot faster and improve end-user experience. Signup Now! http://pubads.g.doubleclick.net/gampad/clk?id=267308311&iu=/4140 _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users