Hi, I am studying a time dependent PDE on a square domain and want to compare exactness of the solution to the number of dofs. I use the Kelly estimator for local refinement. For testing the program, I use a parameter setting allowing for a rotational symmetric analytic solution at every time. Now I observe that the locally refined mesh is by far NOT rotational symmetric, but rather symmetric along the coordinate axes, although the (analytic) Kelly estimator is rotational symmetric for rotational symmetric functions. Now I wonder if anyone has experienced similar problems? Is there an easy way to get around this problem which seemingly comes from the square grid geometry? Right now I am getting correct solutions, but local adaptation works so poorly that there is nearly no advantage compared to uniform meshes.
Thanks, Jochen _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
