I'm using libmesh to solve a bunch of poisson-like problems on the same mesh but with different Dirichlet boundary conditions on (say) one side. I would like to use a DirchletBoundary object to impose these conditions, one iteration at the time. However I would like to perform all calculations using a single system.
I have two questions that probably have related answers: 1) My DirichletBoundary object is defined in my derived system's init_data() member function, before Parent::init_data(). init_data() is called when I call init() on the EquationSystem that owns the system, but how can I reinitialize my DirichletBoundary object? 2) If I call init() on the EquationSystems repeatedly, it is as if the Dirichlet conditions that were imposed at previous iterations are not deleted. For example, if I at one iteration impose a zero Dirichlet condition, the solution at all later iteration will satisfy that condition. How can I remove the existing constraints? Thanks a lot! Jens Eftang ------------------------------------------------------------------------------ For Developers, A Lot Can Happen In A Second. Boundary is the first to Know...and Tell You. Monitor Your Applications in Ultra-Fine Resolution. Try it FREE! http://p.sf.net/sfu/Boundary-d2dvs2 _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
