On Wed, 11 Mar 2009, John Peterson wrote:
On Wed, Mar 11, 2009 at 12:52 PM, Tim Kroeger
but the first TransientLinearImplicitSystem triggers the crash, when projecting its _transient_old_local_solution. variables, both of the same FE-type.All variables are first order Lagrange. Hence, the systems differ only in the system type and the number of variables: No. Type #Vars ----------------------------------------- 0 ExplicitSystem 5 1 LinearImplicitSystem 1 2 LinearImplicitSystem 1 3 ExplicitSystem 1 4 TransientLinearImplicitSystem 2 5 TransientLinearImplicitSystem 2 6 TransientLinearImplicitSystem 1Hmm... the biggest difference with the TransientLinearImplicitSystem would probably be the "old_local_solution" and "older_local_solution" vectors... I'm guessing you guys have already investigated those though?
As said above, old_local_solution is the vector that crashes first. But I don't know what's so different with this vector from current_local_solution.
Wait... there *is* a difference! System::current_local_solution is not projected at all. Rather, only System::solution is projected (employing a temporary parallel vector) and after that localized to System::current_local_solution.
Hence, Roy's suspicion that dofs required via constraints are not explicitly included in the send list would completely explain the problem.
By the way, why doesn't System::project_vector() in case of a parallel vector use a ghosted (rather than a serial) as the temporary?
Best Regards, Tim -- Dr. Tim Kroeger tim.kroe...@mevis.fraunhofer.de Phone +49-421-218-7710 tim.kroe...@cevis.uni-bremen.de Fax +49-421-218-4236 Fraunhofer MEVIS, Institute for Medical Image Computing Universitaetsallee 29, 28359 Bremen, Germany
------------------------------------------------------------------------------ Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are powering Web 2.0 with engaging, cross-platform capabilities. Quickly and easily build your RIAs with Flex Builder, the Eclipse(TM)based development software that enables intelligent coding and step-through debugging. Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com
_______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
