> Thank you for your reply. Where the discontinuities are depends on the > initial_refinement integer. When this was 6 it produced the figure I sent > originally - sorry I should have said this. This is what I found > particularly odd, that for different levels of refinement, the behaviour > is different, as the VectorTools::projection occurs after the refinement > of the grid.
It is an optical illusion. You get to see the overshoots if you generate the plot with data_out.build_patches(2) or higher because are using quadratic elements. The reason why the overshoots are different in x- and y-direction is because your mesh is stretched in one direction (i.e., delta-x is larger than delta-y) and consequently the projection is anisotropic as well. So I believe that all is working as it should be. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
