[Trilinos-Users] Question on Epetra_BlockMap with variable element length

Littlewood, David John djlittl at sandia.gov
Thu Jan 27 16:53:27 MST 2011


Thanks for your nearly instantaneous response!

Your suggestion makes sense, but unfortunately it is not compatible with other parts of the code that operate directly on the underlying arrays of doubles and assume that there will be no entries in the array for points with no neighbors.  This perhaps points to shortcomings in our code that have nothing to do with Epetra, we'll have to review what we're doing there.

For the time being I'm just going to omit points with no neighbors when constructing the Epetra_BlockMap for the storage of neighbor-related data.  This works but has the unfortunate consequence of making the number of entries in the map for bond data different from the number of entries in the map for point-wise data.  Not a huge disaster, but the bookkeeping becomes a bit ugly.

Best Regards,

On 1/27/11 11:23 AM, "Heroux, Mike" <MHeroux at CSBSJU.EDU> wrote:


BlockMap was not designed to handle a block entry size of 0.  I am not sure what would be required to do so.  Might it be possible to keep the ID of the self, then check the ID list and ignore the self ID?

This would allow you to proceed.


On 1/27/11 12:11 PM, "Littlewood, David John" <djlittl at sandia.gov> wrote:


When constructing an Epetra_BlockMap with a variable element length, is it possible to have an element with length zero?

I am using a map to store the list of neighbors in a peridynamics simulation, and it is possible to have zero neighbors.  It appears that one cannot construct an Epetra_BlockMap with an element that has zero length (throws exception).


Trilinos-Users mailing list
Trilinos-Users at software.sandia.gov

-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://software.sandia.gov/pipermail/trilinos-users/attachments/20110127/7c88d122/attachment.html 

More information about the Trilinos-Users mailing list