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

Reply via email to