> Things would be much easier if deal.ii's Triangulation class supported > the notion of (topologically) identical vertices with different > coordinates. Otherwise, combining periodic boundary conditions with mesh > refinement seems to require significant amounts of code on the > application side.
Yes, I agree with this. I've seen this in another library recently and thought it was a very cool feature. I wished we had thought of this 10 years ago :-( > I use prepare_coarsening_and_refinement for grid smoothing and then fix > any discrepancy in refine/coarsen requests between the time boundaries. > This is repeated until the refine/coarsen flags do not change any more. I think this is the way to go. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________
