Yeah, I have actually the same mesh (these are two consecutive steps of a Newton method) so I should use:
Real normh1 = system.calculate_norm(uh - vh,H1); ????? uh and vh are NumericVector<Number> type right? Thanks guys On Mon, Oct 6, 2014 at 3:30 PM, Paul T. Bauman <[email protected]> wrote: > If you have the exact same mesh and discretization for u_h and v_h, you > can use what Roy pointed out: > > system.calculate_norm(u_h - v_h,...) > > > http://libmesh.sourceforge.net/doxygen/classlibMesh_1_1System.php#a3153155431e918a782fabedcaf58e4e6 > > > http://libmesh.sourceforge.net/doxygen/classlibMesh_1_1System.php#a1105a347b46aceaefb97e5e97ce7802d > > If they are different, I would suggest attaching one to the ExactSolution > object > > > http://libmesh.sourceforge.net/doxygen/classlibMesh_1_1ExactSolution.php#a8aaee778ee07a4c872cdba1e4e538329 > > And use that to compute the error. (There may be better ways, but that's > what occurred to me first). > > On Mon, Oct 6, 2014 at 9:10 AM, Rodrigo Broggi <[email protected]> > wrote: > >> It is almost that, actually I would like to compute the h1norm error >> between two approximate functions, not between an exact and an approximate >> ( ||uh - vh||_h1). Does this class work also for this case? >> Il 06/ott/2014 15:06 "Rodrigo Broggi" <[email protected]> ha scritto: >> >> > It is almost that, actually I would like to compute the h1norm error >> > between two approximate functions, not between an exact and an >> approximate >> > ( ||uh - vh||_h1). Do this class works also for this case? >> > Il 06/ott/2014 14:10 "Paul T. Bauman" <[email protected]> ha scritto: >> > >> >> Not sure exactly what you're looking for, but there are facilities to >> >> compute L2/H1/etc. norms of || u - u_h ||: >> >> >> >> >> http://libmesh.sourceforge.net/doxygen/classlibMesh_1_1ExactSolution.php >> >> >> >> http://libmesh.sourceforge.net/adaptivity_ex3.php >> >> >> >> Is this what you had in mind? >> >> >> >> On Mon, Oct 6, 2014 at 5:39 AM, Rodrigo Broggi <[email protected]> >> >> wrote: >> >> >> >>> Hi, >> >>> >> >>> A simple question: is there any facility/method to calculate the h1 >> norm >> >>> of >> >>> a solution as there is to calculate te l2_norm of a vector? >> >>> >> >>> Thanks, >> >>> >> >>> Rodrigo Broggi >> >>> >> >>> >> ------------------------------------------------------------------------------ >> >>> Slashdot TV. Videos for Nerds. Stuff that Matters. >> >>> >> >>> >> http://pubads.g.doubleclick.net/gampad/clk?id=160591471&iu=/4140/ostg.clktrk >> >>> _______________________________________________ >> >>> Libmesh-users mailing list >> >>> [email protected] >> >>> https://lists.sourceforge.net/lists/listinfo/libmesh-users >> >>> >> >> >> >> >> >> ------------------------------------------------------------------------------ >> Slashdot TV. Videos for Nerds. Stuff that Matters. >> >> http://pubads.g.doubleclick.net/gampad/clk?id=160591471&iu=/4140/ostg.clktrk >> _______________________________________________ >> Libmesh-users mailing list >> [email protected] >> https://lists.sourceforge.net/lists/listinfo/libmesh-users >> > > ------------------------------------------------------------------------------ Slashdot TV. Videos for Nerds. Stuff that Matters. http://pubads.g.doubleclick.net/gampad/clk?id=160591471&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
