Re this: So I want to have an EquationSystems that holds two systems, one that is read from file, and one that is not.
I tried a couple of things: Read system from file, then adding another system, then init() and add a system, then init(), and then read from file. None of this works and seems to, when I print_info(), duplicate the number of dofs in whatever system is inited twice. Also, if I add a vector to this system the size of that vector is twice as large as it should be. What does seem to work is to first add a system, and then read from file, and not call init() at all. Is this the right way to do this? I guess it makes sense since read() also initializes ... Best, Jens On 01/25/2013 02:02 PM, Jens Lohne Eftang wrote: > This was not quite accurate ... it should be: > > It seems that if I have an EquationSystems es and do es.read( ... ) > followed by es.init() followed by es.print_info(), the n_dofs in the > read system are counted twice. > > I know that read also does initialization, but I want to add more > systems to es, and hence I need to call es.init() after I've added > these additional systems. > > Any idea? > > Best, > Jens > > On 01/25/2013 02:00 PM, Jens Lohne Eftang wrote: >> All, >> >> It seems that if I have an EquationSystems es and do es.read( ... ) >> followed by es.init() or es.reinit() followed by es.print_info(), the >> n_dofs in the read system are counted twice. >> >> I know that read also does initialization, but I want to add more >> systems to es, and hence I need to call es.init() or es.reinit() >> after I've added these additional systems. >> >> Any idea? >> >> Best, >> Jens > ------------------------------------------------------------------------------ Master Visual Studio, SharePoint, SQL, ASP.NET, C# 2012, HTML5, CSS, MVC, Windows 8 Apps, JavaScript and much more. Keep your skills current with LearnDevNow - 3,200 step-by-step video tutorials by Microsoft MVPs and experts. ON SALE this month only -- learn more at: http://p.sf.net/sfu/learnnow-d2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
