Roy, I was looking at this thread since Derek has asked me to look into this problem a little more and noticed that your original patch to fix the problem was all written within fe_base dealing with constraints. I'm no expert but this doesn't seem like enough when you consider cascading effects that show up when you enforce the level one rule. Instead I was thinking about fixing the problem by essentially applying your interior AMR algorithms inside of the MeshRefinment class across the periodic boundaries. Specifically, I was considering making changes to the make_(?:refinement|coarsening)_compatible functions. Each of these functions contains a loop that runs to maintain the level one rule by adjusting the refinement flags each time the mesh is coarsened and/or refined. All of the right logic is present to force the right refinement patterns throughout the interior mesh elements already. Think about a worst case scenario where you have refinement happening right on one periodic boundary and then "just inside" of the mirrored boundary. If you take some naive approach and only work with flags on the boundary you start to see unrefined patches of elements showing up or other interior pattern cases which again are already handled for interior elements. We just need to make the elements on the periodic boundaries "see" their respective mirrored elements as neighbors for the purpose of this algorithm. This could be done by adding in the appropriate logic right inside of the respective "_maintain_level_one" loops. Does this seem reasonable? I'd be happy to take a crack at this unless you have any other ideas.
Thanks, Cody On Oct 30, 2010 4:14am, Ming Q <[email protected]> wrote: > Hi again, > It seem the periodic BC functions is not working with the current SVN > version. > I resend the test case I have sent before. > Regards, > / Ming Q. ------------------------------------------------------------------------------ Nokia and AT&T present the 2010 Calling All Innovators-North America contest Create new apps & games for the Nokia N8 for consumers in U.S. and Canada $10 million total in prizes - $4M cash, 500 devices, nearly $6M in marketing Develop with Nokia Qt SDK, Web Runtime, or Java and Publish to Ovi Store http://p.sf.net/sfu/nokia-dev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
