On Tue, 10 Mar 2009, Roy Stogner wrote: > On Tue, 10 Mar 2009, Roy Stogner wrote: > >> if I can't find the bug right away. > > Found it.
Great. Now, the application crashes at the "original" crash point. It happens in System::project_vector() for the case of a ghosted vector. This calls DofMap::enforce_constraints_exactly(), which in line 680 of dof_map_constraints.C calls NumericVector::operator(). It hits the libmesh_assert() near the end of PetscVector::map_global_to_local_index(), i.e. the index supplied to operator() is neither a local nor a ghost index. There are 7 systems in the application (all on the same grid). Two ExplicitSystem's and two LinearImplicitSystem's are projected successfully, but the first TransientLinearImplicitSystem triggers the crash, when projecting its _transient_old_local_solution. That system contains two variables, both of the same FE-type. To reproduce it more easily, I wrote the grid to an .xdr file directly before the call to MeshRefinement::refine_and_coarsen_elements(), plus a file containing the refinement flags of all active elements (stored in the native order, i.e. the order that the iterator steps through the elements). Then, I wrote a test program that reads in the grid, initializes a TransientLinearImplicitSystem with two variables on that mesh, reads the refinement flags, and calls MeshRefinement::refine_and_coarsen_elements() and EquationSystems::reinit(). I ran this program on the same number of processors as the main application (that is 8) -- but that program does *not* crash. I think some very odd things must be going on here. Do you have any idea how to track this down further? 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