On Wed, 18 Feb 2009, Tim Kroeger wrote: > What's the current state about the ghosted vectors?
They appear to be working for stationary problems, but there's definitely a bug which corrupts transient results - possibly in the System::project_vector code, possibly in the enforce_constraints_exactly it calls, probably in both since the ghost-specific code is very similar. > Are you waiting for some input from me, or did you just not have > time in the meanwhile? I'm swamped right now. The next step is to reduce ex10 to a simpler case (maybe solving du/dt = 0 in nD with some non-zero ICs?) where I can step through the whole projection in gdb, but it's likely to be weeks before I can get to that. --- Roy ------------------------------------------------------------------------------ Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, CA -OSBC tackles the biggest issue in open source: Open Sourcing the Enterprise -Strategies to boost innovation and cut costs with open source participation -Receive a $600 discount off the registration fee with the source code: SFAD http://p.sf.net/sfu/XcvMzF8H _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel