On Wed, Mar 27, 2013 at 3:32 PM, Roy Stogner <[email protected]>wrote:
> > I appreciate the persistence, and I believe you're seeing a problem, > but this proposed fix just looks wrong. It's essentially doing the > same current->old solution copy of an advance_timestep() without the > user having called advance_timestep(), and without updating > system.time, neither of which should ever happen. > > Any chance you could distill a problem code into something you could > send me? There's got to be something different between your > time+adaptivity loops and mine. > I had essentially started with the transient Navier-Stokes system example, and had modified it to do what I needed. The behavior that I was seeing was that after the call to UnsteadySolver::reinit(), the old_local_nonlinear_solution was zero (as opposed to projecting the solution before refinement), and this was directly used in calculating the dU in euler_solver. I will try to create a small example code and send it your way. Manav ------------------------------------------------------------------------------ Own the Future-Intel® Level Up Game Demo Contest 2013 Rise to greatness in Intel's independent game demo contest. Compete for recognition, cash, and the chance to get your game on Steam. $5K grand prize plus 10 genre and skill prizes. Submit your demo by 6/6/13. http://p.sf.net/sfu/intel_levelupd2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
