Hi all (I think I accidentally sent this to the old mailing list yesterday, so apologies if anyone is receiving it twice),
I'm using a 1D model to try to test out an a posteriori estimator (I figure 1D will make it easier to identify whether the estimator is more effective than a generic one, and also I have some nice quasi-analytic 1D solutions). I realized that there is no mesh regularization implemented for the call to the 1-D version of prepare_coarsening_and_refinement(), which could degrade accuracy, right? I figure it should not be so difficult to implement the regularization step in 1D on my own, but I was wondering if I could get some assistance from someone who has done it or thought about it? Would I do a similar reverse sweep, from highest down to (2nd) lowest level, to the one described in the docs, where at each level I check that if a cell's level after refinement would be 2 higher than that of its neighboring cell, then the neighboring cell is marked for refinement? And I noticed the fix_coarsen_flags() function is called even in 1D.. would it just need to be called once after the aforementioned sweep? That's actually the part I'm concerned about. Thanks very much, Dan _______________________________________________
