On Fri, 30 May 2014, Robert wrote: > Are there any additional step necessary in each timestep to apply > the boundary conditions or what else am I missing?
One little step, one big one: To apply the boundary conditions, you'll want to run EquationSystems::reinit() before each time step. Don't trust the accuracy without testing. If you're not updating System::time to match the *end* of your time step rather than the beginning when reinit gets called, then your functions will get called with the wrong time and you'll have an O(deltat) error from that. --- Roy ------------------------------------------------------------------------------ Time is money. Stop wasting it! Get your web API in 5 minutes. www.restlet.com/download http://p.sf.net/sfu/restlet _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
