My fault, I was treating the vector as if I was using a lagrangian basis. Solved doing L2 projections, sorry for bothering :)
On 05/03/2017 03:16 PM, Giacomo Rosilho de Souza wrote: > Hello, > I'm implementing discontinuous galerkin following your example and I > have a strange visualization problem (maybe not only visualization). The > code seems to work well, I did the convergence tests and they give good > results. But when I modify a vector "by hand", iterating over its > components, I obtain strange outputs. For example, if I iterate over the > vector components and set them to 1, when I visualize it in paraview, > some nodes have value 2. But, if I print the L infinity norm of the > vector I get 1. The same happens if I do v.zero(); and then v.add(1.);. > Instead, the vectors coming from the solution of a linear system are > visualizing right. I'm doing something wrong? Here below theres a link > for a working example. > > Thanks :) > > https://www.dropbox.com/s/c10xhrwfdseoolr/example.C?dl=0 > > > ------------------------------------------------------------------------------ > Check out the vibrant tech community on one of the world's most > engaging tech sites, Slashdot.org! http://sdm.link/slashdot > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
