On Tue, 5 Apr 2011, Derek Gaston wrote:
> On Apr 5, 2011, at 3:42 PM, Roy Stogner wrote: > >> Okay. But the problem isn't that uniformly_refine() needs to call >> clean_refinement_flags() - that would fix the >> uniformly_refine-then-read-then-reinit case, but it would break the >> read-then-uniformly_refine-then-reinit case; those flags are the only >> way reinit() knows it needs to do a projection from the coarse mesh! >> >> I believe the right fix would be a clean_refinement_flags() at the end >> of the step which actually makes that projection unnecessary - the >> EquationSystems::read(). Give that a try? > > Did you mean EquationSystems::init()??? If so... I added a call to > clean_refinement_flags() at the end of init() (I had to make it > public and create a MeshRefinement object in init() ) and that seems > to have worked. > > What do you think? This sounds good. I think it's only half of the fix we need, though: I meant *both* read() and init(). I think either way should have the same problem. read() ends up calling init() if READ_HEADER is set, so we should be okay there, but we probably need to do the clean_refinement_flags() by hand in the READ_DATA && !READ_HEADER case. --- Roy ------------------------------------------------------------------------------ Xperia(TM) PLAY It's a major breakthrough. An authentic gaming smartphone on the nation's most reliable network. And it wants your games. http://p.sf.net/sfu/verizon-sfdev _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
