On Tue, 11 May 2010, David Knezevic wrote: > On 05/11/2010 12:44 PM, Roy Stogner wrote: >> >> So it's asserting that the solution vector hasn't been closed after >> the solve? That mesh_position_get() call would have unclosed the >> solution vector, but the solver should have closed it again. > > I think it's asserting that the solution isn't closed before the solve.
Ah, of course! UnsteadySolver has to fill the old_nonlinear_solution before the first solve; that operator= is where it's failing. > Adding this->solution->close(); in FEMSystem::mesh_position_get() resolves > this. And that's the right solution; I'll commit it to SVN. > So now I'm getting mesh motion, though if GMV is to be believed it's only the > outermost layer of elements that are moving. I'll have a closer look at > what's happening... Thanks! --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
