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
