On Fri, Jan 29, 2016 at 5:05 PM, Roy Stogner <royst...@ices.utexas.edu> wrote:
> > On Fri, 29 Jan 2016, Junchao Zhang wrote: > > I tested the modified code with uniform refinement and >> exact_sol.extra_quadrature_order(1). I printed out both l2 error and h1 >> error and got >> >> time = 0.025, refinement step = 0, elements = 10240, l2_error = >> 0.002885, h1_error = 0.239394 >> Time = 0.025, refinement step = 1, elements = 40960, l2_error = >> 0.002643, h1_error = 0.207076 >> Time = 0.025, refinement step = 2, elements = 163840, l2_error = >> 0.002596, h1_error = 0.210480 >> Time = 0.025, refinement step = 3, elements = 655360, l2_error = >> 0.002575, h1_error = 0.214309 >> Time = 0.025, refinement step = 4, elements = 2621440, l2_error = >> 0.002569, h1_error = 0.215482 >> >> We can see l2_error monotonically decreases, however h1_error does not. >> What is the possible reason? >> > > Both should be decreasing much faster than that. You're refining the > initial mesh, before the initial conditions are projected? > I don't understand the question. The mesh is uniformly refined at https://github.com/libMesh/libmesh/blob/master/examples/adaptivity/adaptivity_ex2/adaptivity_ex2.C#L212 Then the initial condition is computed at https://github.com/libMesh/libmesh/blob/master/examples/adaptivity/adaptivity_ex2/adaptivity_ex2.C#L233 > --- > Roy ------------------------------------------------------------------------------ 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
