Hi Daniel, > Second - maybe I am misunderstanding terminology? Regularization, I > thought, is making sure that a cell is no more than once more refined > than any of its neighbors? I thought that, in 2D and 3D, > prepare_coarsening_and_refinement () does this even with none of the > flags set that you mentioned below..? I am talking about the case when > all flags in the member smooth_grid are false. I looked at Triangulation and found that terminology is a bit odd: making a grid one-irregular, that is, not more than one one-sided refinement of a face is not considered regularization there, it is what we do anyway. Except in 1D, where this is not necessary for grid integrity in our sense. Therefore, this step is not done in 1D.
The regularization flags control additional smoothing of refinement. Some of these do make sense in 1D and you are most welcome to implement them. We's also be grateful if you would provide your results for incorporation in the library. And right, no need to add anisotropic refinement in 1D. Best, Guido _______________________________________________
