David:

I have two distinct distributed meshes representing the same spatial domain and (for now) with equal refinement (number of cells),

Do I understand correctly that these meshes are identical, but that you store them in separate variables? Why not just one mesh that two DoFHandlers build on?


The violated condition was:

    partitioner->in_local_range(global_index) || partitioner- >ghost_indices().is_element(global_index)

Additional information:
    You tried to access element 54 of a distributed vector, but this     element is not stored on the current processor. Note: The range of
     locally owned elements is [135,242], and there are 27 ghost elements
     that this vector can access.

     A common source for this kind of problem is that you are passing a
     'fully distributed' vector into a function that needs read access to
     vector elements that correspond to degrees of freedom on ghost cells
     (or at least to 'locally active' degrees of freedom that are not also
     'locally owned'). You need to pass a vector that has these elements as
     ghost entries.

Before we go into the weeds, is this true in your case? Is the vector you are passing into the FEFieldFunction (where the error originates) ghosted?

Best
 Wolfgang

--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
https://groups.google.com/d/forum/dealii?hl=en
--- You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email 
to dealii+unsubscr...@googlegroups.com.
To view this discussion visit 
https://groups.google.com/d/msgid/dealii/5ac8b3da-c1f6-4d9b-9933-cdb9608c6d26%40colostate.edu.

Reply via email to