Re: [Libmesh-devel] Doing System::update in EquationSystems::read?

[Lorenzo Botti]

Hi Roy, it seems that ES already has a call to update() in the end of read(). This should avoid the user to call  System::update(). But if you are trying to restart a transient system in parallel the problem is the lack of an update() method for such a system. I posted a message on the list about this problem some time ago. Here is it. Hi All, I have finally find out why I couldn't restart properly a transient system in parallel. The problem is that transient system doens't have an update(...) method. So I have added one coping it from the re_update(...) method. It works for me but it will be good if someone could check it and add it to the svn. Another problem is that we should avoid the reading of the RHS_Vector becuase as I said in a previous message Petsc doesn't allowed to set the value of a Petsc_vector an then zero it without calling close(). This now is happening in ImplicitSystem::init(..). I don't have a good solution for this... of course we can call close() before zero() but this is an horrible hack. Any idea? Thanks. Lorenzo

transient_system.C
Description: Binary data

transient_system.h
Description: Binary data

Lorenzo Il giorno 20/giu/08, alle ore 19:14, Roy Stogner ha scritto: In ex10.C (or anywhere we load a restart file), we have to follow up an EquationSystems::read() call with a System::update() for every system in the EquationSystems object.  Would it make sense for the library to do that job itself in ES::read?  It seems to me that users are never going to want their vectors left in an inconsistent state right after a read. (yes, I just got bit by the discovery that my own restart code was only working in serial; why do you ask?) --- Roy ------------------------------------------------------------------------- Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://sourceforge.net/services/buy/index.php _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel

-------------------------------------------------------------------------
Check out the new SourceForge.net Marketplace.
It's the best place to buy or sell services for
just about anything Open Source.
http://sourceforge.net/services/buy/index.php
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel